Academic Editor: Giovanna Zamboni
Background: The relationship between switching rate of multilayer functional network and cognitive ability in mild cognitive impairment (MCI) and Alzheimers’ disease remains unclear. Methods: We followed up MCI patients for one year and analyzed the association of switching rates with cognitive decline. The iterative and ordinal Louvain algorithm tracked the switching of functional networks, while elastic network regression and Bayesian belief networks were used to test the relationship between network switching rate and cognitive performance cross-sectionally and longitudinally. Results: The switching rate of the default mode network positively correlated with better cognitive function, while that of salience and executive control network was negatively associated with memory and executive function. The lower default mode network (DMN) switching rate predicted MCI progression to dementia, while the lower sensorimotor network switching rate heralded in slower cognitive decline. Conclusions: The present study investigated the predictive effect of switching rate on cognitive performance, as well as MCI progression to dementia. The inverse effect from different functional networks may become useful for early diagnosis and revealing the mechanism of neural networks in cognitive decline.
As the potential early stage of Alzheimer’s disease (AD), mild cognitive impairment (MCI) was largely studied about its progression to AD. About 10–12% of MCI patients converted to AD annually, and the rate became 25–30% at three years. Overall, more than 60% of MCI patients eventually developed AD [1, 2]. The judgment of MCI progression in the early stage was a crucial step to attain further medical intervention. Various predictive models have been generated, combing demographics, genetical factors, brain structure, and functional connectivity, amyloid deposition, and tauopathy [3, 4, 5, 6]. The comprehensive analysis of multimodal information has made these models of relatively high accuracy. However, we still need to explore how the cognition related brain activity changes in this critical stage, as well as predictive effect of brain network.
During the past two decades, resting-state functional magnetic resonance imaging (rs-fMRI) has developed rapidly, which has promoted the understanding of the human brain. The neural network was regarded as the basis of human cognition. However, the traditional study of the brain network was based on the static description of brain functional connectivity (FC) from rs-fMRI. In the process of cognitive activity, the neural network reorganized dynamically . Based on the analysis of temporally dynamic changes of the blood oxygen level-dependent (BOLD) signal in the resting state, it is found that the dynamic functional network connection (dFNC) modes fluctuated spontaneously.
The dynamic functional connection of the brain network was tried to predict the change of cognition, which greatly promoted the understanding of brain cognition, brain development, and the neural mechanism of brain diseases . However, the dynamic FC or dFNC analysis still did not reflect real temporal changes of specific networks, as there would be a gap between BOLD signals and real neural network. Moreover, traditional clustering method, such as k-means or distribution-based clustering, did not fully make use of the dynamic temporal information.
Multilayer network analysis provided a novel method to process the complex data of multivariate and multi-scale information . The multi-layer modularization algorithm introduced the time and space of networks, and the decomposed networks had a non-overlapping period and space. The multilayer modular model could quantify the dynamic FC and track the temporal changes of each network. In this way each network had a unique switching rate, which refers to the degree of network switching between different modules in the multilayer network. In other words, the multilayer modularity model and switching rate allowed us to track the temporal network of a specific one and quantify their changes .
Previous studies had used the switching rate as a marker to predict various cognitive performance in healthy people [7, 11, 12], as well as in the severity of mental and neurological diseases with cognitive impairment, such as schizophrenia and epilepsy [13, 14].
It is the first time that the dynamic network switching rate from fMRI is estimated as an effective index to predict cognition in MCI. We assume that the switching rate of networks from fMRI can predict the cognitive changes of MCI and become a new functional neuroimaging biomarker for MCI progression. Given that the network serves as the basis of cognitive activity, we also tend to find critical networks accounting for cognition changes during the MCI stage.
The study was observational and used data from the control group of MCI in the SIMPLE cohort (detailed information is available in ). The cohort recruited MCI patients with age between 50 to 80 years from memory clinics or community. After screened by neurologists, patients with cognitive decline caused by other diseases were excluded (including but not limited to cerebrovascular disease, central nervous system infections, Parkinson’s disease, metabolic encephalopathy, deficiency of folic acid/vitamin B12 and hypothyroidism). The patients were diagnosed according to the National Institute on Aging-Alzheimer’s Association (NIA-AA) workgroups . Participants present mild cognitive problems based on the Mini-Mental State Examination and Clinical Dementia Rating scale (0.5). All MCI cases had medial temporal lobe atrophy (MTA) score above Grade II. Further, we exclude patients with MRI artifacts. The imaging and cognitive data were collected in the same day.
In total, we included 118 MCI patients (mean age, 65.39
After the baseline assessment, patients would be followed up without any
intervention for cognition and had the identical neuropsychological assessment
after 12 months. AD was diagnosed at the 12-month follow-up if patients met the
diagnostic criteria of dementia due to AD (described in , according to the
National Institute on Aging-Alzheimer’s Association (NIA-AA) workgroups ).
The AD criteria included decline from previous levels of functioning and
significant cognitive symptoms which impaired activities of daily living. The
diagnosis was made by a skilled clinician based on the individual circumstances
of the patient and were further were confirmed by
Baseline demographics include sex, age, and education level, listed in Table 1. Neuropsychological battery for multiple cognitive domains was administered to all patients, including Alzheimer’s Disease Assessment Scale-Cognitive subscale (ADAS-Cog) , the Auditory Verbal Learning Test (AVLT) - Huashan version , the Trail Making Test (TMT, including Part A and B) , the Rey-Osterrieth Complex Figure Test (CFT) , the Stroop Color-World Test (SCWT) and Boston Naming Test. In the present study, we had considered only participants who underwent the neuropsychological assessment both at the first time point and after 12 months.
|Stable MCI||Progressing MCI||p|
|APOE e4 carrier||36||8||0.4588|
|AVLT 5 min||4.67
|AVLT 20 min||4.52
|CFT recall 30 min||16.09
|M/F, Male/Female; sMCI, Stable mild cognitive impairment; pMCI, Progressed mild cognitive impairment; MMSE, Mini-mental State Examination; AVLT, Auditory Verbal Learning Test; CFT, Complex Figure Test; SCWT, Stroop Color-Word Test; TMT, the Trail Making Test; BNT, Boston naming test; ADAS-Cog, Alzheimer’s Disease Assessment Scale–Cognitive Subscale. * Comparison between baseline and 12-month.|
MRI images were acquired using a 3.0T uMR-890 MRI scanner (United Imaging, Shanghai, China) with a 64-channel coil and simultaneous multislice imaging techniques. The patients were placed in a supine position quietly with their eyes closed but asked to stay awake during the MRI. The MR sequences included high-resolution 3-dimensional T1-weighted imaging and resting-state fMRI. We collected data in Clinical Neuroscience Center, Ruijin Hospital LuWan Branch, Shanghai Jiao Tong University School of Medicine, Shanghai, China.
High-resolution T1-weighted images were acquired as the templates of further
functional images for coregistration. The sequences had the following parameters:
repetition time (TR) = 7.52 ms, echo time (TE) = 3.4 ms, field of view (FOV) =
The rs-fMRI data were preprocessed using the SPM12 software package (v7487,
Wellcome Department of Imaging Neuroscience, London, United Kingdom) based on
MATLAB 2020a (Mathworks Inc., Sherborn, MA, USA). Briefly, the preprocessing
steps consisted of removal of the first ten volumes of resting-state data for
magnetization stabilization, field bias correction, slice time correction, and
Friston 24-parameter motion correction. Spatial normalization was performed by
Diffeomorphic Anatomical Registration Through Exponentiated Lie algebra (DARTEL)
tool for creating a group specific template and for normalizing functional images
to the common space . Spatial smoothing was performed in the resultant images
using a Gaussian kernel with an 6-mm full-width half-maximum with the modified
MATLAB toolbox “Data Processing & Analysis of Brain Imaging (DPABI, version
3.1)” . We excluded patients who had head motion
Then a group-level independent component analysis (ICA) was used to define the
25 brain nodes of interest (independent component, IC) as described in previous
studies [23, 24]. We performed spatial group independent component analysis (ICA)
implemented in the Group ICA of functional MRI Toolbox (GIFT v4.0a;
http://icatb.sourceforge.net) [25, 26]. Resting-state functional MRI data of all
patients were firstly decomposed into principal components for subject-specific
data reduction, and then spatially decomposed into 25 ICs, each of which
exhibited a unique time course profile, based on the Infomax algorithm . The
resultant data was converted to Z scores, and ICs with Z score
As depicted in Fig. 1, the sliding-window approach was used to explore time-resolved dFNC. Resting-state time-series data were extracted and segmented into a 30-repetition time (TR) window with a size of 21 s, which is convolved with Hamming function to mitigate edge artifacts of the windows and attenuate potentially noisy signals . The window was sliding step-wise by 1 TR along the 660-TR length scan, resulting in a total number of 631 consecutive windows. We chose 30 TRs segmented window length as the grade-off between dynamics of functional connectivity  and 1/f0 wavelength criterion proposed by Zalesky and Breakspear . Here, the pairwise Pearson’s correlation coefficient between 25 ICs was calculated in each 21s time window. We thus obtained 631 correlation matrix, and positive matrix values entered the next multilayer network switching analysis. Although some identified clusters might be related to noise, we did not exclude ICs from ICA for the following reasons as in the previous studies [14, 30]. In the preprocessing step, nuisance covariates, including white matter, cerebrospinal fluid and global signals were regressed out. Meanwhile, the ICs identified from ICA could not be 100% from noise as they may be in the mid brain (i.e., IC 1). In addition, the switching rate represented the dynamic change of each IC, and it was not greatly influenced by other ICs which might be noise. The elastic net regression analysis automatically selected switching rate of ICs and it could further rule out switching rate of noise.
Flow diagram shows data processing and patients grouping for the present study. Abbreviation: MCI, mild cognitive impairment; IC, independent component.
To quantify spatiotemporal network switching, we used an iterative and ordinal
Louvain algorithm to track network function over time , with the code
governing parameters (
Modularity was quantified by Q ranging from 0 (low network segregation) to 1
(high network segregation).
Networks had average modularity (Q) of 0.536
We used elastic net regression to test the association between network switching
rate of 25 ICs and the individual’s cognitive performance at baseline, as well as
its changes (12 months - baseline), adjusted for age, sex, and education, and
APOE genotype. Elastic net enabled data-driven regression analysis by enforcing
sparsity of regression output values and it provided automatic variable selection
by removing all ICs not related to cognitive performance. Network switching
rates, cognition data were normalized into z-scores to ensure all data were
scaled equally. According to previous study, we set the
Where X is the original values of our neuropsychological variables, and
In relationship with MMSE, AVLT, STT, CFT, and Boston naming test, positive
For assessing the association between groups (sMCI and pMCI) and ICs’ switching rate, a Bayesian belief network was introduced. Package Bnlearn (version 4.4, Marco Scutari, Lugano, Switzerland) and Caret (version 6.0-81, Max Kuhn, New London, CT, USA) in RStudio software (the R Project for Statistical Computing, R software, version 3.1.0, R Core Team, Boston, MA, USA) was introduced for building a Bayesian belief network and directed acyclic graph. To quantitively define their relationship, we calculated the mean receiver operating characteristic (ROC) curve for the results from the Bayesian belief network. For multiple comparisions, we used Bonferroni correction to minimize the chance of false positive rate.
The age, sex, years of education, and ApoE genotype were of no
difference between the two groups (Mann-Whitney U-test for continuous data or
Chi-square test for categorical data). The pMCI group had significantly declined
performance mainly in SCWT and ADAS-Cog test (p
We identified 25 ICs from baseline resting-state fMRI using the group ICA. Some ICs are topologically consistent with previous results on the networks by visual inspection [24, 32]: default mode network (DMN, IC 15), salience network (IC 9), auditory network (IC 10, 14), sensorimotor network (SMN, IC 17), visual network (IC 3, 11), executive control network (ECN, IC 5, 13, 16) and cerebellum network (IC 19, 22). The detailed spatial maps of independent components are listed in Supplementary Fig. 1.
The switching rate of each IC was calculated, as shown in Supplementary
Table 1, and normalized further. The regularization parameter in elastic
|R||p||Lambda||MinSE||ICs in the regression model|
|MMSE||0.471||0.0002||0.1153||0.0766||1, 3, 5, 13, 19|
|AVLT immediate||0.502||0.0001||0.1011||0.0625||3, 14|
|AVLT 5 min||0.574||0.0001||0.0931||0.0690||7, 14, 15|
|CFT recall 30 min||0.647||0.0001||0.0879||0.0629||10, 12, 17, 21|
|SCWT_Interfere||0.444||0.0005||0.0414||0.0211||5, 7, 15, 25|
|TMT_B||0.498||0.0001||0.0869||0.0374||1, 6, 19|
|ADAS-word recall||0.823||0.0001||0.0236||0.0313||7, 9, 14, 15, 16|
|IC, independent component; MinSE, minimal mean square error; MMSE, Mini-mental State Examination; AVLT, Auditory Verbal Learning Test; CFT, Complex Figure Test; SCWT, Stroop Color-Word Test; TMT, the Trail Making Test; BNT, Boston naming test; ADAS-Cog, Alzheimer’s Disease Assessment Scale–Cognitive Subscale; p, significance of regression model.|
Besides, the switching rate of the auditory network (IC 14) involved in the AVLT immediate and 5-min recall, and it also positively contributed to ADAS-Cog word recall. DMN (IC 15) also positively contributed to AVLT 5-min recall and ADAS-Cog word recall. The switching rate of SAN (IC 9) and ECN (IC 16) negatively associated with performance in ADAS-Cog word recall. Another component in ECN (IC 5) also negatively associated with MMSE and SCWT-interference (Fig. 2).
Cross-sectional switching rate and cognition. (A) Contribution of switching rate in each IC to cross-sectional cognitive performance. The positive effect suggested that a higher switching rate correlated with better cognition and vice versa. (B) The fitting between original cognitive performance and predicted data from the elastic regression model, with the input of demographics and switching rate. Abbreviation: IC, independent component; MMSE, Mini-mental State Examination; AVLT, Auditory Verbal Learning Test; CFT, Complex Figure Test; SCWT, Stroop Color-Word Test; TMT, the Trail Making Test; BNT, Boston naming test; ADAS-Cog, Alzheimer’s Disease Assessment Scale–Cognitive Subscale.
We compared the switching rate of each IC between sMCI and pMCI groups at
baseline, and IC 15 had a significantly higher switching rate in the sMCI group
|Changes||R||p||Lambda||MinSE||ICs in the regression model|
|AVLT 5 min||0.289||0.028||0.1034||0.0300||17|
|AVLT 20 min||0.554||0.001||0.0745||0.0379||1, 10, 17|
|CFT copy||0.294||0.025||0.0605||0.0158||10, 17|
|BNT||0.330||0.011||0.0842||0.0423||14, 20, 21|
|ADAS-Cog||0.554||0.001||0.0674||0.0376||2, 11, 17|
|sMCI vs pMCI||0.339||0.009||0.1691||0.1090||15|
|No. of IC, Number of independent component; MinSE, minimal mean square error; sMCI, stable mild cognitive impairment; pMCI, progressed mild cognitive impairment; MMSE, Mini-mental state examination; AVLT, Auditory Verbal Learning Test; CFT, Complex Figure Test; SCWT, Stroop Color-Word Test; TMT, the Trail Making Test; BNT, Boston naming test; ADAS-Cog, Alzheimer’s Disease Assessment Scale–Cognitive Subscale; p, significance of regression model.|
Switching rate and longitudinal cognition changes. (A) Contribution of switching rate in each IC to longitudinal cognitive changes. The positive effect suggested that a higher switching rate correlated with better cognition endpoints. (B) The mapping of IC 15 and the ROC curve for MCI grouping by IC 15’s switching rate. (C) The fitting between original cognitive performance changes and predicted changes from the elastic regression model. Abbreviation: IC, independent component; MMSE, Mini-mental State Examination; AVLT, Auditory Verbal Learning Test; CFT, Complex Figure Test; SCWT, Stroop Color-Word Test; TMT, the Trail Making Test; BNT, Boston naming test; ADAS-Cog, Alzheimer’s Disease Assessment Scale–Cognitive Subscale.
To make the results more consolidated, we used the Bayesian belief network and
found that only a lower IC 15 switching rate has an associated effect on the pMCI
group (Coefficients of IC 15: –8.135). Since we found the topological
association between progress and IC 15, we built a classifier just using IC 15
with a random forest model. The area under the curve was 0.820
It has been largely studied that dynamic functional network connectivity serves as a physiological biomarker of cognitive performance. The switching rate of nodes in the network provides novel time-vary connectivity changes in adjacent ICA time windows. Using the data-driven approach of independent component analysis, we observed a set of independently coherent nodes belonging to networks in MCI, which were coincided with functional networks, as previously reported . The relationship between switching rates and behavior has been extensively investigated. In healthy individuals, the whole-brain modularity steadily increased during training for both conditions of the dual n-back task, especially the autonomy of the default mode system. Its integration among task-positive systems was modulated by training . In addition, major depression and bipolar disorder patients decreased network switching rate of key hubs in default mode network . By contrast, higher ‘flexibility’ (switching between multilayer network communities) was suggested to be a feature of schizophrenia during working memory task . There was significantly higher flexibility in the thalamus due to default-mode sensory/motor transitions . The present study is the first one to investigate the predictive effect of resting state switching rate in MCI on cross-sectional cognitive performance, as well as MCI conversion to dementia.
In the cross-sectional part, the compelling finding was that, in MCI patients, more frequent switching in the DMN generally suggested better cognition, while the networks related to executive control (ECN) and attention (SAN) had contrary results. It was well-established that the frontal and parietal regions were crucial aspects of the ECN, which sent rich sensory information not only for movement controls, but also for other cognitive abilities, especially in executive function . A meta-analysis study suggested an increase of functional changes in the frontal and parietal regions of MCI , which was in line with our results. The attention network was associated with working memory and episodic memory encoding. Recent studies have found intra-network and inter-network functional disruptions in the DAN in MCI patients [39, 40]. Mind-body exercise significantly increased the selective attention of MCI patients, and meanwhile decreased functional connectivity in attention network .
However, the relationship between the executive control and attention networks with cognition was not found significant longitudinally. The relationship might be obscured by longitudinal u-shaped trajectory in inter-network FC between ECN and DMN . Another possible reason related to neural compensation, which involved frontal regions and the dorsal attention network. Neural compensation was suggested to be common in mild cognitive impairment and Alzheimer’s disease .
The brain regions of DMN include the precuneus/posterior cingulate cortex (PCC),
medial prefrontal cortex (MPFC) and medial, lateral, and inferior parietal cortex
, which are consistent with our IC 15 mapping. The resting-state DMN changes
in aging, MCI, and AD patients have been largely studied, and its activation and
inter-network connectivity are generally considered biomarker to dementia [45, 46]. The anterior aspect of the DMN was negatively related to cognitive decline
in the older group . In the individuals at risk for Alzheimer’s disease, APOE
We further focused on the switching rate of DMN and its relationship with MCI cognition. The resting brain’s functional organization was suggested to be configured to maintain a balance between network segregation and integration. This functional balance was associated with better memory. Furthermore, brain tending toward stronger segregation versus integration foster different cognitive abilities . Thus, our results suggested DMN switching rate as marker for cognitive abilities in MCI and AD.
Amyloid and tau burden also showed their association with DMN. DMN functional connectivity during rest is altered with increasing amyloid-PET signal levels in aging and AD patients [53, 54]. In the early phase of AD, hyperconnectivity in the DMN was associated with neocortical tau in the positive amyloid individuals, while hypoconnectivity was observed when tau level elevated . As switching of DMN is presumed to be related to increased information load on the related brain regions , our observations closely resemble the trend of DMN changes in MCI and AD. This reduction in DMN switching rate was closely linked to the clinical severity of MCI patients, especially the short-term memory and executive function by data-driven elastic net regression. The difference in IC 15 switching rate between groups was consolidated by elastic net regression for grouping and Bayesian belief network validation.
While the DMN is generally regarded as active when the brain is at wakeful rest, it is negatively correlated with other networks . We also observed inverse switching rate changes in some other resting-state networks. Worse memory and executive function in MCI correlated with increased switching rate in ECN and SAN. It implies the higher entropy or information load on the specific brain regions illustrated in IC 5, 9, 13, and 16 , especially on the bilateral inferior parietal cortex. The association cortex was known to integrate information between a range of different networks, and might act as compensation . The observed increased switching rate in ECN showed some resemblance to the functional connectivity changes in ECN of MCI and AD patients. The mild AD showed significantly increased resting-state functional in ECN compared to healthy controls, as well as increased rsFC in ECN compared to MCI . In the effective connectivity study for MCI and AD, in contrary to connectivity from the DMN to the other resting-state networks, increased connectivity was evident between the memory network and the ECN in the AD and MCI patients .
Longitudinal data in the present study also demonstrated the ominous effects of the higher switching rate of SMN. Several studies of Alzheimer’s disease have also identified changed functional connectivity in the sensorimotor network indicative of compensation or covert biomarker [62, 63]. Amnesic MCI showed increased FC between the supplementary motor area with the superior parietal lobe in the SMN, while the DMN has reduced inter-network connectivity with the SMN . Although statistical association was found between cognitive tests (i.e., AVLT, ADAS-Cog and CFT-copy) and switching rate of IC 17 in SMN, the association was not strong to further show the relationship between SMN and MCI/AD conversion. The potential compensation role of SMN in MCI stages need to be further studied.
It is suggested that the motor cortex is involved in the initial stages of AD,
and its hyperexcitability is a well-defined neurophysiological feature of early
AD [54, 65, 66]. The increased ability of excitation was regarded as its ability
to plastically reorganize itself via recruiting extra cortical circuits in the
SMN, or other alternative circuits , due to its natural distributed network
with multiple representations of the motor maps . The motor cortex
hyperexcitability would be due to an imbalance between excitatory and inhibitory
circuits, probably induced by
There are a few limitations that should be considered in the interpretation of our results. First, the parameter of temporal resolution and length of acquisition in BOLD signal sampling should be carefully checked for reliable results. Previous studies suggested much higher temporal resolution [14, 71] and longer resting-state session (57.6 min), which could merely be possible in clinical practice. Though temporal down-sampling of fMRI networks has similar results in the previous study , it is necessary to re-test the multilayer network modularity detection in a larger sample of MCI for further validation. Second, the length of the time window was discussed and varied in different studies [28, 33]. The reliability caused by different time window needs to be optimized by more accurate computational methods. Finally, due to the moderate sample size (118 MCI) and relatively short follow-up period for patients (1 year), the MCI conversion rate is low (14 in 118) and posed difficulty in constructing the predictive model. Further studies with larger sample size and longer follow-up periods would allow for the examination of MCI conversion.
In summary, the switching rate of resting-state networks serves as a novel biomarker for cognitive performance in MCI patients, as well as a predictive marker for disease progression. Among them, DMN switches as a representation of slow cognitive decline, while switching of SMN heralds worsening. The biomarker-guided predictive models will further permit researchers to identify and stratify MCI populations.
ZH and BL designed the research study. ZH and YD performed the research. BL analyzed the data. ZH and BL wrote the manuscript. All authors contributed to editorial changes in the manuscript. All authors read and approved the final manuscript.
All subjects gave their informed consent for inclusion before they participated in the study. The study was conducted in accordance with the Declaration of Helsinki, and the protocol was approved by the Ethics Committee of Ruijin Hospital, Shanghai Jiao Tong University School of Medicine (No. 2021-337).
We thank all the researchers in the SIMPLE study. We acknowledge the altruism of SIMPLE participants and their families.
This work was supported by Shanghai Rising-Star Program (21QA1405800).
The authors declare no conflict of interest.
Publisher’s Note: IMR Press stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.