Abstract

Resting-state functional magnetic resonance imaging (rs-fMRI) has been widely employed to examine brain functional connectivity (FC) alterations in various neurological disorders. At present, various computational methods have been proposed to estimate connectivity strength between different brain regions, as the edge weight of FC networks. However, little is known about which model is more sensitive to Alzheimer’s disease (AD) progression. This study comparatively characterized topological properties of rs-FC networks constructed with Pearson correlation (PC), dynamic time warping (DTW), and group information guided independent component analysis (GIG-ICA), aimed at investigating the sensitivity and effectivity of these methods in differentiating AD stages. A total of 54 subjects from Alzheimer’s Disease Neuroimaging Initiative (ANDI) database, divided into healthy control (HC), mild cognition impairment (MCI), and AD groups, were included in this study. Network-level (global efficiency and characteristic path length) and nodal (clustering coefficient) metrics were used to capture groupwise difference across HC, MCI, and AD groups. The results showed that almost no significant differences were found according to global efficiency and characteristic path length. However, in terms of clustering coefficient, 52 brain parcels sensitive to AD progression were identified in rs-FC networks built with GIG-ICA, much more than PC (6 parcels) and DTW (3 parcels). This indicates that GIG-ICA is more sensitive to AD progression than PC and DTW. The findings also confirmed that the AD-linked FC alterations mostly appeared in temporal, cingulate, and angular areas, which might contribute to clinical diagnosis of AD. Overall, this study provides insights into the topological properties of rs-FC networks over AD progression, suggesting that FC strength estimation of FC networks cannot be neglected in AD-related graph analysis.

1. Introduction

Alzheimer’s disease (AD) is the most common cause for dementia, which has an insidious preclinical stage in which pathological tau accumulates slowly until clinical symptoms are observable in prodromal stages [1, 2]. For this reason, the neuroscience community is focusing on investigating early signs of AD which could lead to the development of validated biomarkers. Now, abundant functional magnetic resonance imaging (fMRI) studies have reported that cognitive impairments in AD were associated with abnormal functional connectivity (FC) interactions between different brain regions [35]. As spontaneous neuronal activity plays an important role in reflecting the brain’s intrinsic mental state and human behavior [6], the investigation of brain resting-state FC (rs-FC) might facilitate our understanding of the neurophysiological mechanisms underlying AD progression. Therefore, resting-state fMRI (rs-fMRI) has been widely employed to study the alterations in neuronal activity of AD patients, with the measurement of resting brain synchronized activity through low frequency oscillations in the blood oxygen level-dependent (BOLD) signals [7, 8].

Currently, many rs-fMRI studies have showed altered functional connectivity in AD [911]. These studies assumed that FC patterns in the brain are static over the course of rs-fMRI scans. However, there is a growing consensus that the spontaneous fluctuations and correlations of signals between two distinct brain regions change over time, even in a resting state [1215]. It is expected that dynamic properties of time-varying functional networks may offer additional information for understanding AD mechanisms. Conventional methods do not account for temporal variability and are not sensitive to associated connectivity changes in AD. As the brain works as a dynamically integrated network [16, 17], topologic characterization of dynamic FC networks may be better for revealing AD-linked FC alterations, which might be difficult to achieve in the traditional static network analysis of rs-fMRI [18].

At present, dynamic FC patterns have been primarily investigated using sliding-window technique [1921], independent component analysis (ICA) [12, 22, 23], and leading eigenvector dynamics analysis (LEiDA) [24]. The sliding window technique has been repeatedly applied to explore how dynamic FC is affected by neurological disorders. In each windowed segment, between-region statistical association was usually estimated with Pearson correlation [11, 25] and dynamic time warping (DTW) [26, 27]. Pearson correlation directly measures the statistical linear relationships between pairs of BOLD series. DTW is an elastic matching algorithm, and as such can account for lag and shape differences between BOLD series [27]. ICA enables us to obtain dominant components by eliminating high-order statistical correlation of concatenated FC matrix [21]. Specially, group information guided ICA (GIG-ICA) showed promising potential for detecting altered brain dynamic functional connectivity [28, 29]. Differently, LEiDA captures the main orientation of BOLD signal phases over time by calculating the leading eigenvector for dynamic phase-locking matrix, which estimates the phase alignment between each pair of brain regions. These pioneering studies motivated us to comparatively investigate the sensitivity and effectivity of static and dynamic FC estimation methods in detecting alterations in rs-FC networks over AD progressive stages. In this study, after rs-FC networks were constructed using Pearson correlation, DTW, and GIG-ICA, topological analyses were performed. 54 subjects obtained from Alzheimer’s Disease Neuroimaging Initiative (ANDI) database were included, consisting of 18 healthy controls (HCs), 18 mild cognition impairments (MCIs), and 18 ADs. The feasibility and effectivity of Pearson correlation, DTW, and GIG-ICA in exploring AD-related alterations were compared and discussed.

