Skip to main content

ORIGINAL RESEARCH article

Front. Aging Neurosci., 06 January 2023
Sec. Neurocognitive Aging and Behavior
This article is part of the Research Topic Perioperative Brain Homeostasis: From the Bench to the Bed View all 6 articles

Decreased default mode network functional connectivity with visual processing regions as potential biomarkers for delayed neurocognitive recovery: A resting-state fMRI study and machine-learning analysis

Zhaoshun Jiang,&#x;Zhaoshun Jiang1,2Yuxi Cai,&#x;Yuxi Cai1,2Songbin Liu,Songbin Liu1,2Pei Ye,Pei Ye1,2Yifeng YangYifeng Yang3Guangwu LinGuangwu Lin4Shihong LiShihong Li4Yan XuYan Xu1Yangjing ZhengYangjing Zheng1Zhijun Bao,,Zhijun Bao2,5,6Shengdong Nie
Shengdong Nie3*Weidong Gu,
Weidong Gu1,2*
  • 1Department of Anesthesiology, Huadong Hospital Affiliated to Fudan University, Shanghai, China
  • 2Shanghai Key Laboratory of Clinical Geriatric Medicine, Huadong Hospital Affiliated to Fudan University, Shanghai, China
  • 3School of Health Science and Engineering, University of Shanghai for Science and Technology, Shanghai, China
  • 4Department of Radiology, Huadong Hospital Affiliated to Fudan University, Shanghai, China
  • 5Department of Geriatric Medicine, Huadong Hospital Affiliated to Fudan University, Shanghai, China
  • 6Research Center on Aging and Medicine, Fudan University, Shanghai, China

Objectives: The abnormal functional connectivity (FC) pattern of default mode network (DMN) may be key markers for early identification of various cognitive disorders. However, the whole-brain FC changes of DMN in delayed neurocognitive recovery (DNR) are still unclear. Our study was aimed at exploring the whole-brain FC patterns of all regions in DMN and the potential features as biomarkers for the prediction of DNR using machine-learning algorithms.

Methods: Resting-state functional magnetic resonance imaging (fMRI) was conducted before surgery on 74 patients undergoing non-cardiac surgery. Seed-based whole-brain FC with 18 core regions located in the DMN was performed, and FC features that were statistically different between the DNR and non-DNR patients after false discovery correction were extracted. Afterward, based on the extracted FC features, machine-learning algorithms such as support vector machine, logistic regression, decision tree, and random forest were established to recognize DNR. The machine learning experiment procedure mainly included three following steps: feature standardization, parameter adjustment, and performance comparison. Finally, independent testing was conducted to validate the established prediction model. The algorithm performance was evaluated by a permutation test.

Results: We found significantly decreased DMN connectivity with the brain regions involved in visual processing in DNR patients than in non-DNR patients. The best result was obtained from the random forest algorithm based on the 20 decision trees (estimators). The random forest model achieved the accuracy, sensitivity, and specificity of 84.0, 63.1, and 89.5%, respectively. The area under the receiver operating characteristic curve of the classifier reached 86.4%. The feature that contributed the most to the random forest model was the FC between the left retrosplenial cortex/posterior cingulate cortex and left precuneus.

Conclusion: The decreased FC of DMN with regions involved in visual processing might be effective markers for the prediction of DNR and could provide new insights into the neural mechanisms of DNR.

Clinical Trial Registration: : Chinese Clinical Trial Registry, ChiCTR-DCD-15006096.

1. Introduction

Delayed neurocognitive recovery (DNR) is a common complication of the neurological system, whose incidence is about 25.8% after non-cardiac surgery in elderly patients (Moller et al., 1998; Evered et al., 2018). With the development of an aging society, more and more geriatric surgical patients suffer from adverse outcomes caused by cognitive dysfunction after surgery, such as prolonged recovery time and increased incidence of complications (Evered et al., 2017; Chen and Qin, 2021). Hence, it is of crucial importance to understand the underlying pathological mechanisms of DNR and identify the potential biomarkers of DNR for its prevention.

In most brain function studies, resting-state functional magnetic resonance imaging (rs-fMRI) served as a useful approach to describe the connection characteristics of the brain networks associated with certain disease states (Rubinov and Sporns, 2010; Wu et al., 2021). Of particular interest is the default mode network (DMN), which is the most representative, stable, and individual-identifiable network in the resting state, as well as the most studied one in Alzheimer’s disease (AD) and mild cognitive impairment (MCI; Li et al., 2022). The DMN is not only engaged in the integration of cognitive, but also in memory recollection and reconstruction of long-term memories during retrieval (Krajcovicova et al., 2017). The medial prefrontal cortex (mPFC) and posterior cingulate cortex (PCC) are the core hubs of DMN. Besides, the lateral parietal cortex, middle cingulate cortex (MCC), and hippocampal gyrus are also components of DMN. These specific regions exhibit highly organized resting functional activity. However, most of the existing studies exclusively focused on the FC of mPFC and PCC, thus overlooking a lot of information and details about the DMN whole connectivity patterns. Hence, in the present study, we proposed a feasible approach to evaluate changes in whole-brain connectivity across all regions of DMN in patients with and without DNR by defining the ROI mask from the DMN atlas.

