Abstract
Alzheimer’s disease (AD) is characterised by cognitive impairment and progressive brain atrophy. Recent human neuroimaging studies showed that such AD symptoms are linked with anatomical and functional changes seen in the default mode network (DMN), but the key brain region whose atrophy disturbs the neural activity of the DMN and consequently contributes to the symptoms of AD remains unclear. Here, we examined the intrinsic neural timescales (INT) of the regions in the DMN aiming to identify such a crucial brain area. We investigated INT since prior work demonstrated its utility in quantifying local neural dynamics and bridging knowledge gaps between atypical local neuroanatomical features and symptoms in neuropsychiatric disorders. We used resting-state functional MRI data and compared INT between individuals diagnosed with AD and age-/sex-/handedness-matched cognitively normal individuals. To exclude the possibility that a region outside of the DMN has such a critical role, we first performed an exploratory whole-brain analysis and found that only the left angular gyrus, a region within the DMN, exhibited a shorter INT in the AD group compared to those with normal cognition. We also identified an AD-specific decrease in the grey matter volume in the left angular gyrus and revealed that such regional atrophy shortened the INT of the angular gyrus and consequently reduced the overall INT of the DMN. Moreover, we found that the overall shorter DMN INT led to the symptoms of AD, particularly, its impairments of attention control. Taken together, our findings indicate that the left angular gyrus serves as a key brain region whose structural atrophy and resultantly shorter INT destabilise DMN neural dynamics and contribute substantially to AD cognitive decline. Clinically, the outcome of this study indicates that INT of the left angular gyrus could serve as a convenient biomarker for the cognitive symptoms of AD.
Introduction
Alzheimer’s disease (AD) is the most common form of dementia, ostensibly due to progressive brain atrophy. In both human and non-human animal model systems, pathological findings show that cognitive impairment in AD and its neural degeneration are specifically associated with amyloid and tau depositions in myriad brain regions, including the medial and lateral temporal lobe and several parietal regions1-7.
Neural degenerative changes in AD commonly occur in regions within the so-called default mode network (DMN)8,9, including the medial prefrontal cortex, precuneus and angular gyrus4,10-15, and the atrophy seen in some of these DMN regions is specific to AD12-14.
Moreover, AD yields atypical activity in the DMN as assessed by resting-state functional magnetic resonance imaging (rsfMRI)16,17: specifically, AD patients showed significantly weaker functional connectivity within the entire DMN16,18 or only its posterior part19,20, and such atypical functional connectivity in the DMN was associated with amyloid4,19,21 and tau deposition22-24 within the DMN’s brain regions. Given the DMN’s crucial role in cognitive capability in older adults25, these findings led us to hypothesise that the anatomical changes within the DMN would disturb its neural dynamics, which then might lead to the cognitive impairment of AD.
However, it remains unclear whether one or more brain regions in the DMN have such a critical role in the pathophysiology of AD. Anatomically, multiple neuroimaging studies described AD-specific atrophy in the temporal-parietal junction12-14, but whether such focal atrophy in the DMN affects the network’s neural activity in AD was not examined. Several functional MRI studies with AD patients identified atypical neural dynamics in parietal regions within the DMN19,20,26; but these studies did not directly test whether such aberrant neural activation was attributable to local brain atrophy.
Consequently, the DMN region whose focal anatomical change induces network-level neural dysfunction that leads to AD symptoms in humans has not yet been identified.
Here, we examined the “intrinsic neural timescale” (INT) of resting-state neural activity for all brain regions27-29 aiming to identify such trigger neural zones that are crucial for the cognitive decline in AD.
INT, also known as temporal receptive windows30-35 or temporal receptive fields36, represents the time period during which a brain region integrates neural inputs from other regions37. Regions within the DMN and frontal-parietal network typically exhibit longer INT compared to other brain regions, which enables the DMN to incorporate diverse neural inputs from various remote brain areas into their intrinsic signal processing30-33,36,38-40. In contrast, sensory-related cortices, such as the primary visual cortex, typically exhibit a shorter INT, which permits these regions to respond to their inputs with higher temporal fidelity30,38,41,42.
In addition to such functional findings, we previously showed that this temporal index for local neural activity is related to neuroanatomical features: a brain region with a larger grey matter volume (GMV) more likely had a longer INT27. This relationship between INT and GMV has biological relevance insofar as a brain area with a higher neuronal density, as proxied by a larger GMV, tends to have more recursive neural processing, which prolongs the autocorrelation of its intrinsic neural activities, thereby yielding a longer INT.
Taken together, these previous findings indicate the possibility that INT could link local neuroanatomical changes to atypical neural activity and symptomatic behaviours in neurological and neuropsychiatric disorders27.
Given this background, we focused on INT aiming to identify a key DMN region whose atrophy affected its INT, which in turn then disturbed DMN neural dynamics, which then potentially led to AD symptoms. By analysing rsfMRI data obtained from AD patients and matched healthy controls, we tested the following four working hypotheses: (i) one or more DMN regions should exhibit a significantly shorter INT in people with AD; (ii) the overall INT of the entire DMN should also be shortened in AD individuals, and the DMN regions with shorter INT should have atypically smaller GMV; (iii) the smaller GMV of the specific DMN regions should induce the shorter INT of the regions, which results in shortened overall INT of the DMN; (iv) the shorter INT of the DMN regions or overall DMN should be associated with AD symptoms.
Materials and methods
Participants
We used an MRI dataset collected by the Washington University Knight Alzheimer Disease Research Center43, which has public availability as the Open Access Series of Imaging Studies (OASIS-3)44. The data acquisition procedures from people diagnosed with AD and matched healthy controls received ethics approval for human subject research by the Washington University Human Research Protection Office, and all participants or their guardians provided written consent.
Participants completed clinical assessment protocols in accordance with the Uniform Data Set (UDS)45, which included medical history, physical examination, neuropsychological test battery and neurological evaluation. The diagnosis of AD was made based on the Clinical Dementia Rating scale (CDR)46,47: individuals with more than a zero CDR score were diagnosed as having AD (‘AD dementia’), while those with zero CDR score were labelled as non-demented (‘cognitively normal’, CN)43,44,48.
Aiming to reduce within-group heterogeneity, we excluded AD individuals who showed relatively normal scores in the Mini-Mental State Examination (MMSE), instead focusing on AD individuals with MMSE scores ≤2349,50.
We also excluded AD participants whose clinical assessment occurred ≥365 days before the MRI scan, whose head motion during scanning exceeded 3 mm, or who had dementia-relevant cardiovascular diseases or strokes. The control group included age-/sex-/handedness-matched adults. With these criteria, we included 32 individuals with AD and 138 CN individuals (Table 1).
MRI data acquisition
The rsfMRI and anatomical MRI data were collected with a 3T MRI scanner (TIM Trio, Siemens) using a 20-channel head coil. The rsfMRI data were acquired using a multi-band echo-planar imaging sequence in a single run lasting ∼6 min (repetition time = 2.2 s, echo time = 27 ms, 36 slices, interleaved, flip angle 90°, in plane resolution = 4 × 4 mm2, slice thickness = 4 mm); for anatomical MRI, structural images were obtained using T1-weighted sequence (repetition time = 2.4 s, echo time = 3.16 ms, flip angle = 8°, in plane resolution = 1 × 1 mm2, slice thickness = 1 mm). The participants were asked to lie quietly with their eyes open during the rsfMRI scanning.
Preprocessing of rsfMRI data
SPM12 (www.fil.ucl.ac.uk/spm) after discarding the first five volumes. This preprocessing procedure consisted of realignment, unwarping, slice-timing correction and normalisation to a standard template (ICBM 152). We then performed regression analyses to remove the effects of head motion, white matter signals and cerebrospinal fluid signals, and finally conducted band-pass temporal filtering (0.01–0.1 Hz). Note that participants who had high head motion (mean motion >3 mm) were excluded, which resulted in no significant difference in the mean head motion (P > 0.1 in a two-sample t-test) or maximum/mean framewise displacement (P > 0.2 in a two-sample t-test) between the AD and CN groups.
Explanatory whole-brain analysis of INT
Using the preprocessed rsfMRI data, we evaluated INT for all brain regions at a single-participant level following previously used procedures27. First, an autocorrelation function (ACF) of the preprocessed rsfMRI signal (time bin = TR) was estimated, and the area under the curve of the ACF values was calculated within the initial period during which the ACF showed positive values. The point at which the ACF first reached zero was the upper limit of the initial period. After repeating this procedure for every voxel, we smoothed the resulting whole-brain INT map (Gaussian kernel, full-width at half maximum = 8 mm).
Using these INT maps, we estimated the group-mean maps by averaging them within groups and then compared the AD map to the CN map, using a random-effects model with age and sex information included as nuisance variables, setting a statistical threshold at PFWE = 0.05 at the cluster level.
INT and AD symptoms
For the brain areas showing significantly different INT between the AD and CN groups, we examined associations between their INT and the severity of AD symptoms. We calculated Pearson correlation coefficients between the average INT in these regions of interest (ROIs) and the MMSE scores in the AD group. The ROIs were defined as a 4mm-radius sphere centred at the peak coordinates found in the above whole-brain analysis, and an INT for each ROI was given as the average within the ROI. The effects of multiple comparisons between these ROIs were corrected using the Bonferroni correction. For these data and other statistical tests, we report exact probabilities from P = 0.05 to P =.0001 and use the less-than-equal symbol for P values less than 0.001.
INT of DMN
We investigated INT within the entire DMN since we hypothesised a critical role of the DMN in AD and also since the DMN regions emerged in the above-noted exploratory analysis (see Results). We defined the DMN based on a functionally-determined brain parcellation system51.
We first calculated INT within the DMN by averaging INT, within-subjects, for all rsfMRI time-series data across all the regions within the DMN. We next assessed whether INT of the ROIs emerging from the exploratory analysis (see Results) correlated with the overall DMN INT.
Regional GMV and INT
We also calculated the GMV of the ROIs that were detected in the above INT analysis. We estimated the GMV based on T1-weighted MRI images using SPM12 for each participant as follows. First, the MRI images were segmented into grey matter, white matter, and cerebrospinal fluid with the New Segment Toolbox52. Second, using the DARTEL Toolbox53, the segmented grey matter images were aligned, warped to a template space, resampled to 1.5-mm isotropic voxels and registered to a participant-specific template. The grey matter images were then smoothed with a Gaussian kernel (full-width at half maximum = 8 mm). Using these preprocessed T1-weighted images, we then calculated the GMV of the ROI that was detected in the above analysis. Next, we examined a Pearson correlation coefficient between the GMV and INT in the ROIs identified in the INT analysis.
Mediation effect of regional INT
To assess potential causal relationships between the local GMV, local INT and DMN INT, we then conducted a nonparametric mediation analysis. In the analysis, we set the independent variable using the GMV of the ROIs found in the exploratory whole-brain analysis and incorporated DMN INT as a dependent variable. The INT of the ROIs were used as a mediator variable. The statistical significance of the mediation effect was calculated using a bootstrapping method.
As a control, we also performed the same mediation analysis using the GMV and INT of the brain regions that were found in the above-mentioned exploratory whole-brain analysis but showed only a marginally significant correlation between its INT and AD symptoms.
As another control, we also conducted the same mediation analysis using the GMV and INT of the brain regions in the other hemisphere compared to the original ROIs that were identified in the whole-brain explanatory analysis and correlation test with AD symptoms.
Associations with AD symptoms
Next, we examined the causal relationships between the local INT, DMN INT and AD symptoms using a nonparametric mediation analysis. Before this analysis, we confirmed whether the DMN INT correlated with the AD symptoms.
We conducted two types of mediation analysis to test the following two models: (i) the atypically shorter INT of the ROIs found in the exploratory whole-brain analysis shortens DMN INT, which results in the AD symptoms; and (ii) the atypicality of the ROI INT is strong enough to account for the link between DMN INT and AD symptoms.
To test model (i), we set the ROI INT as an independent variable and defined the mediator variable using DMN INT. To test model (ii), we used DMN INT as an independent variable and adopted the ROI INT as a mediator variable. In both cases, the MMSE scores indicating the AD symptoms were set as a dependent variable.
After determining which model worked in the current dataset, we examined the specificity of the DMN in this brain-behaviour relationship. For this positive control test, we conducted another mediation analysis after replacing DMN INT with INT of the other six brain networks defined in the same brain parcellation system51 (i.e., frontal-parietal network, dorsal attention network, ventral attention network, limbic system, sensory-motor network and visual network).
Associations with different cognitive components
We also searched for cognitive components that were specifically affected by the changes in DMN INT. To this end, we conducted a regression analysis between DMN INT and the scores of the following neuropsychological tests45,54: ‘Digit Span Test ‘for the measurement of attention; ‘Trail Making Test Part A ‘and ‘Digit Symbol Test ‘to evaluate processing speed; ‘Trail Making Test Part B ‘to quantify executive function; ‘Logical Memory Test ‘to measure episodic memory;
‘Boston Naming Test ‘and ‘Category Fluency Test ‘to evaluate language function. All the scores were converted to z-scores. The scores of Trail Making Tests were inverted to represent that a higher value indicated better cognitive performance.
Results
Whole-brain analysis of INT
We first performed an exploratory whole-brain analysis of INT (Fig. 1A) to identify brain regions that showed an atypical INT in the AD individuals compared to the age-/sex-/handedness-matched controls (CN; Fig. 1B and Table 1). We found that the AD group yielded significantly shorter INT compared to CN in the left angular gyrus (AG) and right supramarginal gyrus (SMG) (P < 0.001 for the left AG, P = 0.002 for the right SMG, PFWE < 0.05; Fig. 1C and Table 2). No region showed a significantly longer INT in the AD individuals compared to the controls.
We then assessed potential associations between the severity of AD and the atypical INT observed in the left AG and right SMG. We found positive correlations between the MMSE and INT in both the left AG (r = 0.40, P = 0.022; Fig. 1E, upper panel) and the right SMG (r = 0.38, P = 0.033; Fig. 1E, lower panel); however, only the correlation seen in the left AG survived the Bonferroni correction for multiple comparisons. Given these results, we focused only on the left AG in subsequent analyses.
INT of DMN
Next, we examined whether the shorter INT with the AG filled the gap between anatomical changes of the AG and the atypical activity of the DMN. We investigated DMN activity because, as stated in the Introduction, we hypothesised that atypical neural dynamics of the DMN should be closely related to AD symptoms and also because the left AG is a component of the DMN (Fig. 2A).
To address this hypothesis, we first calculated the INT of the DMN and confirmed that DMN INT was significantly shorter in AD compared to matched controls (t168 = 2.33, P = 0.021, η2= 0.03 in a two-sample t-test; Fig. 2B). In addition, we found that DMN INT correlated with AG INT in the AD group (r = 0.80, P ≤ 0.001; Fig. 2C).
GMV of AG
To continue testing our hypotheses, we then calculated the GMV of the left AG and confirmed that AG GMV was significantly reduced in the AD group compared to that in the CN group (t168 = 4.05, P ≤ 0.001, η2= 0.09 in a two-sample t-test; Fig. 2D). Additionally, left AG GMV correlated with left AG INT in the AD group (r = 0.47, P = 0.006; Fig. 2E). This GMV-INT association is consistent with our previous finding27.
Mediation effect of AG INT
Using these aforementioned estimations, we then performed a nonparametric mediation analysis and examined the extent to which left AG INT explained the seemingly significant direct correlation between left AG GMV and DMN INT (r = 0.38, P = 0.03; Fig. 2F).
We found that the total effect of left AG GMV on DMN INT (c = 1.53, P = 0.005; Fig. 2G) was entirely explained by a mediation effect of left AG INT (c’ = 0.05, P = 0.09; ab = 1.48, P = 0.006). That is, this complete mediation model suggests that the smaller GMV of the left AG reduces left AG INT, which results in the shorter DMN INT.
To assess the specificity of the left AG GMV on DMN neural dynamics, we assessed whether GMV and/or INT of the other DMN regions mediated the total DMN INT. Since INT in the right SMG differed between the AD and CN groups (Fig. 1D) and showed a marginally significant correlation with cognitive decline in the AD individuals (Fig. 1E), we reasoned that SMG INT might also exert a mediation effect on the relationship between SMG GMV and DMN INT. However, this analysis failed to show such a significant mediation effect (P = 0.26; Fig. 2H). Moreover, we confirmed that INT of the right AG, another DMN region and counterpart of the left AG, did not have a significant mediation effect on the relationship between GMV of the right AG and DMN INT (P = 0.75). We also found no significant mediation effect for any of the other DMN regions (superior medial prefrontal gyrus, anterior cingulate, bilateral superior frontal gyri, bilateral inferior temporal gyri, bilateral parahippocampal gyri and middle cingulate). These negative findings suggest that the left AG has a specific role in linking its local anatomical change to the overall DMN activity in AD individuals.
Association with AD symptoms
We next investigated whether the atypically shorter INT occurring in both the left AG and overall DMN in AD patients could explain AD symptoms. Since both left AG INT and DMN INT correlated with AD symptoms (Fig. 1E for left AG INT and Fig. 3A for DMN INT), we tested the following two models to understand the interrelationships among left AG INT, DMN INT and AD symptoms, as indexed by the MMSE test. First, we assessed whether atypically shorter left AG INT reduces DMN INT and consequently contributes to AD symptoms, and second, whether the shortened left AG INT is so substantial that it can explain the association between DMN INT and AD symptoms.
We examined these two models using a mediation analysis and only found support for the first model. That is, DMN INT mediated the entire relationship between left AG INT and AD symptoms (P = 0.01; Fig. 3B). By contrast, left AG INT did not account for the significant direct association between reduced DMN INT of AD patients and their symptoms (P = 0.78; Fig. 3C).
To control for the possibility that DMN INT simply reflected a generalised decrease in neural dynamics across all brain regions, we assessed whether the INT of brain networks other than the DMN mediated the direct effect of left AG INT on the MMSE score. We found that, except for the DMN, none of the brain networks exhibited atypical INT that could explain the direct influence of the left AG INT on the MMSE score (mediation effects, P > 0.08). These null results provide further evidence that the DMN has a specific functional role in bridging the gap between atypical neural processing in the left AG and AD symptoms.
Considering this result (Fig. 3B) along with other findings (Fig. 2G), we infer that, in the AD patients, their atypical reduction in GMV of the left AG shortened its INT, shorter AG INT reduced DMN INT, and finally that the shorter DMN INT led to their cognitive impairment.
Specific cognitive components
Finally, we assessed whether DMN INT exerted specific effects on cognitive components in AD (Fig. 4). We found that DMN INT correlated significantly in the AD patients only with attention (R2 = 0.28, P = 0.002) and not with processing speed, executive function, memory, or language (P > 0.01). These results suggest that the alteration of DMN INT in AD has a specific effect on attention-related symptoms.
Discussion
By focusing on INT, as measured by rsfMRI, we identified the left AG as a key region whose atrophy induces its atypical intrinsic neural activity, which then appears to disturb neural dynamics within the overall DMN and consequently appears to contribute to cognitive impairment in AD, especially attention-related symptoms (Fig. 5). Given that INT length indicates the capacity of local information integration, these findings suggest that, in individuals with AD, the atrophy of the left AG affects the DMN’s capability to integrate a wide range of neural inputs to the network, which then appears to diminish efficient cognitive processing and leads to some of the AD symptoms, such as impaired attention control.
The ostensible linkage between shorter DMN INT and AD symptoms is consistent with previous studies on the DMN in healthy humans. That is, several human neuroimaging studies reported that regions within the DMN are more likely to exhibit longer INT than those in sensory-related brain networks29,55-57, which in turn permits the DMN to pool and integrate neural inputs from other regions for complex behavioural and cognitive processes38,41,55,56,58. In particular, the DMN has a central role in orchestrating diverse cognitive functions in older individuals25. Given these previous reports, the current findings on the adverse effects of the atypically shorter INT of the DMN on cognitive capacities seem neurobiologically reasonable.
The crucial role of the AG in AD has been supported by outcomes of previous studies. In rsfMRI studies59,60, the AG showed atypically weak functional connectivity with other brain regions.
Research employing positron emission tomography reported hypometabolism, atrophy, and amyloid and tau deposition in parietal cortical regions, including the AG, in AD26,61-64. The current observation about the importance of the left AG in AD cognitive decline represents a key observation that fills a gap between AG anatomy and function and provides a more integrated view on the biological mechanisms underpinning AD.
Our findings suggest the possibility that the left AG could represent a target area for non-invasive brain stimulation to slow or even modify AD symptoms. Indeed, accumulating evidence shows the beneficial effects of repetitive transcranial magnetic stimulation (rTMS) on AD65: for example, a recent rTMS study demonstrated that stimulation of the precuneus could mitigate the cognitive decline seen in AD patients66. Considering our findings, a brain stimulation method to modify INT of the AG would be a new approach to alleviate the progressive cognitive impairments of the neurodegenerative dementia.
In addition, INT of the left AG could serve as a biomarker for AD, which also is relatively facile to measure. If the current finding were verified in a larger sample, a mere 15-min structural and functional MRI scan would enable evaluation of the degree of the neurodegeneration and resultant dysfunction of the DMN in individuals with AD67.
One limitation of the current study is that we used a cross-sectional dataset with a relatively limited AD sample size. Nevertheless, the mediation analyses that we conducted indicated direct links between local atrophy, local brain dynamics, network neural dynamics and behaviour (Fig. 5). However, longitudinal studies with larger-sized cohorts will be necessary for the direct examination of the causal effects of such alteration of the local brain activity on AD symptoms.
Another limitation concerns the absence of directly investigating pathological links between INT and amyloid and tau depositions, which was far beyond the scope of the employed data set. However, we confirmed the effects of the GMV reduction on the shorter INT in the left AG, and such a decrease in GMV is thought to reflect neuronal atrophy68. Future studies would be necessary to directly examine the relationships between INT and amyloid/tau depositions in individuals with AD.
Data Availability
All data produced in the present study are available upon reasonable request to the authors.
Acknowledgements
This work was supported by Grant-in-aid for Research Activity from Japan Society for Promotion of Sciences (SAM, 22K20879; TM, 22K07512; TW, 19H03535, 21H05679, 23H04217), AMED (TM, 22dk0207061s0101) and JST Moonshot R & D Program (TW, JPMJMS2021).
Footnotes
Competing interests The authors declare that no competing interests exist.