Identification of Cognitive Dysfunction in Patients with T2DM Using Whole Brain Functional Connectivity

Zhenyu Liu,Jiangang Liu,Huijuan Yuan,Taiyuan Liu,Xingwei Cui,Zhenchao Tang,Yang Du,Meiyun Wang*,Yusong Lin*,Jie Tian

1CAS Key Laboratory of Molecular Imaging,Institute of Automation,Chinese Academy of Sciences,Beijing 100190,China

2School of Computer and Information Technology,Beijing Jiaotong University,Beijing 100044,China

3Department of Endocrinology and Metabolism,Henan Provincial People's Hospital&the People's Hospital of Zhengzhou University,Zhengzhou 450003,China

4Department of Radiology,Henan Provincial People's Hospital&the People's Hospital of Zhengzhou University,Zhengzhou 450003,China

5Cooperative Innovation Center for Internet Healthcare&School of Software,Zhengzhou University,Zhengzhou 450003,China

6School of Mechanical,Electrical&Information Engineering,Shandong University(Weihai),Weihai 264209,China

7University of Chinese Academy of Sciences,Beijing 100080,China

8Beijing Advanced Innovation Center for Big Data-Based Precision Medicine,School of Medicine,Beihang University,Beijing 100191,China

9Engineering Research Center of Molecular and Neuro Imaging of Ministry of Education,School of Life Science and Technology,Xidian University,Xi'an 710126,China

Abstract Majority of type 2 diabetes mellitus(T2DM)patients are highly susceptible to several forms of cognitive impairments,particularly dementia.However,the underlying neural mechanism of these cognitive impairments remains unclear.We aimed to investigate the correlation between whole brain resting state functional connections(RSFCs)and the cognitive status in 95 patients with T2DM.We constructed an elastic net model to estimate the Montreal Cognitive Assessment(MoCA)scores,which served as an index of the cognitive status of the patients,and to select the RSFCs for further prediction.Subsequently,we utilized a machine learning technique to evaluate the discriminative ability of the connectivity pattern associated with the selected RSFCs.The estimated and chronological MoCA scores were significantly correlated with R=0.81 and the mean absolute error(MAE)=1.20.Additionally,cognitive impairments of patients with T2DM can be identified using the RSFC pattern with classification accuracy of 90.54%and the area under the receiver operating characteristic(ROC)curve(AUC)of 0.9737.This connectivity pattern not only included the connections between regions within the default mode network(DMN),but also the functional connectivity between the task-positive networks and the DMN,as well as those within the task-positive networks.The results suggest that an RSFC pattern could be regarded as a potential biomarker to identify the cognitive status of patients with T2DM.

KEYWORDS Type 2 diabetes mellitus;Resting state functional connectivity;Elastic net;Support vector machines;MoCA

Introduction

Type 2 diabetes mellitus(T2DM)is typically accompanied by cognitive impairments and is associated with a much higher risk of dementia[1,2].Patients with this disorder may experience a deterioration of memory,attention,information processing speed,and executive function[1,3].An in-depth understanding of the causative neuro-mechanism of cognitive impairment in the early stages of T2DM could help clinicians to identify patients with a high risk of dementia,and to consequently introduce effective interventions to retard and even arrest the progression of subtle cognitive decrements[4].However,the mechanism behind the cognitive impairments in patients with T2DM remains unclear.Recently,an increasing number of neuroimaging studies on patients with T2DM revealed alterations in the gray matter volume and white matter integrity associated with cognitive impairment[5,6].This suggests that changes in the neural substrates of patients with T2DM are likely linked to cognitive dysfunction.