Machine learning can incorporate diversified rs-fMRI variables and identify potential biomarkers to build corresponding classification models (Solomon et al., 2020). As the information obtained from rs-fMRI data may provide objective evidence for the classification of DNR, combining it with machine learning algorithms may help improve the accuracy of early prediction of DNR. Rs-fMRI data mining has various stages, including rs-fMRI data preprocessing and analysis, extracted rs-fMRI features standardization (data preparation), model selection (machine-learning algorithms), the training and evaluation of models, parameter adjustment, and finally prediction (Sayadi et al., 2022).

In the present study, we used all subregions of DMN as the ROIs to determine whether DNR patients show abnormal FC between DMN and other voxels of the whole brain before surgery. Second, we adopted the support vector machine (SVM), logistic regression, decision tree, and random forest algorithms to develop machine-learning classification models to identify DNR patients based on rs-fMRI data. We hypothesized that whole-brain FC changes of DMN could serve as effective biomarkers to predict patients with DNR, which might help us to better understand the neural mechanisms of DNR and provide support for future clinical early screening and intervention in DNR patients.

2. Materials and methods

2.1. Ethics approval and participants

This nested case–control study was approved by the Ethics Committee of Huadong Hospital Affiliated to Fudan University (approval number: 20170020). The trial was registered on the Chinese Clinical Trial Registry1 with the identification number of ChiCTR-DCD-15006096. The principal investigator was Weidong Gu, and the date of registration was 16 March 2015.

Subjects were recruited at the Huadong Hospital Affiliated to Fudan University between September 2017 and February 2019. Written informed consent was obtained from all subjects participating in the trial before enrollment. The inclusion and exclusion criteria of our study population have been previously described (Jiang et al., 2020, 2021). The inclusion criteria were as follows: (1) patients scheduled to undergo non-cardiac surgery; (2) age ≥60 years old; (3) American Society of Anesthesiologists (ASA) classification I-III; and (4) right-handedness. The exclusion criteria were as follows: (1) education duration <6 years; (2) baseline mini-mental state examination (MMSE) score <24 points; (3) preexisting mental or psychiatric disease, cardiac or central nervous system vascular disease, Parkinson’s disease, cardiac or cranial surgery history; (4) major surgery within the past 12 months; (5) taking sedatives or antidepressants during the nearest year; (6) alcohol or drug abuse; and (7) vision and audition impairment or speech issues impeding communication.

2.2. Neurocognitive assessment

The Z-score method recommended by the International study of postoperative cognitive dysfunction (IPSOCD1) was used to identify the DNR (Moller et al., 1998). The comprehensive neurocognitive tests were conducted before surgery and at 7–14 days after the surgery for each subject. The neurocognitive tests consisted of MMSE, digit symbol substitution test (DSST), digit span forward and backward test (DSF/DSB), trail-making test-part A (TMT-A), and verbal fluency test (VFT). The Z score for each test and the composite Z score were calculated (Jiang et al., 2021). A patient was diagnosed with DNR when the Z scores of at least two of the neurocognitive tests or the composite Z score were ≥1.96.

2.3. Rs-fMRI data acquisition and preprocessing

Siemens Skyra 3.0 T MRI scanner was used to obtain all MRI data before surgery. The complete MRI acquisition protocol included three-dimensional (3D) anatomical T1-weighted imaging and fMRI echo-planar imaging. The 3D anatomical T1-weighted imaging parameters were as follows: 176 sagittal slices, repetition time = 1,900 ms, echo time = 3.57 ms, voxel size = 1 × 1 × 1 mm, and flip angle = 9°. The echo-planar imaging sequence parameters were as follows: 33 axial slices, slices thickness = 4 mm with a 0-mm gap, repetition time = 3,000 ms, echo time = 30 ms, voxel size = 3.4 × 3.4 × 4 mm, and flip angle = 90°. During the fMRI imaging, 120 volumes were obtained that lasted 8.5 min.

All rs-fMRI data were processed using Statistical Parametric Mapping version 12 and RESTplus version 1.24, based on MATLAB version R2013b. The preprocessing of rs-fMRI data consisted of eight following steps, where the first step included discarding the first five volumes to avoid the potential noise related to the participants’ adaptation to the scanner. Then, the remaining volumes were subjected to the slice-timing correction, head motion correction, spatial normalization, smoothing with a 6 × 6 × 6 mm Gaussian kernel, low-frequency filtering (0.01–0.08 Hz), linear trend of time course removal, and nuisance covariates regression (motion artifact, white matter signal, and cerebrospinal fluid signal; Jiang et al., 2021).

2.4. Functional connectivity analysis of DMN

A seed-based FC analysis was performed to explore the preoperative whole-brain voxel-wise FC alteration of the DMN in DNR patients. Nine regions in the dorsal DMN and nine regions in the ventral DMN were used as ROIs based on the DMN atlas from the Functional Imaging in Neuropsychiatric Disorders Lab at Stanford University (Shirer et al., 2012). Nine ROIs in the dorsal DMN included mPFC/anterior cingulate cortex, right superior frontal gyrus (SFG), PCC, MCC, left and right angular gyrus, left and right hippocampus, and thalamus. Nine ROIs in the ventral DMN included left and right retrosplenial cortex (RSC)/PCC, left and right middle frontal gyrus, left and right parahippocampal gyrus, left and right middle occipital gyrus, and precuneus.

