Aberrant Brain Regional Homogeneity and Functional Connectivity in Middle-Aged T2DM Patients: A Resting-State Functional MRI Study

Introduction

Type 2 diabetes mellitus (T2DM) is characterized by progressive insulin secretion defects on the basis of insulin resistance (American Diabetes Association, 2015) and is accompanied by hyperinsulinemia, hyperglycemia, increased inflammatory mediators, advanced glycation end products and amyloid deposition. These risk factors have caused T2DM to become closely associated with cognitive dysfunction, which can even progress to dementia (Strachan et al., 2011). Data from the International Diabetes Federation Annual Report¹ showed that as many as 415 million people were living with diabetes in 2015. With a risk of dementia greater than 6–7% (Koekkoek et al., 2015), the report called for effective means to cope with the emerging novel complications. It is highly important to illuminate the neurological mechanisms of cognitive dysfunction induced by T2DM, as this could enable the application of preventive treatment strategies.

Resting-state functional magnetic resonance imaging (rs-fMRI) can provide excellent opportunities to investigate cognitive dysfunction related to T2DM in vivo (Musen et al., 2012; Xia et al., 2013; Chen et al., 2014; Cui et al., 2014, 2015). Among several rs-fMRI indices, regional homogeneity (ReHo) and functional connectivity (FC) are often used to evaluate brain activity synchronization in normal subjects or patients. ReHo measures the similarity of the time series of a given voxel to its nearest neighborly voxels, and it is related to local synchronization (Zang et al., 2004). There were only two studies which investigated ReHo changes in T2DM patients (Cui et al., 2014; Peng et al., 2016). Cui et al. (2014) reported a decrease of ReHo primarily in the occipital lobe that was related to poor cognitive performance in T2DM patients. Peng et al. (2016) revealed abnormal brain activity in T2DM patients with and without microangiopathy using ReHo analysis. It appears that the occipital lobe, temporal lobe, precuneus, frontal gyrus and insula brain regions are susceptible to T2DM (Peng et al., 2016). Regarding FC, it measures the similarity of the time series of two relatively remote brain regions (Biswal et al., 1995). Alterations in the FC of the default mode network (DMN) were recently investigated using independent component analysis (Cui et al., 2015). T2DM patients exhibited increased connectivity in the anterior sub-network of the DMN, but decreased connectivity in the posterior sub-network of the DMN (Cui et al., 2015).

As effective indicators reflecting the intrinsic organization of the resting brain, ReHo and FC have been conjunctively applied in a number of studies of complex brain functional activity and its alterations in diseases and aging (Damoiseaux et al., 2008; Jiang et al., 2015; Cui et al., 2016; Wang et al., 2016). These two indices have proved to be mutually complementary for detecting local and remote brain activity synchronization. However, the aberrant brain function activity in middle-aged T2DM patients has not been elucidated using these two methods together. Therefore, it has been suggested that the combination of ReHo and FC analyses could provide more information about aberrant brain activity synchronization in T2DM patients than either index alone.

Based on the above-mentioned findings, we hypothesized that T2DM relates to aberrant spatially local and remote brain activity synchronization that might be associated with cognitive impairment. In the present study, we first investigated the abnormal brain activity using ReHo analysis. According to the peak Montreal Neurological Institute (MNI) coordinates of brain regions with abnormal ReHo, we further explored the FC between these brain regions and all other brain voxels. Finally, we investigated the relationships of aberrant ReHo and FC with neuropsychological test scores and diabetes-related parameters. Our exploratory study provides complementary information for a deeper understanding of the mechanisms underlying T2DM-associated brain dysfunction.

Materials and Methods

Subjects

The subjects were recruited from inpatients and communities between December 2013 and January 2015. The T2DM patients met the criteria published by the World Health Organization in 1999 (Alberti and Zimmet, 1998), and they were diagnosed with diabetes for more than 2 years. All of the subjects were right handed, between the ages of 45 and 60 years old and had at least 6 years of education. Healthy controls and T2DM patients were group-matched by age, sex, education level and body mass index (BMI). The exclusion criteria for all subjects were a history of stroke, white matter hyperintensity showing confluence or severe changes (Wahlund et al., 2001; Chen et al., 2014), Parkinson’s disease, probable dementia (Mini-Mental State Examination [MMSE] score ≤24; Galea and Woodward, 2005), depression (Hamilton Depression Rating Scale [HAMD] >7; Hamilton, 1960) or any other neuropsychiatric disorder, thyroid dysfunction, severe renal impairment, severe visual or hearing loss and any contraindications to MRI. Of the 64 recruited subjects, 25 T2DM patients and 25 healthy controls were screened for subsequent studies (Table 1). All the participants provided written informed consent, and the study was approved by the local Medical Research Ethics Committee of the Southwest Hospital (Chongqing, China).

TABLE 1

Table 1. Demographic and clinical data comparison results of T2DM patients and healthy controls.

Clinical Data

