Abstract

Slower psychomotor speed is very common in patients with type 1 diabetes mellitus (T1D), but the underlying mechanisms are not clear. We propose that hyperglycemia is associated with slower psychomotor speed via disruption of brain activation. Eighty-five adults (48% women, mean age: 49.0 years, mean duration: 40.8) with childhood onset T1D were recruited for this cross-sectional study. Median response time in seconds (longer = worse performance) and brain activation were measured while performing a psychomotor speed task. Exposure to hyperglycemia, measured as glycosylated hemoglobin A1c, was associated with longer response time and with higher activation in the inferior frontal gyrus and primary sensorimotor and dorsal cingulate cortex. Higher activation in inferior frontal gyrus, primary sensorimotor cortex, thalamus, and cuneus was related to longer response times; in contrast, higher activation in the superior parietal lobe was associated with shorter response times. Associations were independent of small vessel disease in the brain or other organs. In this group of middle-aged adults with T1D, the pathway linking chronic hyperglycemia with slower processing speed appears to include increased brain activation, but not small vessel disease. Activation in the superior parietal lobe may compensate for dysregulation in brain activation in the presence of hyperglycemia.

1. Introduction

Compared to similarly aged nondiabetic adults, persons with type 1 diabetes (T1D) are more likely to exhibit cognitive deficits, especially localized within the psychomotor and information processing domains [15]. Slower psychomotor speed is associated with several negative health outcomes, in persons both with [68] and without diabetes [9, 10]. Deficits in psychomotor speed appear to increase as persons with T1D age and/or live longer with diabetes [1113]. With the rising number of persons with T1D now surviving into older age [14], it is increasingly important to characterize the predictors of slower psychomotor speed to design interventions that reduce its negative impact on health outcomes.

The mechanisms underlying slower psychomotor speed in T1D are not entirely clear. Chronic hyperglycemia has been implicated in cognitive complications in several studies of T1D, including ours [1518]. It has been suggested that chronic hyperglycemia may compromise cognitive function because of its detrimental effects on the microvasculature [1820]. Indeed, small vessel diseases, including brain white matter hyperintensities and retinopathy, have been associated with slower psychomotor speed in persons with T1D [1, 2124]. We have recently shown that white matter hyperintensities partially explain the slower psychomotor speed observed in middle-aged adults with T1D as compared to similarly aged non-T1D adults [25].

One hypothesized pathway linking chronic hyperglycemia, brain small vessel disease, and slower psychomotor speed could be via disrupted neuronal activation. Thickened cerebral cortical capillary basement membranes in diabetics and decreased capillary density are two hallmark signs of small vessel disease [26] which can potentially reduce oxygen delivery and vasoreactivity [27] and therefore alter neuronal activation. Moreover, positive relationships between imaging signs of small vessel disease and higher brain activation have been recently reported in neuroimaging studies of older adults [2830] and in one small study of selected patients with T1D, using retinal imaging [20].

We propose that hyperglycemia affects psychomotor speed via disruption of brain activation and that brain small vessel disease mediates these associations. To test this hypothesis, we measured patterns of brain activation while performing a task of psychomotor speed, concurrently with historic measures of exposure to hyperglycemia and small vessel disease in the brain and other target organs.

2. Methods

2.1. Participants

In 2010–2013, individuals from the Pittsburgh Epidemiology of Diabetes Complications Study (EDC) were invited to a neurovascular ancillary imaging study. All participants had childhood-onset (<17 years of age) T1D and were diagnosed or seen within one year of diagnosis (1950–1980) at Children’s Hospital of Pittsburgh. Following the first clinical assessment (1986–1988), when average participant age and diabetes duration were 28 and 19 years, respectively, biennial examinations were conducted for 10 years, with further examinations at 18 and 25 years [31]. Among the 263 locally resident individuals who participated in the 24-year EDC followup and were invited for this neurovascular imaging study, 149 accepted the invitation, 112 of whom did not have contraindication for brain MRI. Brain MRIs were obtained on participants, and 89 of these had both functional brain MRI and performance measures while performing the task in the scanner. Of these, 85 also had complete cognitive data obtained outside the scanner and are included in this analysis.