The time series of each region were averaged and correlated with every other voxel within the gray-matter mask. The FC correlation maps were converted using Fisher’s r-to-z transform. A two-sample t-test was performed to explore the differences in seed-based FC between the DNR patients and non-DNR patients. The statistical maps were thresholded at p < 0.001 at the voxel level and p < 0.05 at the cluster level by false discovery rate (FDR) corrected for age, sex, and education duration (Chumbley and Friston, 2009).

2.5. Machine learning modeling

The prediction models of DNR were established using the SVM, logistic regression, decision tree, and random forest algorithms based on the FC features. Considering that the eigenvalues of FC features differ greatly, Scikit-learn’s “Preprocessing.StandardScaler” function was used to convert FC features into standard normally distributed data with zero mean and unit variance (Cho et al., 2022). The dataset was randomly divided into a training and a testing set at a ratio of 7:3. To reduce selection bias or overfitting, a 10-fold cross-validation method was employed for internal validation in the training set. In addition, because of the minority of DNR patients in the data set, the “balanced” class weight mode in the machine-learning algorithms was conducted to automatically adjust weights inversely proportional to class frequencies in the input data.

We further adopted the grid search method to identify the optimal parameters of the established prediction models. The main parameter to be adjusted in the LinearSVM algorithm was the penalty coefficient C (Wang L. et al., 2021). The logistic regression algorithm adopted the L2 penalty term, and adjusted the parameters including the inverse of regularization strength C and the maximum number of iterations (max_iter). The maximum depth of the tree (max_depth), the minimum number of samples required to split an internal node (min_samples_split), and the minimum number of samples required to be at a leaf node (min_samples leaf) were the most important parameters of decision tree algorithm. The random forest algorithm was an extension of the bagging method that fits a number of decision tree classifiers on various sub-samples of the dataset and uses averaging to improve the predictive accuracy and control overfitting. In addition to the parameters of the decision tree algorithm, the number of trees in the forest (n_estimators) was also a parameter that needed to be adjusted in the random forest classifier. The machine-learning algorithms can provide the corresponding weight of each variable, thereby identifying the variables significantly influencing the classification model.

The mean and standard deviation (SD) of the classification accuracy, sensitivity, specificity, and area under the receiver operating characteristic curve (AUC) of each model were computed in 100 runs of repeated re-divided dataset randomly. For further model assessment, a permutation test was performed to assess the statistical significance of the classification accuracy. The permutation test was repeated 1,000 times (Tonga and Huangb, 2022). During each time, the classifier reallocated labels of DNR patients and non-DNR patients to the training data randomly and repeated the whole classification process. The p value was acquired after the entire permutation was finished. Finally, Welch’s analysis of variance (ANOVA) and Post-hoc analysis (Games–Howell test with a Tukey correction) were conducted to compare the classification accuracies of different models. The statistical analyses of machine learning algorithms were performed in Python.

2.6. Statistical analysis

IBM SPSS 22.0 and GraphPad Prism were used to perform demographical statistical analyses. Continuous variables were assessed with a two-sample t-test or Mann–Whitney test and were presented as mean ± SD or median (interquartile range, IQR). Categorical variables were assessed with Chi-squared or Fisher exact test. For neurocognitive tests, analysis of covariance (ANCOVA) was used to compare follow-up scores between the two groups, with baseline scores as covariance. In addition, the paired-sample t-test analysis was used to compare the tests before and after surgery in the DNR patients and non-DNR patients. p < 0.05 indicated statistical significance.

3. Results

3.1. Subject characteristics

The flowchart of patient enrollment is provided in Supplementary Figure 1. In the present study, 74 patients completed the rs-fMRI scan and the neurocognitive tests. A total of 16 patients were diagnosed with DNR at 7–14 days postoperatively. Also, the duration of education was significantly shorter in the DNR patients [median (IQR), 6 (6, 9)] compared to the non-DNR patients [median (IQR), 9 (9, 12)] (p = 0.002). No significant differences were detected in the other baseline characteristics, including age, sex, height, weight, body mass index, surgical history, and comorbidities between the two groups. In addition, there were no significant differences in the surgical duration and nature of surgery between the groups (Table 1).

TABLE 1
www.frontiersin.org

Table 1. Baseline characteristics.

3.2. Neurocognitive assessments

We investigated group differences in follow-up scores of neurocognitive tests using ANCOVA, wherein we controlled for baseline scores. Our results showed that the DNR patients exhibited significantly lower MMSE, VFT, DSF/DSB, and DSST follow-up scores and higher TMT-A follow-up scores (seconds) compared to the non-DNR patients (all p < 0.05, Table 2). Paired-samples t-test analysis revealed that the performance of postoperative DSST and TMT-A assessment was the opposite in the two groups. For DNR patients, the DSST score after surgery was significantly worse than before, and the TMT-A time after surgery was significantly longer than before surgery. However, non-DNR patients showed improved performance in DSST and TMT-A post-surgery compared with pre-surgery. These results indicated that the DNR patients had poor performance after non-cardiac surgery in many domains of cognitive function, especially in those assessed by DSST and TMT-A, such as attention, visual–spatial ability, visual scanning, visual search speed, executive function, and visual motor coordination.