Data on height, weight and arterial blood pressure were obtained by measurements on the spot. Medical history and medication use were acquired from medical records, questionnaires and interviews. The T2DM patients were undergoing treatment with one or more of the following medications: hypoglycemic agents (insulin, n = 12; biguanides, n = 14; sulfonylureas, n = 8; glinides, n = 2; α glucosidase inhibitors, n = 5; thiazolidinediones, n = 3); antihypertensives (angiotensin II receptor blockers, n = 3; angiotensin-converting enzyme inhibitors, n = 1; calcium channel blockers, n = 1); and statins (n = 2). Then, BMI was calculated ([weight in kg]/[height in m]²). Blood samples were collected by venipuncture for the measurement of fasting plasma glucose, fasting insulin, fasting C-peptide, glycosylated hemoglobin (HbA_1c), triglycerides, total cholesterol, triglyceride, high density lipoprotein (HDL) cholesterol, low density lipoprotein (LDL) cholesterol, homocysteine, blood urea nitrogen, serum creatine, urinary microalbumin, cystatin C, uric acid, free triiodothyronine (FT3), free thyroxine (FT4) and thyroid stimulating hormone (TSH) levels, with the patient in the fasting state after overnight abrosia. Because of the influence of exogenous insulin (primarily the injection of insulin) on fasting plasma insulin levels, the score for homeostasis model assessment of insulin resistance was not calculated.

Neuropsychological Tests

The cognitive status of the participants was evaluated by a battery of neuropsychological tests covering major cognitive domains. The MMSE and Montreal Cognitive Assessment (MoCA) tests were conducted to evaluate the general level of cognition. The HAMD was used to exclude possible cases of depression. Then, attention, executive function and psychomotor speed were evaluated by the Trail Making Test parts A and B (TMT-A and TMT-B; Bowie and Harvey, 2006). Working memory was evaluated by the Digital Span Test (DST, forwards and backwards; Gong, 1992). Semantic memory was evaluated by the Verbal Fluency Test (VFT; Brucki and Rocha, 2004). Episodic memory was evaluated by the Auditory Verbal Learning Test (AVLT, including immediate recall, short term delayed recall, long term delayed recall, long term delayed recognition and total score; Schmidt, 1996). All of the tests were administered by a trained neuropsychologist who was blinded to the subjects’ diagnoses. Each individual required approximately 60 min to complete the test in a fixed order.

MRI Data Acquisition

MRI data were acquired using 3.0-T MR scanner (Trio, Siemens Medical, Erlangen, Germany) with a 12-channel head coil at approximately 9:00 a.m. after clinical data collection and neuropsychological tests. Earplugs were used to diminish the influence of scanner noise. All of the subjects were instructed to remain awake with their eyes closed and to not think about anything in particular. T2-weighted and fluid attenuated inversion recovery (FLAIR) images were obtained for the assessment of white matter hyperintensity and intracranial structural lesions. High-resolution T1-weighted images were obtained sagittally using a volumetric 3D magnetization prepared rapid gradient-echo sequence: TR 1900 ms, TE 2.52 ms, flip angle 9°, FOV 256 mm × 256 mm, slices = 176, thickness = 1 mm, in-plane matrix = 256 × 256 and voxel size = 1 mm × 1 mm × 1 mm. Resting-state functional images were obtained transversely using an echo planar imaging (EPI) sequence: TR 2000 ms, TE 30 ms, flip angle 90°, FOV 192 mm × 192 mm, slices = 36, thickness = 3 mm, in-plane matrix = 64 × 64 and voxel size = 3 mm × 3 mm × 3 mm; 240 volume acquisition cost 8 min and 8 s.

Data Analysis

The T2-weighted and FLAIR images were separately reviewed by two experienced neurologists blinded to the group status. The two neurologists reached consensus on brain disease assessment through discussion. T1-weighted images were processed by the VBM 8 toolbox, based on Statistical Parametric Mapping 8 software (SPM 8, Wellcome Department of Imaging Neurosciences, University College London, UK²) running on MATLAB 7.14 (Math-Works, Natick, MA, USA³). The process obeyed the following standard protocol (Whitwell, 2009). Intracranial tissue was segmented into gray matter, white matter and cerebrospinal fluid. The volumes of the gray matter, white matter and cerebrospinal fluid were automatically calculated, and the brain parenchyma volume was equal to the sum of the gray matter and white matter volumes. Functional image preprocessing and ReHo calculation were conducted with Data Processing Assistant for Resting-State fMRI (DPARSF v2.3⁴). The first 10 volumes of functional images were removed for their possible heterogeneity. Slice-timing and realignment were performed for head motion using the remaining images. No subjects were excluded from the analyses due to head motion >2.0 mm in any direction of x, y and z or 2.0° of any angular motion. Then, the functional images were coregistered to the high-resolution T1-weighted images. Subsequently, the coregistered images were normalized to the MNI template. Detrending was performed to remove linear trends from the image time series, and the data were filtered at the 0.01–0.08 Hz band. ReHo calculation was performed using the preprocessed images, and the resulting images were smoothed with an isotropic Gaussian kernel of 4 mm full-width half-maximum.

FC calculation was conducted with rs-fMRI data analysis toolkits (REST v1.8⁵). The coordinates of the center of the seed region were based on the inter-group analysis of ReHo, and the radius was 4 mm. Correlation analysis of time course was performed between the spherical seed region and every voxel of the whole brain for each subject. Finally, Fisher’s r-to-z transformation was applied to improve the normality of the ReHo and FC maps (Song et al., 2011).

