Functional Brain Network Learning Based on Spatial Similarity for Brain Disorders Identification

Functional brain network (FBN) measures based on functional magnetic resonance imaging (fMRI) data, has become important biomarkers for early diagnosis and prediction of clinical outcomes in neurological diseases, such as Alzheimer’s diseases (AD) and its prodromal state (i.e., Mild cognitive impairment, MCI). In the past decades, researchers have developed numbers of approaches for FBN estimation, including Pearson’s correction (PC), sparse representation (SR), and so on. Despite their popularity and wide applica-tions in current studies, most of the approaches for FBN estimation only consider the dependency between the measured blood oxygen level dependent (BOLD) time series, but ignore the spatial relationships between pairs of brain regions. In practice, the strength of functional connection between brain regions will decrease as their distance increases. Inspired by this, we proposed a new approach for FBN estimation based on the assumption that the closer brain regions tend to share stronger relationships or similarities. To verify the effectiveness of the proposed method, we conduct experiments on a public dataset to identify the patients with MCIs from health controls (HCs) using the estimated FBNs. Experimental results demonstrate that the proposed approach yields statistically significant improvement in seven performance metrics over using the baseline methods.


Introduction
As the common neurological diseases, Alzheimer's disease (AD) is defined as diagnostic criteria of progressive cognitive and memory impairment [1]. It is reported that increasing populations will suffer from AD over time, which will severely affect their health and daily life [2]. Consequently, it is urgently needed for early detection and correct diagnosis for mild cognitive impairment (i.e., MCI, the early prodromal stage of AD) populations [3].
As a widely used non-invasive technique for measuring brain activities, functional magnetic resonance imaging (fMRI) has been successfully applied to explore early MCI diagnosis [4]. In particular, the popular diagnosis models include Bayesian network (BN) [5], support vector machine (SVM) [6], and graph learning (GL) [7], but, due to the randomness and the asynchronization of the spontaneous brain activities, it is difficult to train these models directly using the fMRI data. In contrast, functional brain network (FBN), estimated by fMRI data, can provide an effective tool to better understand and study the brain structure of neurological diseases (i.e., MCI) [8]. In other words, estimating a high-quality FBN may be helpful for the subsequent neurological identification task [9] [10].
Mathematically, each node in the FBN corresponds to the brain region of interest (ROI), and each edge shows the dependency relationships between pairs of ROIs. To estimate these edges, researchers have proposed numbers of approaches, in which the standard correction-based approaches include Pearson's correction (PC)-based approach [11], (regularized) partial correction-based approach [12], and dynamic causal model [13] respectively. In this paper, we only consider second-order approaches, i.e., PC-based method and (regularized) partial correction-based method, since a recent study has empirically validated that they are more effective than many complex higher-order approaches.
Despite its simplicity and popularity, PC always produces dense FBNs that all nodes are fully-connected by edges, since it is sensitive to both direct and indirect functional connections in the brain. Thus, PC-based FBNs generally contain noisy or uninformative information. To remove the spurious information from the estimated FBNs, a common strategy is thresholding [14]. However, the optimal value of thresholding is currently an open problem. In contrast, partial correction is only sensitive to direct connections that can explain more complex interactions among brain regions. Even so, the estimation of partial correction matrix is usually ill-posed due to the singularity of the sample covariance matrix.
To over this problem, the regularization term is generally contained the partial correction model, which generates a sparse representation (SR)-based method.
Similarly, it is difficult to select the optimal value for the regularization term. In addition, both the above methods do not consider the spatial relationships of the BOLD time series.
In practice, according to previous work, the strength of functional connection between two ROIs will decrease as their distance increases. For instance, people in a society are analogous as a node (or ROI), and the social relationship between people is analogous as edge (or functional connectivity) in the estimated FBN [15]. If two persons are spaced close together, they may have similar rela-tionships or social circles. Thus, we simply assume that the brain regions close to each other tend to share stronger relationships or similarity. Inspired by this assumption, we propose a new approach for FBN estimation based on the spatial distance between ROIs. To verify the effectiveness of our model, we utilize the FBNs estimated by our method to identify the patients with MCI from health controls (HCs). Experiments conducting on a public dataset demonstrate that the proposed model outperforms than baseline methods (i.e., PC and SR). To better understand, we describe the whole pipeline of MCI classification task used in this study in Figure 1 that contains three main steps: 1) Data Acquisition, 2) FBN Estimation, and 3) Classification, which will be shown in Section II, III, and IV respectively.
The rest of this paper is organized as follows: in Section 2, we introduce the data sources and preprocessing step; in Section 3, we review several related works and develop our method for FBN estimation; in Section 4, we describe the experimental setting and verify the proposed method for MCI identification task; in Section V, we conclude this paper.