TABLE 2
www.frontiersin.org

Table 2. Neurocognitive assessments.

3.3. Altered seed-based FC patterns of DMN

Seed-based FC analysis revealed that the preoperative whole-brain resting-state connections of DMN core regions were significantly weakened in the DNR patients (Figure 1; Table 3). In the dorsal DMN, altered whole-brain FC was observed only in the right SFG. The DNR group exhibited decreased FC of the right SFG to the right calcarine and right middle frontal gyrus compared to the non-DNR group. Regarding the altered patterns of the ventral DMN, decreased connections were detected in the DNR group compared to the non-DNR group between the left parahippocampal gyrus and right cuneus. Also, the DNR group exhibited decreased FC of bilateral RSC/PCC to the right superior occipital gyrus/calcarine and left precuneus. In addition, decreased FC of the right RSC/PCC to the right precuneus was also found in the DNR patients. The results indicated that the key ROIs of DMN in DNR patients had altered resting-state functional connections with visual processing-related cortical areas, which mainly manifested as weakened connectivity. Age, sex, and education duration were included as covariates in the seed-based whole-brain FC analysis. The results were corrected using the cluster-based FDR (uncorrected voxel p < 0.001 and corrected cluster p < 0.05).

FIGURE 1
www.frontiersin.org

Figure 1. Altered seed-based whole-brain FC of DMN. (A) Brain regions with decreased connectivity to DMN core ROIs in the DNR patients compared to the non-DNR patients after adjusting for age, sex, and education duration (cluster p < 0.05, FDR-corrected). (B) Visualized nodes with peak T value of such brain regions. CAL, calcarine; CUN, cuneus; DMN, default mode network; DNR, delayed neurocognitive recovery; FC, functional connectivity; FDR, false discovery rate; MFG, middle frontal gyrus; PCC, posterior cingulate cortex; PCUN, precuneus; PHG, parahippocampal gyrus; ROIs, regions of interest; RSC, retrosplenial cortex; SFG, superior frontal gyrus; SOG, superior occipital gyrus; A, anterior; P, posterior; L, left; R, right.

TABLE 3
www.frontiersin.org

Table 3. Regions with decreased connectivity to DMN core ROIs in the DNR patients.

3.4. Machine learning prediction models

The resting-state FC data were standardized before training to ensure that machine learning estimators could correctly model all features. The mean and variance of each feature were converted to 0.00 and 1.00 by standardization. The data distribution before and after preprocessing is shown in Supplementary Figure 2. In order to develop the prediction model, the dataset was randomly divided into a training set containing 70% of the samples and a testing set containing the remaining 30%, whereas the models were validated using the 10-fold cross-validation method.

By the grid search method, the optimal penalty coefficient C of the SVM classifier was identified as 8.88889. The inverse of regularization strength C of 2 and max_iter of 10 were the optimal parameters of the logistic regression classifier. The optimal parameters of the decision tree classifier were max_depth of 9, min_samples_split of 2, and min_samples leaf of 2. And the optimal parameters of the random forest classifier were max_depth of 2, min_samples_split of 5, min_samples leaf of 2, and n_estimators of 20.

Applying the 100 runs of repeated re-divided dataset randomly, the SVM classifier with optimal parameters achieved a classification accuracy of 75.3 ± 9.4%. The logistic regression classifier with optimal parameters achieved a classification accuracy of 78.1 ± 7.8%. The decision tree classifier with optimal parameters achieved a classification accuracy of 77.2 ± 7.8%. The random forest classifier with optimal parameters achieved a classification accuracy of 84.0 ± 6.5%. The AUC of all classification models is shown in Figure 2. The mean values and standard deviation of AUC, accuracy, sensitivity, and specificity of all classification models are summarized in Figure 3.

FIGURE 2
www.frontiersin.org

Figure 2. The ROC curves of all classification models in 100 runs of repeated re-divided dataset randomly. (A) The SVM model; (B) the logistic regression model; (C) the decision tree model; (D) the random forest model. AUC, an area under the curve; ROC, the receiver operating characteristic; SVM, support vector machine.

FIGURE 3
www.frontiersin.org

Figure 3. The classification performance of models. The mean AUC, accuracy, sensitivity, and specificity were calculated based on 100 runs of repeated re-divided dataset randomly. Abbreviations: AUC, an area under the curve; SVM, support vector machine; LR, logistic regression; DT, decision tree; RF, random forest.

The permutation test can be seen that the SVM model (p = 0.022), the logistic regression model (p = 0.023), and the random forest model (p = 0.019) could distinguish the DNR group and the non-DNR group, whereas the decision tree model could not (p = 0.118). Furthermore, Welch’s ANOVA was performed to compare the classification accuracies of different models. It revealed significant differences among the SVM, logistic regression, and random forest models in the classification accuracies (Welch F = 34.002, p < 0.001). Post-hoc tests showed that the accuracy of the random forest model was higher than the SVM model (p < 0.001, 95% CI 0.060 to 0.114) and the logistic regression model (p < 0.001, 95% CI 0.035 to 0.084), while there was no significant difference in the accuracy between the SVM model and the logistic regression model (p = 0.068, 95% CI −0.002 to 0.056). The results indicated that the random forest classifier with optimal parameters could achieve the best prediction performance for DNR. The random forest algorithm also provided the corresponding weight of each feature, thereby identifying the features that influence the predictions (Corradi et al., 2018). Figure 4 shows the corresponding weights of each variable calculated by the random forest algorithm to reveal the contributions of these variables to the model.

