Research Paper Advance Articles
Differential impact of glycemic control and comorbid conditions on the neurophysiology underlying task switching in older adults with type 2 diabetes
- 1 Institute for Human Neuroscience, Boys Town National Research Hospital, Boys Town, NE 68010, USA
- 2 Department of Psychology, University of Nebraska Omaha, Omaha, NE 68182, USA
- 3 Department of Internal Medicine, Division of Diabetes, Endocrinology and Metabolism, University of Nebraska Medical Center, Omaha, NE 68198, USA
Received: April 11, 2022 Accepted: June 4, 2022 Published: June 17, 2022https://doi.org/10.18632/aging.204129
How to Cite
Copyright: © 2022 Embury et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Type 2 diabetes is known to negatively affect higher order cognition and the brain, but the underlying mechanisms are not fully understood. In particular, glycemic control and common comorbidities are both thought to contribute to alterations in cortical neurophysiology in type 2 diabetes, but their specific impact remains unknown. The current study probed the dynamics underlying cognitive control in older participants with type 2 diabetes, with and without additional comorbid conditions (i.e., cardiovascular disease, nephropathy, peripheral neuropathy, retinopathy), using a task switching paradigm and a dynamic functional brain mapping method based on magnetoencephalography (MEG). We hypothesized that neural dynamics would be differentially impacted by the level of glycemic control (i.e., diabetes itself) and the burden of additional comorbid conditions. Supporting this hypothesis, our findings indicated separable, but widespread alterations across frontal, parietal, temporal and cerebellum regions in neural task-switch costs in type 2 diabetes that were differentially attributable to glycemic control and the presence of comorbid conditions. These effects were spatially non-overlapping and the effects were not statistically related to one another. Further, several of the effects that were related to the presence of comorbidities were associated with behavioral performance, indicating progressive deficits in brain function with extended disease. These findings provide insight on the underlying neuropathology and may inform future treatment plans to curtail the neural impact of type 2 diabetes.
Type 2 diabetes is a chronic metabolic disorder characterized by progressive insulin resistance that is often accompanied by a collection of comorbid conditions, many of which are also progressive. Common comorbid conditions include retinopathy, micro- and macro-vascular disease, peripheral neuropathy, and nephropathy. Previous studies have also found significant cognitive and neural deficits in people with type 2 diabetes. Specific cognitive domains affected include memory, attention, and executive functioning, with effect sizes ranging from 0.25–0.5 [1, 2], although moderate to large effect sizes have been found in people over 65 years-old , suggesting these decrements increase with aging. Studies focusing on the neurological impact of type 2 diabetes have shown widespread impact, including deficits in grey matter volume and white matter integrity [1, 4–7], with some studies linking these declines to cognitive outcomes and an increased risk for dementia [6, 8]. Likewise, functional imaging studies have shown alterations in resting state functional connectivity and task-based neural responses across widespread networks, including the default mode and frontal-parietal networks [9–12]. Increased numbers of white matter hyperintensities and greater cerebrovascular damage have also been noted in people with type 2 diabetes [1, 10, 13].
Executive functioning is one cognitive domain particularly affected in type 2 diabetes [1, 2]. One task known to tap executive function is task switching, which requires participants to flexibly adjust goal sets and behavior to match contextual cues . Behaviorally, switching rule sets requires longer reaction times relative to stay or no switch trial types, revealing a consistent behavioral switch cost . These paradigms have generally been found to elicit activity across the frontal-parietal and cingulo-opercular networks . Previous EEG studies have also shown greater alpha, beta and theta recruitment in switch, relative to same set trials, reflecting neural switch costs [17, 18]. An MEG study has also found a role for gamma in switch tasks, where neural switch costs were greater across frontal-parietal and cingulate regions . However, such neural dynamics have not yet been probed in the context of type 2 diabetes.
While studies have shown that type 2 diabetes affects executive function, whether the driving force behind these deficits is primarily related to dysglycemia or the increased number of comorbid conditions in those with type 2 diabetes remains to be elucidated. Studies have found mediating effects of high body mass index, age, microvascular damage, and hypertension, as well as glycated hemoglobin (HbA1c) and other metrics of glycemic control on the cognitive decrements associated with type 2 diabetes [1, 12, 20]. Insulin resistance, a hallmark of the disease, has a direct impact on the brain. Insulin receptors are found across the cortex and subcortical structures including the hippocampus, hypothalamus, and amygdala [21, 22], and insulin signaling has been implicated in several cognitive processes, most notably episodic memory performance . Both glycemic clamp and insulin administration studies show immediate improvement in cognition with normalized glucose levels, especially within the memory domain [23–25].
In the current study, we used a task switching paradigm  and MEG to probe the neural dynamics and circuitry serving this essential component of executive function in a cohort of older adults with type 2 diabetes, about half of whom had common comorbidities. We hypothesized that both glycemic control and comorbid conditions would alter oscillatory activity across relevant networks, including occipital, parietal, and prefrontal regions, suggesting an accelerated aging phenotype with worse glycemic control and the presence of comorbidities. We further posited that the effects of glycemic dysregulation and comorbidities would be separable and that some effects would relate to behavioral switch costs.
Demographic, behavioral and disease status results
Four participants had to be excluded due to MEG technical issues or artifactual MEG data. Two additional participants were excluded due to behavioral performance at or near chance levels on the task switching paradigm (see Figure 1A). This resulted in a sample of 48 participants (25 with significant comorbidities, 23 without). See demographics and blood panel values in Table 1. Groups did not differ by age (p = 0.497), sex (p = 0.951), race (p = 0.845), education (p = 0.132), handedness (p = 0.502), BMI (p = 0.180), and HbA1c (p = 0.205). Of those with comorbidities, 11 had nephropathy (GFR 45–60, ACR 41.45–375), five had retinopathy (mild), 14 had peripheral neuropathy (generally mild, limited to specific extremities), and six had cardiovascular disease (three stents, three bypasses), with nine (36%) of these participants having two or more of these conditions. Participants performed well, with an average accuracy of 95.5 ± 3.6%, average reaction time of 1527.67 ± 216.02 ms, and an average behavioral switch cost, or the difference in reaction times between conditions, of 79.27 ± 96.37 ms. Accuracy did not differ between groups (t46 = −0.88, p = 0.381), but there was a trend for behavioral switch costs differences by group (t46 = 1.79, p = 0.080, see Figure 1B), as those without comorbidities were less affected by the switch.
Figure 1. Task switch paradigm and behavioral switch costs. (A) A fixation cross was presented for 2700 ± 100 ms, followed by a number (possible values: 1–4, 6–9) surrounded by a square or diamond, indicating the rule set for that trial, presented for 2500 ms. For trials with a number within a square, participants had to respond by button press as to whether the number was below (index finger) or above (middle) 5. For trials with a number within a diamond, participants had to answer whether the number was odd (index finger) or even (middle) by button press. Trials were pseudorandomized such that 50% of trials were the same rule set as the previous trial (No Switch) and 50% were a different rule set from the previous trial (Switch). All analyses were relative to the difference between these conditions, or switch costs. (B) Accuracy did not differ by group (not shown), but there was a trend for group differences in behavioral switch costs (t46 = 1.79, p = 0.080), where those with type 2 diabetes and additional comorbidities (shown in green) had longer reaction times relative to those without comorbidities (shown in blue). # denotes 0.050 < p < 0.100.
Table 1. Task switch: demographics and laboratory tests.
|Type 2 diabetes with comorbidities||Type 2 diabetes without comorbidities|
|Age (years)||64.0 ± 4.8||62.7 ± 5.1|
|Sex||10 Males; 15 Females||8 Males; 15 Females|
|Handedness||24 Right; 1 Left||21 Right; 2 Left|
|Disease duration (years)||14.1 ± 2.0||10.0 ± 7.5|
|HbA1c (mmol/mol)||59 ± 11.7||55 ± 9.2|
|HbA1c (%)||7.52 ± 1.07||7.16 ± 0.84|
|Creatinine (mg/dL)||0.84 ± 0.22||0.83 ± 0.23|
|Glucose (mg/dL)||119.52 ± 22.75||121.91 ± 41.01|
|Albumin/Creatinine (ugAL/mgCR)||45.08 ± 81.41||17.17 ± 13.40|
|Thyroid-stimulating hormone (mcIU/mL)||2.72 ± 1.72||1.77 ± 0.86|
|B12 (pg/mL)||451.45 ± 361.81||343.08 ± 222.00|
|Values depict means ± SD.|
Sensor level results
To identify significant oscillatory responses for dynamic imaging, t-tests between the active and baseline periods per time-frequency bin were computed across all participants and both conditions, which were then corrected for multiple comparisons using nonparametric permutation testing. Significant time-frequency windows were found in alpha (8–12 Hz) from 350–1050 ms and in gamma (40–52 Hz, 62–74 Hz) from 150–650 ms (p < 0.001; Figure 2). To study early and late dynamics of task switch processing, alpha was imaged from 350–700 and 700–1050, while both gamma frequency bins were imaged from 150–400 and 400–650 ms. To identify the anatomical origins of these oscillatory responses, images per response window were grand averaged across conditions and participants and this revealed that alpha responses originated in parieto-occipital regions throughout the task period, while both low (40–52 Hz) and high (62–74 Hz) gamma were strongest in bilateral occipital cortices (Figure 2).
Figure 2. Task switch spectrograms. Left: Grand-averaged time-frequency spectrograms are shown, derived from representative parieto-occipital MEG sensors. Time is shown on the x-axis in seconds, frequency is shown on the y-axis in Hz. Spectrograms are shown from 4–30 Hz (bottom), 30–60 Hz (middle), and 60–90 Hz (top). The colors reflect power increases (red) and decreases (blue) relative to the baseline, with the scale bar shown to the right of each time-frequency plot. Time-frequency windows for source imaging (beamforming) were derived from statistical analysis of the sensor-level data across all participants (ps < 0.001). A clear alpha (8–12 Hz) decrease can be seen throughout the task period, as well as a lower gamma (40–52 Hz) and a higher gamma (62–74 Hz) synchronization during the task period. Right: Grand average images from each time-frequency bin show the alpha dynamics in parietal and occipital regions, and gamma dynamics largely restricted to bilateral occipital cortices.
Source level results
Neural switch cost differences by patient groups
In order to examine neural switch cost effects, whole brain subtraction maps were computed (switch – no-switch) and these maps were compared using independent-sample t-tests to identify group effects (patients with comorbidities vs. patients without comorbidities). As shown in Figure 3, our results revealed high gamma (62–74 Hz) differences from 150–400 ms in the left cerebellum (t45 = 3.32, p = 0.002), left inferior frontal (t45 = 3.78, p < 0.001), left superior frontal (t45 = 3.70, p < 0.001), and right inferior parietal (t45 = 3.76, p < 0.001) cortices. In the later 400–650 ms window, high gamma group differences were found in the left precentral (t45 = 4.55, p < 0.001), left occipital (t45 = 3.53, p < 0.001), left cerebellum (t45 = 3.68, p < 0.001), and right parieto-occipital (t45 = −3.03, p = 0.004) regions. In all cases except the right parieto-occipital, the group with comorbidities exhibited stronger high gamma oscillations during task performance (Figure 3). No group differences were detected in low gamma. As per alpha, a significant group difference in neural switch cost was found in the right parietal lobe during the 350–700 ms time window (t44 = −3.42, p = 0.001) such that those with comorbidities had greater switch costs in this region, while no group differences were found in the later alpha time window. Notably, none of these effects were associated with HbA1c level (ps from 0.436 to 0.797).
Figure 3. Neural switch costs differ by comorbidity status. Left: Significant group effects on neural switch costs were found in both time bins of high gamma, including the left cerebellum, left inferior (not shown) and superior frontal regions and the right parietal cortex in the early high gamma window, and left precentral, left occipital, left cerebellum, and right parieto-occipital regions in late high gamma activity. Right: Greater neural switch costs in early bin parietal alpha activity were found in type 2 diabetes patients with additional comorbidities, relative to those without comorbidities. These effects show the impact of comorbidities on neural switch costs across the spectrum. Images are thresholded from p < 0.01 to p < 0.001, with a minimum cluster size of 20 4 mm3 voxels.
Neuro-behavioral switch cost associations
Peak values were extracted from each of the significant group difference clusters and correlations with behavior were computed. Significant correlations across both groups were found between neural and behavioral switch costs. Specifically, in the early high gamma window, a significant relationship between the left cerebellar peak and behavioral switch cost was apparent (r45 = 0.34, p = 0.020; Figure 4A), such that as the behavioral switch cost became larger, so too did the neural switch costs in this region. In the later high gamma window, a significant inverse relationship between the right parieto-occipital peak and behavioral switch cost emerged (r45 = −0.35, p = 0.016), as increased behavioral switch cost was associated with reduced neural switch costs (Figure 4B). Because of the differing trajectories of these responses across time, we probed the relationship between neural switch costs in these two regions and found a significant reciprocal effect, whereby increased switch costs in one region tracked with decreased switch costs in the other. To determine whether this reciprocal pattern was related to group, a repeated-measures ANCOVA was computed to determine the relationship between neural and behavioral switch costs and group identity. In this model, a significant interaction between neural switch costs at each peak and behavioral switch cost was again apparent (F1,43 = 4.31, p = 0.044). Further, controlling for the effect of behavior, there was a significant interaction of peak by group, such that each group tended to have greater neural switch costs in one peak over the other to complete the task (F1,43 = 5.52, p = 0.024, see Figure 4C). In this case, participants with type 2 diabetes and comorbidities tended to have higher switch costs in the early bin left cerebellum peak and lower switch costs in the late bin right parieto-occipital peak, whereas participants with type 2 diabetes and no comorbidities tended to have higher switch costs in the late bin right parieto-occipital peak and lower switch costs in the early bin left cerebellum peak. No significant interaction of group by behavioral switch costs (F1,43 = 0.12, p = 0.735), nor peak by group by behavioral switch cost (F1,43 = 0.04, p = 0.851) was found. These findings suggest a differential impact of comorbidities on neural oscillatory responses serving higher order cognition in type 2 diabetes, but group membership did not specifically predict behavioral outcomes.
Figure 4. Neural and behavioral switch cost correlations. Significant correlations between neural and behavioral switch costs were apparent in the early high gamma left cerebellar peak ((A) r45 = 0.34, p = 0.020) and the late high gamma right parieto-occipital peak ((B) r45 = −0.35, p = 0.016). In the left cerebellum, greater neural switch costs were associated with greater behavioral switch costs, while in the right parieto-occipital lower neural switch costs were related to greater behavioral switch costs. This distinct pattern of neural responses leading to differential behavioral outcomes was further investigated, and this analysis showed group-specific recruitment of each region ((C) F1,43 = 5.52, p = 0.024), whereby those with comorbid conditions (shown in green) had greater switch costs in left cerebellum, while those without comorbidities (shown in blue) had greater switch costs in the right parieto-occipital region. *p < 0.05, **p < 0.01, ***p < 0.001.
Neural switch cost correlations with clinical metrics
To investigate the relationship between switch costs and glycemic control, we ran whole brain correlations on the neural switch cost maps. Several areas correlated with HbA1c level in the late high gamma window including the left parietal (r45 = 0.48, p < 0.001), left lateral occipital (r45 = 0.48, p < 0.001), right inferior occipital (r45 = 0.46, p = 0.001), and left inferior temporal (r45 = 0.40, p = 0.005) regions. These effects were such that as HbA1c increased, high gamma switch costs also increased in each region (Figure 5A–5D). A significant correlation with HbA1c was also apparent in the early low gamma window in the left occipital cortex (r45 = −0.52, p < 0.001), such that as HbA1c levels increased, low gamma switch costs decreased (Figure 5E). None of these effects differed by comorbidity distinction (ps from 0.386 to 0.965). These results suggest largely distinct effects of glycemic dysregulation and the presence of comorbid conditions.
Figure 5. Glycemic control and neural switch cost associations. Significant correlations between neural switch cost dynamics and glycemic control level were apparent in the late high gamma and late low gamma windows. Specifically, gamma responses in the left lateral occipital ((A) r45 = 0.48, p < 0.001), left parietal ((B) r45 = 0.48, p < 0.001), right inferior occipital ((C) r45 = 0.46, p = 0.001), and left inferior temporal ((D) r45 = 0.40, p = 0.005) regions exhibited greater neural switch costs in the context of higher HbA1c levels. Interestingly, greater neural switch costs were associated with lower HbA1c levels in the left occipital region’s late low gamma response ((E) r45 = −0.53, p < 0.001). These results show a differential impact of glycemic control level on neural switch costs, with separable effects from those seen in differences by comorbidity status. *p < 0.05, **p < 0.01, ***p < 0.001.
We found that glycemic dysregulation and the presence of comorbid conditions have largely distinct effects on the neurophysiological activity underlying task switching in a sample of older adults with type 2 diabetes. The impact of comorbid conditions was largely an increase in behavioral switch costs, as well as an increase in the strength of oscillations in the high gamma range across multiple brain regions, including the left superior and inferior frontal regions, left precentral gyrus, left cerebellum, and bilateral occipital areas. Comorbidities were also associated with aberrations in right parietal alpha oscillations. Such accentuated neural switch costs across frontal-parietal network regions suggests that comorbidities may affect hub regions involved in cognitive control, as well as cerebellar activity that has been linked to error monitoring and conflict resolution [26, 27], all of which are subprocesses contributing to task switching performance. Similarly, inhibitory control has also been shown to be diminished in type 2 diabetes [28, 29], so it follows that the underlying neural correlates would be impacted. Regarding glycemic control, increased HbA1c levels were associated with larger neural switch costs (i.e., stronger neural oscillations) in the high gamma range within superior parietal, lateral occipital, inferior temporal, and ventral occipital cortices, as well as smaller neural switch costs in medial occipital cortex in the lower gamma range. Thus, the brain regions most strongly affected by glycemic dysregulation were largely independent of those associated with significant comorbidities. Below, we discuss the implications of these findings for understanding the neural aberrations often seen in those with type 2 diabetes, especially regarding how these relate to decrements in task switching.
Some of our key findings were that comorbid conditions were associated with increases in the strength of high gamma oscillations across the left superior and inferior frontal regions, left precentral gyrus, left cerebellum, and bilateral occipital areas, as well as right parietal alpha. The frequency specificity of these findings likely reflects the underlying processes directly affected by comorbid conditions associated with type 2 diabetes. Specifically, frontal and occipital gamma responses have previously been implicated in tasks probing executive functioning and cognitive control [30–32]. These gamma responses have also been implicated in a wide range of functions including attention, working and long-term memory, and feature binding [33, 34]. Parietal alpha has an important role in the dorsal attention network, where it has been linked to the top-down direction of attention toward task goals . While the widespread nature of effects also suggests a more global impact of comorbid conditions in type 2 diabetes, the relatively large number of findings in the gamma band may suggest effects on specific mechanisms underlying inhibitory processes, such as gamma-aminobutyric acid (GABA) signaling. GABA is one of the key receptors known to be downregulated in the insulin resistant brain . Gamma dynamics are thought to regulated by local GABA-ergic inhibitory circuitry [37, 38]. Indeed, in previous studies in people with type 2 diabetes, GABA levels in occipital and prefrontal cortices were shown to be significantly altered, with some effects directly associated with cognitive task performance [39–42]. These alterations in GABA could be a critical mechanism underlying many of our findings in the current study, as many effects in our task switching paradigm were in the gamma range. Future studies should further examine GABA-gamma links in this population using combined MEG and MRI-based GABA spectroscopy.
We also showed significant associations between neural and behavioral switch costs. In this case, greater neural switch costs in the left cerebellum’s early gamma response was associated with greater behavioral switch costs, while smaller neural switch costs in the right parieto-occipital cortex’s late gamma response was associated with greater behavioral switch costs. Further, these relative switch costs were separable by group, reflecting the unique impact of comorbidities in type 2 diabetes. Of note, participants with comorbidities generally exhibited more aberrant activity in the left cerebellum, where greater neural switch costs led to greater behavioral switch costs. On the other hand, participants without comorbid conditions tended toward the more compensatory pattern, where increased right parieto-occipital activity led to less behavioral switch costs. These patterns of findings suggest that the additional comorbidities lead to an inefficient system underlying cognition and behavior. Previous studies have also shown that the presence of comorbidities, particularly those affecting vascular integrity, have strong detrimental effects on the brain and cognition [1, 3, 43, 44], albeit without the frequency specificity of neural effects we were able to discern in the current study.
Importantly, the peak neural differences between those with and without comorbidities did not spatially overlap with those that varied based on level of glycemic control, as measured by HbA1c. Further, comorbidity status associated effects did not relate to HbA1c level and glycemic control findings did not differ by comorbidity status. Instead, neural switch costs in several other brain regions were found to significantly correlate with HbA1c, including the left occipital cortex in the late low gamma window and left lateral occipital, left parietal, left inferior temporal, and right inferior occipital regions in late high gamma window. Higher HbA1c levels were associated with larger neural switch costs in the high gamma peaks, and weaker neural switch costs in the low gamma range. Higher HbA1c has been previously shown to have worse cognitive outcomes [45, 46], although some studies suggest that the relationship is weak, only accounting for about 10% of the variance in measured decrements . While it is generally true that worse glycemic control generally leads to greater incidence of comorbidities, from our findings, the effects of glycemic control and the presence of comorbid conditions differentially impact the brain and behavior. Further, neural switch cost differences by comorbidity status related more closely to behavior, where in general, greater neural switch costs lead to worse behavior in those with comorbid conditions. Taken together, the relative impact of both worse glycemic control and the presence of comorbidities seems to be widespread, suggesting that better disease management may prevent worse cognitive and neural outcomes. These findings suggest that it is important to separate the effects specific to the disease processes from those that are due to comorbid conditions, particularly since the effects relative to comorbidity status showed specific relationships with behavior.
While this study advances the field’s understanding of the impact of type 2 diabetes and associated comorbid conditions on higher order cognition and the underlying neurophysiology, some limitations need to be acknowledged. These include the characterization of comorbid conditions (yes/no) and the use of HbA1c as a coarse measure of glycemic control. Future studies should examine these variables more in depth to further clarify their causal links to the neural deficits observed. In closing, we found distinct effects of comorbid conditions and glycemic control level on neural switch costs during an executive function task in people with type 2 diabetes. The dynamics of these response differences were largely found within the gamma band, with aberrations in higher frequency dynamics linked to both comorbid conditions and worse glycemic control, although in different brain regions. Neural switch cost differences were related to behavioral switch costs, with the directionality of effects separable by comorbidity status.
These findings add to the growing body of literature illustrating the cognitive decrements and neural dysfunction associated with type 2 diabetes, but are novel in that they identify the separable neurophysiological signatures of comorbidity status and glycemic control level. Future studies should expand on the nature of these effects, including examining whole brain connectivity dynamics, cross-frequency interactions, and structure-function relationships, which may further elucidate the complex nature of how type 2 diabetes affects the brain. Examining individual differences in disease management, comorbid conditions, age, BMI, and other relevant factors and their independent and synergistic effects on the brain may also inform future studies, as the heterogeneity of pathological features and expression may impact cognitive and neural outcomes.
Fifty-four participants with type 2 diabetes (25 without and 29 with significant comorbid conditions) were recruited from the Diabetes Clinic at the University of Nebraska Medical Center (UNMC; age range: 55–73 years, 33 females). Participants had to be formally diagnosed with type 2 diabetes for at least 1 year prior to entering the study. Participants with additional complications and comorbidities were considered as such if any of the following conditions were met: (1) nephropathy, as characterized by glomerular filter rate of 45–60 and/or microalbuminuria of 30–300, (2) retinopathy, excluding proliferative retinopathy, (3) peripheral neuropathy, confirmed by documented symptoms or monofilament testing, and/or (4) cardiovascular disease, except recent major cardiovascular events (heart attack, stroke) in the previous 6 months from enrollment. Exclusionary criteria included: (1) any medical diagnosis affecting brain function (e.g., psychiatric and/or neurological disease), (2) known brain neoplasm or lesion, (3) history of cerebrovascular events (i.e., CVA, stroke, TIA) based on previous diagnosis and chart review, (4) history of significant head trauma, seizures or epilepsy, (5) current substance use disorder within the past 6 months, (6) pregnancy or lactation, (7) hospitalization within the previous 90 days, (8) any type of cancer diagnosis or treatment in the past 5 years, except skin cancer, (9) uncontrolled hypertension, with blood pressures greater than 140/90 or 160/100 if currently on medication treatment, (10) body mass index of greater than or equal to 40, (11) liver disease (AST or ALT > 3× normal), (12) any untreated thyroid or B12 deficiencies, (13) history of brain radiation, (14) treatment with antipsychotics, antidepressants, and related medications known to affect brain function, with the exception of as-needed antidepressants following a 24 hour washout period, and (15) ferromagnetic implants.
To study the long-term effects of type 2 diabetes on the brain in a more controlled manner, relative euglycemia (70 to 200 mg/dl) in these participants at the time of study was a prerequisite for participation, measured by point-of-care device prior to study commencement. Written informed consent was obtained from each participant following the guidelines of UNMC’s Institutional Review Board, who approved the study protocol, in accordance with the Declaration of Helsinki.
Prior to MEG, participants underwent a panel of blood tests according to the standards of care described by the American Diabetes Association . Furthermore, demographic and medical history data were collected. For full lab results and general characteristics of the patient groups, see Table 1. Participants then completed cognitive tasks while undergoing MEG, to probe the neurophysiological changes related to essential cognitive processes in this sample.
Task switch paradigm
During the MEG session, participants were seated in a nonmagnetic chair and instructed to fixate on a crosshair presented centrally for 2600–2800 ms. Following fixation, participants were shown a single centrally-presented integer from 1–4 or 6–9 for 2500 ms, surrounded by either a square or a diamond shape which cued the participant to the appropriate rule set for the trial. If the number was surrounded by a square, the participant had to respond by button press whether the number was less than (index) or greater than (middle) five. If the number was surrounded by a diamond, the participant had to respond by button press whether the number was odd (index) or even (middle). Each rule set occurred 50% of the time, with a pseudo-randomized trial order to have about 50% of all trials repeat the previous rule set (No-Switch) and the remaining trials change rule set from the previous trial (Switch). This ensured participants could not anticipate whether any particular trial would switch or not from the previous rule set. Each trial was 5200 +/− 100 ms, with a total of 200 trials (see Figure 1A). This task has been previously validated by our group . Behavioral metrics were collected concurrently with MEG data, including reaction time and accuracy. Behavioral switch costs, or the difference (i.e., increase) in reaction time between Switch and No-Switch conditions, were computed to characterize the cognitive control component of the task more directly.
MEG methods and analyses
MEG acquisition and analysis methodology followed standardized pipelines, corresponding to MEG studies previously published by our group . Briefly, MEG data were recorded using a 306-sensor Elekta Neuromag MEG system (Elekta/MEGIN, Helsinki, Finland), sampled at 1 kHz with an acquisition bandwidth of 0.1 – 330 Hz. Each participant’s data were corrected for head motion and subjected to noise reduction using the signal space separation method with a temporal extension . Each participant’s MEG data were then coregistered with structural T1-weighted MRI data. Blink and cardiac artifacts were removed by signal-space projection .
The continuous time series was divided into 5000 ms epochs, with a 400 ms baseline directly preceding stimulus onset (−400 to 0 ms). Epochs were rejected using a fixed threshold method, supplemented with visual inspection. The distribution of amplitude and gradient was examined per participant across all trials, and the epochs with highest extrema amplitude and/or gradient values relative to that individual’s distribution were rejected with fixed amplitude and gradient thresholds. Cutoffs were determined individually so as to minimize bias due to head size, sensor proximity, and related variables that influence the signal amplitude, since magnetic field strength decreases exponentially with distance from the source. When appropriate, artifact-free trials in some participants were randomly omitted to balance the total number of trials across groups and/or conditions, which ensures results are not confounded by signal-to-noise ratio differences arising from imbalances in the number of trials. Only correctly answered trials were included in final analyses and participants with poor performance (i.e., near chance accuracy) were excluded from final analyses. The resulting artifact-free trials were transformed in the time-frequency domain using complex demodulation, and the resulting spectral power estimations per sensor were averaged over trials to generate time-frequency plots of mean spectral density. Sensor-level data were normalized by dividing the power of each time-frequency bin by the mean baseline power. Statistical analysis of the sensor-level spectrograms comparing the active window relative to the baseline period across the array of gradiometers during the task period were then computed to determine time-frequency windows to be examined in source-space.
Each data point in the spectrogram was initially evaluated using a mass univariate approach based on the general linear model (GLM), and then corrected for multiple comparisons in stage two. First, t-tests were conducted on each data point against the mean baseline value at that frequency and the output spectrograms of t-values were thresholded at p < 0.05 to define time-frequency bins containing potentially significant oscillatory deviations across all participants. In stage two, time-frequency bins that survived the threshold were clustered with other significant temporally and/or spectrally neighboring bins, and a cluster value was derived by summing the t-values of all data points in the cluster. Nonparametric permutation testing was then used to derive a distribution of cluster-values and the significance level of the observed clusters were tested directly using this distribution [51, 52]. From these analyses, time-frequency windows with significant oscillatory responses across all participants and conditions were then imaged at a 4.0 × 4.0 × 4.0 mm resolution using a linearly-constrained minimum variance vector beamformer [53, 54], as implemented in the Brain Electrical Source Analysis (BESA) software (v 7.1). This beamformer method uses spatial filters in the time-frequency domain to calculate source power for the entire brain volume. Following convention, we computed noise-normalized source power per voxel in each participant using active (i.e., task) and passive (i.e., baseline) periods of equal duration and bandwidth. As mentioned above, the passive period was defined as the 400 ms pre-stimulus baseline period (i.e., −400 to 0 ms). The resultant images are referred to as pseudo-t maps, with units that reflect noise-normalized power differences at each voxel. Each participant’s functional MEG images were transformed into standardized space using the transform that was previously applied to the structural images and spatially-resampled.
Oscillatory responses that were longer than the baseline in duration (>400 ms) were broken down into non-overlapping time windows of 400 ms or less and imaged separately using the same baseline. These images were then averaged across all relevant windows comprising the oscillatory response to derive a single map per oscillatory response per participant. These whole brain maps were then averaged across all participants to determine the anatomical basis of the oscillatory responses identified in the sensor-level analyses. Next, the whole brain maps of activity were subtracted (switch – no-switch) per individual to compute neural switch costs, or the increased amount of recruitment required to switch rule sets, per oscillatory response. T-tests for differences between patients with and without comorbid conditions were computed for each time-frequency window and all output statistical maps of the neural switch cost were thresholded at p < 0.01, using a spatial extent threshold (i.e., cluster restriction; k = 20 4 mm3 voxels) based on the theory of Gaussian random fields . Significant peaks from these t-tests were extracted for post-hoc analyses that focused on correlations with behavioral switch costs. A follow-up analysis of covariance (ANCOVA) was also computed to test the interaction between neural and behavioral switch costs by group, to determine whether the directionality of effects was driven by the presence or absence of comorbid conditions. Further, to examine the specific effects of glycemic control level on task switching in the brain, a whole brain correlational approach was used, with HbA1c as our measure of glycemic control.
ACR: Albumin Creatinine Ratio; ALT: Alanine aminotransferase; AST: Aspartate Transaminase; BESA: Brain Electrical Source Analysis; GABA: Gamma-Amino Butyric Acid; MEG: Magneto Encephalo Graphy; UNMC: University of Nebraska Medical Center.
C.M.E., A.T.D., C.V.D., and T.W.W. contributed to experimental design. C.M.E., G.H.L., A.T.D., and C.V.D., collected data and assisted in participant recruitment. C.M.E., C.V.D., and T.W.W. contributed to data analysis and interpretation. C.M.E. drafted the article. C.M.E., A.T.D., C.V.D., and T.W.W. contributed to article revisions. All authors take responsibility for the contents of this article.
We thank our participants for volunteering their time to the current study.
Conflicts of Interest
Authors C.M.E., G.H.L., A.T.D., and T.W.W. have no conflicts of interest to declare. C.V.D. is a consultant with NovoNordisk, Bayer, and AstraZeneca, and receives research funding from Theracos and Sanofi.
This work was supported by the National Institutes of Health [R01-MH103220 (TWW), R01-MH116782 (TWW), R01-MH118013 (TWW), P20-GM144641 (TWW), and an award from the National Science Foundation (#1539067; TWW)]. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
- 1. McCrimmon RJ, Ryan CM, Frier BM. Diabetes and cognitive dysfunction. Lancet. 2012; 379:2291–9. https://doi.org/10.1016/S0140-6736(12)60360-2 [PubMed]
- 2. Palta P, Schneider AL, Biessels GJ, Touradji P, Hill-Briggs F. Magnitude of cognitive dysfunction in adults with type 2 diabetes: a meta-analysis of six cognitive domains and the most frequently reported neuropsychological testswithin domains. J Int Neuropsychol Soc. 2014; 20:278–91. https://doi.org/10.1017/S1355617713001483 [PubMed]
- 3. Biessels GJ, Deary IJ, Ryan CM. Cognition and diabetes: a lifespan perspective. Lancet Neurol. 2008; 7:184–90. https://doi.org/10.1016/S1474-4422(08)70021-8 [PubMed]
- 4. Raffield LM, Cox AJ, Freedman BI, Hugenschmidt CE, Hsu FC, Wagner BC, Xu J, Maldjian JA, Bowden DW. Analysis of the relationships between type 2 diabetes status, glycemic control, and neuroimaging measures in the Diabetes Heart Study Mind. Acta Diabetol. 2016; 53:439–47. https://doi.org/10.1007/s00592-015-0815-z [PubMed]
- 5. Moran C, Beare R, Wang W, Callisaya M, Srikanth V, and Alzheimer's Disease Neuroimaging Initiative (ADNI). Type 2 diabetes mellitus, brain atrophy, and cognitive decline. Neurology. 2019; 92:e823–30. https://doi.org/10.1212/WNL.0000000000006955 [PubMed]
- 6. Moulton CD, Costafreda SG, Horton P, Ismail K, Fu CH. Meta-analyses of structural regional cerebral effects in type 1 and type 2 diabetes. Brain Imaging Behav. 2015; 9:651–62. https://doi.org/10.1007/s11682-014-9348-2 [PubMed]
- 7. van Bussel FCG, Backes WH, Hofman PAM, van Oostenbrugge RJ, van Boxtel MPJ, Verhey FRJ, Steinbusch HWM, Schram MT, Stehouwer CDA, Wildberger JE, Jansen JFA. Cerebral Pathology and Cognition in Diabetes: The Merits of Multiparametric Neuroimaging. Front Neurosci. 2017; 11:188. https://doi.org/10.3389/fnins.2017.00188 [PubMed]
Roberts RO, Knopman DS, Przybelski SA, Mielke MM, Kantarci K, Preboske GM, Senjem ML, Pankratz VS, Geda YE, Boeve BF, Ivnik RJ, Rocca WA, Petersen RC, Jack CR
Jr. Association of type 2 diabetes with brain atrophy and cognitive impairment. Neurology. 2014; 82:1132–41. https://doi.org/10.1212/WNL.0000000000000269 [PubMed]
- 9. Zhang Y, Lu S, Liu C, Zhang H, Zhou X, Ni C, Qin W, Zhang Q. Altered brain activation and functional connectivity in working memory related networks in patients with type 2 diabetes: An ICA-based analysis. Sci Rep. 2016; 6:23767. https://doi.org/10.1038/srep23767 [PubMed]
- 10. Biessels GJ, Reijmer YD. Brain changes underlying cognitive dysfunction in diabetes: what can we learn from MRI? Diabetes. 2014; 63:2244–52. https://doi.org/10.2337/db14-0348 [PubMed]
- 11. Chen Y, Zhou Z, Liang Y, Tan X, Li Y, Qin C, Feng Y, Ma X, Mo Z, Xia J, Zhang H, Qiu S, Shen D. Classification of type 2 diabetes mellitus with or without cognitive impairment from healthy controls using high-order functional connectivity. Hum Brain Mapp. 2021; 42:4671–84. https://doi.org/10.1002/hbm.25575 [PubMed]
- 12. Mansur RB, Cha DS, Woldeyohannes HO, Soczynska JK, Zugman A, Brietzke E, McIntyre RS. Diabetes mellitus and disturbances in brain connectivity: a bidirectional relationship? Neuromolecular Med. 2014; 16:658–68. https://doi.org/10.1007/s12017-014-8316-8 [PubMed]
- 13. Biessels GJ, Nobili F, Teunissen CE, Simó R, Scheltens P. Understanding multifactorial brain changes in type 2 diabetes: a biomarker perspective. Lancet Neurol. 2020; 19:699–710. https://doi.org/10.1016/S1474-4422(20)30139-3 [PubMed]
- 14. Monsell S. Task switching. Trends Cogn Sci. 2003; 7:134–40. https://doi.org/10.1016/s1364-6613(03)00028-7 [PubMed]
- 15. Rogers R, Monsell S. Costs of a predictible switch between simple cognitive tasks. J Exp Psychol Gen. 1995; 124:207–31. https://doi.org/10.1037/0096-3418.104.22.168
- 16. Dosenbach NU, Fair DA, Cohen AL, Schlaggar BL, Petersen SE. A dual-networks architecture of top-down control. Trends Cogn Sci. 2008; 12:99–105. https://doi.org/10.1016/j.tics.2008.01.001 [PubMed]
- 17. Capizzi M, Ambrosini E, Arbula S, Mazzonetto I, Vallesi A. Electrophysiological Evidence for Domain-General Processes in Task-Switching. Front Hum Neurosci. 2016; 10:124. https://doi.org/10.3389/fnhum.2016.00124 [PubMed]
- 18. Cooper PS, Wong ASW, McKewen M, Michie PT, Karayanidis F. Frontoparietal theta oscillations during proactive control are associated with goal-updating and reduced behavioral variability. Biol Psychol. 2017; 129:253–64. https://doi.org/10.1016/j.biopsycho.2017.09.008 [PubMed]
- 19. Proskovec AL, Wiesman AI, Wilson TW. The strength of alpha and gamma oscillations predicts behavioral switch costs. Neuroimage. 2019; 188:274–81. https://doi.org/10.1016/j.neuroimage.2018.12.016 [PubMed]
- 20. Moheet A, Mangia S, Seaquist ER. Impact of diabetes on cognitive function and brain structure. Ann N Y Acad Sci. 2015; 1353:60–71. https://doi.org/10.1111/nyas.12807 [PubMed]
- 21. Kullmann S, Heni M, Hallschmid M, Fritsche A, Preissl H, Häring HU. Brain Insulin Resistance at the Crossroads of Metabolic and Cognitive Disorders in Humans. Physiol Rev. 2016; 96:1169–209. https://doi.org/10.1152/physrev.00032.2015 [PubMed]
- 22. Kullmann S, Heni M, Fritsche A, Preissl H. Insulin action in the human brain: evidence from neuroimaging studies. J Neuroendocrinol. 2015; 27:419–23. https://doi.org/10.1111/jne.12254 [PubMed]
- 23. Backeström A, Papadopoulos K, Eriksson S, Olsson T, Andersson M, Blennow K, Zetterberg H, Nyberg L, Rolandsson O. Acute hyperglycaemia leads to altered frontal lobe brain activity and reduced working memory in type 2 diabetes. PLoS One. 2021; 16:e0247753. https://doi.org/10.1371/journal.pone.0247753 [PubMed]
- 24. Ott V, Benedict C, Schultes B, Born J, Hallschmid M. Intranasal administration of insulin to the brain impacts cognitive function and peripheral metabolism. Diabetes Obes Metab. 2012; 14:214–21. https://doi.org/10.1111/j.1463-1326.2011.01490.x [PubMed]
- 25. Alagiakrishnan K, Sankaralingam S, Ghosh M, Mereu L, Senior P. Antidiabetic drugs and their potential role in treating mild cognitive impairment and Alzheimer's disease. Discov Med. 2013; 16:277–86. [PubMed]
- 26. Dosenbach NU, Fair DA, Miezin FM, Cohen AL, Wenger KK, Dosenbach RA, Fox MD, Snyder AZ, Vincent JL, Raichle ME, Schlaggar BL, Petersen SE. Distinct brain networks for adaptive and stable task control in humans. Proc Natl Acad Sci U S A. 2007; 104:11073–8. https://doi.org/10.1073/pnas.0704320104 [PubMed]
- 27. Schweizer TA, Oriet C, Meiran N, Alexander MP, Cusimano M, Stuss DT. The cerebellum mediates conflict resolution. J Cogn Neurosci. 2007; 19:1974–82. https://doi.org/10.1162/jocn.2007.19.12.1974 [PubMed]
- 28. Eckstrand KL, Mummareddy N, Kang H, Cowan R, Zhou M, Zald D, Silver HJ, Niswender KD, Avison MJ. An insulin resistance associated neural correlate of impulsivity in type 2 diabetes mellitus. PLoS One. 2017; 12:e0189113. https://doi.org/10.1371/journal.pone.0189113 [PubMed]
- 29. Ishizawa KT, Kumano H, Sato A, Sakura H, Iwamoto Y. Decreased response inhibition in middle-aged male patients with type 2 diabetes. Biopsychosoc Med. 2010; 4:1. https://doi.org/10.1186/1751-0759-4-1 [PubMed]
- 30. Keizer AW, Verment RS, Hommel B. Enhancing cognitive control through neurofeedback: a role of gamma-band activity in managing episodic retrieval. Neuroimage. 2010; 49:3404–13. https://doi.org/10.1016/j.neuroimage.2009.11.023 [PubMed]
- 31. Cho RY, Konecky RO, Carter CS. Impairments in frontal cortical gamma synchrony and cognitive control in schizophrenia. Proc Natl Acad Sci U S A. 2006; 103:19878–83. https://doi.org/10.1073/pnas.0609440103 [PubMed]
- 32. Minzenberg MJ, Firl AJ, Yoon JH, Gomes GC, Reinking C, Carter CS. Gamma oscillatory power is impaired during cognitive control independent of medication status in first-episode schizophrenia. Neuropsychopharmacology. 2010; 35:2590–9. https://doi.org/10.1038/npp.2010.150 [PubMed]
- 33. Jensen O, Kaiser J, Lachaux JP. Human gamma-frequency oscillations associated with attention and memory. Trends Neurosci. 2007; 30:317–24. https://doi.org/10.1016/j.tins.2007.05.001 [PubMed]
- 34. Wilson TW, Heinrichs-Graham E, Proskovec AL, McDermott TJ. Neuroimaging with magnetoencephalography: A dynamic view of brain pathophysiology. Transl Res. 2016; 175:17–36. https://doi.org/10.1016/j.trsl.2016.01.007 [PubMed]
- 35. Corbetta M, Shulman GL. Control of goal-directed and stimulus-driven attention in the brain. Nat Rev Neurosci. 2002; 3:201–15. https://doi.org/10.1038/nrn755 [PubMed]
- 36. Arnold SE, Arvanitakis Z, Macauley-Rambach SL, Koenig AM, Wang HY, Ahima RS, Craft S, Gandy S, Buettner C, Stoeckel LE, Holtzman DM, Nathan DM. Brain insulin resistance in type 2 diabetes and Alzheimer disease: concepts and conundrums. Nat Rev Neurol. 2018; 14:168–81. https://doi.org/10.1038/nrneurol.2017.185 [PubMed]
- 37. Kujala J, Jung J, Bouvard S, Lecaignard F, Lothe A, Bouet R, Ciumas C, Ryvlin P, Jerbi K. Gamma oscillations in V1 are correlated with GABA(A) receptor density: A multi-modal MEG and Flumazenil-PET study. Sci Rep. 2015; 5:16347. https://doi.org/10.1038/srep16347 [PubMed]
- 38. Uhlhaas PJ, Pipa G, Neuenschwander S, Wibral M, Singer W. A new look at gamma? High- (>60 Hz) γ-band activity in cortical networks: function, mechanisms and impairment. Prog Biophys Mol Biol. 2011; 105:14–28. https://doi.org/10.1016/j.pbiomolbio.2010.10.004 [PubMed]
- 39. Sanches M, Abuhaiba SI, d'Almeida OC, Quendera B, Gomes L, Moreno C, Guelho D, Castelo-Branco M. Diabetic brain or retina? Visual psychophysical performance in diabetic patients in relation to GABA levels in occipital cortex. Metab Brain Dis. 2017; 32:913–21. https://doi.org/10.1007/s11011-017-9986-3 [PubMed]
- 40. Thielen JW, Gancheva S, Hong D, Rohani Rankouhi S, Chen B, Apostolopoulou M, Anadol-Schmitz E, Roden M, Norris DG, Tendolkar I. Higher GABA concentration in the medial prefrontal cortex of Type 2 diabetes patients is associated with episodic memory dysfunction. Hum Brain Mapp. 2019; 40:4287–95. https://doi.org/10.1002/hbm.24702 [PubMed]
- 41. d'Almeida OC, Violante IR, Quendera B, Moreno C, Gomes L, Castelo-Branco M. The neurometabolic profiles of GABA and Glutamate as revealed by proton magnetic resonance spectroscopy in type 1 and type 2 diabetes. PLoS One. 2020; 15:e0240907. https://doi.org/10.1371/journal.pone.0240907 [PubMed]
- 42. van Bussel FCG, Backes WH, Hofman PAM, Puts NAJ, Edden RAE, van Boxtel MPJ, Schram MT, Stehouwer CDA, Wildberger JE, Jansen JFA. Increased GABA concentrations in type 2 diabetes mellitus are related to lower cognitive functioning. Medicine (Baltimore). 2016; 95:e4803. https://doi.org/10.1097/MD.0000000000004803 [PubMed]
- 43. Moran C, Beare R, Phan T, Starkstein S, Bruce D, Romina M, Srikanth V. Neuroimaging and its Relevance to Understanding Pathways Linking Diabetes and Cognitive Dysfunction. J Alzheimers Dis. 2017; 59:405–19. https://doi.org/10.3233/JAD-161166 [PubMed]
- 44. Kodl CT, Seaquist ER. Cognitive dysfunction and diabetes mellitus. Endocr Rev. 2008; 29:494–511. https://doi.org/10.1210/er.2007-0034 [PubMed]
- 45. West RK, Ravona-Springer R, Schmeidler J, Leroith D, Koifman K, Guerrero-Berroa E, Preiss R, Hoffman H, Silverman JM, Heymann A, Schnaider-Beeri M. The association of duration of type 2 diabetes with cognitive performance is modulated by long-term glycemic control. Am J Geriatr Psychiatry. 2014; 22:1055–9. https://doi.org/10.1016/j.jagp.2014.01.010 [PubMed]
- 46. Cukierman-Yaffe T, Gerstein HC, Williamson JD, Lazar RM, Lovato L, Miller ME, Coker LH, Murray A, Sullivan MD, Marcovina SM, Launer LJ, and Action to Control Cardiovascular Risk in Diabetes-Memory in Diabetes (ACCORD-MIND) Investigators. Relationship between baseline glycemic control and cognitive function in individuals with type 2 diabetes and other cardiovascular risk factors: the action to control cardiovascular risk in diabetes-memory in diabetes (ACCORD-MIND) trial. Diabetes Care. 2009; 32:221–6. https://doi.org/10.2337/dc08-1153 [PubMed]
- 47. Geijselaers SLC, Sep SJS, Stehouwer CDA, Biessels GJ. Glucose regulation, cognition, and brain MRI in type 2 diabetes: a systematic review. Lancet Diabetes Endocrinol. 2015; 3:75–89. https://doi.org/10.1016/S2213-8587(14)70148-2 [PubMed]
- 48. American Diabetes Association. 4. Comprehensive Medical Evaluation and Assessment of Comorbidities: Standards of Medical Care in Diabetes-2021. Diabetes Care. 2021 (Suppl 1); 44:S40–52. https://doi.org/10.2337/dc21-S004 [PubMed]
- 49. Taulu S, Simola J. Spatiotemporal signal space separation method for rejecting nearby interference in MEG measurements. Phys Med Biol. 2006; 51:1759–68. https://doi.org/10.1088/0031-9155/51/7/008 [PubMed]
- 50. Uusitalo MA, Ilmoniemi RJ. Signal-space projection method for separating MEG or EEG into components. Med Biol Eng Comput. 1997; 35:135–40. https://doi.org/10.1007/BF02534144 [PubMed]
- 51. Ernst MD. Permutation Methods: A Basis for Exact Inference. Statist Sci. 2004; 19:676–85. https://doi.org/10.1214/088342304000000396
- 52. Maris E, Oostenveld R. Nonparametric statistical testing of EEG- and MEG-data. J Neurosci Methods. 2007; 164:177–90. https://doi.org/10.1016/j.jneumeth.2007.03.024 [PubMed]
- 53. Gross J, Kujala J, Hamalainen M, Timmermann L, Schnitzler A, Salmelin R. Dynamic imaging of coherent sources: Studying neural interactions in the human brain. Proc Natl Acad Sci U S A. 2001; 98:694–9. https://doi.org/10.1073/pnas.98.2.694 [PubMed]
- 54. Hillebrand A, Singh KD, Holliday IE, Furlong PL, Barnes GR. A new approach to neuroimaging with magnetoencephalography. Hum Brain Mapp. 2005; 25:199–211. https://doi.org/10.1002/hbm.20102 [PubMed]
- 55. Worsley KJ, Marrett S, Neelin P, Vandal AC, Friston KJ, Evans AC. A unified statistical approach for determining significant signals in images of cerebral activation. Hum Brain Mapp. 1996; 4:58–73. https://doi.org/10.1002/(SICI)1097-0193(1996)4:1<58::AID-HBM4>3.0.CO;2-O [PubMed]