Skip to main content

ORIGINAL RESEARCH article

Front. Comput. Neurosci., 02 June 2015
This article is part of the Research Topic Application of Nonlinear Analysis to the Study of Complex Systems in Neuroscience and Behavioral Research View all 34 articles

Detection of subjects and brain regions related to Alzheimer's disease using 3D MRI scans based on eigenbrain and machine learning

\r\nYudong Zhang*Yudong Zhang1*Zhengchao DongZhengchao Dong2Preetha PhillipsPreetha Phillips3Shuihua Wang,Shuihua Wang1,4Genlin Ji,Genlin Ji1,5Jiquan YangJiquan Yang5Ti-Fei Yuan*Ti-Fei Yuan6*
  • 1School of Computer Science and Technology, Nanjing Normal University, Nanjing, China
  • 2Division of Translational Imaging and MRI Unit, New York State Psychiatric Institute, Columbia University, New York, NY, USA
  • 3School of Natural Sciences and Mathematics, Shepherd University, Shepherdstown, WV, USA
  • 4School of Electronic Science and Engineering, Nanjing University, Nanjing, China
  • 5Jiangsu Key Laboratory of 3D Printing Equipment and Manufacturing, Nanjing, China
  • 6School of Psychology, Nanjing Normal University, Nanjing, China

Purpose: Early diagnosis or detection of Alzheimer's disease (AD) from the normal elder control (NC) is very important. However, the computer-aided diagnosis (CAD) was not widely used, and the classification performance did not reach the standard of practical use. We proposed a novel CAD system for MR brain images based on eigenbrains and machine learning with two goals: accurate detection of both AD subjects and AD-related brain regions.

Method: First, we used maximum inter-class variance (ICV) to select key slices from 3D volumetric data. Second, we generated an eigenbrain set for each subject. Third, the most important eigenbrain (MIE) was obtained by Welch's t-test (WTT). Finally, kernel support-vector-machines with different kernels that were trained by particle swarm optimization, were used to make an accurate prediction of AD subjects. Coefficients of MIE with values higher than 0.98 quantile were highlighted to obtain the discriminant regions that distinguish AD from NC.

Results: The experiments showed that the proposed method can predict AD subjects with a competitive performance with existing methods, especially the accuracy of the polynomial kernel (92.36 ± 0.94) was better than the linear kernel of 91.47 ± 1.02 and the radial basis function (RBF) kernel of 86.71 ± 1.93. The proposed eigenbrain-based CAD system detected 30 AD-related brain regions (Anterior Cingulate, Caudate Nucleus, Cerebellum, Cingulate Gyrus, Claustrum, Inferior Frontal Gyrus, Inferior Parietal Lobule, Insula, Lateral Ventricle, Lentiform Nucleus, Lingual Gyrus, Medial Frontal Gyrus, Middle Frontal Gyrus, Middle Occipital Gyrus, Middle Temporal Gyrus, Paracentral Lobule, Parahippocampal Gyrus, Postcentral Gyrus, Posterial Cingulate, Precentral Gyrus, Precuneus, Subcallosal Gyrus, Sub-Gyral, Superior Frontal Gyrus, Superior Parietal Lobule, Superior Temporal Gyrus, Supramarginal Gyrus, Thalamus, Transverse Temporal Gyrus, and Uncus). The results were coherent with existing literatures.

Conclusion: The eigenbrain method was effective in AD subject prediction and discriminant brain-region detection in MRI scanning.

Introduction

Alzheimer's disease (AD) is not a normal part of aging. It is a type of dementia that causes problems with memory, thinking, and behavior. Symptoms usually develop slowly and worsen over time. Symptoms may become severe enough to interfere with daily life, and lead to death (Hahn et al., 2013). There is no cure for this disease. In 2006, 26.6 million people worldwide suffered from this disease. AD is predicted to affect 1 in 85 people globally by 2050, and at least 43% of prevalent cases need high level of care (Brookmeyer et al., 2007). As the world is evolving into an aging society, the burdens and impacts caused by AD on families and the society has also increased significantly. In the US, healthcare on people with AD currently costs roughly $100 billion per year and is predicted to cost $1 trillion per year by 2050 (Miller et al., 2012).

Early and accurate detection of AD is beneficial for the management of the disease (Han et al., 2011). Presently, a multitude of neurologists and medical researchers have been dedicating considerable time and energy toward this goal, and promising results have been continually springing up (Xinyun et al., 2011). Magnetic resonance imaging (MRI) is an imaging technique that produces high quality images of the anatomical structures of the human body, especially in the brain, and provides rich information for clinical diagnosis and biomedical research (Shamonin et al., 2014). The diagnostic values of MRI are greatly enhanced by the automated and accurate classification of the MR images (Goh et al., 2014; Zhang et al., 2015a,b). It already plays an important role in detecting AD subjects from normal elder controls (NC) (Angelini et al., 2012; Smal et al., 2012; Nambakhsh et al., 2013; Hamy et al., 2014; Jeurissen et al., 2014).

In earlier cases, most diagnosis work was done to measure manually or semi-manually a priori region of interest (ROI) of magnetic resonance (MR) images, based on the fact that AD patients suffer more cerebral atrophy compared to NCs (Kubota et al., 2006; Anagnostopoulos et al., 2013). Most of these ROI-based analyses focused on the shrinkage of hippocampus and cortex, and enlarged ventricles (Pennanen et al., 2004). Somehow, the ROI-based methods suffer from some limitations. First, the methods focus on the ROIs need prior knowledge. Second, the accuracy of early detection depends heavily on the experiences of the examiners. Third, the mutual information among the voxels is difficult to operate (Xinyun et al., 2011; Lee et al., 2013). Finally, there is no evidence that other regions (except hippocampus and entorhinal cortex) did not provide any information related to AD. Also, the auto-segmentation of ROI is not feasible in practice, and examiners tend to segment the brain manually.

On the other hand, multivariate approaches that consider all the voxels in a scan as one observation offer an alternative method to ROI-based methods. The advantages of multivariate approaches are that they are data driven, which means that the analyses are fully based on the data without any prior knowledge and that the interactions among voxels and error effects are assessed statistically. However, multivariate approaches suffer from either the curse of dimension problem or the small sample size problems or the lack of the capability, to make statistical inferences about regionally specific changes (Álvarez et al., 2009b).

The Eigenbrain was an excellent multivariate approach that solves both the curse of dimensionality and the problems in small sample size. It was proposed by Alvarez et al. (2009a) and Lopez et al. (2009), and was applied on Single Photon Emission Computed Tomography (SPECT) images. In their research, the eigenbrain approach was shown to efficiently reduce the feature space from ~5 × 105 to only ~102, and therefore, was able to achieve excellent classification accuracy. In this study, we make a tentative test of applying eigenbrains in MRI scans for AD detection.

Support vector machine (SVM) has been arguably regarded as one of the most excellent classification methods in machine learning (Zhang and Wu, 2012a). Original SVMs are linear classifiers, and do not perform well on nonlinear data. Hence, we introduced in the kernel SVMs (KSVMs), which extends original linear SVMs to nonlinear SVM classifiers by applying the kernel function to replace the dot product form in the original SVMs (Gomes et al., 2012). Compared with the original plain SVM, the KSVMs allows one to fit the maximum-margin hyperplane in a transformed feature space (Garcia et al., 2010). The transformation may be nonlinear and the transformed space is high dimensional; thus although the classifier is a hyperplane in the high-dimensional feature space, it may be nonlinear in the original input space (Hable, 2012).