FIGURE 4
www.frontiersin.org

Figure 4. The feature weights of the random forest prediction models.

4. Discussion

The results of the present study showed that the preoperative resting-state whole-brain FC of DMN features could be used to predict DNR by establishing a machine-learning model. The key findings are as follows: (1) the FC between the key ROIs of DMN and other brain regions, especially in visual processing-related cortical areas, in DNR patients before surgery was more weakened compared to non-DNR patients; (2) the random forest machine learning model showed that whole-brain FC features of DMN could be used to predict postoperative DNR with high accuracy.

The human brain can be viewed as an integrated network consisting of brain regions that are anatomically separated but functionally connected (van den Heuvel and Hulshoff Pol, 2010). Hence, compared with the activity of local brain regions, FC analysis can better reflect the actual operation mode of the whole brain. The FC approach has obtained good results in fMRI data classification studies (Wang et al., 2019; Dai et al., 2022). Recently, an increasing number of rs-fMRI studies have reported that resting-state FC before surgery has a primary role in predicting cognitive dysfunction after surgery. Our previous study found that patients with preoperative decreased FC between the bilateral MCC and left calcarine were more susceptible to DNR following non-cardiac surgery by logistic regression (Jiang et al., 2020). By exploring the whole-brain FC of key hubs in cognitive-related brain networks, we found that the altered whole-brain FC of DMN and central executive network could predict DNR in elderly patients using machine learning algorithms (Jiang et al., 2021). In addition, Wu et al. (2021) established a machine learning DNR classification model based on the multi-order brain FC network features, achieving good prediction performance. These rs-fMRI studies examined preoperative FC changes of brain networks in patients with cognitive decline after surgery, revealing some consistent findings, such as the importance of DMN connectivity.

The basal activity of the DMN is quite consistent during the resting state. The DMN has an essential role in episodic memory, processing of cognition and emotion, and so on (Shi et al., 2020). Functional abnormalities of DMN have been well-documented in various diseases related to cognitive impairment, such as AD and MCI (Zhang et al., 2020; Wang S. M. et al., 2021). Our preliminary findings suggested that the whole-brain FC of certain brain regions in the DMN was altered in DNR patients. Therefore, in the present study, we further explored the FC patterns of all components of DMN in DNR patients based on prior knowledge. We also found that the DMN connectivity was significantly decreased in DNR patients. Interestingly, the brain regions with altered DMN connectivity were predominantly in visual processing areas. Similar changes were detected in AD and MCI in previous studies. Krajcovicova et al. (2017) found that the patients with MCI/AD had disrupted PCC connectivity with areas of the ventral visual pathway, which suggested that the abnormal DMN connectivity patterns may provide a possible underlying mechanism of impaired visual processing in MCI/AD.

Consistent with cognitive disorders like AD and MCI, and in addition to memory deficits, impairment in other cognitive domains is present in DNR. Many studies used neuropsychological tests to demonstrate visual cognitive dysfunction in DNR patients (Lage et al., 2021; Ji and Li, 2022; Zhao et al., 2022). We also found that the DNR patients showed poor performance postoperatively in visual attention, visual memory, visual–spatial scanning, and visual search speed assessed by DSST and TMT-A. However, the current diagnosis methods of DNR, which are mainly based on assessments of neuropsychological tests, may easily be impacted by human factors. As the information obtained from rs-fMRI data could provide objective evidence for the brain activity characteristics of DNR, combining it with the manual assessment of DNR diagnosis may reduce the diagnosis errors of DNR and identify the potential biomarkers. Using different FC analysis methods, we demonstrated that the DNR patients showed abnormal connectivity in regions related to visual processing, especially decreased connection to the key regions in DMN. The altered MCC spontaneous activity in DNR patients was detected based on the Amplitude of Low-frequency Fluctuation (ALFF) analysis. In our previous study, MCC, a subregion of DMN, showed a significantly decreased FC to the calcarine in DNR patients by ALFF-based FC analysis (Jiang et al., 2020). In addition, we found that the other key region of DMN, the lateral parietal cortex, also exhibited a significantly decreased FC to the calcarine in DNR patients by seed-based FC analysis (Jiang et al., 2021). The calcarine lies in the center of the primary visual cortex and contributes to visual processing, such as visual attention and visual memory (Pratte and Tong, 2014; Bergmann et al., 2016). A large body of studies has suggested that abnormal calcarine FC might be associated with the occurrence and development of multiple cognitive disorders (Moon and Jeong, 2017; Zuo et al., 2017; Zhuang et al., 2019). Among the brain regions with abnormal FC to DMN identified in the present study, in addition to the calcarine, the cuneus, precuneus, and superior occipital gyrus, which are involved in visual processing, were identified. We hypothesized that the visual cognitive impairment in DNR patients was mainly caused by the inhibition of FC between the DMN and visual processing areas. The association of DMN and visual processing could be worthy of further study in future DNR research.