Resting state functional connections(RSFCs)are typically used in neuroimaging studies to measure the correlation between the fMRI time series of different brain regions without any external disturbance,they also can reflect on certain intrinsic mechanisms in the human brain[7].The whole brain RSFCs are a cumulation of the connections of all paired regions of our brain and are relatively rich in information relevant to the intrinsic interactions among the regions of the brain that are induced by the spontaneous neural activities.Therefore,an analysis of the connectivity patterns based on the whole brain RSFCs,in comparison with the techniques based on individual connectivity,could provide us a much more comprehensive understanding on the neural mechanism of certain cognitive disorders.Whole brain RSFCs have previously been used in a number of studies addressing cognitive disorders.Previous studies have reported that patients with Alzheimer’s disease(AD)or mild cognitive impairment(MCI)had abnormal connectivity patterns[8,9].Whole brain RSFCs have also been used to successfully select biomarkers in recent fMRI studies,such as age estimation[10]and identifi cation of psychiatric disorders[11].Thus,it can be speculated that some of connectivity patterns could be regarded as potential biomarkers to either evaluate or identify cognitive impairment in patients with T2DM.

Recently,a number of RSFC-based studies investigated differences between patients with T2DM and normal controls using their neural substrate and some have reported abnormal RSFCs in their default mode network(DMN)[12-14].Yang et al.had particularly[14]reported abnormalities in the connectivity within the DMN along with those among other RSFC networks in patients with T2DM related cognitive impairment.Another resting state study demonstrated a decrease in RSFCs associated with the attention network of patients with T2DM and cognitive impairment[15].In addition to RSFC,Cui et al.[16]reported that the amplitude of low frequency fluctuation and regional homogeneity,which were both calculated from the resting state brain activity,changed in regions of DMN and other brain regions,such as the visual and auditory network of patients with T2DM and cognitive impairment.These findings suggested that the whole brain RSFCs provided more information about the mechanism of T2DM-related cognitive impairment.However,these studies focused on the difference in individual connectivity at a group level between normal controls and patients with T2DM,and did not directly examine the relationship between the whole brain RSFC patterns and cognitive status in patients with T2DM.Thus,these neural substrate findings cannot be used in clinical practice for quantitatively evaluate the progression of cognitive impairment in patients with T2DM.

To bridge this gap, the present study used a multi-variate pattern analysis (MVPA), which has been typically used in neuroimaging studies and demonstrated to be an effective method to analyze fMRI data [10,11,17] and examine the correlations between cognitive status of patients with T2DM and whole brain RSFCs. We particularly aimed to identify potential biomarkers that could be used to detect patients with T2DM, who have a high risk of presenting cognitive impairment.

We first constructed a predictive model with whole brain RSFCs serving as predictors and the Montreal Cognitive Assessment(MoCA),which can be used to measure the degree of general cognition of all participants(regardless of cognitive impairment)serving as the dependent variable.However,whole brain RSFCs include a large number of predictors,which is significantly larger than the number of participants(i.e.,the number of measurement).Furthermore,there are a number of correlations among features of whole brain RSFCs,making it difficult to estimate a model while utilizing an ordinary regressing strategy such as general linear modeling.To address these problems,we used an elastic net(E-Net)to construct a regression model to reduce the dimensions of the features.The E-Net can both successfully select an appropriate feature and estimate a model by balancing the goodness offit and model complexity and compensating for the correlations among the whole brain RSFC features[18].

The surviving predictors after estimating the E-Net model,namely the selected features of whole brain RSFCs,are the potential biomarkers to evaluate the MoCA scores of participants.To further validate the roles of these selected features,a machine learning method was used to evaluate their discriminative ability.This step can not only evaluate the diagnostic valuation of these features,but also help in revealing the neural mechanism behind cognitive impairment in patients with T2DM.

Results and discussion

Characteristics of patients

The characteristics of the enrolled patients are shown in Table 1.The average age of all patients was 54.51 years and 68.4%patients were men(65 of 95).Moreover,the average MoCA score of the patients was 25.85,which was similar to normal cognition.Functional MRI data of the 95 patients enrolled,were preprocessed and we performed a functional connectivity analysis to detect whole brain RSFCs for all the patients.Subsequently,functional connectivity between two of the 90 brain regions were obtained for every patient.Thus,we identified 4005 features for each patient with T2DM for further analysis.

Cognition estimation and selection of key RSFCs