The aim of our study was to develop a novel classification system based on eigenbrain and machine learning, in order to grow a computer-aided diagnosis (CAD) system for the early detection of AD subjects and AD-related brain regions. Our goal was not to replace clinicians, but to provide an assisting tool. The rest of the paper was organized as follows: the next section reviewed relates literatures from two aspects: the extracted features and the classification methods. Section The Proposed Method describes the methodology of the proposed CAD. Section Experiments and Results contains the experiments and results. Section Discussion analyzes the reason behind the experiment results. Finally, Section Conclusion and Future Research is devoted to conclusion and future research. For ease in reading, the acronyms and their meanings of this study are listed in Table 12 in the appendix.

The contributions of the paper fell within the following five aspects: (i) We generalized the Eigenbrain to MR images, and proved its effectiveness; (ii) We proposed a hybrid eigenbrain-based CAD system that can not only detect AD from NC, but also detect brain regions that related to AD. (iii) We proved the proposed method had classification accuracy comparable to state-of-the-art methods, and the detected brain regions were in line with 16 existing literatures. (iv) We used inter-class variance (ICV) and Welch's t-test (WTT) to reduce redundant data; (v) We found POL kernel is better than linear and RBF kernel for this study.

Literature Review

In common convention, the automatic classification consisted of two stages: feature extraction and classifier construction. We reviewed over ten literatures, and analyzed themthrough the two stages.

Features of MR Images

Scholars have proposed numerous methods to extract various features1. Chaplot et al. (2006) used the approximation coefficients obtained by discrete wavelet transform (DWT). Maitra and Chatterjee (2006) employed the Slantlet transform, which is an improved version of DWT. Their feature vector of each image was created by considering the magnitudes of Slantlet transform outputs corresponding to six spatial positions that were chosen according to a specific logic. El-Dahshan et al. (2010) extracted the approximation and detail coefficients of 3-level DWT. Plant et al. (2010) used brain region cluster (BRC). They suggested to use information gain (IG) to rate the interestingness of a voxel, and applied clustering algorithm to identify groups of adjacent voxels with a high discriminatory power. Zhang et al. (2011) exclusively used the approximation coefficients of 3-level decomposition, and used PCA to reduce the features. Ramasamy and Anandhakumar (2011) used fast Fourier transform (FFT) as features. Saritha et al. (2013) proposed a novel feature of wavelet-entropy, and employed spider-web plots to further reduce features. Zhang et al. (2013) employed digital wavelet transform to extract features then used principal component analysis (PCA) to reduce the feature space. Savio and Grana (2013) proposed to use deformation-based morphometry (DBM) techniques, and proposed five features as Jacobian map, modulated GM (MGM), trace of Jacobian matrix (TJM), magnitude of the displacement field, and Geodesic Anisotropy (GEODAN). In addition, they suggested the use of Pearson's correlation (PEC), Bhattacharyya distance (BD), and WTT to measure the significance of voxel site. Das et al. (2013) suggested to use Ripplet transform, followed by PCA to reduce features. Kalbkhani et al. (2013) modeled the detail coefficients of 2-level DWT by generalizing autoregressive conditional heteroscedasticity (GARCH) statistical model, and the parameters of GARCH model were considered as the primary feature vector. Zhang et al. (2014) used an undersampling (US) technique on the volumetric image, followed by singular value decomposition (SVD) to select features. El-Dahshan et al. (2014) proposed to add a preprocessing technique that used pulse-coupled neural network (PCNN) for image segmentation. Zhou et al. (2015) used wavelet-entropy as the feature space. Zhang et al. (2015a) used discrete wavelet packet transform (DWPT), and harnessed Tsallis entropy to obtain features from DWPT coefficients. Yang et al. (2015) selected wavelet-energy as the features.

From the literature used, the DWT based features were proven to be efficient. In this study, we suggested using a novel feature of eigenbrain, which was used for SPECT images but was never been used in MR images.

Classification Model in MRI

There are numerous classification models, but only a few of them are suitable for MR images. Chaplot et al. (2006) employed the self-organizing map (SOM) neural network and SVM. Maitra and Chatterjee (2006) used the common artificial neural network (ANN). El-Dahshan et al. (2010) used ANN and K-nearest neighbor (KNN) classifiers. Plant et al. (2010) used SVM, Bayes statistics, and voting feature intervals (VFI) to derive the quantitative index of pattern matching. Zhang et al. (2011) suggested to use ANN. The weights of ANN were trained by scaled-conjugate-gradient method. Ramasamy and Anandhakumar (2011) proposed to use Expectation and Maximization Gaussian Mixture Model algorithm (EM-GMM). Saritha et al. (2013) used the probabilistic neural network (PNN). Zhang et al. (2013) constructed a kernel SVM with RBF kernel, using particle swarm optimization (PSO) to optimize the parameters C and sigma. Savio and Grana (2013) chose SVM, and used grid search for tuning parameters. Das et al. (2013) used least-square SVM, and their 5 × 5 CV showed high classification accuracy. Kalbkhani et al. (2013) tested the KNN and SVM models. Zhang et al. (2014) proposed to combine KSVM and decision tree, and their method was dubbed KSVM-DT. El-Dahshan et al. (2014) used feed forward back-propagation neural network (FFBPNN). Zhou et al. (2015) used a Naive Bayes classifier (NBC) as classification method. Zhang et al. (2015a) used a generalized eigenvalue proximal SVM (GEPSVM) with RBF kernel. Yang et al. (2015) used SVM as the classifier, and employed biogeography-based optimization (BBO) to train the classifier.

After reviewing the latest literatures that were related to classifiers, we found that SVMs had significant advantages of high accuracy, elegant mathematical tractability, and direct geometric interpretation, compared with other classification methods (Collins and Pape, 2011). In addition, it did not need a large number of training samples to avoid overfitting (Li et al., 2010). Kernel technique further enhanced the performance of SVM. Therefore, KSVM was harnessed in this study.

The Proposed Method

Preprocessing on Volumetric Data

For each individual, all available 3 or 4 volumetric 3D MR brain images were motion-corrected, and coregistered to form an averaged 3D image. Then, those 3D images were spatially normalized to the Talairach coordinate space and brain-masked. CDR was interpreted as the target (label). It is a numeric scale quantifying the severity of symptoms of dementia (Williams et al., 2013). The patient's cognitive and functional performances were assessed in six areas: memory, orientation, judgment and problem solving, community affairs, home and hobbies, and personal care. In this study, we chose two types of CDR, i.e., the subjects with CDR of 0 were considered as NC and subjects with CDR of 1 were considered as AD (Marcus et al., 2007).

Calculating eigenbrains on the entire brain was difficult. Instead, we proposed a simplified method that selected several key slices that capture structures indicative of AD from NC. The procedure was as follows: we established the ICV v as

v(k)=μAD(Slice=k)μNC(Slice=k)2    (1)

where k was the index of key slice, μAD and μNC represented the mean of gray-level values of the kth slice of AD subjects and NC subjects, respectively, ||.||2 represented the l2-norm. Then, we selected the key-slices of ICV larger than 50% of maximum ICV, with 10× undersampling factor (i.e., every 10 slices).

In addition, the slice direction can be chosen as either axial, sagittal, or coronal. Usually coronal direction will give a clearer view than the other two directions. Figure 1 showed that the coronal slice had an advantage over other directions in that it can cover three of the most important tissues within one slice. Those tissues were seen as indicative of AD. These tissues are the cerebral cortex, the ventricle, and the hippocampus. If we used axial or sagittal slice, then we may need to record two or even more slices to cover those tissues. Therefore, we chose the coronal direction for key slice selection, with the aim of recording only one slice.

FIGURE 1
www.frontiersin.org

Figure 1. Difference between (A) a healthy brain and (B) an AD brain. The labeled three regions are (i) cerebral cortex (ii) ventricle, and (iii) hippocampus.

Eigenbrain