In recent years, machine learning has been widely used by cognitive dysfunction disease prediction studies due to its individualized classification ability (Mufti et al., 2019; Hu et al., 2021). The combination of brain network research based on fMRI and machine learning has become a topic of great interest. Yet, the data obtained by fMRI are highly dimensional, sample-limited, and non-linear. Traditional modeling methods or linear dimensionality reduction cannot meet the requirements of the model (Tan et al., 2022). Machine-learning algorithms such as SVM and random forest had advantages in dealing with small samples and non-linear features. In addition, since the fMRI original data contains a large amount of redundant information, selecting the proper feature subset is essential to avoid reducing the model’s performance (Shuai et al., 2017). The present study uses the two-sample t-test after FDR correction to extract features with intergroup differences. Based on the FC features of DMN selected by feature screening, we built a prediction DNR model using the random forest algorithm, achieving great performance.

There are some limitations in the present study. First, the postoperative FC patterns of DMN in DNR patients could not be identified because most subjects refused to perform another fMRI scan postoperatively. Second, the amount of experimental data is limited. A sufficient amount of data ensures the reliability and stability of machine learning. Although we have conducted a strict division of training and testing data and performed the SVM and random forest algorithm, which are suitable for small samples, the stability and reliability of the prediction models still need to be verified on more datasets, especially multi-center data.

In conclusion, the present study demonstrated that the decreased FC between the DMN and brain regions involved in visual processing was detected in elderly patients who developed DNR after non-cardiac surgery. Also, the characteristic FC patterns could be effective biomarkers for the prediction of DNR.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving human participants were reviewed and approved by Ethics Committee of Huadong Hospital Affiliated to Fudan University. The patients/participants provided their written informed consent to participate in this study.

Author contributions

ZJ and YC made substantial contributions to the design, acquisition, analysis, and interpretation of study data, and drafted the manuscript. SLiu and PY made substantial contributions to the acquisition and analysis of study data, and provided critical revisions to the manuscript. YY made substantial contributions to the analysis and interpretation of study data, and provided critical revisions to the manuscript. GL made substantial contributions to the design and interpretation of study data, and provided critical revisions to the manuscript. SLi, YX, and YZ made substantial contributions to the acquisition of study data, and provided critical revisions to the manuscript. ZB made substantial contributions to the design of the study, and provided critical revisions to the manuscript. SN and WG made substantial contributions to the conception, design, acquisition, analysis, and interpretation of study data, drafted the manuscript, and provided critical revisions to the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by the project of National Natural Science Foundation of China (82271286), Science and Technology Commission of Shanghai Municipality (20Y11900200), the Shanghai Municipal Health Commission (2020YJZX0119), Hospital research project of Huadong Hospital (2019jc016), the National Key R&D Program of China (2018YFC2002000), and the Shanghai Municipal Key Clinical Specialty (shslczdzk02801).

Acknowledgments

The authors would like to acknowledge Zhichao Jin (Department of Health Statistics, Second Military Medical University, Shanghai, China) for his assistance in statistical consultation.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2022.1109485/full#supplementary-material

Footnotes

References

Bergmann, J., Genç, E., Kohler, A., Singer, W., and Pearson, J. (2016). Neural anatomy of primary visual cortex limits visual working memory. Cereb. Cortex 26, 43–50. doi: 10.1093/cercor/bhu168

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., and Qin, J. (2021). Modified frailty index independently predicts postoperative delirium and delayed neurocognitive recovery after elective total joint Arthroplasty. J. Arthroplast. 36, 449–453. doi: 10.1016/j.arth.2020.07.074

PubMed Abstract | CrossRef Full Text | Google Scholar

Cho, C., Kwon, K., and Wu, C. (2022). On weather data-based prediction of gamma exposure rates using gradient boosting learning for environmental radiation monitoring. Sensors (Basel) 22:7062. doi: 10.3390/s22187062

PubMed Abstract | CrossRef Full Text | Google Scholar

Chumbley, J. R., and Friston, K. J. (2009). False discovery rate revisited: FDR and topological inference using Gaussian random fields. NeuroImage 44, 62–70. doi: 10.1016/j.neuroimage.2008.05.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Corradi, J. P., Thompson, S., Mather, J. F., Waszynski, C. M., and Dicks, R. S. (2018). Prediction of incident delirium using a random forest classifier. J. Med. Syst. 42:261. doi: 10.1007/s10916-018-1109-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, P., Xiong, T., Zhou, X., Ou, Y., Li, Y., Kui, X., et al. (2022). The alterations of brain functional connectivity networks in major depressive disorder detected by machine learning through multisite rs-fMRI data. Behav. Brain Res. 435:114058. doi: 10.1016/j.bbr.2022.114058

PubMed Abstract | CrossRef Full Text | Google Scholar

Evered, L., Scott, D. A., and Silbert, B. (2017). Cognitive decline associated with anesthesia and surgery in the elderly: does this contribute to dementia prevalence? Curr. Opin. Psychiatry 30, 220–226. doi: 10.1097/yco.0000000000000321

CrossRef Full Text | Google Scholar