To explore the relationship between cognition of patients with T2DM and whole brain RSFCs,we first constructed a cognition estimation model to reveal the association between theMoCA score and whole brain RSFCs.Considering the limited sample size in comparison with the number of RSFCs,we aimed to reduce the latter used in the cognition prediction model to avoid a potential overfitting.Since each RSFC was not directly correlated with a MoCA score,we utilized an E-net to obtain the most valuable RSFCs to estimate the MoCA.E-net was a very useful tool as it could reduce the dimensionality of features,and concurrently isolate key features.We used a 10-fold cross-validation via minimum criteria to select the tuning parameter in the E-net model (Figure 1). The best MoCA estimated model was detected as α= 0.9 and λ=0.424.The results of the MoCA estimation revealed an observable association between the MoCA score and whole brain RSFCs.It showed that,the estimated MoCA scores and real MoCA scores was significantly correlated with R=0.81(Figure 2)and P=0.002.The results indicated that there was no better prediction in the 500 permutations(Figure 2).Additionally,the MAE between the estimated and real MoCA score was 1.20(Figure 2),suggesting that the cognitive decline of patients with T2DM that was indexed by MoCA scores could be predicted by a combination of some RSFCs(i.e.,connectivity pattern)with a relatively good performance.

Table 1 Characteristics of the T2DM patients examined in this study

We selected a connectivity pattern consisting of 23 RSFCs using a non-zero coefficient in the best MoCA scores estimation model corresponding to greater contributions to the estimation of the MoCA scores.We did not include any clinical characteristics in the best model.The selected 23 RSFCs,as key features of the MoCA estimation,are illustrated in Figure 3,and detailed information about the same,such as their correlations with MoCA as well as their contributions to the classification of cognition is listed in Table 2.As indicated in Table 2,the selected connectivity pattern mainly consisted of RSFCs within the DMN,followed by between the DMN and other resting state brain networks such as the audio network(AN),the visual network(VN),effective control network(ECN),and the motion network(MN).Thesefindings suggest that the cognitive impairment of patients with T2DM may be associated with the abnormality of connectivity with and between these different resting state networks.

Among the selected 23 RSFCs,there were 3 between the brain regions within DMN(Figure 3).The first was the functional connection between the right caudate and the right hippocampus.The second was the connectivity between the right superior frontal gyrus(orbital)and the right superior frontal gyrus.The third was the connectivity between the right rectus and the left superior frontal gyrus(orbital).DMN is considered to be a major contributor to cognitive function[19,20],and the abnormal activity of the DMN has been demonstrated to be related to some mental disorder such as MCI[21],AD[22]or other mental disorders[23].Particularly,Yang et al.[14]found that cognitive impairment in patients with T2DM presented a significantly decreased RSFC strength within the DMN than the normal controls;however,such an abnormality was not observed when comparing patients with T2DM who have normal cognition and normal controls.This evidence along with our findings suggested that the cognitive impairment of patient with T2DM might be related to the normality of connections within the DMN.

Figure 1 Feature selection with E-net

Among the 23 selected key RSFCs,there were 9 between regions of DMN and regions of task-positive networks(Figure 3,Table 2).First,there were 5 connections between the DMN and the ECN,namely the ones between the left hippocampus and left inferior frontal gyrus(opercular),the left precentral gyrus and left inferior frontal gyrus(orbital),the right putamen and left inferior frontal gyrus(opercular),the left inferior parietal gyrus and left inferior frontal gyrus(opercular),and between the right parahippocampal gyrus and left middle frontal gyrus.Second,there were two connections between the DMN and the AN,namely the ones between the right anterior cingulum cortex and left superior temporal gyrus(pole),and between the right precuneus and right heschl.Lastly,there were two connections between the DMN and the MoN,namely the ones between the right superior frontal gyrus(medial)and right middle cingulum cortex,and between the left superior frontal gyrus(medial)and right paracentral lobule.

The present study also revealed the connections among task-positive networks(Figure 3,Table 2).There were 3 connections between the AN and the ECN,3 connections between the AN and the VN,and 2 connections between the VN and ECN.Additionally,there were 2 connections within the ECN and 1 connection within the AN.The results suggested that hippocampus could be a crucial central hub in the selected RSFCs,which was consistent with findings of previous studies focusing on the insulin resistance in youth[24-26].Although the dysconnectivity patterns in hippocampal and striatal reward regions in youth susceptible to diabetes relates directly to the degree of insulin resistance,our finding further suggested thatRSFCsrelated with hippocampuswere significantly correlated with the cognitive impairments in patients with T2DM.