All study procedures were approved by the Institutional Review Board at University of Pittsburgh Medical Center.

2.2. Measures of Interest

Prevalence of diabetes-related complications and other health and lifestyle factors were assessed periodically from 1986–1988 to 2004–2006 via questionnaire, physical exam, and medical records [31].

2.2.1. Hyperglycemia

Glycosylated hemoglobin A1c (HbA1c) was measured at time of brain MRI and also at repeated intervals since study entry, using automated high-performance liquid chromatography as previously described [32]. A measure assessing degree and duration of hyperglycemia is computed as “HbA1c months” (for details, see [32]). HbA1c is a well-established marker of mean glucose concentrations during the prior 2-3 months [33].

2.2.2. Brain MRI Protocol

Brain MRI scanning was conducted at the Magnetic Resonance Research Center at the University of Pittsburgh using a 3-Tesla Siemens Trio TIM scanner, with a 12-channel head coil. Structural scans were acquired prior to functional scans. Details on the protocol of acquisition are described in greater detail elsewhere [34]. A fully automated method for quantifying and localizing white matter hyperintensities on MR images was used. The algorithm uses a fully automated fuzzy connectedness algorithm. Seeds are generated from the histogram of the FLAIR image and thus user input is not needed compared to other algorithms. Each seed has its own parameters (providing more specific segmentation) which are calculated automatically (fuzzy adjacency and affinity) [35]. Magnetization-prepared rapid gradient echo (MPRAGE) T1-weighted images were acquired in the axial plane: TR = 2300 ms; TE = 3.43 ms; TI = 900 ms: flip angle = 9°; slice thickness = 1 mm; FOV = 256 mm × 224 mm; voxel size = 1 mm × 1 mm; matrix size 256 × 240; number of slices = 176. T2-weighted and FLAIR images were acquired as previously described [25].

2.2.3. Functional MRI Paradigm

In this study, the computerized version of the Digit Symbol Substitution Tests (DSST) [36, 37] was adapted for the scanning procedure to test psychomotor speed, following a protocol described previously [38, 39]. The paradigm tests motor speed, response selection [40], shifting of attention [41], and perceptual speed [42, 43]. Psychomotor speed is used to actively maintain and manipulate information over a brief period of time and to allocate attentional resources among competing subtasks [44, 45]. Frontal and parietal lobes, thalamus, precuneus, and cingulate cortex are expected to be involved in performing these tasks [38, 39, 46, 47]. For this protocol, the subject sees one number-symbol matching pair (cue) on a computer screen. After the cue disappears, an answer key (probe) appears containing a grid of four number-symbol matching pairs. The subject is instructed to push the right index finger button if the probe contains a number-symbol that matches the cue or to push the left index finger button if the probe does not contain a number-symbol that matches the cue. Instructions are to respond “as fast as you possibly can.” Participants were instructed on the task outside the magnet for as long as needed to familiarize them with the task (usually 5–10 min). The DSST was presented using a block design with 8 trials/block, in which a block of the experimental DSST condition is alternated with control condition for a total of 10 blocks: 5 DSST blocks alternating with 5 control blocks for a total of 9 min and 20 s. Each block lasted for 56 s including 8 s of instructions at the start of each block to remind the subjects of the task. The 8 s of instruction and the first trial of each block were excluded from the analyses. The different task conditions (matching versus nonmatching) were randomized 1 : 1 across the block. The nonmatching condition was designed as a control condition to account for the non-task-specific brain activation, frontal eye field and visual cortex for eye movements and motor cortex for the index fingers movement, which are likely to be elicited during DSST performance in addition to the task-specific activation of the executive control function network. Responses occurring after the probe disappeared (after 3.6 s) were termed as “no response” and were omitted from the calculation of accuracy and response time. Median response time during the task was computed as the interval of time between appearance of the probe and the subject’s response. Accuracy was recorded as number of correct responses divided by total number of responses [38, 39].

2.2.4. Small Vessel Disease in the Brain

Total white matter hyperintensities (WMHs) volume was calculated using a previously published automated method [35, 48] and normalized for brain volume. Burden of WMH was also measured using Fazekas ratings as previously described [25].