AD has different physical structures from NC. Revisit Figure 1 which indicated the AD subjects had severe atrophy of the cerebral cortex (region i), severely enlarged ventricles (region ii), and extreme shrinkage of hippocampus (region iii). Therefore, eigenbrain tried to capture those different characteristic changes of anatomical structures between AD and NC.

Eigenbrain is carried out by PCA, which is a statistical procedure that uses an orthogonal transformation to convert a set of observations of possibly correlated variables into a set of values of linearly uncorrelated variables called principal components (PC). For 2D images the PCs are extended naturally to the 2D eigenbrains.

Suppose X is a given data matrix with size of N × A, where N represents the number of samples and A number of attributes (For a 256 × 256 image, we need to vectorize it to a 1 × 65536 vector, hence A = 65536). First, we normalized the dataset matrix X, so that each sample in the normalized matrix Z was mean-centered and unit-variance scaled, by subtracting its mean value and dividing the difference by its standard deviation.

ZXμ(X)σ(X)    (2)

Next, we estimated the covariance matrix C with size of A × A by

C1N1ZTZ    (3)

Here we used N − 1 instead of N in order to produce an unbiased estimator of the variance (See Bessel's correction (Russell and Cohn, 2012) for details).

Third, we perform the eigendecomposition of C:

C=UU1    (4)

where U is an A × (N − 1) matrix, whose columns are the eigenvectors of covariance matrix C, matrix Λ is an (N − 1) × (N − 1) diagonal matrix whose diagonal elements are eigenvalues of C, each corresponding to an eigenvector of A. It is common to sort the eigenvalue matrix Λ and eigenvector matrix U in order of decreasing eigenvalue u1 > u2 > … > uN. To view the ith eigenbrain u(i), the ith column of U was reshaped to an image. Suppose the ith column of U contains 65536 elements, then the reshaped image was 256 × 256.

u(i)=reshape (U(:,i))    (5)

Note that in our situation (N ~ 102 and A ~ 104, where ~ denotes the same order of magnitude), the computation burdens of eigendecomposition of equation (4) are enormous. It can be accelerated by replacing C in equation (3) with C′, since N << A.

C1N1ZZT    (6)

The size of C′ is N × N, which can significantly reduce the computation burden. Using Matlab, the eigenbrain can be done by a simple “PCA” command without considering these issues. The flowchart of calculating eigenbrain is shown in Figure 2.

FIGURE 2
www.frontiersin.org

Figure 2. Flowchart of calculating eigenbrain.

The eigenvalues represent the distribution of energy of the source data among each of the eigenbrains, where the eigenbrains form a basis for original data.

To further select an eigenbrain that is the most statistically significant, we employ the two-sample location test. Saritha et al. (2013) selected the Student's t-test which assumes both the means and variances of the two data are equal. The assumption of equal variances was not necessary and can be dropped; while the assumption of equal means is essential to select significantly important eigenbrains. Therefore, we used WTT that is an adaption of the Student's t-test and checks nothing except the two populations that have equal means.

The null hypothesis is that the eigenvalues of AD and NC have equal means, without assuming they have equal variances. The alternative hypothesis is they have unequal means. WTT was carried out at the 95% confidence interval. The eigenvalues of the selected most important eigenbrain (MIE) were used as input features for following classification.

Region Detection

We proposed a visual interpretation method of Eigenbrain to detect regions that can distinguish AD and NC, which is not reported in literatures of Alvarez et al. (2009a) and Lopez et al. (2009). The interpretation in a four-stage process is listed in Table 1.

TABLE 1
www.frontiersin.org

Table 1. Four-stage region detection method.

Classifier

SVM was used as the classifier. In addition, sequential minimal optimization (SMO) is chosen to train SVM due to simple and fast speed (Zhang and Wu, 2012b). Traditional linear SVMs cannot separate intricately distributed data. In order to generalize SVMs to create nonlinear hyperplane, the kernel trick is applied. The KSVMs allows us to fit the maximum-margin hyperplane in a transformed feature space (Liu et al., 2014). The transformation may be nonlinear and the transformed space is a higher dimensional space. Though the classifier is a hyperplane in the higher-dimensional feature space, it may be nonlinear in the original input space.

The radial basis function (RBF) kernel is one of the most widely used kernels with the form as Zhang and Wu (2012b).

κ (xm,xn)=exp (xmxn2σ2)    (7)

where κ is the kernel function, σ the scaling factor, and xm and xn are vectors in the input space.

Another commonly used kernel is polynomial (POL) kernel defined as

κ (xm,xn)=(xmTxn+c)d    (8)

where d is the degree of polynomial, and c a soft margin constant trading off the influence of higher-order vs. lower-order terms in the polynomial.

Based on the two kernels, we tested RBF-KSVM and POL-KSVM for our models. To obtain the best parameter of kernels (the scaling factor σ of RBF, or the degree d and soft margin constant c of POL), PSO was employed since it has been used successfully to tune parameters of KSVM in various problems (Aich and Banerjee, 2014; Khazaee and Zadeh, 2014; Xue et al., 2014).

K-fold CV was employed, and K was assigned with a value of 10 considering the best compromise between computational cost and reliable estimates, i.e., the dataset is randomly divided into 10 mutually exclusively subsets of approximately equal size, in which 10 − 1 = 9 subsets were used as training set and the last subset was used as the validation set. The procedure that was mentioned above was repeated 10 times, so each subset was used once for validation. The K results from the K folds were combined together, to yield a single estimation of the whole dataset.

The K-fold CV repeated 50 times, i.e., we carried out a 50 × 10-fold CV. For each time, we used four measures: accuracy, sensitivity, specificity, and precision (Table 2), to assess the performance. Here TP, FP, TN, and FN represented the instance number of true positive, false positive, true negative, and false negative, respectively. We considered a correctly identified AD case as a true positive, following the common convention. Summarizing the 50 repetitions, we reported the final measures of both the mean and standard deviation (SD) of the four measures.

TABLE 2
www.frontiersin.org

Table 2. Assessment of classification performance.

Implementation

The purpose of the proposed method is two-fold: (i) to find discriminant voxels that distinguish AD from NC; and (ii) to develop a CAD system and report its performance. The pseudocode is listed in Table 3.

TABLE 3
www.frontiersin.org

Table 3. Pseudocode of proposed method.

Experiments and Results

The programs were in-house developed using Matlab 2014a, and ran on IBM laptop with 3 GHz Intel i3 dual-processor and 8 GB RAM. Readers could repeat our results on any machine where Matlab is available.

Data Source

We downloaded the dataset from Open Access Series of Imaging Studies (OASIS) (Ardekani et al., 2013, 2014). We chose the cross-sectional dataset corresponding to MRI scans of individuals at a single time point (Bin Tufail et al., 2012). The OASIS dataset consists of 416 subjects aged 18–96, who are all right-handed. We excluded subjects under 60 years old and those with missing records and then picked 126 subjects (98 NCs and 28 ADs) from the rest of the subjects. The demographic statuses of the included subjects were summarized in Table 4. Here SES, CDR, and MMSE represent socioeconomic status, clinical dementia rating, and mini-mental state examination, respectively.

TABLE 4
www.frontiersin.org

Table 4. Subject demographics status.

Preprocessing

Table 5 shows an example of the combination of 3 individual scans of a subject. The resolution is 1 × 1 × 1.25 mm. The preprocessing performed motion-correction on the 3D MR images, registered them to form a combined image in the native acquisition space, and resampled to 1 × 1 × 1 mm. Afterwards, the combined image was spatially normalized to the Talairach coordinate space, and brain-extracted (Table 5).

TABLE 5
www.frontiersin.org

Table 5. Preprocessing of a specified subject.

Key-slice Selection by ICV

The curve of ICV against slice index was shown in Figure 3A. We selected 10 coronal slices (60, 70, 80, 90, 100, 110, 120, 130, 140, and 150). Their corresponding ICVs were all higher than 50% of the maximum. Figures 3B,C showed the axial and sagittal view of the 10 key-slices. Table 6 showed the comparison between NC and AD in the selected 10 key-slices.

FIGURE 3
www.frontiersin.org

Figure 3. Key-Slice selection (The red lines correspond to key-slices). (A) The curve of ICV against coronal slice index. (B) axial view of key-slices. (C) sagittal view of key-slices.

TABLE 6
www.frontiersin.org

Table 6. Difference between NC and AD on key-slices.

Eigenbrains

Table 7 showed the eigenbrain results obtained by running PCA on the slices of all subjects. For each slice, we had a set of 125 eigenbrains in total. Due to the page limit, we selected and listed the first 6 eigenbrains. The eigenbrains were sorted in the order of decreasing eigenvalues. In general, the eigenbrains in the previous columns were more important than in latter columns.

TABLE 7
www.frontiersin.org

Table 7. Eigenbrain results.

Most Important Eigenbrain

WTT was conducted to give quantified proof of why the first eigenbrain was MIE. We performed WTT for the first six eigenbrains of all key-slices between eigenvalues to characterize those that were AD and those that were NC. The results were shown in Table 8, and p-values less than 0.05 were marked in bold. Only the first eigenvalues of all slices were less than 0.05; therefore, the first eigenbrain was indeed the MIE, and we assigned the eigenvalues of MIE of all 10 key-slices (namely, 10 × 1 = 10 features) of each subject to classification.

TABLE 8
www.frontiersin.org

Table 8. WTT of the first six eigenvalues of 10 key-slices.

Classification Comparison

The two classes in order were AD and NC, following common convention. Here we designed three tasks. The first did not use the kernel technique, i.e., the basic linear SVM; the second used RBF-KSVM; and the third used POL-KSVM. The kernel parameters and error penalty were optimized by PSO method. The classification results were listed in Table 9, in addition with the results of state-of-the-art methods.

TABLE 9
www.frontiersin.org

Table 9. Comparison of classification results.

Region Detection

We carried out the region detection procedure from MIE as Section Region Detection described. Table 10 showed the result, in which the green points represented the discriminant voxels.

TABLE 10
www.frontiersin.org

Table 10. Discriminant voxels.

Here we reported the discriminative regions interpreted by eigenbrain in Table 11, where BA represented Brodmann area.

TABLE 11
www.frontiersin.org

Table 11. Regions found by Eigenbrain.

Discussion

It is clearly observed in Table 6 that the selected coronal slices are significant in detecting AD from NC. In particular, the AD subjects show the cerebrospinal fluid (CSF) in the areas occupied by brain matter in the NC subjects. We conclude that 10× is reasonable because of following three reasons: (1) The 10× key-slice undersampling (i.e., select only one slice from 10 consecutive slices) yields a coarser brain while still capturing most tissues in the brain (Compare Table 6 with Figure 1). (2) It is very hard to define a fitness (optimization) function to find the optimal undersampling factor. (3) The classification system has a good accuracy in distinguishing AD from NC, and it detects correct AD-related brain regions (See Tables 9, 11). As there are spatial redundancy for neighboring coronal slices, the undersampling could reduce this redundancy to a rather small degree.

Overall, the eigenbrains in Table 7 capture both similarities and differences of structural features between AD and NC. The first eigenbrain capture the significant feature of AD from NC, and the second and following eigenbrains capture general brain structure. Revisiting the hippocampus part in the first eigenbrain of all key-slices, it is easily perceived that the body lateral ventricles area of AD are highlighted, which is indeed a distinct attribute between AD and NC. Our experiment extends the eigenbrain on SPECT images by Alvarez et al. (2009a) and Lopez et al. (2009) and shows that eigenbrain is also suitable for MRI scans.

The p-values in Table 8 show that the first eigenvalue λ1 are all less than 0.05 for all key-slices. It indicates that mean values of λ1 of AD and NC are significantly different. Hence, the most dominating eigenvalue characterizing AD and NC is the one corresponding to the first eigenbrain. For other eigenvalues, merely 1 of 10 p-values is less than 0.05, which indicates that those eigenbrains are not dominating features indicative of AD from NC. Therefore, the first eigenvalue is MIE and was selected.

Classification results in Table 9 compare the proposed three classifiers with state-of-the-art methods, in which Zhang's results (Table 7 in Zhang et al., 2014) are calculated through a single K-fold CV experiment. Plant's results (Task 1 in Table 3 Plant et al., 2010) offer the means together with 95% confidence intervals. Savio's results (Table 5 Savio and Grana, 2013) give the means with SD. For the proposed methods, it is unexpected that the POL-KSVM produces better classification accuracy of 92.36 ± 0.94 than linear SVM of 91.47 ± 1.02 and RBF-KSVM of 86.71 ± 1.93, because RBF was reported as the most widely used kernel. Our results are better than or comparable to other approaches to AD prediction from MR brain images of NC, e.g., US + SVD-PCA + SVM-DT of 90% (Zhang et al., 2014), BRC + IG + SVM of 90% (Plant et al., 2010), BRC + IG + Bayes of 92% (Plant et al., 2010), MGM + PEC + SVM of 92.07% (Savio and Grana, 2013), GEODAN + BD + SVM of 92.09% (Savio and Grana, 2013), and TJM + WTT + SVM of 92.83% (Savio and Grana, 2013). There were many other methods (Gray et al., 2012; Arbizu et al., 2013; Chaves et al., 2013; Dukart et al., 2013; Cohen and Klunk, 2014) proposed for detecting AD from NC, however, they treated images from other modalities (such as SPECT and PET). Therefore, it is not appropriate to compare the proposed methods with them. We will test our methods on SPECT and PET images in the future.

Table 11 shows that eigenbrains interpret the discriminant voxels involving the following regions reported in existing literatures: Anterior Cingulate (BA-24, BA-32) (Schultz et al., 2014), Caudate Nucleus (Head, body, and tail) (Möller et al., 2015), Cerebellum (Colloby et al., 2014), Cingulate Gyrus (BA-23, BA-24, BA-31) (Yu et al., 2014), Claustrum (De Reuck et al., 2014), Inferior Frontal Gyrus (BA-47) (Eliasova et al., 2014), Inferior Parietal Lobule (BA-40) (Wang et al., 2015), Insula (BA-13) (He et al., 2015), Lateral Ventricle (Voevodskaya et al., 2014), Lentiform Nucleus (Möller et al., 2015), Lingual gyrus (Lehmann et al., 2013), Medial Frontal Gyrus (BA-10, BA-11, BA-25, BA-6) (Kang et al., 2013), Middle Frontal Gyrus (BA-11) (Schultz et al., 2014), Middle Occipital Gyrus (Lehmann et al., 2013), Middle Temporal Gyrus (Aubry et al., 2015), Paracentral Lobule (BA-3, BA-4, BA-5, BA-6, BA-7) (Kang et al., 2013), Parahippocampal Gyrus (Amygdala, BA-28, BA-35, Hippocampus) (Eskildsen et al., 2015), Postcentral Gyrus (BA-5) (Kang et al., 2013), Posterior Cingulate (Shinohara et al., 2014), Precentral Gyrus (BA-4) (Kang et al., 2013), Precuneus (BA-7, BA-31) (Kang et al., 2013), Subcallosal Gyrus (BA-25, BA-34, BA-47) (Paakki et al., 2010), Sub-Gyral (BA-40, Corpus Callosum, Hippocampus) (Streitburger et al., 2012), Superior Frontal Gyrus (Chen et al., 2014), Superior Parietal Lobule (Quiroz et al., 2013), Superior Temporal Gyrus (BA-38) (Paakki et al., 2010), Supramarginal Gyrus (Quiroz et al., 2013), Thalamus (Medial Geniculum Body, Pulvinar, Ventral Lateral Nucleus) (He et al., 2015), Transverse Temporal Gyrus (BA-41) (Kim et al., 2012), and Uncus (BA-28) (Bangen et al., 2014).

Nevertheless, some regions reported to be associated with AD are not interpreted by Eigenbrain, such as subthalamic nucleus (De Reuck et al., 2014). The reason may lie in three aspects. First, the quantile of our method is assigned with a value of 0.98, which is considered high. Reducing the quantile value may include more regions. Second, some literature used other advanced imaging modalities, such as MRSI and fMRI for metabolism detection and function analysis. Third, the key-slice selection procedure may miss important regions.

From another point of view, Table 11 demonstrates the power of the eigenbrain. Our study uses only one feature (eigenbrain) on 10 key-slices of a simple 3D structural MR image, nevertheless, our findings cover 30 related regions reported by over twenty literatures, which used various feature extraction methods and advanced imaging technologies.

The contributions of the paper fall within the following five aspects: (i) We generalize the Eigenbrain to MR images, and prove its effectiveness; (ii) We propose a hybrid eigenbrain-based CAD system that can not only detect AD from NC, but also detect brain regions that related to AD. (iii) We prove the proposed method has a classification accuracy comparable to state-of-the-art methods, and the detected brain regions are in line with 16 existing literatures. (iv) We use ICV and WTT to reduce redundant data; (v) we find POL kernel is better than linear and RBF kernel for this study.

In conclusion, the advantages of eigenbrain are three-fold: (i) it reaches very high classification accuracy, which was better than or competitive with state-of-the-art methods (Plant et al., 2010; Savio and Grana, 2013; Zhang et al., 2014); (ii) it can directly find discriminant voxels/regions within the whole brain; (iii) it can be combined with other features, in order to increase the classification performance. On the other hand, the disadvantages of eigenbrain also exist: (i) it is essentially two-dimensional, which does not reduce the redundancy along the slice direction; (ii) it needs preprocessing of spatial registration, which costs large amount of computation resources.

To the policy-makers, this study suggests the eigenbrain technique can achieve comparable results to traditional methods. It may offer a ray of hope for AD diagnosis with unconventional means with the combination of eigenbrain and machine learning. This preclinical study suggests that hospitals and medical laboratories enroll more computer scientists and engineers, with the aim of developing efficient AD diagnosis and region detection systems.

Conclusion and Future Research

We presented an automated and accurate classification method that was based on eigenbrains and machine learning, in order to detect AD subjects and AD-related brain regions using 3D MR images. The results showed the proposed POL-KSVM method achieved 92.36% accuracy, which was competitive with state-of-the-art methods.

In the future, we will focus our research in the following aspects: (i) We shall generalize the eigenbrain to three dimensional, so the procedure of key-slice selection can be removed; (ii) We shall test other kernels for SVM, and try to replace KSVM with other advanced pattern recognition tools. (iii) Eigenbrain can be used in combination with DWT-based features and others, and an increase in classification accuracy is expected.

Conflict of Interest Statement

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.

Acknowledgments

This work was supported by NSFC (610011024, 61273243, 51407095), Program of Natural Science Research of Jiangsu Higher Education Institutions (13KJB460011, 14KJB520021), Jiangsu Key Laboratory of 3D Printing Equipment and Manufacturing (BM2013006), Key Supporting Science and Technology Program (Industry) of Jiangsu Province (BE2012201, BE2014009-3, BE2013012-2), Special Funds for Scientific and Technological Achievement Transformation Project in Jiangsu Province (BA2013058), and Nanjing Normal University Research Foundation for Talented Scholars (2013119XGQ0061, 2014119XGQ0080). The authors acknowledge their gratitude to the OASIS dataset that came from NIH grants P50AG05681, P01 AG03991, R01 AG021910, P50 MH071616, U24 RR021382 and R01 MH56584.

Supplementary Material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fncom.2015.00066/abstract

Footnotes

1. ^Some abbreviations are modified to avoid conflict within this paper.

References

Aich, U., and Banerjee, S. (2014). Modeling of EDM responses by support vector machine regression with parameters selected by particle swarm optimization. Appl. Math. Model. 38, 2800–2818. doi: 10.1016/j.apm.2013.10.073

CrossRef Full Text | Google Scholar

Alvarez, I., Gorriz, J. M., Ramirez, J., Salas-Gonzalez, D., Lopez, M., Puntonet, C. G., et al. (2009a). Alzheimer's diagnosis using eigenbrains and support vector machines. Electron. Lett. 45, 342–343. doi: 10.1049/el.2009.3415

CrossRef Full Text | Google Scholar

Álvarez, I., Górriz, J. M., Ramírez, J., Salas-Gonzalez, D., López, M., Segovia, F., et al. (2009b). “Alzheimer's diagnosis using eigenbrains and support vector machines,” in Bio-Inspired Systems: Computational and Ambient Intelligence, Vol. 5517, eds J. Cabestany, F. Sandoval, A. Prieto, and J. Corchado (Berlin: Springer), 973–980.

Google Scholar

Anagnostopoulos, C. N., Giannoukos, I., Spenger, C., Simmons, A., Mecocci, P., Soininen, H., et al. (2013). “Classification models for Alzheimer's disease Detection,” in Engineering Applications of Neural Networks, Vol. 384(Pt II), eds L. Iliadis, H. Papadopoulos, and C. Jayne (Berlin; Heidelberg: Springer), 193–202. doi: 10.1007/978-3-642-41016-1_21

CrossRef Full Text | Google Scholar

Angelini, E. D., Delon, J., Bah, A. B., Capelle, L., and Mandonnet, E. (2012). Differential MRI analysis for quantification of low grade glioma growth. Med. Image Anal. 16, 114–126. doi: 10.1016/j.media.2011.05.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Arbizu, J., Prieto, E., Martinez-Lage, P., Marti-Climent, J. M., Garcia-Granero, M., Lamet, I., et al. (2013). Automated analysis of FDG PET as a tool for single-subject probabilistic prediction and detection of Alzheimer's disease dementia. Eur. J. Nucl. Med. Mol. Imaging 40, 1394–1405. doi: 10.1007/s00259-013-2458-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Ardekani, B. A., Bachman, A. H., Figarsky, K., and Sidtis, J. J. (2014). Corpus callosum shape changes in early Alzheimer's disease: an MRI study using the OASIS brain database. Brain Struct. Funct. 219, 343–352. doi: 10.1007/s00429-013-0503-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Ardekani, B. A., Figarsky, K., and Sidtis, J. J. (2013). Sexual dimorphism in the human corpus callosum: an MRI study using the OASIS brain database. Cereb. Cortex 23, 2514–2520. doi: 10.1093/cercor/bhs253

PubMed Abstract | CrossRef Full Text | Google Scholar

Aubry, S., Shin, W., Crary, J. F., Lefort, R., Qureshi, Y. H., Lefebvre, C., et al. (2015). Assembly and interrogation of Alzheimer's disease genetic networks reveal novel regulators of progression. PLoS ONE 10:25. doi: 10.1371/journal.pone.0120352

CrossRef Full Text | Google Scholar

Bangen, K. J., Nation, D. A., Clark, L. R., Harmell, A. L., Wierenga, C. E., Dev, S. I., et al. (2014). Interactive effects of vascular risk burden and advanced age on cerebral blood flow. Front. Aging Neurosci. 6:159. doi: 10.3389/fnagi.2014.00159

PubMed Abstract | CrossRef Full Text | Google Scholar

Bin Tufail, A., Abidi, A., Siddiqui, A. M., and Younis, M. S. (2012). “Multiclass classification of initial stages of Alzheimer's disease using structural MRI phase images,” in Proceedings of the IEEE International Conference in Control System, Computing and Engineering (ICCSCE) (Penang: IEEE), 317–321. doi: 10.1109/ICCSCE.2012.6487163

CrossRef Full Text | Google Scholar

Brookmeyer, R., Johnson, E., Ziegler-Graham, K., and Arrighi, H. M. (2007). Forecasting the global burden of Alzheimer's disease. Alzheimers Dement. 3, 186–191. doi: 10.1016/j.jalz.2007.04.381

PubMed Abstract | CrossRef Full Text | Google Scholar

Chaplot, S., Patnaik, L. M., and Jagannathan, N. R. (2006). Classification of magnetic resonance brain images using wavelets as input to support vector machine and neural network. Biomed. Signal Process. Control 1, 86–92. doi: 10.1016/j.bspc.2006.05.002

CrossRef Full Text | Google Scholar

Chaves, R., Ramirez, J., Gorriz, J. M., and Alzheimer's Dis, N. (2013). Integrating discretization and association rule-based classification for Alzheimer's disease diagnosis. Expert Syst. Appl. 40, 1571–1578. doi: 10.1016/j.eswa.2012.09.003

CrossRef Full Text | Google Scholar

Chen, Y., Liu, Z., Zhang, J., Xu, K., Zhang, S., Wei, D., et al. (2014). Altered brain activation patterns under different working memory loads in patients with Type 2 diabetes. Diabetes Care 37, 3157–3163. doi: 10.2337/dc14-1683

PubMed Abstract | CrossRef Full Text | Google Scholar

Cohen, A. D., and Klunk, W. E. (2014). Early detection of Alzheimer's disease using PiB and FDG PET. Neurobiol. Dis. 72, 117–122. doi: 10.1016/j.nbd.2014.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Collins, M. P., and Pape, S. E. (2011). The potential of support vector machine as the diagnostic tool for schizophrenia: a systematic literature review of neuroimaging studies. Eur. Psychiatry 26, 117–122. doi: 10.1016/S0924-9338(11)73068-1

CrossRef Full Text | Google Scholar

Colloby, S. J., O'Brien, J. T., and Taylor, J. P. (2014). Patterns of cerebellar volume loss in dementia with Lewy bodies and Alzheimer's disease: A VBM-DARTEL study. Psychiatry Res. 223, 187–191. doi: 10.1016/j.pscychresns.2014.06.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Das, S., Chowdhury, M., and Kundu, M. K. (2013). Brain MR image classification using multiscale geometric analysis of ripplet. Prog. Electromagn. Res. 137, 1–17. doi: 10.2528/PIER13010105

CrossRef Full Text | Google Scholar

De Reuck, J. L., Deramecourt, V., Auger, F., Durieux, N., Cordonnier, C., Devos, D., et al. (2014). Iron deposits in post-mortem brains of patients with neurodegenerative and cerebrovascular diseases: a semi-quantitative 7.0 T magnetic resonance imaging study. Eur. J. Neurol. 21, 1026–1031. doi: 10.1111/ene.12432

PubMed Abstract | CrossRef Full Text | Google Scholar

Dukart, J., Mueller, K., Barthel, H., Villringer, A., Sabri, O., Schroeter, M. L., et al. (2013). Meta-analysis based SVM classification enables accurate detection of Alzheimer's disease across different clinical centers using FDG-PET and MRI. Psychiatry Res. 212, 230–236. doi: 10.1016/j.pscychresns.2012.04.007

PubMed Abstract | CrossRef Full Text | Google Scholar

El-Dahshan, E. S. A., Hosny, T., and Salem, A. B. M. (2010). Hybrid intelligent techniques for MRI brain images classification. Digit. Signal Process. 20, 433–441. doi: 10.1016/j.dsp.2009.07.002

CrossRef Full Text | Google Scholar

El-Dahshan, E. S. A., Mohsen, H. M., Revett, K., and Salem, A. B. M. (2014). Computer-aided diagnosis of human brain tumor through MRI: a survey and a new algorithm. Expert Syst. Appl. 41, 5526–5545. doi: 10.1016/j.eswa.2014.01.021

CrossRef Full Text | Google Scholar

Eliasova, I., Anderkova, L., Marecek, R., and Rektorova, I. (2014). Non-invasive brain stimulation of the right inferior frontal gyrus may improve attention in early Alzheimer's disease: a pilot study. J. Neurol. Sci. 346, 318–322. doi: 10.1016/j.jns.2014.08.036

PubMed Abstract | CrossRef Full Text | Google Scholar

Eskildsen, S. F., Coupé, P., Fonov, V. S., Pruessner, J. C., and Collins, D. L. (2015). Structural imaging biomarkers of Alzheimer's disease: predicting disease progression. Neurobiol. Aging 36(Suppl. 1), S23–S31. doi: 10.1016/j.neurobiolaging.2014.04.034

PubMed Abstract | CrossRef Full Text | Google Scholar

Garcia, S., Fernandez, A., Luengo, J., and Herrera, F. (2010). Advanced nonparametric tests for multiple comparisons in the design of experiments in computational intelligence and data mining: experimental analysis of power. Inf. Sci. 180, 2044–2064. doi: 10.1016/j.ins.2009.12.010

CrossRef Full Text | Google Scholar

Goh, S., Dong, Z., Zhang, Y., DiMauro, S., and Peterson, B. S. (2014). Mitochondrial dysfunction as a neurobiological subtype of autism spectrum disorder: evidence from brain imaging. JAMA Psychiatry 71, 665–671. doi: 10.1001/jamapsychiatry.2014.179

PubMed Abstract | CrossRef Full Text | Google Scholar

Gomes, T. A. F., Prudêncio, R. B. C., Soares, C., Rossi, A. L. D., and Carvalho, A. (2012). Combining meta-learning and search techniques to select parameters for support vector machines. Neurocomputing 75, 3–13. doi: 10.1016/j.neucom.2011.07.005

CrossRef Full Text | Google Scholar

Gray, K. R., Wolz, R., Heckemann, R. A., Aljabar, P., Hammers, A., Rueckert, D., et al. (2012). Multi-region analysis of longitudinal FDG-PET for the classification of Alzheimer's disease. Neuroimage 60, 221–229. doi: 10.1016/j.neuroimage.2011.12.071

PubMed Abstract | CrossRef Full Text | Google Scholar

Hable, R. (2012). Asymptotic normality of support vector machine variants and other regularized kernel methods. J. Multivar. Anal. 106, 92–117. doi: 10.1016/j.jmva.2011.11.004

CrossRef Full Text | Google Scholar

Hahn, K., Myers, N., Prigarin, S., Rodenacker, K., Kurz, A., Förstl, H., et al. (2013). Selectively and progressively disrupted structural connectivity of functional brain networks in Alzheimer's disease—Revealed by a novel framework to analyze edge distributions of networks detecting disruptions with strong statistical evidence. Neuroimage 81, 96–109. doi: 10.1016/j.neuroimage.2013.05.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Hamy, V., Dikaios, N., Punwani, S., Melbourne, A., Latifoltojar, A., Makanyanga, J., et al. (2014). Respiratory motion correction in dynamic MRI using robust data decomposition registration – Application to DCE-MRI. Med. Image Anal. 18, 301–313. doi: 10.1016/j.media.2013.10.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, J. W., Kim, T. H., Lee, S. B., Park, J. H., Lee, J. J., Huh, Y., et al. (2011). 327 Diagnostic Stability of Mild Cognitive Impairment Subtype. Asian J. Psychiatry 4(Suppl. 1), S65–S66. doi: 10.1016/s1876-2018(11)60250-5

CrossRef Full Text

He, W., Liu, D., Radua, J., Li, G., Han, B., and Sun, Z. (2015). Meta-analytic comparison between PIB-PET and FDG-PET results in Alzheimer's disease and MCI. Cell Biochem. Biophys. 71, 17–26. doi: 10.1007/s12013-014-0138-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Jeurissen, B., Leemans, A., and Sijbers, J. (2014). Automated correction of improperly rotated diffusion gradient orientations in diffusion weighted MRI. Med. Image Anal. 18, 953–962. doi: 10.1016/j.media.2014.05.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Kalbkhani, H., Shayesteh, M. G., and Zali-Vargahan, B. (2013). Robust algorithm for brain magnetic resonance image (MRI) classification based on GARCH variances series. Biomed. Signal Process. Control 8, 909–919. doi: 10.1016/j.bspc.2013.09.001

CrossRef Full Text | Google Scholar

Kang, K., Yoon, U., Lee, J. M., and Lee, H. W. (2013). Idiopathic normal-pressure hydrocephalus, cortical thinning, and the cerebrospinal fluid tap test. J. Neurol. Sci. 334, 55–62. doi: 10.1016/j.jns.2013.07.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Khazaee, A., and Zadeh, A. E. (2014). ECG beat classification using particle swarm optimization and support vector machine. Front. Comput. Sci. 8, 217–231. doi: 10.1007/s11704-014-2398-1

CrossRef Full Text | Google Scholar

Kim, J. S., Lee, S. H., Park, G., Kim, S., Bae, S. M., Kim, D. W., et al. (2012). Clinical implications of quantitative electroencephalography and current source density in patients with Alzheimer's disease. Brain Topogr. 25, 461–474. doi: 10.1007/s10548-012-0234-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Kubota, T., Ushijima, Y., and Nishimura, T. (2006). A region-of-interest (ROI) template for three-dimensional stereotactic surface projection (3D-SSP) images: initial application to analysis of Alzheimer disease and mild cognitive impairment. Int. Congr. Ser. 1290, 128–134. doi: 10.1016/j.ics.2005.11.104

CrossRef Full Text | Google Scholar

Lee, W., Park, B., and Han, K. (2013). Classification of diffusion tensor images for the early detection of Alzheimer's disease. Comput. Biol. Med. 43, 1313–1320. doi: 10.1016/j.compbiomed.2013.07.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Lehmann, M., Ghosh, P. M., Madison, C., Laforce, R., Corbetta-Rastelli, C., Weiner, M. W., et al. (2013). Diverging patterns of amyloid deposition and hypometabolism in clinical variants of probable Alzheimer's disease. Brain 136, 844–858. doi: 10.1093/brain/aws327

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, D., Yang, W., and Wang, S. (2010). Classification of foreign fibers in cotton lint using machine vision and multi-class support vector machine. Comput. Electron. Agric. 74, 274–279. doi: 10.1016/j.compag.2010.09.002

CrossRef Full Text | Google Scholar

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 Inform. 18, 984–990. doi: 10.1109/JBHI.2013.2285378

PubMed Abstract | CrossRef Full Text | Google Scholar

Lopez, M., Ramirez, J., Gorriz, J. M., Alvarez, I., Salas-Gonzalez, D., Segovia, F., et al. (2009). “Automatic system for Alzheimer's disease diagnosis using eigenbrains and bayesian classification rules,” Bio-Inspired Systems: Computational and Ambient Intelligence, Vol. 5517, eds J. Cabestany, A. Prieto, F. Sandoval, and J. M. Corchado (Berlin: Springer-Verlag Berlin), 949–956.

Google Scholar

Maitra, M., and Chatterjee, A. (2006). A Slantlet transform based intelligent system for magnetic resonance brain image classification. Biomed. Signal Process. Control 1, 299–306. doi: 10.1016/j.bspc.2006.12.001

CrossRef Full Text | Google Scholar

Marcus, D. S., Wang, T. H., Parker, J., Csernansky, J. G., Morris, J. C., and Buckner, R. L. (2007). Open Access Series of Imaging Studies (OASIS): cross-sectional MRI data in young, middle aged, nondemented, and demented older adults. J. Cogn. Neurosci. 19, 1498–1507. doi: 10.1162/jocn.2007.19.9.1498

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, V., Erlien, S., and Piersol, J. (2012). Identifying Dementia in MRI Scans using Machine Learning. Stanford, CA: Stanford University.

Google Scholar

Möller, C., Dieleman, N., van der Flier, W. M., Versteeg, A., Pijnenburg, Y., Scheltens, P., et al. (2015). More atrophy of deep gray matter structures in frontotemporal dementia compared to Alzheimer's disease. J. Alzheimers Dis. 44, 635–647. doi: 10.3233/JAD-141230

PubMed Abstract | CrossRef Full Text

Nambakhsh, C. M. S., Yuan, J., Punithakumar, K., Goela, A., Rajchl, M., Peters, T. M., et al. (2013). Left ventricle segmentation in MRI via convex relaxed distribution matching. Med. Image Anal. 17, 1010–1024. doi: 10.1016/j.media.2013.05.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Paakki, J. J., Rahko, J., Long, X., Moilanen, I., Tervonen, O., Nikkinen, J., et al. (2010). Alterations in regional homogeneity of resting-state brain activity in autism spectrum disorders. Brain Res. 1321, 169–179. doi: 10.1016/j.brainres.2009.12.081

PubMed Abstract | CrossRef Full Text | Google Scholar

Pennanen, C., Kivipelto, M., Tuomainen, S., Hartikainen, P., Hänninen, T., Laakso, M. P., et al. (2004). Hippocampus and entorhinal cortex in mild cognitive impairment and early AD. Neurobiol. Aging 25, 303–310. doi: 10.1016/S0197-4580(03)00084-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Plant, C., Teipel, S. J., Oswald, A., Böhm, C., Meindl, T., Mourao-Miranda, J., et al. (2010). Automated detection of brain atrophy patterns based on MRI for the prediction of Alzheimer's disease. Neuroimage 50, 162–174. doi: 10.1016/j.neuroimage.2009.11.046

PubMed Abstract | CrossRef Full Text | Google Scholar

Quiroz, Y. T., Stern, C. E., Reiman, E. M., Brickhouse, M., Ruiz, A., Sperling, R. A., et al. (2013). Cortical atrophy in presymptomatic Alzheimer's disease presenilin 1 mutation carriers. J. Neurol. Neurosurg. Psychiatry 84, 556–561. doi: 10.1136/jnnp-2012-303299

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramasamy, R., and Anandhakumar, P. (2011). “Brain tissue classification of MR images using fast fourier transform based expectation- maximization gaussian mixture model,” in Advances in Computing and Information Technology, Vol. 198, D. C. Wyld, M. Wozniak, N. Chaki, N. Meghanathan, and D. Nagamalai (Berlin, Springer-Verlag Berlin), 387–398.

Google Scholar

Russell, J., and Cohn, R. (2012). Bessel's Correction.

Saritha, M., Joseph, K. P., and Mathew, A. T. (2013). Classification of MRI brain images using combined wavelet entropy based spider web plots and probabilistic neural network. Pattern Recognit. Lett. 34, 2151–2156. doi: 10.1016/j.patrec.2013.08.017

CrossRef Full Text | Google Scholar

Savio, A., and Grana, M. (2013). Deformation based feature selection for computer aided diagnosis of Alzheimer's Disease. Expert Syst. Appl. 40, 1619–1628. doi: 10.1016/j.eswa.2012.09.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Schultz, S. A., Larson, J., Oh, J., Koscik, R., Dowling, M. N., Gallagher, C. L., et al. (2014). Participation in cognitively-stimulating activities is associated with brain structure and cognitive function in preclinical Alzheimer's disease. Brain Imaging Behav. doi: 10.1007/s11682-014-9329-5. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Shamonin, D. P., Bron, E. E., Lelieveldt, B. P. F., Smits, M., Klein, S., and Staring, M. (2014). Fast Parallel Image Registration on CPU and GPU for Diagnostic Classification of Alzheimer's Disease. Front. Neuroinform. 7:50. doi: 10.3389/fninf.2013.00050

PubMed Abstract | CrossRef Full Text | Google Scholar

Shinohara, M., Fujioka, S., Murray, M. E., Wojtas, A., Baker, M., Rovelet-Lecrux, A., et al. (2014). Regional distribution of synaptic markers and APP correlate with distinct clinicopathological features in sporadic and familial Alzheimer's disease. Brain 137, 1533–1549. doi: 10.1093/brain/awu046

PubMed Abstract | CrossRef Full Text | Google Scholar

Smal, I., Carranza-Herrezuelo, N., Klein, S., Wielopolski, P., Moelker, A., Springeling, T., et al. (2012). Reversible jump MCMC methods for fully automatic motion analysis in tagged MRI. Med. Image Anal. 16, 301–324. doi: 10.1016/j.media.2011.08.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Streitburger, D. P., Möller, H. E., Tittgemeyer, M., Hund-Georgiadis, M., Schroeter, M. L., and Mueller, K. (2012). Investigating structural brain changes of dehydration using voxel-based morphometry. PLoS ONE 7:e44195. doi: 10.1371/journal.pone.0044195

PubMed Abstract | CrossRef Full Text | Google Scholar

Voevodskaya, O., Simmons, A., Nordenskjold, R., Kullberg, J., Ahlstrom, H., Lind, L., et al. (2014). The effects of intracranial volume adjustment approaches on multiple regional MRI volumes in healthy aging and Alzheimer's disease. Front. Aging Neurosci. 6:264. doi: 10.3389/fnagi.2014.00264

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Z., Xia, M., Dai, Z., Liang, X., Song, H., He, Y., et al. (2015). Differentially disrupted functional connectivity of the subregions of the inferior parietal lobule in Alzheimer's disease. Brain Struct. Funct. 220, 745–762. doi: 10.1007/s00429-013-0681-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Williams, M. M., Storandt, M., Roe, C. M., and Morris, J. C. (2013). Progression of Alzheimer's disease as measured by clinical dementia rating sum of boxes scores. Alzheimers Dement. 9(1, Suppl.), S39-S44. doi: 10.1016/j.jalz.2012.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Xinyun, C., Wenlu, Y., and Xudong, H. (2011). “ICA-based classification of MCI vs HC. Natural Computation (ICNC),” Seventh International Conference, Vol. 3 (Shanghai: IEEE), 1658–1662. doi: 10.1109/ICNC.2011.6022275

CrossRef Full Text | Google Scholar

Xue, Z. H., Du, P. J., and Su, H. J. (2014). Harmonic analysis for hyperspectral image classification integrated with PSO optimized SVM. J. Select. Topics Appl. Earth Obs. Remote Sens IEEE 7, 2131–2146. doi: 10.1109/JSTARS.2014.2307091

CrossRef Full Text

Yang, G., Zhang, Y., Yang, J., Ji, G., Dong, Z., Wang, S., et al. (2015). Automated classification of brain images using wavelet-energy and biogeography-based optimization. Multimed. Tools Appl. 1–17. doi: 10.1007/s11042-015-2649-7

CrossRef Full Text | Google Scholar

Yu, Q., Peng, Y., Mishra, V., Ouyang, A., Li, H., Zhang, H., et al. (2014). Microstructure, length, and connection of limbic tracts in normal human brain development. Front. Aging Neurosci. 6:228. doi: 10.3389/fnagi.2014.00228

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Dong, Z., Wang, S., Ji, G., and Yang, J. (2015a). Preclinical Diagnosis of Magnetic Resonance (MR) Brain Images via Discrete Wavelet Packet Transform with Tsallis Entropy and Generalized Eigenvalue Proximal Support Vector Machine (GEPSVM). Entropy 17, 1795–1813. doi: 10.3390/e170a41795

CrossRef Full Text | Google Scholar

Zhang, Y., Dong, Z., Wu, L., and Wang, S. (2011). A hybrid method for MRI brain image classification. Expert Syst. Appl. 38, 10049–10053. doi: 10.1016/j.eswa.2011.02.012

CrossRef Full Text | Google Scholar

Zhang, Y., Wang, S., and Dong, Z. (2014). Classification of Alzheimer disease based on structural magnetic resonance imaging by kernel support vector machine decision tree. Prog. Electromagn. Res. 144, 171–184. doi: 10.2528/PIER13121310

CrossRef Full Text | Google Scholar

Zhang, Y., Wang, S., Ji, G., and Dong, Z. (2013). An MR brain images classifier system via particle swarm optimization and kernel support vector machine. Scientific World Journal 2013:130134. doi: 10.1155/2013/130134

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Wang, S., Ji, G., and Dong, Z. (2015b). Exponential wavelet iterative shrinkage thresholding algorithm with random shift for compressed sensing magnetic resonance imaging. IEEJ Trans. Electr. Electron. Eng. 10, 116–117. doi: 10.1002/tee.22059

CrossRef Full Text | Google Scholar

Zhang, Y., and Wu, L. (2012a). Classification of fruits using computer vision and a multiclass support vector machine. Sensors 12, 12489–12505. doi: 10.3390/s120912489

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., and Wu, L. (2012b). An MR brain images classifier via principal component analysis and kernel support vector machine. Prog. Electromagn. Res. 130, 369–388. doi: 10.2528/PIER12061410

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, X., Wang, S., Xu, W., Ji, G., Phillips, P., Sun, P., et al. (2015). “Detection of pathological brain in MRI scanning based on wavelet-entropy and naive bayes classifier,” in Bioinformatics and Biomedical Engineering, Vol. 9043, eds F. Ortuño and I. Rojas (Granada: Springer International Publishing), 201–209. doi: 10.1007/978-3-319-16483-0_20

CrossRef Full Text | Google Scholar

Keywords: Alzheimer's disease, Welch's t-test, magnetic resonance imaging, machine learning, machine vision, eigenbrain, support vector machine, particle swarm optimization

Citation: Zhang Y, Dong Z, Phillips P, Wang S, Ji G, Yang J and Yuan T-F (2015) Detection of subjects and brain regions related to Alzheimer's disease using 3D MRI scans based on eigenbrain and machine learning. Front. Comput. Neurosci. 9:66. doi: 10.3389/fncom.2015.00066

Received: 12 February 2015; Accepted: 17 May 2015;
Published: 02 June 2015.

Edited by:

Tobias Alecio Mattei, Brain and Spine Center - InvisionHealth - Kenmore Mercy Hospital, USA

Reviewed by:

Fahad Sultan, University Tübingen, Germany
Petia D. Koprinkova-Hristova, Bulgarian Academy of Sciences, Bulgaria

Copyright © 2015 Zhang, Dong, Phillips, Wang, Ji, Yang and Yuan. 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) or licensor 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: Yudong Zhang, School of Computer Science and Technology, Nanjing Normal University, 1 Wenyuan, Nanjing, Jiangsu 210023, China, zhangyudong@njnu.edu.cn;
Ti-Fei Yuan, School of Psychology, Nanjing Normal University, 22 Ninghai Rd., Nanjing, Jiangsu 210008, China, ytf0707@126.com

Disclaimer: 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.