Cognition classification and performance evaluation

We conducted a correlation analysis between the selected features and MoCA scores to investigate the relationship of these selected RSFCs to the cognitive performance.Seven of the selected 23 RSFCs demonstrated significant positive correlation with MoCA scores(r＞0,and P＜0.05 after Bonferroni correction for the number of features),and two of the selected 23 RSFCs showed a significant negative correlation with MoCA scores(r＞0,and P＜0.05 after Bonferroni correction for the number of features).Additionally,there were three RSFCs that showed marginally significant negative correlations with MoCA scores.Detailed information regarding the correlation analysis is shown in Figure 4.Moreover,these RSFCs were presented in Figure 5.We can see that five of the nine RSFCs with significant or marginally significant correlations with MoCA scores involved in brain regions of DMN,suggesting the importance of DMN.To investigate if there were differences between the patients with normal and impaired cognition,a two-sampled t-test of the selected 23 RSFCs were performed between the two groups(Figure 6).Most of these RSFCs(18/23)showed significant differences between these two groups(P＜0.05),while all of the 5 RSFCs without significant differences were either between regions of DMN and regions of task-positive networks or among the regions of task-positive networks.

Figure 2 Results of MoCA estimation and cognition classification

Classification of the results of cognition indicated that the SVM classifier achieved a performance with an AUC of 0.9737 and the accuracy of classification was 90.54%,based on the selected 23 RSFCs while using a 10-fold cross validation(Figure 2),and there was an absence of better classification results that were detected using the non-parametric permutations with a P value of P=0.002(Figure 2).Additionally,the sensitivity and specificity of the classification model were 92.11%and 88.89%,respectively.The contribution of each selected RSFC to the classification of cognition was analyzed with its weight in the SVM model.Detailed information about the contributions of the 23 RSFCs is provided in Table 2.

The relatively high sensitivity and specificity of the classi fication model ensured that the 23 RSFCs selected via E-net analysis could be used as the potential biomarkers to either evaluate or identify the cognitive impairment of T2DM patients.It was noted that,as revealed via correlation analysis,more than half of the 23 RSFCs(14 out of 23)were not significantly correlated with the MoCA score.Therefore,the singlevariable method(e.g.,Pearson correlation) ‘underestimates”the role of these RSFCs in evaluating the cognitive decline of patients with T2DM.Furthermore,due to the limited number of brain regions and connectivity between these regions,the connectivity patterns do not individually match with a large number of cognitive processes.Therefore,according to a network theory,an interaction among different regions or networks may be more relevant to a certain cognitive process than any particular individual region or network[22,27,28].If this idea is correct,the RSFC pattern that we identified using a MVPA may provide a relatively complete profile of the neural mechanism of cognitive impairment in patients with T2DM,and reveal potential biomarkers to evaluate or identify the associated cognitive impairment.

Limitations

Figure 3 The RSFCs selected from the MoCA estimation model

The present study has a number of limitations.First,it was a cross-sectional investigation that only focused on the cognitive impairment in patients with T2DM.However,patients with T2DM and cognitive impairment are at a significant risk of developing AD and MCI;therefore,further longitudinal studies are required to elucidate the changes in the whole brain RSFCs with regard to this cognitive impairment,particularly as it progresses into either AD or MCI.Additionally,although the sample size in the present study was large for an fMRI study,it was relatively small for a clinical study,and this is also a reason to explain the absence of independent validation.Future studies with a larger population sample or multicenter imaging data may further confirm the potentiality of the RSFCs in predicting the cognitive impairment in patients with T2DM,and develop an outperformed prediction model with independent validation.Lastly,possible effects of different hypoglycemic agents may impact the results.Future studies should take this into consideration and explore its potential effects.

Conclusion