2.2.5. Small Vessel Disease in Other Target Organs

Eye Disease. Proliferative retinopathy was assessed by three-standard field stereo fundus photography at each EDC examination by the Madison Reading Center in Wisconsin and graded according to the ETDRS Classifications. A history of laser therapy for proliferative retinopathy is used as a proxy in the absence of photos.

Neuropathy. Distal symmetric polyneuropathy is determined by medical history and clinical examination using the DCCT protocol (i.e., the presence of two or more of the following: symptoms (consistent with DSP), reduced deep tendon reflexes, and signs of sensory loss) and it is further “confirmed” by determining age specific abnormal vibratory thresholds (Vibratron II) of the great toe. Cardiac autonomic neuropathy is assessed by heart variation during deep breathing and heart rate and BP response to standing.

Renal Disease. Overt nephropathy (ON), defined as the presence of renal failure (serum creatinine >5 mg/dL and/or ESRD, i.e., dialysis and/or renal transplant) or an albumin excretion rate >200 μg/min in at least 2 of 3 timed urine samples or (in the absence of any urine collections) serum creatinine >2 mg/dL and microalbuminuria, defined as an albumin excretion rate between 20 and 200 μg/min in at least 2 or 3 urine samples. If the adequacy of timed urine samples is questionable in terms of creatinine excretion, a previously validated urinary albumin to creatinine ratio (mg/gm) was used (i.e., >0.30 ON; 0.03–0.29 microalbuminuria).

2.2.6. Cognitive Function

Participants were evaluated using a comprehensive battery of tests [17]: estimated verbal intelligence (North American Adult Reading Test [NAART]); psychomotor efficiency (Digit Symbol Substitution Test [DSST] and Grooved Pegboard); executive function (Trail Making Test, Part B [TMTB], Verbal Fluency [VF: animal naming and F-A-S], Stroop Color-Word, and Letter/Number Sequencing [LN Sequence]); learning and working memory (Rey Auditory Verbal Learning Tests [RVLT], Four Word Short-Term Memory [4WSTM] and Rey Osterrieth Complex Figure Delayed task [ROCF-delay]); and visuoconstruction skills (Rey Osterrieth Complex Figure Copy task [ROCF-copy]). Details on each of these tests can be found elsewhere (see Strauss et al. [49]). Individual raw test scores that were ≥1.5 SD worse than demographically appropriate published norms [4951] were classified as “impaired” and impairment indicator variables were created for each test [17].

2.3. Statistical Analysis

Associations between sample characteristics and median response time while performing the task in the scanner were tested using Spearman or Pearson correlation coefficients depending on whether the distribution of the variables of interest was normal. The association of HbA1c months with median reaction time was tested in multivariable linear regression models after adjustment for markers of small vessel disease (brain WMHs, neuropathy, nephropathy, and retinopathy, and cardiac autonomic neuropathy). Since the markers of small vessel disease are highly correlated with each other, they entered the model separately. Small vessel disease in the brain was considered a potential mediator of the association between HbA1c and functional activation. Small vessel complications outside the brain can represent peripheral biomarkers of underlying brain damage; and were therefore considered confounders of the association between hyperglycemia and functional activation. The change in the regression coefficient of HbA1c months predicting median reaction time after adjusting for the above variables was measured; a change of >10% was considered as an indication that small vessel complications substantially attenuated the main association [52]. Models were further adjusted for the sample characteristics associated with median reaction time in univariate models.