Statistical Analysis

Demographic/clinical data and neuropsychological test scores were compared between the T2DM patients and healthy controls using SPSS software (version 20.0; SPSS Inc., Chicago, IL, USA). We used the χ² test for proportions, the independent samples t-test for normally distributed continuous data and the Mann-Whitney U test for non-normally distributed data. The distribution of the data was analyzed using the Kolmogorov-Smirnov test. A p-value < 0.05 was considered statistically significant.

ReHo and FC map analyses were conducted with REST software. Jarque-Bera goodness-of-fit test was applied to test the distribution of ReHo and FC maps within 61 × 73 × 61 brain mask (70,831 voxels). The results showed that both ReHo and FC maps were normally distributed (100% voxel with p > 0.05). Intra-group patterns were analyzed using the one-sample t-test, with ReHo maps compared to “1” (Chao-Gan and Yu-Feng, 2010) and FC maps compared to “0”. The inter-group differences on the ReHo and FC maps were analyzed using the two-sample t-test. Age, sex, education and BMI were imported as covariates. Gray matter maps obtained from VBM analysis were also used as covariates to exclude the possible effects of gray matter atrophy on the results. The voxel level was p < 0.005. The outcome maps of ReHo and FC were then corrected with AlphaSim (cluster p < 0.01, cluster size >16 voxels, full width at half maximum = 4 mm, number of Monte Carlo simulations = 1000, cluster connection radius: rmm = 5.00) to detect the brain regions that survived the comparison and to determine their peak MNI coordinates. Finally, the z scores of abnormal brain regions were extracted for subsequent correlation analyses.

To explore the relationships among ReHo values, FC z scores, neuropsychological performances and diabetes-related parameters in T2DM patients, partial correlation analyses were conducted with SPSS software (age, sex, education and BMI were included as control variables).

Results

Demographic and Clinical Data

The subjects’ demographic and clinical characteristics are reported in Table 1. The two groups did not show significant differences in terms of age, sex, education level, BMI, systolic/diastolic blood pressure, total cholesterol, triglycerides, LDL cholesterol, homocysteine, blood urea nitrogen, urinary microalbumin, cystatin C, uric acid, FT4 or TSH (all p > 0.05). However, the T2DM patients exhibited significantly higher HbA_1c and fasting plasma glucose than healthy controls, but lower HDL cholesterol and serum creatine than healthy controls (all p < 0.05). The FT3 levels of the T2DM patients were significantly higher than those of the healthy controls (p < 0.05) but within the normal range. No subject was excluded for severe white matter hyperintensity or intracranial lesions.

Neuropsychological Test Results

The T2DM patients scored poorer on the TMT-A, TMT-B and DST forward tests (all p < 0.05) and had no significant decreases in the other neuropsychological tests (all p > 0.05; Table 2).

TABLE 2

Table 2. Comparison of neuropsychological test results of T2DM patients and healthy controls.

Comparison of Brain Volume Across Groups

The gray matter, white matter and brain parenchyma volumes (the total volumes of gray and white matter) did not differ significantly between the two groups (Table 1).

ReHo Analysis

The T2DM patients exhibited significantly higher ReHo values than the healthy controls in the bilateral anterior cingulate gyrus (ACG) and significantly lower ReHo values in the right fusiform gyrus (FFG), right precentral gyrus (PreCG) and right medial orbital part of the superior frontal gyrus (SFG; Table 3; Figures 1, 2A).

TABLE 3

Table 3. Brain regions with different regional homogeneity (ReHo) values between T2DM patients and healthy control groups.

FIGURE 1

Figure 1. Regional homogeneity (ReHo) and functional connectivity (FC) maps of intra-group and inter-group comparisons. Voxel p < 0.005 (AlphaSim correction, cluster p < 0.01, cluster size >16 voxels). Color scale denotes the t value. ACG, anterior cingulate gyrus; FFG, fusiform gyrus; PreCG, precentral gyrus; ORBsupmed, superior frontal gyrus, medial orbital; INS, insula; ANG, angular gyrus; HIP, hippocampus; SMA, supplementary motor area; PCUN, precuneus; MTG, middle temporal gyrus; SFGdor, superior frontal gyrus, dorsolateral; ORBsup, superior frontal gyrus, orbital part; MFG, middle frontal gyrus; R, Right; L, Left.

FIGURE 2

Figure 2. Brain regions with significant differences in ReHo (A) and FC (B–E) between Type 2 diabetes mellitus (T2DM) patients and healthy controls. Voxel p < 0.005 (AlphaSim correction, cluster p < 0.01, cluster size >16 voxels). ACG, anterior cingulate gyrus; FFG, fusiform gyrus; PreCG, precentral gyrus; ORBsupmed, superior frontal gyrus, medial orbital; INS, insula; ANG, angular gyrus; HIP, hippocampus; SMA, supplementary motor area; PCUN, precuneus; MTG, middle temporal gyrus; SFGdor, superior frontal gyrus, dorsolateral; ORBsup, superior frontal gyrus, orbital part; MFG, middle frontal gyrus; R, Right; L, Left.