Data Acquisition
In this section, ADNI (http://adni.loni.ucla.edu) dataset was utilized to identify the participants with MCI from HCs based on the estimated FBNs to demonstrate the effectiveness of the proposed method. Concretely, for MCI vs. NC classification task, the preprocessed dataset with 137 subjects (68 MCIs and 69 HCs) was utilized in this study, in which all participants were scanned by 3.0T Philip scanners. Further, to remain the signal stabilization, the first three volumes of each participant were removed from the fMRI time series/BOLD signals. Then, the remaining volumes were processed by Data Processing Assistant for Resting-State fMRI (DPARSF) toolbox [16]. The preprocessing step can be summarized into the following steps: 1) head motion correction (i.e., subjects with head motion larger than 2 mm or 2˚ were excluded), 2) Friston 24-parameters nuisance regression model (i.e., the influence of the ventricle, white matter signals, and the high-order effect of head motion were removed ); 3) registration to standard Montreal Neurological Institute (MNI) space; 4) spatially smoothing and temporal band-pass filtering (0.01 -0.1 Hz). Finally, according to the automated anatomical labeling (AAL) atlas [17], the brain was divided into 116 ROIs, and all BOLD time series of the whole ROI was placed into 137 116 Figure 1. The classification pipeline used in our study.

Pearson's Correlation
PC, as the most popular and simplest method to estimate FBNs, can be defined as follows mathematically: where , 1, , is the mean time series extracted from ith brain region, T is the number of time points in each series, m is the total number of ROIs, is the connection weight between ith ROI and jth ROI. In particular, we redefine i.e., the new i x has been centralized by i i x x − and normalized by which corresponds to the optimal solution of the following form: ( Equivalently, Equation (2) can be further transformed into the following matrix model: where S is the connectivity matrix between ROIs, is the data matrix, and F ⋅ represents the F-norm of a matrix. In fact, FBN estimated by PC is very dense that contains noises or redundant connections. Thus, the threshold strategy is utilized to remove these redundant connections. However, selecting the optimal value of the threshold is still difficult.

Sparse Representation
As an alternative PC, SR is one of the commonly-used approaches to calculate the partial correction, which can estimate more reliable relationships between two ROIs by regressing out the confounding effect from other ROIs. The mathematical model of SR can be defined below: Accordingly, it can be further rewritten by the following matrix form: . It can be observed that the data fitting term of Equation (5) can obtain the sparse FBNs, since the data fitting term estimates the partial corrections between two ROIs, and the regularization parameter plays an important role to control the sparsity of the estimated FBNs. However, this method ignores the spatial relationship of pairs of time series.

Model
The proposed method can be described by the following optimization model: Equally, Equation (6) can be rewritten by the following matrix form in mathematic:  (7) is to prevent trivial solutions (i.e., 0, , 1, , ij S i j m = ∀ =  ). To our best knowledge, the brain is not a fully connected network. To estimate the sparse FBNs, we can utilize the following step: 1) deciding the number of edges of estimated FBNs based on utilizing k-nearest neighbor (k-NN) [18], that is, the brain region i and j are connected by an edge if i is among k nearest neighbors of j or j is among k nearest neighbors of i; 2) choosing the weights to edges between ROIs based on Equation (7) reasonably. In particular, the brain region i and j are spatially connected close together, they will share the stronger similarity. In contrast, they will obtain the small weight. In Table 1, we summarize the main algorithm for solving the proposed FBN estimation methods.

Experimental Setting
In our paper, we utilize PC and SR as the baseline methods to compare our method. Note that, the proposed method constrains the edge weights to be nonnegative. Thus, besides PC and SR, we also select PC+ and SR+ as the baseline methods based on the original PC and SR. In particular, PC+ and SR+ only keep the positive edges and turn the negative edges into zero in PC and SR. Since the threshold, regularization parameter or the number of nearest neighbor involved in both five methods (i.e., PC, PC+, SR, SR+ and ours) play a significant impact on the final classification accuracy, we set the threshold corresponding to different sparsity level in the set of [0, 10% •••, 90%, 99%] for PC-based method, in which each percentage represents the proportion of edges that are removed, set the regularization parameter in the range of [2 −5 , 2 −4 , •••, 2 4 , 2 5 ] for SR-based Step 1 Construct the connectivity matrix based on k-NN Step 2 For 1, , N ε =  Calculate S by solving Equation (7), and initial objective function value of Equation (7): 0