Here we used an MVPA method to investigate the correlations between the RSFCs and cognitive impairments in patients with T2DM.We observed that the T2DM associated cognitive decline could be predicted via an RSFC pattern.This connectivity pattern included the connections between regions within the DMN,along with those between the task-positive networks and the DMN,as well as those within the taskpositive networks.Although majority of these connections do not individually demonstrate significant correlation with the cognition of patients with T2DM,all of them together are significantly correlated with cognitive impairment.These results suggest that this RSFC pattern may play important roles in the T2DM related cognitive decline;therefore,can be regarded as potential biomarkers to evaluate or identify this cognitive impairment.

Materials and methods

Patients

We recruited 95 patients with T2DM from the Henan Provincial People’s Hospital.We used the newest criteria of the American Diabetes Association[29]to define the disease.All these patients were between 43 and 75 years of age,had been suffering from the disease for longer than a year,were routinely treated with hypoglycemic agents,and were closely self-monitored.Our exclusion criteria were as follows:(1)a history of brain lesions such as stroke or brain tumors;(2)a history of alcohol or substance abuse;(3)a psychiatric or neurological disease,including major depression or any otherpsychopathology;(4)a history of hypertension;(5)a history of hypoglycemic episodes;and(6)contraindications to MRI.All patients underwent a MoCA test to examine their general cognition[30];an education-adjusted MoCA score was detected for every patient[31].The Ethics Committee of Henan Provincial People’s Hospital approved the present study in accordance with the Helsinki declaration,and written informed consent was obtained from all patients prior to their participation in our study.The characteristics of the patients are summarized in Table 1.

Table 2 RSFCs selected as the key features from the MoCA estimation model for cognition classification

Figure 4 Correlations between the selected RSFCs and MoCA scores

Data acquisition and imaging processing

Figure 5 RSFCs demonstrate significant or marginally significant correlations with MoCA scores

All MRI images were captured using a Magnetom Trio 3.0 T scanner(Siemens,Erlangen,Germany)at the Radiology Department of Henan Provincial People’s Hospital utilizing a 12-channel receive-only head coil.The patients were instructed simply to rest with their eyes closed,to relax,but not fall asleep,and to keep still.Functional resting state images were acquired using a gradient echo T2-weighted pulse sequence with TR=2000 ms,TE=30 ms,matrix=64×64,field of view(FOV)=240 mm×240 mm,thickness=4 mm,and flip angle=90°.The functional resting state scan lasted for 7 min,and 210 volumes were collected.