FC Analysis

The results of the FC group comparisons are presented in Table 4 and Figures 1, 2B–E. The T2DM patients exhibited weaker negative connectivity relative to the healthy controls between the right PreCG and bilateral precuneus. In addition, negative to positive transformations of connectivity were observed in several brain regions, such as the FC between the right FFG and right middle temporal gyrus (MTG). Weaker positive connectivity was observed between the right orbital part and right medial orbit of the SFG. There were also brain regions that experienced positive to negative transformations of connectivity, such as the FC between the ACG and the right precuneus. In addition, the T2DM patients exhibited greater negative connectivity between the ACG and the left angular gyrus (ANG).

TABLE 4

Table 4. Brain regions with different functional connectivity (FC) between T2DM patients and healthy control groups.

Correlation Analysis

Better DST forward scores were correlated with higher ReHo values of the right PreCG (ρ = 0.527, p = 0.014; Figure 3A) and stronger negative FC between the right FFG and the MTG (ρ = −0.437, p = 0.048; Figure 3B). Higher levels of HbA_1c were associated with a stronger positive FC between the right orbital part and the right medial orbit of the SFG (HbA_1c [%], ρ = 0.476, p = 0.029; HbA_1c [mmol/mol], ρ = 0.484, p = 0.026; Figures 3C,D). No significant correlations were observed between the FC of other brain regions and clinical data or cognitive test scores.

FIGURE 3

Figure 3. The relationships of aberrant ReHo/FC with neuropsychological performances/diabetes-related parameters after adjustment for age, sex, education and body mass index (BMI). (A) ReHo values of PreCG.R vs. DST forwards. (B) z scores of FFG.R-MTG.R vs. DST forwards. (C) z scores of ORBsupmed.R-ORBsup.R vs. HbA_1c (%). (D) z scores of ORBsupmed.R-ORBsup.R vs. HbA_1c (mmol/mol). DST, Digital Span Test; PreCG, precentral gyrus; FFG, fusiform gyrus; MTG, middle temporal gyrus; ORBsupmed, superior frontal gyrus, medial orbital; ORBsup, superior frontal gyrus, orbital part; R, Right; L, Left.

Discussion

It has been proposed that the normal brain development experiences a “local to remote” FC organization pattern with important implications for revealing the neural substrate of cognition (Fair et al., 2009). Previous studies have also reported that local alteration of functional synchronization can diffuse to remote synchronization in diseases such as schizophrenia (Jiang et al., 2015), subjective tinnitus (Chen et al., 2015b) and relapsing-remitting multiple sclerosis (Wu et al., 2016). In the present study, the T2DM patients exhibited abnormal local brain activity synchronization in the bilateral ACG, right FFG, PreCG and medial orbital part of the SFG, as detected by ReHo analysis. Using FC analysis, we further found aberrant brain activity synchronization between the frontal and parietal lobes. These abnormalities suggested that the T2DM patients were subjected to disease attacks, which could lead to brain dysfunction from the local to global scales.

Although the T2DM patients’ general cognition assessed by MMSE and MoCA was not significantly decreased, they indeed suffered from cognitive impairment according to the current rs-fMRI study and neuropsychological test results. Compared to healthy controls, the T2DM patients had decreased ReHo values, and the decreased ReHo values in the right PreCG were associated with lower DST forward scores. These decreases could be harmful to brain function, according to previous studies (Cui et al., 2014; Peng et al., 2016). Furthermore, evidence from a recent functional MRI study showed that the precuneus, which is the central node of the posterior DMN, exhibited decreased connectivity in T2DM patients by local brain independent component analysis (Cui et al., 2015). In the present study, the changed FC between the ACG and the right precuneus might suggest the disassociation of anterior and posterior components of the DMN, which is a major contributor to normal cognitive function (Cui et al., 2015).

More specifically, our findings provide information about the domains of cognitive impairment. The aberrant FC observed between the ACG and right ANG was consistent with previous studies (Zhou et al., 2010; Mormino et al., 2011; Hoogenboom et al., 2014). Decreased FC between the ACG and left hippocampus has been observed in elderly T2DM patients (Zhou et al., 2010) and is also apparent in middle-aged patients. These changed brain areas have been related to episodic memory (Zhou et al., 2010; Mormino et al., 2011), suggesting that episodic memory might be a vulnerable cognitive domain in T2DM patients. Moreover, our findings suggest dysfunction of working memory in T2DM patients. Brain regions, including the supplementary motor area (SMA), parietal areas and hippocampus, which were involved in aberrant FC in the present study, were significantly activated in the spatial n-back task experiment (Huang et al., 2015). Further, with the FC between the right FFG and MTG deviating from the baseline, T2DM patients scored poorer on DST forward tests. This evidence could support the notion that working memory is impaired in T2DM patients. However, the diabetic brain might attempt to maintain cognitive function to overcome the attack by T2DM on memory. The ACG has been reported to have increased neural intensity in the resting-state (Cui et al., 2014) or during a working memory task (He et al., 2015) in T2DM patients. Higher ReHo in the ACG might represent compensatory neural activities for maintaining normal cognition (He et al., 2015). In addition, poor performance on the TMT-A and TMT-B indicated that T2DM patients exhibited executive, attention and psychomotor speed dysfunctions, consistent with previous studies (Chen et al., 2014, 2015a; Cui et al., 2014, 2015). The PreCG is the site of the primary motor cortex, and it is reportedly atrophic in T2DM patients (Chen et al., 2012). Thus, deactivation of the PreCG might account for the slower drawing speed on the TMT-A and TMT-B tests.

