Abstract
Alzheimer’s disease (AD) is associated with abnormal resting-state network (RSN) architecture of the default mode network (DMN), the dorsal attention network (DAN), the executive control network (CON), the salience network (SAL), and the sensory-motor network (SMN). However, little is known about the disrupted intra- and inter-network architecture in mild cognitive impairment (MCI). Here, we employed a priori defined regions of interest to investigate the intra- and inter-network functional connectivity profiles of these RSNs in longitudinal participants, including normal controls (n = 23), participants with early MCI (n = 26), and participants with late MCI (n = 19). We found longitudinal alterations of functional connectivity within the DMN, where they were correlated with variation in cognitive ability. The SAL as well as the interaction between the DMN and the SAL were disrupted in MCI. Furthermore, our results demonstrate that longitudinal alterations of functional connectivity are more profound in earlier stages as opposed to later stages of the disease. The increased severity of cognitive impairment is associated with increasingly altered RSN connectivity patterns, suggesting that disruptions in functional connectivity may contribute to cognitive dysfunction and may represent a potential biomarker of impaired cognitive ability in MCI. Earlier prevention and treatment may help to delay disease progression to AD.
Keywords
INTRODUCTION
Alzheimer’s disease (AD) is the most common neurodegenerative disorder, with symptoms typically including executive dysfunction, cognitive decline, and language deficits [1–3]. Mild cognitive impairment (MCI) is a transitional state between healthy aging and the dementia that eventually occurs in AD [4, 5]. The prodromal phase of AD before dementia may last for a decade or more. Therefore, MCI has received extensive investigative attention because of the high risk of converting to AD (annual conversion rate 10% –15% in MCI, compared to 1% –2% in the healthy elderly population) [6–9]. To define an incipient stage for disease detection, the Alzheimer’s Disease Neuroimaging Initiative (ADNI, ADNI GO and ADNI 2) has sub-classified MCI into early MCI (EMCI) and late MCI (LMCI). EMCI refers to performance between 1.0 and 1.5 standard deviations (SD) below the normative mean on a standard memory test, while LMCI refers to performance more than 1.5 SD below the normative mean on the same test[10, 11].
Resting-state functional magnetic resonance imaging (rs-fMRI) has been increasingly used to explore brain function by measuring correlations between spontaneous blood oxygen level-dependent (BOLD) signal of topographically-separated regions during “rest” (absence of task), since the first study by Biswal and colleagues [12]. Correlated spatially-distributed brain regions have been classified into several resting-state networks (RSNs), which reflect functional integration and segregation of brain activity [13, 14]. Convergent evidence suggests that the default mode network (DMN), the executive control network (CON), the salience network (SAL), the dorsal attention network (DAN), and the sensory-motor network (SMN) are among the most reproducible RSNs [15–19]. These RSNs not only reflect human brain structure and function but also provide potential clinical value as sensitive markers of diseaseprocesses [20–22].
Many studies have reported functional abnormalities and structural atrophy of DMN, in support of the hypothesis that this RSN is preferentially affected by AD [21, 23–25] and MCI [17, 27]. Previous longitudinal studies of MCI and AD revealed longitudinal functional connectivity patterns of DMN deficits using independent component analysis (ICA) [28, 29]. In addition, DAN, associated with top-down attentional processing and externally-directed cognition [30], was reported to exhibit reduced functional connectivity in patients with MCI [17, 31]. A recent study revealed that AD pathology is present in sensory- and motor-related regions such as auditory and visual association cortex and motor pathways [32]. More importantly, the human brain may support information integration across spatially segregated brain regions (interactions between RSNs) in a manner that confers resilience against pathological attack [33, 34]. Several studies indicate that dysfunction of core networks of the brain such as DMN, CON, and SAL results in abnormalities in other networks with which these core networks interact [35–37]. Therefore, studying the organization of intra- and inter-network connectivity in the context of cognitive changes observed in individuals at high risk for AD may shed light on the neurological basis of cognitive decline in AD [17, 38].
The sub-classification into early and late MCI provides an opportunity to investigate alterations related to disease severity in early stages of the disease. We hypothesized that both intra-network functional connectivity and the interactions between RSNs would be impaired in MCI compared to the normal control (NC) group. We also hypothesized that longitudinal alterations in brain functional connectivity would covary with the severity of cognitive impairment. To test these hypotheses, we performed a region of interest (ROI) correlation analysis on a large, clinically-graded longitude sample. We constructed connected brain networks for 23 NC subjects, 26 patients with EMCI, and 19 patients with LMCI using a priori defined brain regions representing the five core RSNs [38]. We tested for disease-related alterations in functional connectivity with three different levels of analysis: ‘nodal integration’ (each node’s connectivity with the rest of the nodes); ‘RSN profile’ (intra- and inter-network connectivity patterns for each network); and ‘large-scale connectivity’ (pairwise connectivity between every ROI). We also tested for continuous relationships between brain network measurements and cognitive ability(Fig. 1).
MATERIALS AND METHODS
Standard protocol approvals, registrations, and patient consents
The dataset used in the current study was obtained from the ADNI database (http://adni.loni.usc.edu/). According to ADNI protocol, written consent was obtained from all subjects participating in the study, and the study was approved by the institutional review board at each participating site.
Subjects
Participants were recruited through ADNI GO and ADNI 2 from over 50 centers in the USA and Canada. This study utilized ADNI data through June 2014, comprising a total of 221 participants with 765 scans. Individuals in the present study were classified as subjects with NC, EMCI and LMCI according to clinical and behavioral measures provided by ADNI GO and ADNI 2 (http://adni.loni.usc.edu/). ADNI was launched in 2003 as a public-private partnership, led by Principal Investigator Michael W. Weiner, MD. The primary goal of ADNI has been to test whether serial magnetic resonance imaging (MRI), positron emission tomography (PET), clinical measures, neuropsychological assessments and other biological markers can be combined to measure the progression of MCI and AD.
The present study included 68 participants: 23 NC subjects (12 males, mean age 74.0 ± 5.4 years at first scan, mean interval between scans 6.5 ± 1.9 months); 26 patients with EMCI (10 males, mean age at first scan 71.2 ± 7.9 years, mean scan interval 6.2 ± 1.3 months); and 19 patients with LMCI (13 males, mean age at first scan 73.2 ± 8.5, mean scan interval 7.2 ± 2.5 months). Each subject underwent a battery of neuropsychological tests and functional tests, including the Alzheimer’s Disease Assessment Scale (ADAS), the Clinical Dementia Rating Scale (CDR), the Mini-Mental State Examination (MMSE), and the Functional Activities Questionnaires (FAQ). Details can be found http://adni.loni.usc.edu/methods/documents/. Demographic, psychological characteristics and statistics related to data preprocessing are summarized in Table 1.
The inclusion criteria were as follows: 1) at least two visits during the time of study; 2) CDR of 0 and MMSE scores of 26 or higher for the NC group; 3) CDR of 0.5 and MMSE scores of 24 or higher for the MCI group; 4) no change in categorical disease status at the follow-up visit. Additionally, because apolipoprotein E (APOE) ɛ4 allele is associated with an increased risk for AD and many studies have revealed that the APOE ɛ4 allele has potential effects on brain activity [20] and functional connectivity [39], we only included subjects with APOE genotype information in order to correct for the genotype effect on functional connectivity. The following exclusion criteria were applied: 1) obvious distortion on the raw imaging files; 2) translation or rotation in any axis of head motion larger than the size of one voxel during the scanning. For full ADNI inclusion/exclusion criteria, see http://www.adni-info.org.
Data acquisition
According to ADNI protocol (http://adni.loni.usc.edu/), the resting-state fMRI data was acquired on 3.0 Tesla Philips scanners (varied models/systems) while subjects were instructed to keep their eyes open. Images were acquired using an echo planar imaging (EPI) sequence with a repetition time (TR) ranging from 2250 ms to 3090 ms; echo time (TE) = 30 ms; flip angle (FA) = 80; slices = 48; field of view (FOV) = 212 mm RL, 198.75 mm AP, 159 mm FH.The imaging resolution in the X and Y dimensions ranged from 2.29 mm to 3.31 mm, and the slice thickness was 3.3 mm. For each subject, 140 volumes were acquired with a total scan duration of 420 seconds. Some individuals received extended scans with 200 volumes for a total duration of 600 seconds. For consistency, only the first 140 volumes of these scans were analyzed.
Data preprocessing
The resting-state fMRI data were pre-processed using the in-house Brainnetome fMRI Toolkit (http://www.brainnetome.org/brainnetometool) based on Statistical Parametric Mapping (SPM8, http://www.fil.ion.ucl.ac.uk/spm). Pre-processing steps included (1) slice-timing correction; (2) realignment to the first volume for head motion correction; (3) spatial normalization to a standard EPI template with re-slicing to 2 mm cubic voxels; (4) regression of nuisance parameters including linear drift, six motion parameters, and the mean time series of all voxels within the white matter and cerebrospinal fluid; (5) temporal filtering (0.01–0.08 Hz) to reduce high frequency and low frequency noise; (6) spatial smoothing with a Gaussian kernel (full-width-at-half-maximum [FWHM] 6 mm) [40–43].
The difference in gray matter volume between patients and normal controls may impact on functional connectivity [44, 45]. In order to control for this confounding variable, we performed voxel-based morphometry (VBM) analysis using the VBM8 toolbox (http://dbm.neuro.uni-jena.de/vbm/). The structural images were preprocessed through standard steps [46]. First, the subjects’ T1-weighted images were aligned with the ICBM template and were tissue-segmented according to the tissue probability map provided in VBM8. Second, the segmented GM was modulated by the non-linear normalization parameters to correct for individual brain size. Finally, the modulated and warped GM was smoothed by convolving with an isotropic Gaussian kernel (8 mm FWHM). The resulting gray matter volume was used as a covariate in the statistical analysis.
Regions of interest
In the present study, 36 voxels’ MNI coordinates taken from the study by Brier and colleagues [38] were used to generate regions of interest (ROIs) (for details see Supplementary Table 1). These ROIs were derived by maximizing the topographic concordance between results obtained by ROI connectivity mapping and by spatial independent component analysis [38]. Then, 36 6-mm-radius spherical ROIs centered on these coordinates were produced (Supplementary Figure 1). We extracted the mean time series within each ROI from the smoothed normalized fMRI time series for each subject. These ROIs represent seven brain networks corresponding to primary auditory, primary visual, somatomotor, DMN, DAN, SAL, and CON. Primary auditory, primary visual, and somatomotor areas were further integrated into a single sensory-motor network (SMN) for consistency with previous studies [38, 47], yielding five RSNs for subsequent analysis: DMN, DAN, SAL, CON, and SMN.
Estimation of regional functional connectivity
The Pearson’s correlation coefficients between mean time series of each pair of ROIs were computed for each subject. The correlation coefficients were converted to η values using an exponential function related to the connectivity “distance” between the two connected ROIs [48]. This step insured that all of the correlations were positive, which was necessary for the analysis of the changes between baseline and follow-up. Specifically, the conversion was η
ij
= e-ξd
ij
, where ξ = 2 is a positive constant [49], d
ij
= (1 - r
ij
)/(1 + r
ij
) is a hyperbolic correlation measure [50] denoting the distance between two nodes, and r
ij
is the Pearson’s correlation coefficient. Using these measures, we proceeded at three levels of analysis: ‘nodal integration’: The total connectivity degree, Γ, of node i was calculated using [49] as the measurement of nodal integration, i.e., each node’s connectivity with the rest of the brain. This measure is equivalent to the weighted degree of the node in graph theory. ‘RSN profile’: Intra- and inter- network analyses were performed as previously described [38]. The intra-network composite score, c, was defined by averaging the transformed correlation coefficients of all ROI pairs in a particular RSN as cN = 〈 ηij 〉, where i and j refer to ROIs in a particular network N and< >denotes the mean across ROI pairs. The inter-network composite score, c, was calculated by averaging the transformed correlation coefficients of all ROI pairs belonging to different RSNs as cX,Y = 〈 ηij 〉, where X and Y denote separate networks, i refers to an ROI in X, and j refers to an ROI in Y [38]. ‘large-scale connectivity’: Treating correlated brain regions as networks achieves data reduction and reduces the effect of sampling error across pairs of nodes, but this step may also obscure focal phenomena [38]. To address this concern, we further investigated the functional connectivity between all the possible pairs of ROIs. In other words, for the selected 36 ROIs, we evaluated the disrupted connectivity profile for each of the 630 ROI pairs among the three participant groups [43].
Statistical analysis
For each level of analysis, one-way analysis of variance (ANOVA) was performed for a group-wise comparison between the NC, EMCI, and LMCI groups after regressing out age, gender, APOE genotype and interactions between them as well as gray matter volume effects, both at the baseline visit and the follow up and longitudinal changes between them. A permutation test was used to control for Type I error. Specifically, for each connectivity measure, the test statistic was computed first for the actual experimental groups. Then, 10,000 permutations of group membership were performed, and the statistics were recalculated for the permuted groups to yield an empirical distribution of the test statistic. The permutation p value was determined by computing the proportion of the 10,000 permutations for which the permuted test statistic was greater than the original test statistic. Post hoc Student’s t-test were also performed.
To determine whether connectivity measures were correlated with cognitive ability, Pearson’s correlation coefficients were calculated between cognitive scores and the connectivity measures in MCI patients. Bootstrapping methods were applied to obtain the 95% confidence interval (CI) of the correlation coefficients. Briefly, the MCI population was resampled with replacement 10,000 times and correlation coefficients were calculated for each resample to generate the 95% CI. If the CI does not contain zero, this is an indication that the observed correlation is statistically significant.
RESULTS
Cognitive and neuropsychological tests
According to the inclusion and exclusion criteria, 68 subjects were included in the study. There were no significant differences in gender(= 0.139) or age (p = 0.400 for the first visit and 0.392 for the second) among the NC, EMCI, and LMCI groups. The averaged MMSE score was 29.2 ± 1.0 for NC group, 28.3 ± 1.2 for EMCI group and 27.6 ± 1.2 for LMCI at the baseline visit (ANOVA, p = 0.002). Compared to normal controls, MCI patients exhibited deficits in modified ADAS-Cog both at baseline (p < 0.001) and follow-up (p < 0.001) (mean ADAS-Cog 7.7 ± 4.3 for NC, 11.6 ± 6.1 for EMCI, and 17.4 ± 6.5 for LMCI at baseline; Table 1).
It has been reported that small amounts of head-motion during data scanning can have a substantial impact on estimations of functional connectivity and functional networks in the resting brain [51–54]. Therefore, we excluded subjects with translations or rotation larger than the size of one voxel in each scan. We also evaluated group differences in head motion among the three groups according to the criteria of Van Dijk and colleagues [53]. The results showed that the three groups had no significant differences in head motion (Table 1).
Nodal integration
For each ROI, we tested for group differences in nodal functional connectivity strength at each of the two time visits as well as for longitudinal changes between visits. Two regions of DAN (laIPS and raIPS) and three regions of SMN (rMC, SMA and lA1) showed significantly lower nodal strength in MCI comparing to NC in the original visit (Supplementary Figures 2A–E). Widely distributed regions (PCC, riTmp, laIPS, dmPFC, rACC, SMA, rV1, and lA1) showed significant correlations with ADAS-Cog (Supplementary Figures 2F-M and Supplementary Table 2). Nodal strength of two regions (riTmp and rPut) were affected in MCI patients at the follow-up visit (Supplementary Figures 3A, B). Additionally, three regions (lLP, lV1, and rV1) showed significant correlations between nodal strength and ADAS-Cog scores at the follow-up visit (Supplementary Figures 3C–E and Supplementary Table 2). The changes in nodal strength between visits showed significant group differences in two regions of the SAL network, rPut (permutation test, p = 0.022) and rIns (permutation test, p = 0.030) (Figs. 2A, B).Post hoc analysis further showed that the main effect of group was between NC and LMCI and between EMCI and LMCI in rPut, while the main effect was between NC and EMCI in rIns. Furthermore, the changes in nodal strength in mdThal (r = 0.312, 95% CI: 0.073, 0.312) and lV1 (r = 0.298, 95% CI: 0.027, 0.534) were strongly positively related to the changes in ADAS-Cog score (Figs. 2A, C; Table 2).
RSN profile: Intra-network connectivity
One-way ANOVA was performed on each individual RSN for each of the two time visits as well for the longitudinal change between visits. SMN, showed statistically significant group differences at the baseline visit (permutation test p = 0.023, Supplementary Figure 4A). There was a trend towards group differences in DAN (p = 0.052) between the three groups. The lower composite z-scores of SMN (r = –0.389, 95% CI: –0.612, –0.145) and CON (r = –0.308 95% CI: –0.499, –0.137) were associated with greater ADAS-Cog scores (Supplementary Figures 4B, C; Supplementary Table 3). No RSN in the follow-up was found abnormal either one-way ANOVA or correlation with ADAS-Cog at statistical level of 0.05. However, change in composite score of the DMN (r = 0.331, 95% CI: 0.066, 0.56) was related to change in ADAS-Cog from follow-up to baseline (Fig. 3A).
RSN profile: Inter-network connectivity
On the original visit, we found that the connectivity strength between the DAN and SMN showed significant group differences, with a difference between NC and LMCI driving the overall effect (Supplementary Figures 4D, E). In addition, the inter-network correlations of DMN-SAL, DMN-SMN, and DAN-SMN were significantly correlated with ADAS-Cog scores (Supplementary Figures 4D, F, H; Supplementary Table 3). There were no significant abnormalities in inter-network correlation at the follow-up visit. To quantify the longitudinal effects, we evaluated differences in longitudinal alterations across RSNs between the three groups. The change in correlation strength of the DMN-SAL pair (permutation test, p = 0.041) was identified to be affected with the main effect between the NC and LMCI groups(Figs. 3B, C).
Large-scale connectivity
At the baseline visit, abnormal functional connectivities were distributed widely between brain regions. The disrupted intra-network ROI connectivities were mainly concentrated in DAN and SMN, and the disrupted inter-network ROI connectivities mainly involved the DAN-SMN pair (Supplementary Figure 5A, Supplementary Table 4). Post hoc analysis revealed that the main effects were between NC and LMCI, and between EMCI and LMCI (Supplementary Figure 5C). Moreover, connectivity strength was associated with severity of cognitive impairment in MCI. The correlation between connectivity strength and ADAS-Cog score mostly involved PCC, riTmp, dmPFC, and SMA (Supplementary Figure 5B and Supplementary Table 4).
At the follow-up visit, the large-scale connectivity results showed disrupted connectivity distributed throughout the brain but most particularly involving riTmp and rPut (Supplementary Figure 6A and Supplementary Table 5). The group differences were mainly between NC and LMCI, and between EMCI and LMCI (Supplementary Figure 6C). Connectivity correlations with ADAS-Cog principally involved lLP, lV1, and rV1 (Supplementary Figure 6B and Supplementary Table 5).
In order to investigate the disrupted functional connectivity in MCI patients, we evaluated the longitudinal alterations of functional connectivity in MCI compared to NC. Our results showed abnormalities in intra-network ROI connectivity (mainly involving SAL and SMN) and abnormalities in inter-network ROI connectivity (mainly involving the DMN-SAL and DAN-SMN pairs) (Fig. 4A and Supplementary Table 6). The longitudinal alterations of functional connectivity in EMCI were greater than those observed in LMCI (Fig. 4C). Few ROI pairs showed correlations between changes in connectivity strength and changes in ADAS-Cog score (Fig. 4B and Supplementary Table 6).
DISCUSSION
The main innovation of the present study is the analysis of intra- and inter-network connectivity patterns at three parallel levels of analysis, from nodal to large-scale connectivity, in a large, clinically graded, longitudinal sample including early and late MCI. This study is also novel in elucidating that longitudinal change in intra- and inter-network connectivity is associated with change in symptomatic impairment measured in the MCI group. These alterations in intra- and inter-network brain organization complement the aberrant brain network architecture observed in subjects at high risk for AD.
Disrupted intra-network connectivity within individual RSNs
The human brain at rest is organized as a complex system with different functional modules termed RSNs— such as DMN, DAN, CON, SAL, and SMN— which are both specialized for processing distinct forms of information and integrated to fulfill the functions of daily life. These vital networks have been hypothesized to be vulnerable to diseases such as MCI/AD [17, 55]. The present study supports this hypothesis by showing alterations within the DMN in a longitudinal analysis of MCI patients. The change in connectivity strength within the DMN was correlated with the change in cognitive ability measured by ADAS-Cog (Fig. 3A). The connectivity strengths of PCC and of riTmp, a core hub region of DMN [56], demonstrated significant associations with cognitive ability (Supplementary Figure 2F, G). Another important region in the DMN, mdThal, showed significant correlations between changes in nodal degree and changes ADAS-Cog score (Fig. 2C). The DMN is associated with self-referential, episodic memory [23, 58], possibly representing the core network targeted in MCI/AD subjects [17, 60]. Our findings reveal that the regional interactions within the DMN become disrupted in early cognitive impairment, confirming the important status of this network for understanding the pathophysiological process of AD.
With disease progression, other representative RSNs become compromised, notably SAL. Our results demonstrate that the rIns and rPut, important regions of the salience network, show significant group differences in nodal strength (Fig. 2B). SAL is involved in the important function of detecting and adapting to salient external stimuli and internal events [61, 62]. Specifically, the insula, an integral hub region in SAL, detects stimulus signals and initiates control signals to coordinate brain activity [36, 63]. Our findings are supported by recent studies of the insula, which revealed regional atrophy and accumulation of Aβ in MCI [64], as well as gray matter loss and functional connectivity reduction in AD [65]. The impaired functional connectivity of this region may result in impaired functional connectivity throughout SAL and other brain networks.
Dysfunctions in sensory or motor networks are often taken as normal signs of aging, when in fact this dysfunction may represent an early sign of AD. AD patients suffer from sensory, cognitive, and motor impairments. A recent study revealed impaired connectivity between the thalamus and the rest of the visual system in EMCI and LMCI, suggesting reduced functional integrity of the visual system [66]. A recent review also emphasized the importance of sensory and motor regions to AD pathology [32]. In the present study, intriguingly, we observe significant group differences in connectivity strength of SMN at all three levels of analysis (Supplementary Figures 2, 4, and 5)). Moreover, differences in SMN (including important regions such as SMA, rV1 and lA1) were associated with alterations in ADAS-Cog score both at baseline (Supplementary Figures 2 and 4)and longitudinally between visits (Fig. 2C). These findings are supported by previous evidence that neurofibrillary tangles and Aβ plaques in AD may result in dysfunction of visual cortex [67, 68]. Moreover, SMN participates with supramarginal gyrus in modulating action recognition and episodic memory [69, 70], impairment of which is one of the key features of AD [71, 72]. Taken together, these results suggest that SMN abnormalities may facilitate the early diagnosis of AD.
Disrupted inter-network connectivity between RSNs
The human brain is a complex hierarchical system that coordinates a variety of specific functions, depending on a balance between local processing and global integration of information. Efficient behavior involves the coordinated activity of brain networks, and this coordination is the biological cornerstone for daily human cognition [33]. This coordination may represent the dynamic interaction of networks as reflected by the connectivity between RSNs [18, 73]. It has been proposed that abnormal cognition results from a disruption in the balance between different brain systems in AD/MCI [74, 75].
Previous cross-sectional studies have demonstrated inter-network disconnections between SAL, CON, and DMN associated with cognitive decline from normal aging to MCI and eventually to AD [65, 76]. In the present study, there was evidence of longitudinal change in DMN-SAL connectivity strength in MCI (Figs. 3, 4). Moreover, inter-network correlation strength (DAN-SMN, DMN-SAL and DMN-SMN pairs) was progressively reduced with diseaseprogression (Supplementary Figures 4 and 5). These findings are supported by previous studies which have provided convergent evidence that not only functional connectivity within RSNs but also between RSNs was impaired in patients with MCI/AD [17, 78]. Notably, two of these three RSN pairs involved the DMN. As discussed above, DMN consists of a set of core hub region that play a key role in modulating daily cognition [59]. It has been proposed that the dysfunction of one network may initiate disordered outputs to other networks, which are initially unaffected by the disease [38, 79]. This is in line with a recent study which supports the hypothesis of the transneuronal spread model that some toxic agents travel between networks [80]. Convergence evidence has revealed that the salience network can causally influence activity of the DMN [81, 82] and mediate engagement and disengagement of corresponding brain regions in response to stimuli [36], a process which is thought to be crucial for cognitive control and daily cognition [35, 84]. Therefore, we speculate that the SAL and SMN are necessary for the efficient regulation of activity in the DMN. A disturbance of this regulation may lead to inefficient cognitive control in MCI. Collectively, our results support the hypothesis that disruption between RSNs begins in networks that are affected early in the disease (i.e., DMN) and subsequently spreads to other networks. This hypothesis needs to be tested in future longitudinal studies.
Disruption and compensatory mechanism coexist in patients with MCI
Accumulating functional imaging evidence suggests that increased recruitment of brain regions may be a compensatory response in patients with cognitive impairment [85–87]. In the present study, we observed significant correlations between changes in connectivity strength within DMN and changes in ADAS-Cog scores in patients with MCI (Fig. 3A), which indicates that reduction of functional connectivity was correlated with cognitive impairment in MCI and AD [43]. This disrupted functional connectivity is a potential biomarker of cognitive impairment in MCI [17]. Additionally, the change in connectivity strength from follow-up to baseline was greater in EMCI patients than in LMCI patients, possibly because of a recruitment of other brain regions to compensate for functional deficits in LMCI. These conclusions are supported by previous cross-sectional studies in which compensatory increases in brain connectivity were identified in AD/MCI subjects [87–89]. Patterns of disruption and compensation may co-exist during the prodromal stage of AD. Taken together, these results strengthen the conclusion that interventions that either prevent disruptions or strengthen compensatory responses could delay conversion to AD.
Methodological considerations
One methodological consideration in our analysis was that the Pearson’s correlation coefficient was converted to the correlation index using an exponential transformation [49]. The advantage of this method is that negative correlations were transformed to positive values, to avoid the offsetting of positive and negative correlations when investigating mean effects. However, negative correlations exist intrinsically between different networks [90], and the amplitude of negative correlations in different networks varies as a function of frequency [91], indicating an elusive biological role for negatively-correlated brain activity in need of further investigation.
Our study addressed potential confounds of genotype and brain atrophy on the functional connectivity analysis. We elected to use of APOE genotype as an inclusion criterion. The APOE ɛ4 allele is associated with increased risk for AD and has potential effects on both brain activity [20] and functional connectivity [39]. Therefore, it was important to correct for the effect of genotype on functional connectivity. As for brain atrophy, although we corrected for total brain volume, localized atrophy in early AD could still confound results. If a volume of interest contained fewer gray matter voxels, then the fMRI signal would be predicted to be noisier. Therefore, in a reanalysis we corrected for the gray matter of individual volumes of interest (Supplementary Figure 7), showing convergent results with the original analysis.
Several other methodological limitations existed in the present study. 1) Only 140 fMRI volumes were used for each subject, which is a relatively short scan length and may result in relatively low signal-to-noise (SNR) in the fMRI analysis. To address this issue, we conducted a supplementary analysis demonstrating that the temporal signal-to-noise (tSNR) [92] in the selected regions is acceptable (Supplementary Figure 8). 2) The correlation analysis between functional connectivity and clinical variables (ADAS-Cog) may be susceptible to Type I error. Our use of an uncorrected p-value threshold of 0.05 combined with a permutation test (in the ANOVA analysis) and a bootstrap procedure (in the correlation analysis) was designed to limit Type I error, and we note that many results far exceeded this nominal threshold. 3) The EMCI patients and the LMCI patients were combined into one group to investigate the relationship between cognitive performance and functional connectivity, so between group differences could be contributing to the relationships that we report. Sub-group analysis revealed trend-level but not statistically significant correlations within each group separately, possibly due to the relatively small sample sizes and the limited ability of ADAS-Cog to detect subtle cognitive changes. 4) We included only two scans for each subject and limited the analysis to patients who did not change diagnostic groups between baseline and follow up visits. Future analyses could expand this area of analysis, for example, by looking explicitly at intra- and inter-network connectivity in subjects who progressed to AD at follow up. 5) Finally, the present study was an ROI-based analysis. Because the regions of interest were defined a priori, group differences may well exist in other parts of the brain. Further investigations using finer brain parcellations should be conducted in the future.
Footnotes
ACKNOWLEDGMENTS
This work was partially supported by the Natural Science Foundation of China (Nos. 81571062, 81270020, 81471120, 81371544, 61305143); the Natural Science Foundation of Beijing (No. 7152096); the Youth Innovation Promotion Association CAS (No. 2014119); the Beijing Nova Program (No. Z1511000003150112); and the Open Project Program of the National Laboratory of Pattern Recognition (No. 201407344).
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.; 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 & 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 (
). 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 SouthernCalifornia.