Functional imaging data were analyzed with Statistical Parametric Mapping 5 (SPM5) using general linear models [53]. For each subject, all Epipolar-Plane Image (EPI) volumes were realigned to the first volume of the time series. A standard SPM high pass filter (1/128 hz) was applied. Smoothing was done at FWHM = 8 mm; physiological noise correction was performed. Each condition was modeled using a delayed boxcar function convolved with the SPM hemodynamic response function. Movement parameters derived from realignment were added as covariates of no interest to correct for confounding effects induced by head movement. Contrasts of interest (DSST > control and control > DSST) were first estimated for each subject individually averaging activation across runs. The resulting statistical maps with the contrast of interest for each subject were our outcome maps for the group analysis. More specifically, the single-subject statistical maps with the contrast of interest were difference of beta coefficients from the regression of task minus control. Three group analyses were conducted: first, the main effect of task across all participants was tested in relation to WMH, HbA1c, and response time; this was done by taking a one-sample -test across all subjects on the beta maps task-control which generates a -statistic map. The relationship between HbA1c and WMH and activation (during condition) was also investigated by testing the null hypothesis that the slope WMH = 0 for our first model and slope median response time = 0 for our second model using a multiple linear regression model and testing the coefficient of interest using a -test. The statistical maps representing the associations between fMRI activation and WMH used a threshold (, uncorrected). To correct for multiple comparisons, permutation methods for peak-cluster level error correction were applied (AlphaSim). Specifically, we compute (via simulations) the probability of a random field of noise producing a cluster of a certain size given some threshold (). We found that a minimum cluster size of 130 voxels was needed to attain significance of , corrected.

3. Results

The characteristics of the participants included in this analysis (Table 1) are similar to those previously reported for this cohort [17, 25, 5456].

3.1. Associations between Sample Characteristics and Response Time While Performing the DSST in the Scanner

Length of response time while performing the DSST in the scanner (e.g., longer time = worse performance) was directly associated with age, T1D duration, HbA1c months, burden of WMH, and higher prevalence of small vessel disease and inversely related to accuracy (Table 1). Of note, associations with HbA1c or glucose levels on the day of MRI were not significant (Table 1). Length of response time while performing the DSST in the scanner was also inversely associated with scores on cognitive tests performed outside the scanner, including the pencil and paper DSST and the composite score from the psychomotor speed domain (Table 1). Longer response time was also associated with higher probability of having clinically relevant cognitive impairment (Table 1). These associations remained significant after adjustment for age or duration (not shown).

3.2. Associations between Hyperglycemia and Response Time While Performing the DSST in the Scanner

In multivariable regression models (Table 2), higher HbA1c months remained associated with longer response time independent of adjustment for age, WMH, or other markers of small vessel diseases. Of note, the size of the regression coefficient was attenuated by >10% after adjustment for proliferative retinopathy while it remained similar after adjustment for other markers of small vessel disease (Table 2).

3.3. Main Effect

Brain activation while performing the DSST in the scanner was significantly greater than brain activation while performing the control condition in regions previously shown to be related to psychomotor speed in general and to this task in particular (Figure 1 and Table 3): dorsolateral prefrontal cortex which extended rostrally to include the inferior frontal gyrus, caudally to include pre- and postcentral gyri and the parietal cortex, and subcortically to include thalamus and dorsal cingulate gyrus. Analyses to test the opposite direction of association (main effect of control condition minus DSST) displayed a network of areas that extended from cortical temporal areas toward the posterior cingulate gyrus and medial frontal gyrus medially.

3.4. Associations between Hyperglycemia, Brain Activation While Performing the DSST in the Scanner, and Response Time

There was a significant and direct association between HbA1c months and brain activation while performing the task (Figure 2), which was localized in some of the regions related to the main effect: inferior frontal, precentral, and dorsal cingulate gyri in the right hemisphere (Table 4). Analyses were similar after adjustment for WMH and complications (not shown).

There was a significant and direct association between brain activation and length of response time while performing the DSST (Figure 3), localized in the inferior frontal gyrus, precentral gyrus in the right hemisphere, the precentral gyrus in the left hemisphere, and the thalamus, postcentral gyrus, and cuneus bilaterally (Table 5). Additionally, there was a significant and inverse association between length of response time and brain activation while performing the DSST (Figure 3), localized in both the medial (precuneus) and lateral areas of the superior parietal lobule (Table 5).

3.5. Contribution of Brain Small Vessel Disease

Higher WMH burden was not associated with brain activation while performing the task in the scanner and the interaction of WMH burden and response time on brain activation while performing the task was also not significant (not shown).

4. Discussion

The results of this functional neuroimaging study suggest an explanatory pathway of slower psychomotor speed in middle-aged persons with childhood-onset T1D. We found that chronic hyperglycemia was related to higher activation in selected brain regions that were engaged in performing the DSST and to worse task performance. Moreover, we found a direct association between activation in some of these brain regions and performance. Contrary to our expectations, brain small vessel disease did not contribute to these relationships.