Emotion and behavior were also affected by T2DM because certain changed brain regions were involved in these functions, except for cognition. T2DM patients often have treatment-related, recurring periods of low blood glucose levels, which are associated with negative mood states (Anderson et al., 2001). The SFG, SMA, insula and FFG are emotion-associated brain regions vulnerable to blood glucose fluctuations, as concluded from a mood induction experiment (Kohn et al., 2015). Although we excluded subjects with depression, the fMRI results indicated that the remaining patients might be at risk of developing depressive symptoms. Moreover, the insula performs an integrative function in the olfacto-gustatory system (Kurth et al., 2010), and it plays an important role in flavor identification (Rolls, 2006). The insula morphology and energy variability might underlie dysfunction of food intake regulation (Jauch-Chara et al., 2015). Therefore, the altered FC between the ACG and left insula could reflect the ingestion behavior response to blood glucose level changes, which contribute to the typical symptoms of T2DM patients.

According to the results, we attempted to depict the patterns of FC alterations in middle-aged T2DM patients. Except for the stronger negative correlation between the ACG and left ANG, the FC of most abnormal brain regions manifested inverse time course correlations compared to healthy controls, including inverse negative and positive correlations. We inferred that the weakened negative and positive correlations would finally develop into inverse connectivity patterns as the disease progresses. Weaker negative correlations between the right PreCG and bilateral precuneus have also been seen in Alzheimer’s disease (Wang et al., 2007), in agreement with our findings. This finding suggested that these two diseases might share similar mechanisms. Few studies have reported these change patterns in T2DM patients. Either inverse or weakened connectivity patterns tend to be associated with poorer cognitive performance, which might reflect the detrimental effects of T2DM on brain function. Cognitive training can hopefully strengthen the positive or negative correlations of FC as well as reverse the aberrant, reversed positive or negative correlations (Cao et al., 2016). It provides a rewarding method to intervene in cognitive impairment related to T2DM. Moreover, the pattern of FC might be a meaningful imaging marker to detect brain dysfunction and to monitor the results of interventions in T2DM patients.

T2DM is a metabolic disorder and is often accompanied by abnormalities in certain clinical variables. It is difficult to eliminate the influences of these clinical variables on the brain’s functional activity. In the current study, HbA_1c levels were positively correlated with the FC between the right orbital part and the medial orbital of the SFG. Our findings were similar to those of a previous study (Hoogenboom et al., 2014). In this study, a positive correlation was found between HbA_1c and uncinate fasciculus fractional anisotropy in both the T2DM patients and healthy controls. However, higher uncinate fasciculus fractional anisotropy was associated with increased information processing speed (Hoogenboom et al., 2014), thus appearing to contradict the idea that T2DM patients may benefit from blood glucose control to prevent cognitive decline (Kravitz et al., 2013; Biessels et al., 2014). However, HbA_1c only reflects the ambient glucose concentration of the previous 2 to 3 months (O’Keeffe et al., 2016). The duration of FC change may be longer than 2 to 3 months. Therefore, HbA_1c might not be the best predictor of the brain’s condition, or the positive relationship might be attributable to a compensatory mechanism in middle-aged patients (He et al., 2015). In summary, this question must be studied further.

Our study has several limitations. First, the relatively small population size might have restricted the ability to investigate the full extent of ReHo and FC. Therefore, the results should be interpreted cautiously. Second, we tried to achieve inter/intra-group homogeneity except for T2DM-related parameters, but inter-group differences existed in terms of HDL cholesterol, serum creatine and FT3 levels, and there was heterogeneity in the drug therapy of the T2DM patients. The influences of these factors on brain function must be studied further. Third, we adopted a traditional method to analyze the FC, which assumed that the connections remained static in the period of MRI data collection. However, dynamic analysis, a method that has recently emerged, revealed the transient connectivity of the brain, which could improve the temporal resolution of fMRI and capture spontaneously jumping changes in thoughts over time (Hutchison et al., 2013).

Conclusion

This study investigated the spatially local and remote brain activity synchronization in middle-aged T2DM patients. We found that aberrant ReHo and FC were mainly distributed in the insula, MTG and frontal and parietal lobes, which are related to cognitive impairment. It is noteworthy that the patterns of ReHo and FC alterations might serve as imaging biomarkers for the diagnosis and evaluation of the effects of diabetes on the brain. Our findings shed light on the mechanisms underlying T2DM-associated brain dysfunction.

Author Contributions