FBN Estimation
In this section, we take one of the participants from ADNI dataset as an example to visualize the FBNs estimated by five different methods (i.e., PC, PC+, SR, SR+ and ours). In Figure 2, we show the experimental results, in which the thresholds or regularization parameters involved in PC-and SR-based methods are Journal of Applied Mathematics and Physics determined based on their best classification accuracy. Specifically, the thresholds used in PC and PC+ are 50% and 40%, while the values of the regularization parameter or nearest neighbor used in SR and SR+ are 2 2 , 2 5 and 22, respectively. It can be visualized in Figure 2 that PC-based FBNs produces dense FBN easily. The reason for this is that the full correction is utilized to model the network adjacency matrix, and thus generates false functional connections. Although FBN estimated by PC+ improves the sparsity, it is still denser than SR and SR+-based FBN. That is because the original PC contains more redundant connections. In contrast, due to the introduction of regularization term, the estimated FBN by SR and SR+ are sparse. Compared with the baseline methods, the FBN estimated by the proposed method looks very clean, and its topological structure is similar to that of SR and SR+-based FBNs. We can argue that the spatial relationship between ROIs plays an important role in FBN estimation.

Classification Results
In Figure 3, we report the quantitative results achieved by five different methods under four accepted p-values in MCI classification tasks. According to the experimental results, we could have three main observations. First, the proposed method generally achieves better classification performance, compared with the baseline methods (i.e., PC, PC+, SR and SR+). For instance, in light of the average ACC values under four accepted p-values (i.e., 0.001, 0.005, 0.01, 0.05), the proposed method achieves the values of 80.29%, 81.75%, 87.59% and 88.32%, which largely outperforms the best results of baseline methods. These results suggest that the proposed method could perform better classification performance. Second, in terms of SEN values under four p-values, the proposed method also obtains the highest classification results, indicating that the proposed method may have a practically meaningful advantage for timely diagnosis of MCI population. For the remaining metrics, the final results also show the effectiveness of the proposed method.

Sensitivity to Network Modelling Parameters
In practice, the final classification accuracies are affected not only by p-values, but also by different network modelling parameters (e.g., thresholds and regula-rization parameters). In order to investigate the influence of parameters on the final classification accuracy, we conduct MCI classification experiments, and report their results with different values of thresholds or regularization parameters under four accepted p-values in Figure 4. It can be observed that the results based on the baseline methods will significantly influence the parametric values.  Note that, the horizontal axis for PC-based, SR-based, and our methods represents different thresholds, the regularization parameters, and for SR-based methods show while the vertical axis for five methods represents the average classification. Figure 5. The discriminative features selected by t-test (p < 0.001) between MCI vs. NC classification task based on AAL atlas for brain parcellation. In particular, this figure is created by a specific function (i.e., Circular Graph) based on Matlab software).
In contrast, the proposed method can obtain the higher classification results, compared with the competing methods in most cases, indicating the results of our method is reliable.

Discriminative Features
To visualize which features (i.e., functional connections) contribute more to the final classification task, we select the top 10 most discriminative features for identifying MCI population based on t-test for feature selection, and show the results in Figure 5. In particular, the thickness of each arc shown in Figure 5 represents the discriminative power that is inversely proportional to corresponding to p-value. In addition, according to these features, we further find several of them correspond to the brain regions, such as the left postcentral region [19], the right precuneus region [20], are reported as potential biomarkers for MCI vs. NC classification task. The results further prove the effectiveness of the proposed method.

Conclusion
In this paper, PC and SR are utilized as the baseline methods to estimate FBNs.
PC is sensitive to the direct and indirect relationship between pairs of ROIs, and thus always produces false functional connections. In contrast, SR can construct sparse FBN based on the obtained BOLD time series, due to the introduction of the regularization term. Even so, this method also ignores the spatial relationship between pairs of time series. Inspired by this, we propose a novel approach based on the spatial distance, and conduct MCI vs. NC classification task to verify the effectiveness of our method. Experimental results show that the proposed method can obtain better accuracy than the baseline methods. In the future, we will propose more novel FBN estimation methods corresponding to practical questions.