Understanding the mechanisms underlying slowing of psychomotor speed in persons with T1D is urgently needed. Psychomotor slowing reflects deficits in attention, planning, and problem-solving skills. It is a potent predictor of poorer diabetes management in adults with type 2 diabetes [68] and is associated with increased risk and fear of falls and greater psychological distress [6]. Psychomotor slowing is also a well-known risk factor for developing disability, dementia, and mortality in older persons without T1D [9, 10]; similar trends may be observed in future years for middle-aged persons with T1D, with rising personal, societal, and public health costs [57]. Both the increasing incidence rate [58] and longer life expectancy of persons with T1D [14] contribute to an increasing number of individuals with T1D who are aging. Therefore, addressing psychomotor slowing in these patients may help reduce the risk of developing these negative health outcomes.

The neurodegenerative effects of exposure to chronic hyperglycemia have been consistently shown. Chronic hyperglycemia appears to have a stronger impact on poor cognitive function for adults with T1D [17, 24] whereas hypoglycemia appears to play a more important role in the relationships with cognitive deficits in younger persons with T1D [59]. Consistent with the literature in adults with T1D, we found that chronic hyperglycemia was associated with poorer performance of a task of psychomotor speed while in the scanner, and this association was robust to adjustment for other covariates.

There appeared to be spatial colocalization of the regions that were more active in relationship to chronic hyperglycemia with the regions that were more active in the presence of slower performance, specifically the right inferior frontal and precentral gyri. One interpretation for these findings is that these regions are vulnerable to hyperglycemia and respond to this condition through hyperactivation, which in turn would negatively interfere with performance. Other studies have found an association between hyperglycemia and higher activation in these areas. A prior fMRI study [60] of visual spatial working memory showed that young persons with T1D (, mean age: 20 years) had higher brain activation compared to similarly aged persons without T1D (, age = 21 years) and this pattern of hyperactivation was localized in the right inferior frontal gyrus, cerebellum, basal ganglia, and substantia nigra. Another small study of newly diagnosed T2D [61] (mean age = 41, ) showed hyperactivation while performing a working memory task in a network of prefrontoparietal regions that included the inferior frontal gyrus, as compared to controls (mean age: 40 years, ). While the overlap of these studies’ results with ours (i.e., hyperactivation in the inferior prefrontal gyrus) is of interest, direct comparisons cannot be done because these prior studies did not directly examine within-group relationships between chronic hyperglycemia, microvascular complications, brain activation, and performance.

Why would chronic hyperglycemia cause higher brain activation? Activation as measured in the scanner is the result of several interrelated physiological processes related to exposure to a mental challenge. Blood flow, neuronal activation, and ability of the neurons to extract oxygen all contribute to the intensity of the activation signal. A mechanism implicating small vessel disease has been suggested to explain the effect of chronic hyperglycemia on lower cognitive function and higher brain activation. Chronic hyperglycemia could disrupt the integrity of the blood brain barrier through several pathways [19], including basal membrane thickening of the small vessels, which would be manifested radiologically as WMH and could compromise efficient oxygen delivery and lead to changes in brain activation while performing a task in the scanner. Although we found a direct relationship between WMH and poorer performance, WMH did not appear to mediate the association between hyperglycemia and performance in the scanner nor did we find an association between WMH and brain activation. A positive relationship between WMH lesion load and higher brain activation in the frontal lobe has been found in some studies of older adults [28, 29]. Unfortunately, evidence in favor of (or against) these associations has been indirect for persons with T1D. We are aware of one prior functional neuroimaging study of 24 adults (mean age: 40 years) with T1D [20] showing that those with small vessel disease, measured as presence of proliferative retinopathy, had higher task-related activation in the orbital frontal gyrus. It is possible that, in our cohort, WMH is not severe enough to interfere with the task we used in the scanner and/or that this task is too easy. Indeed, most participants performed the task with high accuracy. Another possibility is that the spatial distribution of WMH is not affecting the areas critical for performing this task. It would be informative to examine the patterns of activation in relationship with WMH while performing other tasks of psychomotor speed that are titrated for difficulty, as well as to examine the spatial distribution of WMH.