DL contributed to performing the experiments, data analysis and writing of the manuscript. SD contributed to performing the experiments and revising the manuscript. JZ designed the experiment and revised the manuscript. CZ contributed to the data collection. ML contributed to the data analysis and manuscript revision. XY contributed to the manuscript revision. PW and JW are the guarantors of this work and, as such, had full access to all of the data in the study, and they accept responsibility for the integrity of the data and the accuracy of the data analysis.

Funding

The study was supported by the National Natural Science Foundation of China (81471647) and the Innovation Fund for Younger Investigators of Southwest Hospital of the Third Military Medical University (SWH2013QN09).

Conflict of Interest Statement

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.

Acknowledgments

The authors thank all of the volunteers who participated in the study and the staffs of the Department of Laboratory Medicine and Department of Nuclear Medicine, the Southwest Hospital of the Third Military Medical University, Chongqing, China for their selfless and valuable assistance.

Footnotes

1. ^ www.idf.org/publications/annual-report
2. ^ www.fil.ion.ucl.ac.uk/spm
3. ^ www.mathworks.com
4. ^ www.restfmri.net/forum/DPARSF
5. ^ www.restfmri.net

References

Alberti, K. G. M. M., and Zimmet, P. Z. (1998). Definition, diagnosis and classification of diabetes mellitus and its complications. Part 1: diagnosis and classification of diabetes mellitus provisional report of a WHO consultation. Diabet. Med. 15, 539–553. doi: 10.1002/(SICI)1096-9136(199807)15:7<539::AID-DIA668>3.0.CO;2-S

PubMed Abstract | CrossRef Full Text | Google Scholar

American Diabetes Association. (2015). Classification and diagnosis of diabetes. Diabetes Care 38, S8–S16. doi: 10.2337/dc15-s005

CrossRef Full Text | Google Scholar

Anderson, R. J., Freedland, K. E., Clouse, R. E., and Lustman, P. J. (2001). The prevalence of comorbid depression in adults with diabetes: a meta-analysis. Diabetes Care 24, 1069–1078. doi: 10.2337/diacare.24.6.1069

PubMed Abstract | CrossRef Full Text | Google Scholar

Biessels, G. J., Strachan, M. W., Visseren, F. L., Kappelle, L. J., and Whitmer, R. A. (2014). Dementia and cognitive decline in type 2 diabetes and prediabetic stages: towards targeted interventions. Lancet Diabetes Endocrinol. 2, 246–255. doi: 10.1016/S2213-8587(13)70088-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Biswal, B., Yetkin, F. Z., Haughton, V. M., and Hyde, J. S. (1995). Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magn. Reson. Med. 34, 537–541. doi: 10.1002/mrm.1910340409

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowie, C. R., and Harvey, P. D. (2006). Administration and interpretation of the trail making test. Nat. Protoc. 1, 2277–2281. doi: 10.1038/nprot.2006.390

PubMed Abstract | CrossRef Full Text | Google Scholar

Brucki, S. M. D., and Rocha, M. S. G. (2004). Category fluency test: effects of age, gender and education on total scores, clustering and switching in Brazilian Portuguese-speaking subjects. Braz. J. Med. Biol. Res. 37, 1771–1777. doi: 10.1590/s0100-879x2004001200002

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, W., Cao, X., Hou, C., Li, T., Cheng, Y., Jiang, L., et al. (2016). Effects of cognitive training on resting-state functional connectivity of default mode, salience and central executive networks. Front. Aging Neurosci. 8:70. doi: 10.3389/fnagi.2016.00070

PubMed Abstract | CrossRef Full Text | Google Scholar

Chao-Gan, Y., and Yu-Feng, Z. (2010). DPARSF: a MATLAB toolbox for “pipeline” data analysis of resting-state fMRI. Front. Syst. Neurosci. 4:13. doi: 10.3389/fnsys.2010.00013

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y.-C., Jiao, Y., Cui, Y., Shang, S.-A., Ding, J., Feng, Y., et al. (2014). Aberrant brain functional connectivity related to insulin resistance in type 2 diabetes: a resting-state fMRI study. Diabetes Care 37, 1689–1696. doi: 10.2337/dc13-2127

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Z., Li, L., Sun, J., and Ma, L. (2012). Mapping the brain in type II diabetes: voxel-based morphometry using DARTEL. Eur. J. Radiol. 81, 1870–1876. doi: 10.1016/j.ejrad.2011.04.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y.-C., Liu, Z., Zhang, J., Tian, G., Li, L., Zhang, S., et al. (2015a). Selectively disrupted functional connectivity networks in type 2 diabetes mellitus. Front. Aging Neurosci. 7:233. doi: 10.3389/fnagi.2015.00233

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y.-C., Zhang, J., Li, X.-W., Xia, W., Feng, X., Qian, C., et al. (2015b). Altered intra- and interregional synchronization in resting-state cerebral networks associated with chronic tinnitus. Neural Plast. 2015:475382. doi: 10.1155/2015/475382

PubMed Abstract | CrossRef Full Text | Google Scholar