2. Materials and Methods

2.1. Data

Data used in our study were obtained from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) database (http://www.adni.loni.usc.edu). Detailed descriptions of the demographic and clinical features of the participants are provided in Table 1. The subjects were examined on a 3T MRI scanner (Magnetom Trio Tim, Siemens, Erlangen, Germany). MRI session included a high-resolution structural T1-weighted image (sagittal MPRAGE;  ms,  ms, , and isometric voxel  mm3) and a 10-min resting-state fMRI ( ms,  ms, matrix, and  mm3). Participants were instructed to rest with eyes opened and stay awake in rs-fMRI scanning.

The entire processing flowchart in this study is shown in Figure 1. The tissues of gray matter, white matter, and cerebrospinal fluid were first segmented from T1-weighted MRI images. Gray matter was then coregistered into Montreal Neurosciences Institute (MNI) space. After original fMRI images were smoothed and coregistered into MNI space, whole-brain parcellation atlas in MNI space was employed to parcellate whole-brain into 132 parcels (Table 2), including cerebrum and cerebellum. ROI-specific time serial was obtained by averaging BOLD signal timeseries corresponding to the voxels within ROIs. At last, FC adjacent matrix computed with Pearson correlation, DTW, and GIG-ICA, respectively.

2.2. rs-FC Networks

In this work, CONN (functional connectivity toolbox, https://web.conn-toolbox.org/) and SPM12 (Statistical Parametric Mapping, https://www.fil.ion.ucl.ac.uk/spm/) were used to perform T1-weighted MRI and rs-fMRI preprocessing, including spatial coregistration, functional realignment and unwarp, slice-timing correction, outlier identification, segmentation and normalization, and smoothing [30, 31]. T1-weighted images were normalized into standard MNI space and segmented into gray matter, white matter, and CSF tissue classes using SPM12 unified segmentation and normalization procedure. rs-fMRI datasets were realigned using SPM12 realign and unwarp procedure, where all scans were coregistered and resampled to a reference image (the first scan) using b-spline interpolation. This procedure successfully addressed potential susceptibility distortion-by-motion interactions by estimating the derivatives of the deformation field with respect to head movement and resampling the functional data to match the deformation field of the reference image. Temporal misalignment between different slices of the functional data, induced by the sequential nature of fMRI acquisition protocol, was corrected using SPM12 slice-timing correction procedure, where the functional data was time-shifted and resampled using Sinc interpolation to match the time in the middle of each acquisition time. Last, functional data was smoothed using spatial convolution with a Gaussian kernel of 8 mm full width half maximum, in order to increase BOLD signal-to-noise ratio (SNR) and reduce the influence of residual variability in functional and gyral anatomy across subjects. Then, functional data was coregistered to the structural data using SPM12 intermodality coregistration procedure with a normalized mutual information cost function. This procedure estimated an optimal affine transformation between the reference functional image (mean BOLD signal) and the reference structural image (T1-weighted volume) that maximizes the mutual information between the two-modal imaging.

MNI-registered gray matter was divided into 132 regions (including 91 cortical regions, 15 subcortical regions, and 26 cerebellar regions) based on FSL Harvard-Oxford Atlas maximum likelihood cortical atlas and the automated anatomical labeling (AAL) template (Table 2) [32]. Then, whole-brain rs-FC networks were built, with gray matter parcels defined as nodes and the FC strength between the nodes considered as edge weights. Thus, weighted matrices were obtained for further analyses. In this study, for comparison, the edges of rs-FC network were estimated using three methods, including Pearson correlation, DTW, and GIG-ICA (https://www.nitrc.org/projects/gig-ica/). Static FC strength was estimated over entire rs-fMRI scan in terms of Pearson correlation coefficient and DTW similarity [27]. After FC matrix was obtained from covariance matrix with a window width of 20 TRs [33], the first (dominant) component extracted by GIG-ICA from the concatenated dynamic FC matrices would be used for topological analysis. As group information captured by standard ICA on the group level is used as guidance to compute individual subject specific independent components (ICs), GIG-ICA is able to obtain subject-specific ICs with better spatial correspondence, higher spatial, and temporal accuracy [28, 34].

2.3. Topological Metrics of rs-FC Networks

For topological analysis, we focused on global efficiency, characteristic path length, and clustering coefficient, which were calculated based on the estimated FC adjacency matrices using BCT toolbox [35]. Our goal is to use these metrics to investigate the FC disruption caused by AD progression. The metrics are summarized as follows.

Global efficiency is inversely proportional to the average shortest path length. This metric measures how globally efficient a network is in the sense of connecting distant nodes together [35, 36].

where is the shortest path length between node and node . is the number of nodes in a network.

Characteristic path length indicates a lower capacity to integrate information using shortest path routing, revealing the integration of a network structure [35].

The clustering coefficient indicates the extent of local interconnectivity or cliquishness in a network. It is defined as the ratio of the number of triangles a given node belongs to over the total number of triangles it could belong to [37]. where is the number of triangles around node and is the node degree.

2.4. Statistical Analysis

Once the graph measures were extracted, we statistically evaluated groupwise difference. After whole-brain static and dynamic FC networks of each subject was built, groupwise differences of topologic metrics were assessed using one-way analysis of covariance (ANCOVA) tests with age and gender as covariates [22, 38]. As an extension of analysis of variance (ANOVA), ANCOVA provides a way of statistically controlling the effect of covariables. In this study, value less than 0.05 (uncorrected) was considered statistically significant.

3. Results

3.1. rs-FC Adjacent Matrices

Figure 2 shows examples of whole-brain FC estimation from HC, MCI, and AD subjects. In this work, rs-FC adjacent matrices were estimated using Pearson correlation, DTW, and GIG-ICA. Figures 2(d)–2(g) denote group-level FC states (GSs), and Figure 2(c) is the dominant GSs which contribute the most information across the entire time-varying connectivity patterns. Figure 2(a) was computed with Pearson correlation, and Figure 2(b) was obtained with DTW.

3.2. Topological Measures

As shown in Figure 3, no significant differences were found in FC networks constructed with DTW, in terms of the both network-level metrics. In FC networks built with Pearson correlation, only MCI and AD groups were significantly identified according to the metric of characteristic path length. In GIG-ICA dominant component networks, significant differences between HC, AD, and MCI groups were found in terms of global efficiency. Overall, significant groupwise differences could hardly be detected in terms of characteristic path length. The mean and standard deviation values of the network-level topological metrics are reported in Table 3.

Figure 4 shows ANCOVA results regarding the nodal metric of clustering coefficient. The box plots of ANCOVA are depicted in Figure 4, and the metric values () are reported in Tables 46. In the rs-FC networks built using Pearson correlation, significant differences across HC, MCI, and AD groups were found for 6 brain parcels, including right angular gyrus (41), left angular gyrus (42), left superior division of lateral occipital cortex (44), right frontal operculum cortex (76), right central opercular cortex (78), and left planum polare (83) (Tables 2 and 7). Only three brain parcels were significantly identified in DTW-constructed networks, namely, left occipital pole (91), right cerebellum crus2 (110), and right cerebellum 7 (118) (Tables 2 and 7). In GIG-ICA dominant component networks, we found that total 52 brain parcels (Tables 2 and 7) exhibit significantly different clustering coefficient across HC, MCI, and AD groups. Figure 5 shows the 3D distribution of these significantly identified brain parcels (red dots).

4. Discussion

At present, there are a few computational models to estimate FC strength between pairs of brain parcels, which were used to investigate FC alterations in various neurological disorders. In this study, we performed topological analyses on rs-FC networks of 18 HCs, 18 MCIs, and 18 Ads, aiming to finding a FC strength computational method sensitive to AD progression. Groupwise differences in global and nodal metrics were examined using ANCOVA tests with gender and age as covariates. Our results showed that almost no significant differences were observed across HC, MCI, and AD groups in terms of global metrics, and that much more affected brain regions by AD were revealed according clustering coefficient extracted from GIG-ICA dominant rs-FC networks. These findings provide new insights into understanding of the macroscopic network mechanisms underlying AD progression.

The human brain is a complex and interconnected network characterized by an efficient small-world network with high local clustering [39, 40]. Compared with healthy controls, MCI and AD patients exhibit a disruption of the global integration of brain networks [5]. However, in this study, almost no significant differences (, uncorrected) among HC, MCI, and AD groups were observed in terms of network-level metrics derived from rs-FC networks, including global efficiency and characteristic path length (Figure 3 and Table 3). In other words, global topologic metrics extracted from rs-FC networks are not sensitive to AD progression. Interestingly, in GIG-ICA dominant rs-FC networks, global efficiency of MCI group significantly decreased compared with HC and AD groups (, uncorrected). By comparison, significant intergroup differences have been identified in structural connectivity (SC) networks [3, 41, 42]. This may due to different nature of FC and SC brain networks [43, 44]. FC reflects neuronal synchronization between brain regions and is vulnerable to physiologic and psychologic conditions. SC serves as the substrate for FC, with relative stability.

As a measure of local interconnectivity, clustering coefficient was estimated to examine FC disruption caused by AD. For rs-FC networks constructed with Pearson correlation and DTW, significant difference across HC, MCI, and AD groups were observed only in 6 (Table 4) and 3 (Table 5) regions. However, statistically significant differences among HC, MCI, and AD groups were observed in total 52 nodes of GIG-ICA dominant rs-FC networks (Table 6), mainly distributed in temporal, cingulate, and angular areas. [45] reported that AD patients exhibit the posterior medial temporal neurodegeneration that was associated with episodic memory disturbance. In line with the previous finding, our results also showed AD patients have low clustering coefficient in the temporal gyrus and cingulate region [4648]. However, some brain parcels displayed higher clustering coefficient over AD progression (Figure 4) such as right anterior division of inferior temporal gyrus (27), right anterior division of parahippocampal gyrus (62), right cerebellum crus2 (110), and right cerebellum 7 (118). This phenomenon can be interpreted as a compensatory mechanism to maintain normal cognitive function under AD pathology [12, 49, 50], which deserves further study. From the ANCOVA results, compared with Pearson correlation and DTW, we could conclude that GIG-ICA was more sensitive to AD progression. This may be attributed to that dynamic rs-FC and the derived FC topological metrics might help reveal changes in macroscopic neural interaction patterns underlying AD progression [23]. Pearson correlation and DTW are here applied to measure static FC. Dynamic FC reveals temporal patterns of FC and provides additional information to static FC, confirming that it might be an efficient way to facilitate our understanding of the underlying neurophysiological mechanisms of AD.

Several limitations should be addressed in the future work. First, an empirically validated fixed sliding window of 20 TRs was selected in GIG-ICA. The selection of sliding window size is still inconclusive and may impact the evaluation of rs-FC connectivity strength. FC changes with separate windows of various window lengths will be evaluated in the future study. Second, the sample sizes included in the current study are relatively small, and the replication study with more participants needs to be performed in the future to verify our findings. Third, the whole brain was divided coarsely into 132 regions based on FSL Harvard-Oxford atlas and AAL template for functional brain network construction. Different parcellation strategies are required to validate our findings in future studies. Finally, the association between the disrupted brain functional networks and the neurological conditions was not explored. Future longitudinal studies are necessary to explore potential altered dynamic FC as a clinical biomarker.

5. Conclusion

rs-FC has been often used to identify abnormal brain connectivity patterns in AD cohorts and to understand mechanisms of abnormal brain function, attempting to comprehensively explore the potential utility of rs-FC as a biomarker for AD progression. In summary, this study revealed that more affected regions can be significantly identified in terms of nodal clustering coefficient extracted from GIG-ICA dominant networks. As GIG-ICA dominant component was extracted from the concatenated time-varying rs-FC matrices, it may be concluded that dynamic FC analysis has the potential to provide more reliable and significant scientific findings than static FC studies. Our results also indicated the potential role of clustering coefficient in determination of early-stage AD patients, providing a clue to AD diagnosis and monitoring in clinical applications. In the next step, joint analysis between structural and functional networks could be a better way to reveal brain connectivity alterations over AD progression.

Data Availability

The data used to support the findings of this study can be requested from the corresponding author.

Conflicts of Interest

There are no conflicts of interest.

Acknowledgments

The research was supported by the Natural Science Foundation of Zhejiang Province (LY20E070005 and LY17E070007). Data collection and sharing for this project was funded by the Alzheimer’s Disease Neuroimaging Initiative (ADNI) (National Institutes of Health Grant U01 AG024904) and DOD ADNI (Department of Defense award number W81XWH-12-2-0012). ADNI is funded by the National Institute on Aging, the National Institute of Biomedical Imaging and Bioengineering, and through generous contributions from the following: AbbVie; Alzheimer’s Association; Alzheimer’s Drug Discovery Foundation; Araclon Biotech; BioClinica, Inc.; Biogen; Bristol-Myers Squibb Company; CereSpir, Inc.; Cogstate; Eisai Inc.; Elan Pharmaceuticals, Inc.; Eli Lilly and Company; EuroImmun; F. Hoffmann-La Roche Ltd and its affiliated company Genentech, Inc.; Fujirebio; GE Healthcare; IXICO Ltd.; Janssen Alzheimer Immunotherapy Research and Development, LLC.; Johnson & Johnson Pharmaceutical Research & Development LLC.; Lumosity; Lundbeck; Merck & Co., Inc.; Meso Scale Diagnostics, LLC.; NeuroRx Research; Neurotrack Technologies; Novartis Pharmaceuticals Corporation; Pfizer Inc.; Piramal Imaging; Servier; Takeda Pharmaceutical Company; and Transition Therapeutics. The Canadian Institutes of Health Research is providing funds to support ADNI clinical sites in Canada. Private sector contributions are facilitated by the Foundation for the National Institutes of Health (http://www.fnih.org/). The grantee organization is the Northern California Institute for Research and Education, and the study is coordinated by the Alzheimer’s Therapeutic Research Institute at the University of Southern California. ADNI data are disseminated by the Laboratory for Neuro Imaging at the University of Southern California.