Abstract
Mild cognitive impairment (MCI) is a heterogeneous condition associated with a high risk of progressing to Alzheimer’s disease (AD). Although functional brain network alterations have been observed in progressive MCI (pMCI), the underlying pathological mechanisms of network alterations remain unclear. In the present study, we evaluated neuropsychological, imaging, and cerebrospinal fluid (CSF) data at baseline across a cohort of: 21 pMCI patients, 33 stable MCI (sMCI) patients, and 29 normal controls. Fast eigenvector centrality mapping (fECM) based on resting-state functional MRI (rsfMRI) was used to investigate brain network organization differences among these groups, and we further assessed its relation to cognition and AD-related pathology. Our results demonstrated that pMCI had decreased eigenvector centrality (EC) in left temporal pole and parahippocampal gyrus, and increased EC in left middle frontal gyrus compared to sMCI. In addition, compared to normal controls, patients with pMCI showed decreased EC in right hippocampus and bilateral parahippocampal gyrus, and sMCI had decreased EC in right middle frontal gyrus and superior parietal lobule. Correlation analysis showed that EC in the left temporal pole was related to Wechsler Memory Scale-Revised Logical Memory (WMS-LM) delay score (r = 0.467, p = 0.044) and total tau (t-tau) level in CSF (r = –0.509, p = 0.026) in pMCI. Our findings implicate EC changes of different brain network nodes in the prognosis of pMCI and sMCI. Importantly, the association between decreased EC of brain network node and pathological changes may provide a deeper understanding of the underlying pathophysiology of pMCI.
Keywords
INTRODUCTION
Alzheimer’s disease (AD) is the most common type of dementia in the elderly, accompanied by progressive and irreversible cognitive decline. No effective treatment currently exists and this morbidity impacts the afflicted individuals and their families alike due to the accompanying financial and emotional stresses associated with caregiving. Recently, research has focused on improving our understanding of the pre-dementia stage, in the hopes of facilitating early intervention and inhibiting progression. Mild cognitive impairment (MCI) refers to the symptomatic pre-dementia phase of AD [1]. While longitudinal studies have demonstrated that the conversion rate from MCI to AD is approximately 8.3 to 28% [2, 3], more than half of all MCI patients will not progress to dementia over a 10-year follow-up period [4, 5]. Owing to the heterogeneous nature of MCI, a better understanding of the pathophysiological mechanisms of progressive MCI (pMCI) versus stable MCI (sMCI) may facilitate earlier diagnosis and treatment.
Extracellular amyloid-β (Aβ) and intercellular phosphorylated tau (p-tau) protein depositions are recognized as the characteristic pathological alterations in AD, which can be detected at a very early stage [6]. By using 11C-PIB PET, some studies found that pMCI had increased Aβ deposition in the posterior cingulate, frontal, temporal, and parietal cortices relative to sMCI [7–9]. In cerebrospinal fluid (CSF), decreased Aβ levels and elevated total-tau (t-tau) and p-tau levels were observed in pMCI relative to sMCI, and could also predict conversion from MCI to AD with an accuracy of 70–94% [10, 11].
Recently, resting-state functional magnetic resonance imaging (rsfMRI) has used to explore brain network functional alterations in patients with pMCI versus sMCI. Based on a priori subnetworks, several studies demonstrated that patients with pMCI showed disrupted functional connectivity in the default mode network (DMN) compared with sMCI [12–15]. Using whole-brain network analysis, another study found that the functional connectivity strengths of the left angular gyrus and middle occipital gyrus were significantly lower in pMCI relative to sMCI groups [16]. Although these studies indicate brain network functional disruptions in pMCI, the underlying pathological mechanisms of network disruptions remain unclear. Previous studies have demonstrated a correlation between brain network functional changes and pathological alterations in AD, MCI, and the normal elderly [17–20], thus we speculate that network alterations in pMCI could be related to pathological abnormalities.
In the present study, we investigated the whole-brain network organization of patients with pMCI and sMCI using fast eigenvector centrality mapping (fECM). This analytical method could detect each brain region serving as a hub, which has greater connectivity with other regions of the brain. The fECM permits the efficient exploration of global organization of the whole-brain network without a priori definitions of subnetworks and brain regions [21]. Second, we explored whether changes in brain network organization of pMCI and sMCI patients were related to cognition, and AD-related pathological markers such as amyloid and tau proteins in CSF.
MATERIAL AND METHODS
Study design
Data used in the preparation of this article were obtained from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) database (http://adni.loni.usc.edu/). The primary goal of ADNI has been to test whether serial MRI, positron emission tomography, other biological markers, and clinical and neuropsychological assessments can be combined to measure the progression of MCI and early AD. Participants in ADNI included cognitively normal older individuals, people with early or late MCI, and people with early AD The follow-up duration of each group is specified in the protocols for ADNI-1, ADNI-2, and ADNI-GO. Subjects originally recruited for ADNI-1 and ADNI-GO had the option to be followed in ADNI-2. The research protocol was approved by each local institutional review board, and written informed consent was obtained from each participant. For up-to-date information, see http://www.adni-info.org.
Subjects
All MCI patients and normal controls were obtained from the ADNI-GO and ADNI-2 database. ADNI criteria for MCI were: 1) subjective memory complaints, either self-reported, reported by a study partner, or reported by a clinician; 2) objective memory loss defined as scoring below an education-adjusted cut-off score on delayed recall of the WMS-R logical memory test; 3) a global Clinical Dementia Rating score of 0.5; 4) a Mini–Mental State Examination (MMSE) score of equal to, or higher than, 24 out of 30; and 5) general cognitive and functional performance sufficiently preserved such that a diagnosis of dementia could not be made by the site physician at the time of screening [22]. The pMCI was defined as diagnostic status changing from MCI to AD according to the NINCDS-ADRDA criteria [23] within four years. The sMCI was defined as diagnostic status remaining as MCI for a minimum of three years. Demographic, neuropsychological, CSF, and imaging data at baseline of pMCI and sMCI patients were collected. Criteria for normal controls were defined as: a MMSE score of equal to or higher than 24 out of 30; a Clinical Dementia Rating score of 0; and having no report of any cognition complaint. According to the criteria described above, 21 pMCI patients, 38 sMCI patients, and 29 normal controls were included in this study. Five sMCI patients were excluded due to poor imaging data quality, as detailed described in imaging preprocessing (below). Finally, the rest of the subjects (21 pMCI, 33 sMCI, and 29 normal controls) were used for further analysis.
Demographic, neuropsychological, and CSF data acquisition
Procedures for assessing demographic status, neuropsychological tests, and CSF biomarkers in the ADNI cohort were previously described [22, 24]. In the present study, data are limited to baseline measures. Demographic status was assessed using age, sex, education level, and APOEɛ4 status. The neuropsychological assessments were comprehensive and involved in multiple cognition domains, including: MMSE; Wechsler Memory Scale-Revised Logical Memory II (WMS-LM) immediate and delayed memory; Trails A; semantic verbal fluency (SVF); and the Boston naming test (BNT).
CSF biomarkers included: amyloid-beta 1–42 (Aβ1–42), total tau (t-tau), and phosphorylated tau at position 181 (p-tau181), which were all measured using the multiplex xMAPLuminex platform (Luminex) with Innogenetics (INNOBIA AlzBio3) immunoassay kit-based reagents. The samples for CSF analyses included 19 out of 21 pMCI patients, 27 out of 33 sMCI patients, and 25 out of 29 normal controls.
Imaging data acquisition
Both baseline structural and resting-state functional image data were downloaded from the ADNI database. All subjects were scanned using a 3.0-Tesla Philips MRI scanner. The structural images were acquired using a 3D MPRAGE T1-weighted sequence with the following parameters: repetition time (TR) = 2300 ms; echo time (TE) = 2.98 ms; inversion time (TI) = 900 ms; 170 sagittal slices; within plane FOV = 256×240 mm2; voxel size = 1.1×1.1×1.2 mm3; flip angle = 9°; bandwidth = 240 Hz/pix. Resting-state functional MRI (rsfMRI) images were obtained using an echo-planar imaging sequence with the following parameters: 140 time points; TR = 3000 ms; TE = 30 ms; number of slices = 48; slice thickness = 3.3 mm; spatial resolution = 3.31×3.31×3.31 mm3; flip angle = 80°; and matrix = 64×64.
Imaging preprocessing
Before image preprocessing, all structural and resting-state functional images were visually inspected. Four sMCI patients (one with an infarct in the left centrum semiovale and three with poor rsfMRI image quality) were excluded from the study. Image data were preprocessed using the Data Processing Assistant for Resting-state fMRI Advanced Edition (DPARSFA, Yan and Zang; http://rfmri.org/DPARSF). Firstly, the first 10 image volumes of rsfMRI images were discarded whilst awaiting signal equilibrium and the subject’s adaptation to the scanning noise, followed by slice timing and head motion correction (six-parameter rigid body). Imaging data with more than 3.0 mm maximum displacement in any of the x, y, or z directions or 3.0° of any angular motion were excluded (including one sMCI patient). Subsequently, the T1-weighted image was co-registered to the mean rsfMRI image based on rigid-body transformation, and spatially normalized to the Montreal Neurological Institute (MNI) space, then re-sampled into 3×3×3 mm3 cubic voxels. The rsfMRI images were spatially smoothed with a Gaussian kernel of 6×6×6 mm3 full width at half maximum. Finally, linear trend adjustment and temporal filtering (0.01 Hz<f<0.08 Hz) were performed. In order to remove any residual effects of motion and other non-neuronal factors, six head motion parameters, white matter signals, and cerebrospinal fluid signals were regressed out. In addition, the framewise displacement (FD) value was computed for each subject in order to correct for the head motion artifacts.
ECM
Pre-processed rsfMRI data were further performed with the fast ECM (fECM) tool [25], which yielded a voxel-wise measure of relevance to the functional brain network. The fECM was transformed by the Fisher’s z test. Specially, the human brain is organized as complex networks with small-world properties. Graph theory-based network analysis can be used to investigate the global organization of the whole-brain networks. Eigenvector centrality (EC) is a whole-brain network analytical measure based on graph theory by weighting every nodes based on their degree (both the number and the quality) of connection with the whole brain [26]. Therefore, brain regions with greater EC values mean greater connectedness with the rest of brain and vice versa. The fECM tool is faster and computationally more efficient than traditional ECM calculation methods, because it computes matrix-vector products without having to compute or store the connectivity matrix [25]. Recently, several studies have reported that fECM is an effective method to evaluate whole-brain network alterations by its successful application in AD and Parkinson’s disease [21, 27].
Statistical analysis
Statistical analysis of demographic, neuropsychological, and CSF data was performed using SPSS 21.0 statistical software. Normal distribution was tested using the Kolmogorov-Smirnov test. For continuous variables, we first performed a group-level analysis for pMCI, sMCI, and normal control groups using one-way analysis of variance (ANOVA) and a Kruskal-Wallis test. If group-level test results were significant, post-hoc pairwise comparisons were additionally performed (Bonferroni multiple comparison correction in parametric tests, Dunn multiple comparison tests in nonparametric tests). Binary data, such as sex and APOEɛ4, were compared by a chi-square test.
The statistical analysis of fECM was performed using REST software (http://www.restfmri.net). ANOVA analysis was performed to identify brain network organization alterations among the three groups with age, gender, education, APOE genotype, and head motion (FD_value) as covariates. The threshold was set at single voxel p < 0.05 and cluster size >86 voxels, corresponding to a corrected p < 0.05, using AlphaSim for multiple comparison corrections (http://afni.nimh.nih.gov/pub/dist/doc/manual/AlphaSim.pdf) within the gray matter mask (excluding the cerebellum due to scanning range). Areas showing significant group effects by ANOVA were used as mask for two-sample t-test post-hoc analyses. The threshold here was set at single voxel p < 0.05 and cluster size >27 voxels, corresponding to a correctedp P < 0.05 using AlphaSim for multiple comparison corrections. Signals from the clusters that showed significant differences between pMCI and sMCI patients were extracted to further test the correlation with cognition and CSF biomarkers using Pearson’s correlations.
RESULTS
Demographic, neuropsychological, and CSF data
There was no significant difference in age, sex, education level, and the percentage of APOEɛ4 allele between the three groups (p > 0.05). The pMCI and sMCI had no significant difference in neuropsychological scores (p > 0.05). In CSF biomarkers, pMCI had higher t-tau level than sMCI (p = 0.023). In addition, there was no difference in head motion (p > 0.05) resulting from the motion correction (Table 1).
EC differences between the three groups
The direct comparison between pMCI and sMCI revealed that pMCI showed significantly decreased EC in the left temporal pole and the parahippocampal gyrus and increased EC in left middle frontal gyrus. When compared with normal control, pMCI had decreased EC in the right hippocampus and bilateral parahippocampal gyrus, sMCI showed decreased EC in the right middle frontal gyrus and the superior parietal lobule (Fig. 1; Table 2).
Correlation of EC differences with cognition and CSF biomarkers in pMCI and sMCI
We studied the relationships between EC in left temporal pole, parahippocampal gyrus, and middle frontal gyrus with cognition and CSF biomarkers, given that EC in these regions significantly differed between pMCI and sMCI. In the pMCI group, significant correlations were only found between EC in the left temporal pole with WMS-LM delay score (r = 0.467, p = 0.044) and t-tau level (r = –0.509, p = 0.026). In the sMCI group, no association between EC alterations with cognition and CSF biomarkers was observed. Across MCI groups, we only found that EC in the left temporal pole was associated with t-tau level (r = –0.376, p = 0.009), and EC in left parahippocampal gyrus was related to t-tau level (r = –0.364, p = 0.013). Besides, after adjusting for age, sex, education, and APOE genotype, correlations between EC alterations and t-tau level remained statistically significant (Table 3, Fig. 2).
DISCUSSION
To our knowledge, this is the first study to explore the correlation of brain network alterations with cognition and CSF biomarkers in pMCI and sMCI patients. The key finding was that EC in the left temporal pole was associated with memory function and CSF t-tau level in pMCI patients.
We found that patients with pMCI showed decreased EC in the left temporal pole and the parahippocampal gyrus compared to patients with sMCI, and decreased EC in the right hippocampus and bilateral parahippocampal gyrus relative to normal control. This finding is consistent with the “disconnection hypothesis” of AD [28]. Regions with decreased EC are recognized as important nodes of the hippocampal network [29], which plays an important role in declarative memory processing and is vulnerable in AD for the convergence of amyloid deposition, brain atrophy, functional disconnection, and hypometabolism [30]. Previous studies have also demonstrated the role of structural and perfusion alterations of hippocampus, parahippocampal gyrus, and temporal pole in the progress of AD. By example, morphology studies have reported that volumes of hippocampus, parahippocampal gyrus, and temporal pole could predict conversion from MCI to AD[11, 32]. Diffusion-tensor imaging studies revealed that pMCI showed disrupted white matter integrity in the hippocampus and the parahippocampal gyrus [33], as well as in the fornix and the posterior cingulum (the information transmission pathways from the hippocampal network to the frontal and parietal lobes) relative to sMCI [34, 35]. Several studies also found that hypoperfusion in the medial and inferior temporal cortices could discriminate pMCI patients from sMCI patients [36, 37]. Thus by combining our findings with previously reported results, we propose that the structural and functional alterations in those regions within hippocampal network are related to progression from MCI to AD. In addition, we found increased EC in left middle frontal gyrus in pMCI compared to sMCI. Binnewijzend et al. [21] also found greater EC in AD compared to normal elderly using fECM, which supported our finding. Increased brain network organization in pMCI may indicate a compensation mechanism for the decreased EC in regions within hippocampal network. It should be noted that previous studies on functional alterations of brain network between pMCI and sMCI were mainly within DMN [12–14]. The inconsistencies with our findings may be attributed to different analytical methods in rsfMRI, or the heterogeneity of MCIpopulation.
Furthermore, in the pMCI group, we found that EC in the left temporal pole was positively associated with WMS-LM delay scores, which mainly reflect the episodic memory performance. This result identified the role of temporal pole in memory function. This observation is largely congruent with one animal study that has reported that the temporal pole is related to memory function during an auditory delayed matching-to-sample task [38]. Importantly, supporting evidence was also from AD-related imaging studies. Specially, functional studies have demonstrated that dysfunction of the temporal pole was associated with episodic memory deficits in AD and MCI patients [39, 40]. In addition, structural [41, 42] and metabolic (using proton spectroscopy [43]) neurodegeneration of the temporal pole also significantly correlated with episodic memory impairment in AD and MCI patients. Collectively, memory is the earliest cognitive domain affected by AD, and measures of memory delayed recall are the best neuropsychological predictors of conversion from MCI to AD [44]. However, although there was no statistically significant difference in the WMS-LM delay scores between pMCI and sMCI groups, decreased EC in left temporal pole may reflect the earliest imaging-derived endophenotype underlying memory decline in patients with pMCI. In addition, our result highlighted hub organization of the left temporal pole which may contribute to memory function in pMCI patients.
In the present study, the pMCI group showed relatively higher t-tau level in CSF than the sMCI group, thus confirming previously reported results [11, 45] and suggesting that increased CSF level of t-tau was associated with an increased risk of progression from MCI to AD [46]. In addition, we further observed the significantly negative correlation between EC in the left temporal pole and CSF t-tau in pMCI group, which may have important implications in the pathogenesis of the pMCI. Particularly, the CSF concentration of t-tau is a sensitive biomarker of incipient neurofibrillary tangles in the brain, which are recognized as the pathological hallmark of AD. Neuropathologic studies in patients with AD also revealed heavy accumulation of neurofibrillary tangles in the temporal pole [47], which may disrupt synaptic plasticity and cause synaptic dysfunction [48]. Moreover, functional connectivity from rsfMRI is thought to reflect the co-varying neuronal signals in spatially distinct brain regions and is recognized as a biomarker of synaptic dysfunction [49]. Thus, we speculated that the association between decreased EC in the left temporal pole and the increased CSF t-tau level indicated that EC could also be a useful imaging marker of AD-related pathological abnormalities in pMCI patients. In addition, previous studies investigating the association of disruption of functional connectivity with CSF biomarkers in AD, MCI patients, and normal elderly also consolidated our findings. Specially, in AD, decreased functional connectivity of dorsal posterior cingulate cortex with retrosplenial cortex in DMN was related to decreased CSF Aβ1–42 level [19]. Li et al. [17] further reported ratio of Aβ42/P-tau181p was positively correlated with functional connectivity of left precuneus/cuneus within the DMN in MCI, implicating that functional connectivity within DMN is affected by pathological changes at early stage in AD. In a large number of older adults with normal cognition, Wang et al. [20] revealed decreased CSF Aβ42 and increased CSF ptau181 were independently associated with reduced functional connectivity in posterior cingulate and medial temporal regions within DMN. Combining our findings and those of previous studies, we propose that functional connectivity changes may reflect AD-related pathological changes of the pMCI patients. However, we did not observe the correlation of EC alterations with CSF Aβ1–42 concentrations in this study. One possible explanation for this is that disrupted EC of pMCI patients in this study mainly occurred in temporal areas, whereas Aβ plaques may start very early in the neocortex with relatively little in temporal structures[50, 51].
Additionally, we found that sMCI group manifested decreased EC in right middle frontal gyrus and superior parietal lobule relative to normal control. These two regions are important components of right frontal–parietal network, the functions of which include visuospatial integration, coordination, and execution [52]. Previous SPECT studies revealed that sMCI showed cerebral hypoperfusion in frontal and parietal regions relative to normal controls [36, 53]. These observations may suggest that sMCI could disrupt functions of regions mainly within the frontal–parietal network rather than hippocampal network regions, thus leading to the cognitive deficits.
It should be noted that there are inherent limitations associated with this study. First, the sample size of our study was relatively small. Although the ADNI database contains a large amount of data, the number of MCI patients with both resting-state functional MRI data and at least 3 years follow-up were relatively small. Second, due to the relatively short follow-up duration of this study, it should be noted that some MCI patients who were diagnosed as sMCI in the study may later progress to AD. Therefore, larger cohorts and longer-term follow-up of MCI patients will be necessary to validate these findings in future studies. Finally, our study lacked precise proportions of early and late MCI, amnesic and non-amnesic MCI, or single and multi-domain MCI, and the heterogeneity of this population leads to different pathologies. Therefore, future studies that incorporate sufficient numbers of the different types of MCI, that each progress to AD are needed.
Conclusions
In summary, our findings indicate that there is a different intrinsic network organization between pMCI and sMCI groups. Importantly, the association between decreased EC of brain network node and pathological changes may provide improved understanding of the underlying pathophysiology of pMCI.
Footnotes
ACKNOWLEDGMENTS
Data collection and sharing for this project were funded by the ADNI (National Institutes of Health [NIH] Grant U01 AG024904). The ADNI is funded by the National Institute on Aging and the National Institute of Biomedical Imaging and Bioengineering and through generous contributions from the following: Abbott, Alzheimer’s Association, Alzheimer’s Drug Discovery Foundation, Amorfix Life Sciences, AstraZeneca, Bayer Healthcare, BioClinica, Biogen Idec, Bristol-Myers Squibb, Eisai, Elan Pharmaceuticals, Eli Lilly, F. Hoffmann-La Roche and its affiliated company Genentech, GE Healthcare, Innogenetics, IXICO, Janssen Alzheimer Immunotherapy Research & Development, Johnson & Johnson Pharmaceutical Research & Development, Medpace, Merck, MesoScale Diagnostics, Novartis Pharmaceuticals, Pfizer, Servier, Synarc, and Takeda Pharmaceutical Company. 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 (
). The grantee organization is the Northern California Institute for Research and Education, and the study is coordinated by the Alzheimer’s Disease Cooperative Study at the University of California, San Diego. ADNI data are disseminated by the Laboratory for Neuro Imaging at the University of California, Los Angeles.
This study was also supported by the 12th Five-year Plan for National Science and Technology Supporting Program of China (Grant NO. 2012BAI10B04) and the National Natural Science Foundation of China (Grant Nos.81371519 and 81301190).