Cui, Y., Jiao, Y., Chen, H.-J., Ding, J., Luo, B., Peng, C.-Y., et al. (2015). Aberrant functional connectivity of default-mode network in type 2 diabetes patients. Eur. Radiol. 25, 3238–3246. doi: 10.1007/s00330-015-3746-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Cui, Y., Jiao, Y., Chen, Y.-C., Wang, K., Gao, B., Wen, S., et al. (2014). Altered spontaneous brain activity in type 2 diabetes: a resting-state functional MRI study. Diabetes 63, 749–760. doi: 10.2337/db13-0519

PubMed Abstract | CrossRef Full Text | Google Scholar

Cui, L. B., Liu, K., Li, C., Wang, L. X., Guo, F., Tian, P., et al. (2016). Putamen-related regional and network functional deficits in first-episode schizophrenia with auditory verbal hallucinations. Schizophr. Res. 173, 13–22. doi: 10.1016/j.schres.2016.02.039

PubMed Abstract | CrossRef Full Text | Google Scholar

Damoiseaux, J. S., Beckmann, C. F., Sanz Arigita, E. J., Barkhof, F., Scheltens, Ph., Stam, C. J., et al. (2008). Reduced resting-state brain activity in the “default network” in normal aging. Cereb. Cortex 18, 1856–1864. doi: 10.1093/cercor/bhm207

PubMed Abstract | CrossRef Full Text | Google Scholar

Fair, D. A., Cohen, A. L., Power, J. D., Dosenbach, N. U. F., Church, J. A., Miezin, F. M., et al. (2009). Functional brain networks develop from a “local to distributed” organization. PLoS Comput. Biol. 5:e1000381. doi: 10.1371/journal.pcbi.1000381

PubMed Abstract | CrossRef Full Text | Google Scholar

Galea, M., and Woodward, M. (2005). Mini-Mental State Examination (MMSE). Aust. J. Physiother. 51:198. doi: 10.1016/S0004-9514(05)70034-9

CrossRef Full Text | Google Scholar

Gong, Y. (1992). Wechsler Adult Intelligence Scale-Revised in China Version. Changsha, Hunan/China: Hunan Medical College.

Hamilton, M. (1960). A rating scale for depression. J. Neurol. Neurosurg. Psychiatry 23, 56–62. doi: 10.1136/jnnp.23.1.56

PubMed Abstract | CrossRef Full Text | Google Scholar

He, X.-S., Wang, Z.-X., Zhu, Y.-Z., Wang, N., Hu, X., Zhang, D.-R., et al. (2015). Hyperactivation of working memory-related brain circuits in newly diagnosed middle-aged type 2 diabetics. Acta Diabetol. 52, 133–142. doi: 10.1007/s00592-014-0618-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoogenboom, W. S., Marder, T. J., Flores, V. L., Huisman, S., Eaton, H. P., Schneiderman, J. S., et al. (2014). Cerebral white matter integrity and resting-state functional connectivity in middle-aged patients with type 2 diabetes. Diabetes 63, 728–738. doi: 10.2337/db13-1219

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, R.-R., Jia, B.-H., Xie, L., Ma, S. H., Yin, J. J., Sun, Z. B., et al. (2015). Spatial working memory impairment in primary onset middle-age type 2 diabetes mellitus: an ethology and BOLD-fMRI study. J. Magn. Reson. Imaging 43, 75–87. doi: 10.1002/jmri.24967

PubMed Abstract | CrossRef Full Text | Google Scholar

Hutchison, R. M., Womelsdorf, T., Allen, E. A., Bandettini, P. A., Calhoun, V. D., Corbetta, M., et al. (2013). Dynamic functional connectivity: promise, issues and interpretations. Neuroimage 80, 360–378. doi: 10.1016/j.neuroimage.2013.05.079

PubMed Abstract | CrossRef Full Text | Google Scholar

Jauch-Chara, K., Binkofski, F., Loebig, M., Reetz, K., Jahn, G., Melchert, U. H., et al. (2015). Blunted brain energy consumption relates to insula atrophy and impaired glucose tolerance in obesity. Diabetes 64, 2082–2091. doi: 10.2337/db14-0421

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, L., Xu, Y., Zhu, X.-T., Yang, Z., Li, H. J., and Zuo, X. N. (2015). Local-to-remote cortical connectivity in early- and adulthood-onset schizophrenia. Transl. Psychiatry 5:e566. doi: 10.1038/tp.2015.59

PubMed Abstract | CrossRef Full Text | Google Scholar

Koekkoek, P. S., Kappelle, L. J., van den Berg, E., Rutten, G. E. H. M., and Biessels, G. J. (2015). Cognitive function in patients with diabetes mellitus: guidance for daily care. Lancet Neurol. 14, 329–340. doi: 10.1016/s1474-4422(14)70249-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Kohn, N., Toygar, T., Weidenfeld, C., Berthold-Losleben, M., Chechko, N., Orfanos, S., et al. (2015). In a sweet mood? Effects of experimental modulation of blood glucose levels on mood-induction during fMRI. Neuroimage 113, 246–256. doi: 10.1016/j.neuroimage.2015.03.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Kravitz, E., Schmeidler, J., and Schnaider Beeri, M. (2013). Type 2 diabetes and cognitive compromise: potential roles of diabetes-related therapies. Endocrinol. Metab. Clin. North Am. 42, 489–501. doi: 10.1016/j.ecl.2013.05.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Kurth, F., Zilles, K., Fox, P. T., Laird, A. R., and Eickhoff, S. B. (2010). A link between the systems: functional differentiation and integration within the human insula revealed by meta-analysis. Brain Struct. Funct. 214, 519–534. doi: 10.1007/s00429-010-0255-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Mormino, E. C., Smiljic, A., Hayenga, A. O., Onami, S. H., Greicius, M. D., Rabinovici, G. D., et al. (2011). Relationships between β-amyloid and functional connectivity in different components of the default mode network in aging. Cereb. Cortex 21, 2399–2407. doi: 10.1093/cercor/bhr025