Evered, L., Silbert, B., Knopman, D. S., Scott, D. A., DeKosky, S. T., Rasmussen, L. S., et al. (2018). Recommendations for the nomenclature of cognitive change associated with anaesthesia and surgery-2018. Br. J. Anaesth. 121, 1005–1012. doi: 10.1016/j.bja.2017.11.087

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, M., Shu, X., Yu, G., Wu, X., Välimäki, M., and Feng, H. (2021). A risk prediction model based on machine learning for cognitive impairment among Chinese community-dwelling elderly people with normal cognition: development and validation study. J. Med. Internet Res. 23:e20298. doi: 10.2196/20298

PubMed Abstract | CrossRef Full Text | Google Scholar

Ji, L., and Li, F. (2022). Potential markers of neurocognitive disorders after cardiac surgery: a bibliometric and visual analysis. Front. Aging Neurosci. 14:868158. doi: 10.3389/fnagi.2022.868158

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, Z., Cai, Y., Zhang, X., Lv, Y., Zhang, M., Li, S., et al. (2021). Predicting delayed neurocognitive recovery after non-cardiac surgery using resting-state brain network patterns combined with machine learning. Front. Aging Neurosci. 13:715517. doi: 10.3389/fnagi.2021.715517

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, Z., Zhang, X., Lv, Y., Zheng, X., Zhang, H., Zhang, X., et al. (2020). Preoperative altered spontaneous brain activity and functional connectivity were independent risk factors for delayed neurocognitive recovery in older adults undergoing noncardiac surgery. Neural Plast. 2020, 9796419–9796410. doi: 10.1155/2020/9796419

PubMed Abstract | CrossRef Full Text | Google Scholar

Krajcovicova, L., Barton, M., Elfmarkova-Nemcova, N., Mikl, M., Marecek, R., and Rektorova, I. (2017). Changes in connectivity of the posterior default network node during visual processing in mild cognitive impairment: staged decline between normal aging and Alzheimer's disease. J. Neural Transm. 124, 1607–1619. doi: 10.1007/s00702-017-1789-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Lage, C., González-Suárez, A., Alcalde-Hierro, M. P., Sampedro-González, M. I., Villanueva-Eguaras, M., Sánchez-Crespo, M. R., et al. (2021). Major surgery affects memory in individuals with cerebral amyloid-β pathology. J. Alzheimers Dis. 79, 863–874. doi: 10.3233/jad-191229

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Liu, A., Fu, X., McKeown, M. J., Wang, Z. J., and Chen, X. (2022). Atlas-guided parcellation: individualized functionally-homogenous parcellation in cerebral cortex. Comput. Biol. Med. 150:106078. doi: 10.1016/j.compbiomed.2022.106078

PubMed Abstract | CrossRef Full Text | Google Scholar

Moller, J. T., Cluitmans, P., Rasmussen, L. S., Houx, P., Rasmussen, H., Canet, J., et al. (1998). Long-term postoperative cognitive dysfunction in the elderly ISPOCD1 study. ISPOCD investigators. International study of post-operative cognitive dysfunction. Lancet 351, 857–861. doi: 10.1016/s0140-6736(97)07382-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Moon, C. M., and Jeong, G. W. (2017). Functional and morphological alterations associated with working memory dysfunction in patients with generalized anxiety disorder. Acta Radiol. 58, 344–352. doi: 10.1177/0284185116649794

PubMed Abstract | CrossRef Full Text | Google Scholar

Mufti, H. N., Hirsch, G. M., Abidi, S. R., and Abidi, S. S. R. (2019). Exploiting machine learning algorithms and methods for the prediction of agitated delirium after cardiac surgery: models development and validation study. JMIR Med. Inform. 7:e14993. doi: 10.2196/14993

PubMed Abstract | CrossRef Full Text | Google Scholar

Pratte, M. S., and Tong, F. (2014). Spatial specificity of working memory representations in the early visual cortex. J. Vis. 14:22. doi: 10.1167/14.3.22

PubMed Abstract | CrossRef Full Text | Google Scholar

Rubinov, M., and Sporns, O. (2010). Complex network measures of brain connectivity: uses and interpretations. NeuroImage 52, 1059–1069. doi: 10.1016/j.neuroimage.2009.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Sayadi, M., Varadarajan, V., Sadoughi, F., Chopannejad, S., and Langarizadeh, M. (2022). A machine learning model for detection of coronary artery disease using noninvasive clinical parameters. Life (Basel) 12:1933. doi: 10.3390/life12111933

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, Y., Li, J., Feng, Z., Xie, H., Duan, J., Chen, F., et al. (2020). Abnormal functional connectivity strength in first-episode, drug-naïve adult patients with major depressive disorder. Prog. Neuro-Psychopharmacol. Biol. Psychiatry 97:109759. doi: 10.1016/j.pnpbp.2019.109759

PubMed Abstract | CrossRef Full Text | Google Scholar

Shirer, W. R., Ryali, S., Rykhlevskaia, E., Menon, V., and Greicius, M. D. (2012). Decoding subject-driven cognitive states with whole-brain connectivity patterns. Cereb. Cortex 22, 158–165. doi: 10.1093/cercor/bhr099

PubMed Abstract | CrossRef Full Text | Google Scholar