The data was preprocessed via the Data Processing Assistant for Resting-State toolbox(DPARSF;http://www.restfmri.net/forum/DPARSF). We discarded thefirst ten volumes to equilibrate the magnetic field.All remaining volumes were then realigned to adjust for head motions using a least-squares minimization technique.Any patient with head motion＞2.0 mm of translation or＞2.0°of rotation in any direction were excluded.We then further processed the image data with the spatial normalization based on Montreal NeurologicalInstitute (MNI)space [32],and resampled to 3×3×3 mm3.Smoothing(full width at half maximum[FWHM]=4 mm),detrending,and filtering(0.01-0.08 Hz)were performed in order.Regression analysis including six head motion parameters and mean time series of global,white matter,and cerebrospinal fluid signals were conducted to remove the possible effects on the results.

Weperformed functionalconnectivity analysisusing Resting-State fMRI Data Analysis Toolkit(REST).The registered images were divided into 90 regions according to the automated anatomical labeling atlas[33].The atlas divided the cerebrum into 90 regions(45 in each hemisphere).The representative time series of each region used to analyze functional connectivity was estimated by averaging the fMRI time series over all voxels in the region.We evaluated the functional connectivity between each pair of regions with a Pearson correlation coefficient.A Fisher’s transform was applied to improve the normality of the correlations.Therefore,a resting-state functional network captured via a 90×90 symmetric matrix was detected for each patient.We removed 90 diagonal elements and the upper triangle elements of the matrix were used as features for further analysis,i.e.,the feature space was spanned by the (90×89)/2=4005-dimensional feature vectors.

MoCA estimation and feature selection

The MoCA scores were estimated with a linear model based on the 4005 RSFCs and the clinical characteristics of patients including their age,sex,fasting glucose,glycated hemoglobin(HbA1c),total cholesterol,body mass index(BMI),and disease duration.The linear model is defined as follows:

Figure 6 RSFCs comparison between normal cognition group and cognition impairment group

Where,n was the number of features used in the model,here n=4012;y was MoCA score of the patients;xi(i=1,2,...,n)was the predicting parameter in the model,such as the RSFC or its clinical characteristics;βi(i=0,1,2,...,n)was coefficient of each parameter,and ε was the error term.E-Net was used to estimate the coefficients of the model[34].Here,we simultaneously estimated the coefficients and selected the feature,minimizing the following cost function:

Where,N was the number of patients;xijwas the jth feature of the ith patients;yiwas the MoCA score of the ith patients;and λ and α were regularization parameters.Subsequently,features with a greater contribution to the MoCA estimation could be selected.

We used glmnet[35]to estimate MoCA.We selected λ and α and estimated MoCA using a 10-fold cross validation.The candidate values for the λ and α parameters were both restricted to a range between 0 and 1 in steps of 0.1,in which the optimal λ and α were selected via minimum criteria.The model goodness criteria that we employed were the MAE and the correlation between the estimated and real MoCA scores.The features of the best estimation model were further used to classify cognition.

Cognition classification and performance evaluation

To classify cognition,36 patients with a MoCA score＞26 were taken as the normal cognition group,while 38 patients with a MoCA score＜26 were taken as the cognition impairment group.We classified cognition for these two groups with the features selected from the MoCA estimation model.When we obtained the data for features that highly correlated with MoCA,we used linear support vector machines(SVM)to perform cognition classification.The toolbox LibSVM[36]was used to construct the classification model based on the selected features.The results were based on the best parameter setting.We then investigated the discrimination performance of the SVM model with a 10-fold cross validation.The classification results were further interpreted using the classification accuracy,specificity,and sensitivity based on the cross-validation.Moreover,the ROC curve and AUC were also calculated.

Evaluation of the selected RSFCs

We performed a correlation analysis between the selected RSFCs and MoCAscores to investigate the relationship of these selected RSFCs to the cognitive performance. A linear regression was identified between the RSFC and MoCA score of each included patient, so we can detect the correlation between them.

We also analyzed the contribution of each RSFC to the cognitive classification.The weight value of the RSFC in the SVM model was considered as contribution to classify cognition.The weight value of each feature was calculated through summing the coefficients allocated to the RSFC across the 10-fold cross validation,and then normalized by dividing the maximal coefficients across all the selected features.

Permutation tests

The framework of permutation tests to assess predictive performance has been used in several previous studies[10,11].We performed 500 permutations to estimate the probability of detecting identical MoCA estimation/classification performance.Specifically,the MoCA(for MoCA estimation)/cognition labels(for the classification of cognition)of the patients were permuted 500 times randomly.Furthermore,the P value of the probability for detecting the estimation/classification accuracy was defined as follows:

Where,N=500 here represents the number of permutations;to estimate MoCA,NBetterpredictionswas the number of permutations with larger correlations(compared to the result based on non-permuted MoCA scores)between the permuted and estimated MoCA scores for the estimation;and to classify cognition,NBetterpredictionswas the number of permutations with higher classification accuracy(compared to the result based on non-permuted labels).

Data availability

Please contact with the corresponding authors for data access.The MRI image data will be available upon request after approval by the Ethics Committee of Henan Provincial People’s Hospital,China.

Authors’contributions

HY,TL,and MW collected the data,MW,YL,and JT contributed to the study design and revised the manuscript for intellectual content,XC,ZT,and YD performed the analysis,ZL and JL wrote the manuscript.JT is the guarantor of this manuscript,has full access to all the data in the study,and takes responsibility for the integrity of the data and the accuracy of the data analysis.All authors have read and approved the final manuscript.

Competing interests

The authors have declared no competing interests.

Acknowledgments

This paper is supported by the National Natural Science Foundation of China(Grant Nos.81772012,81227901,61673051,81641168,31470047,81271565,81527805,and 61231004),the National Key R&D Program of China(Grant No.2017YFA0205200),the Youth Innovation Promotion Association,Chinese Academy of Sciences(GrantNo.2019136).