PubMed Abstract | CrossRef Full Text | Google Scholar

Musen, G., Jacobson, A. M., Bolo, N. R., Simonson, D. C., Shenton, M. E., McCartney, R. L., et al. (2012). Resting-state brain functional connectivity is altered in type 2 diabetes. Diabetes 61, 2375–2379. doi: 10.2337/db11-1669

PubMed Abstract | CrossRef Full Text | Google Scholar

O’Keeffe, D. T., Maraka, S., and Rizza, R. A. (2016). HbA_1c in the evaluation of diabetes mellitus. JAMA 315, 605–606. doi: 10.1001/jama.2015.16561

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, J., Qu, H., Peng, J., Luo, T.-Y., Lv, Fa. J., Chen, L., et al. (2016). Abnormal spontaneous brain activity in type 2 diabetes with and without microangiopathy revealed by regional homogeneity. Eur. J. Radiol. 85, 607–615. doi: 10.1016/j.ejrad.2015.12.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Rolls, E. T. (2006). Brain mechanisms underlying flavour and appetite. Philos. Trans. R. Soc. Lond. B Biol. Sci. 361, 1123–1136. doi: 10.1098/rstb.2006.1852

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmidt, M. (1996). Rey Auditory Verbal Learning Test: A Handbook. Los Angeles, CA: Western Psychological Services.

Google Scholar

Song, X.-W., Dong, Z.-Y., Long, X.-Y., Li, S.-F., Zuo, X.-N., Zhu, C.-Z., et al. (2011). REST: a toolkit for resting-state functional magnetic resonance imaging data processing. PLoS One 6:e25031. doi: 10.1371/journal.pone.0025031

PubMed Abstract | CrossRef Full Text | Google Scholar

Strachan, M. W. J., Reynolds, R. M., Marioni, R. E., and Price, J. F. (2011). Cognitive function, dementia and type 2 diabetes mellitus in the elderly. Nat. Rev. Endocrinol. 7, 108–114. doi: 10.1038/nrendo.2010.228

PubMed Abstract | CrossRef Full Text | Google Scholar

Wahlund, L. O., Barkhof, F., Fazekas, F., Bronge, L., Augustin, M., Sjögren, M., et al. (2001). A new rating scale for age-related white matter changes applicable to MRI and CT. Stroke 32, 1318–1322. doi: 10.1161/01.STR.32.6.1318

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, K., Liang, M., Wang, L., Tian, L., Zhang, X., Li, K., et al. (2007). Altered functional connectivity in early Alzheimer’s disease: a resting-state fMRI study. Hum. Brain Mapp. 28, 967–978. doi: 10.1002/hbm.20324

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, P., Li, R., Yu, J., Huang, Z., and Li, J. (2016). Frequency-dependent brain regional homogeneity alterations in patients with mild cognitive impairment during working memory state relative to resting state. Front. Aging Neurosci. 8:60. doi: 10.3389/fnagi.2016.00060

PubMed Abstract | CrossRef Full Text | Google Scholar

Whitwell, J. L. (2009). Voxel-based morphometry: an automated technique for assessing structural changes in the brain. J. Neurosci. 29, 9661–9664. doi: 10.1523/JNEUROSCI.2160-09.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, L., Zhang, Y., Zhou, F., Gao, L., He, L., Zeng, X., et al. (2016). Altered intra- and interregional synchronization in relapsing-remitting multiple sclerosis: a resting-state fMRI study. Neuropsychiatr. Dis. Treat. 12, 853–862. doi: 10.2147/NDT.s98962

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, W., Wang, S., Sun, Z., Bai, F., Zhou, Y., Yang, Y., et al. (2013). Altered baseline brain activity in type 2 diabetes: a resting-state fMRI study. Psychoneuroendocrinology 38, 2493–2501. doi: 10.1016/j.psyneuen.2013.05.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Zang, Y., Jiang, T., Lu, Y., He, Y., and Tian, L. (2004). Regional homogeneity approach to fMRI data analysis. Neuroimage 22, 394–400. doi: 10.1016/j.neuroimage.2003.12.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, H., Lu, W., Shi, Y., Bai, F., Chang, J., Yuan, Y., et al. (2010). Impairments in cognition and resting-state connectivity of the hippocampus in elderly subjects with type 2 diabetes. Neurosci. Lett. 473, 5–10. doi: 10.1016/j.neulet.2009.12.057

PubMed Abstract | CrossRef Full Text | Google Scholar