Shuai, M., Changle, Z., Na, G., Yan, W., Yanwu, Y., Xin, G., et al. (2017). A study of feature extraction for Alzheimer's disease based on resting-state fMRI. Annu. Int. Conf. IEEE Eng. Med. Biol. Soc. 2017, 517–520. doi: 10.1109/embc.2017.8036875

CrossRef Full Text | Google Scholar

Solomon, S. C., Saxena, R. C., Neradilek, M. B., Hau, V., Fong, C. T., Lang, J. D., et al. (2020). Forecasting a crisis: machine-learning models predict occurrence of intraoperative Bradycardia associated with hypotension. Anesth. Analg. 130, 1201–1210. doi: 10.1213/ane.0000000000004636

CrossRef Full Text | Google Scholar

Tan, G., Huang, B., Cui, Z., Dou, H., Zheng, S., and Zhou, T. (2022). A noise-immune reinforcement learning method for early diagnosis of neuropsychiatric systemic lupus erythematosus. Math. Biosci. Eng. 19, 2219–2239. doi: 10.3934/mbe.2022104

PubMed Abstract | CrossRef Full Text | Google Scholar

Tonga, Y., and Huangb, X. (2022). Altered voxel-wise degree centrality associated with patients with comitant exotropia. Neuroreport 33, 526–533. doi: 10.1097/wnr.0000000000001813

PubMed Abstract | CrossRef Full Text | Google Scholar

van den Heuvel, M. P., and Hulshoff Pol, H. E. (2010). Exploring the brain network: a review on resting-state fMRI functional connectivity. Eur. Neuropsychopharmacol. 20, 519–534. doi: 10.1016/j.euroneuro.2010.03.008

CrossRef Full Text | Google Scholar

Wang, S. M., Kim, N. Y., Um, Y. H., Kang, D. W., Na, H. R., Lee, C. U., et al. (2021). Default mode network dissociation linking cerebral beta amyloid retention and depression in cognitively normal older adults. Neuropsychopharmacology 46, 2180–2187. doi: 10.1038/s41386-021-01072-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, C., Xiao, Z., and Wu, J. (2019). Functional connectivity-based classification of autism and control using SVM-RFECV on rs-fMRI data. Phys. Med. 65, 99–105. doi: 10.1016/j.ejmp.2019.08.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Zhang, J., Shan, G., Liang, J., Jin, W., Li, Y., et al. (2021). Establishment of a lung cancer discriminative model based on an optimized support vector machine algorithm and study of key targets of Wogonin in lung cancer. Front. Pharmacol. 12:728937. doi: 10.3389/fphar.2021.728937

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, G., Jiang, Z., Cai, Y., Zhang, X., Lv, Y., Li, S., et al. (2021). Multi-order brain functional connectivity network-based machine learning method for recognition of delayed neurocognitive recovery in older adults undergoing non-cardiac surgery. Front. Neurosci. 15:707944. doi: 10.3389/fnins.2021.707944

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X. Y., Wang, Y. F., Zheng, L. J., Zhang, H., Lin, L., Lu, G. M., et al. (2020). Impacts of AD-related ABCA7 and CLU variants on default mode network connectivity in healthy middle-age adults. Front. Mol. Neurosci. 13:145. doi: 10.3389/fnmol.2020.00145

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, H., Han, Q., Shi, C., and Feng, Y. (2022). The effect of opioid-sparing anesthesia regimen on short-term cognitive function after thoracoscopic surgery: a prospective cohort study. Perioper. Med. 11:45. doi: 10.1186/s13741-022-00278-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhuang, L., Liu, X., Shi, Y., Liu, X., and Luo, B. (2019). Genetic variants of PICALM rs541458 modulate brain spontaneous activity in older adults with amnestic mild cognitive impairment. Front. Neurol. 10:494. doi: 10.3389/fneur.2019.00494

PubMed Abstract | CrossRef Full Text | Google Scholar

Zuo, L., Wang, S., Yuan, J., Gu, H., Zhou, Y., and Jiang, T. (2017). Altered intra-and interregional synchronization in the absence of the corpus callosum: a resting-state fMRI study. Neurol. Sci. 38, 1279–1286. doi: 10.1007/s10072-017-2953-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: whole-brain functional connectivity, default mode network, delayed neurocognitive recovery, machine learning, visual processing

Citation: Jiang Z, Cai Y, Liu S, Ye P, Yang Y, Lin G, Li S, Xu Y, Zheng Y, Bao Z, Nie S and Gu W (2023) Decreased default mode network functional connectivity with visual processing regions as potential biomarkers for delayed neurocognitive recovery: A resting-state fMRI study and machine-learning analysis. Front. Aging Neurosci. 14:1109485. doi: 10.3389/fnagi.2022.1109485

Received: 27 November 2022; Accepted: 16 December 2022;
Published: 06 January 2023.

Edited by:

Hong Hai Zhang, Zhejiang University School of Medicine, China

Reviewed by:

Yao Lu, First Affiliated Hospital of Anhui Medical University, China
Weidong Yao, Wannan Medical College, China
Kai Mo, Southern Medical University, China

Copyright © 2023 Jiang, Cai, Liu, Ye, Yang, Lin, Li, Xu, Zheng, Bao, Nie and Gu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Weidong Gu, hdmz0800@163.com; Shengdong Nie, nsd4647@163.com

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.