One other potential pathway linking hyperglycemia to higher activation could be via greater production of the excitatory neurotransmitter glutamate. It is known that hyperglycemia is related to higher levels of glutamate and a relationship between hyperglycemia, higher glutamate, and higher brain activation in the frontal lobe has been shown in patients with T1D [62]. Chronic exposure to glutamate is known to cause neurotoxicity; hence, higher activation related to higher glutamate could also explain the relationship between higher activation and lower performance.

The association between hyperglycemia and higher activation could also be due to the direct toxic effects of higher levels of glucose at time of testing. Prior studies examining the effect of transient, acute changes in glucose levels on brain activation, both within and outside physiological ranges, have mostly focused on hypoglycemia and have yielded conflicting results. Some studies have found hypoglycemia is related to lower brain activation [63, 64] while others have found that acute, transient hypoglycemia is related to higher activation [65, 66]. We did not find an association between glucose or HbA1c levels on day of testing with brain activation or performance; thus, it is more likely that the effects observed in our study would be due to longer-term exposure to hyperglycemia.

It remains unclear whether brain activation and performance are monotonically related [67, 68]. Higher activation has been found in the presence of higher performance, possibly to compensate for other deficits, and also in relationship with lower performance in persons with Alzheimer’s disease, as a sign of dedifferentiation. In general, higher activation in the presence of lower performance has been postulated to represent signs of dysregulation of neuronal activation. In our study, we found a dual pattern of activation in relationship to performance, consisting of inverse associations in cortico-subcortical networks that included those regions related to chronic hyperglycemia (e.g., inferior frontal gyrus, precentral gyrus), as well as direct associations with activation in the left superior parietal lobule. The direct association between performance and superior parietal activation is consistent with what we have previously found in a group of older adults with a range of chronic health conditions [38, 39]. In this prior study of older adults without T1D, we also found that activation in the superior parietal cortex was higher for those adults who maintained higher performance in the presence of higher WMH burden. Thus, higher activation in the superior parietal lobe may be compensating for abnormalities occurring in other tissues (e.g., presence of WMH in the study of older adults) and/or for the dysregulation of activation in regions vulnerable to hyperglycemia (e.g., hyperactivation of inferior frontal gyrus and precentral gyrus in the presence of hyperglycemia in this study). If confirmed in other studies, it would be very informative to examine the factors most related to integrity of the superior parietal lobe and whether improvement of activation in the superior parietal lobe would also improve psychomotor speed.

Strengths of our study include the large sample size and the extensive characterization of cognitive function using pencil and paper tasks, as well as information on complications obtained over prior years. The performance measure that we chose in this study, length of response time in the scanner, appeared to effectively capture the severity of impairment in psychomotor speed and the presence of cognitive impairment as shown by the strong correlations with the score on the pencil and paper tests. Another strength of this study is that accuracy and reaction time of the participants were measured while recording their actual brain signal, thus allowing correlation of performance with brain activation. Results of this study should be considered with caution, especially with regard to generalization of results. It is likely that this group of patients, who have lived most of their lives with T1D into their fifth and sixth decade of life, may harbor unique characteristics compared to younger persons with T1D. Moreover, we did not test all the potential pathways linking hyperglycemia with brain activation. A limitation inherent in all task-related fMRI studies is that results are task specific. It is possible that a different task of psychomotor speed would reveal a different spatial distribution of hyperactivation.

5. Conclusions

Additional studies are needed to assess whether patterns of brain activation would be helpful markers of cognitive complications in persons with T1D. It would also be important to know whether better glycemic control improves cognitive function via changes in the patterns of brain activation. This would not only have tremendous implications to further our understanding of the mechanisms underlying psychomotor speed but can also help understand whether the brains of persons with T1D retain plasticity and reserve.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to thank these individuals with type 1 diabetes for their willingness to participate in this study. This study was funded by the National Institutes of Health (NIH)/National Institute of Diabetes and Digestive and Kidney Disease (NIDDK) Grants RO1 DK089028 and R01 DK034818-25.