Abstract
Aging is accompanied by profound changes in the brain's dopamine system that affect cognitive function. Evidence of powerful individual differences in cognitive aging has sharpened focus on identifying biological factors underlying relative preservation versus vulnerability to decline. Dopamine represents a key target in these efforts. Alterations of dopamine receptors and dopamine synthesis are seen in aging, with receptors generally showing reduction and synthesis demonstrating increases. Using the PET tracer 6-[18F]fluoro-l-m-tyrosine, we found strong support for upregulated striatal dopamine synthesis capacity in healthy older adult humans free of amyloid pathology, relative to young people. We next used fMRI to define the functional impact of elevated synthesis capacity on cognitive flexibility, a core component of executive function. We found clear evidence in young adults that low levels of synthesis capacity were suboptimal, associated with diminished cognitive flexibility and altered frontoparietal activation relative to young adults with highest synthesis values. Critically, these relationships between dopamine, performance, and activation were transformed in older adults with higher synthesis capacity. Variability in synthesis capacity was related to intrinsic frontoparietal functional connectivity across groups, suggesting that striatal dopamine synthesis influences the tuning of networks underlying cognitive flexibility. Together, these findings define striatal dopamine's association with cognitive flexibility and its neural underpinnings in young adults, and reveal the alteration in dopamine-related neural processes in aging.
SIGNIFICANCE STATEMENT Few studies have combined measurement of brain dopamine with examination of the neural basis of cognition in youth and aging to delineate the underlying mechanisms of these associations. Combining in vivo PET imaging of dopamine synthesis capacity, fMRI, and a sensitive measure of cognitive flexibility, we reveal three core findings. First, we find evidence supporting older adults' capacity to upregulate dopamine synthesis. Second, we define relationships between dopamine, cognition, and frontoparietal activity in young adults indicating high levels of synthesis capacity are optimal. Third, we demonstrate alteration of these relationships in older adults, suggesting neurochemical modulation of cognitive flexibility changes with age.
Introduction
Healthy aging is associated with prominent declines in multiple components of the brain's dopamine system, which have been linked to disrupted cognition (Wang et al., 1998; Bäckman et al., 2000; Mozley et al., 2001; Erixon-Lindroth et al., 2005). Even in healthy young adults, individual differences in dopamine profoundly influence executive functions, including cognitive flexibility (Cools et al., 2008; Stelzel et al., 2010; Samanez-Larkin et al., 2013). Dopamine therefore represents a prime target for investigating the neurochemical basis of variability in executive function in youth, and may represent a driver of powerful individual differences in cognitive aging. However, these relationships are rarely studied directly in humans because of difficulty assessing endogenous dopamine function. The present study measured dopamine synthesis capacity with the PET tracer 6-[18F]fluoro-l-m-tyrosine (FMT), and examined its relationship with age, cognitive flexibility, and frontoparietal activation during an fMRI study of task switching.
Task switching, an experimental paradigm requiring the rapid adjustment of behavior according to changing task rules, is a well-known measure of cognitive flexibility that is sensitive to brain dopamine (e.g., Samanez-Larkin et al., 2013). Previous task switching studies implicate dorsal striatum and specifically dorsal caudate nucleus (DCA) (for review, see Klanker et al., 2013), which receives rich dopaminergic input, and is highly connected with cortex via established cortico-striato-thalamo-cortical loops (Alexander et al., 1986). fMRI has further specified a role of lateral prefrontal cortex approximating BA 44/9, and posterior parietal cortex approximating BA 7 in task switching (Braver et al., 2003). Linking dopamine, task switching, fMRI activation, a D2 receptor polymorphism that elevates receptor density is associated with both increased switch costs and prefrontal inefficiency measured as higher BOLD activation relative to controls (Stelzel et al., 2010).
Although multiple components of the dopamine system decline in older adults, including loss of dopamine producing neurons in the substantia nigra (Fearnley and Lees, 1991), and reduced receptor and transporter density (Seeman et al., 1987; Volkow et al., 1994; Frey et al., 1996; Ishibashi et al., 2009), parallel decline in dopamine synthesis is not evident. Indeed, there is evidence that synthesis may be higher in older adults relative to young. FMT studies in humans and nonhuman primates show age-related increases in synthesis capacity (DeJesus et al., 2001; Braskie et al., 2008) in line with findings from postmortem tissue (Kish et al., 1995). However, little is known about the impact of upregulated synthesis on cognitive performance and underlying neural activity, although it has been linked to working memory disruption (Klostermann et al., 2012; but see Landau et al., 2009).
Cognitive flexibility is acutely sensitive to dopaminergic function as revealed by endogenous PET and pharmacological measures (Samanez-Larkin et al., 2013), with subtle augmentation having sometimes dramatic effects on cognitive function. The inverted-U-shaped model describes dopamine's influence on cognition, which proposes that shifts above or below an “optimal” dopamine tipping-point disrupt performance (for review of the inverted-U-shaped model of dopamine, see Cools and Robbins, 2004; Cools and D'Esposito, 2011). Although to date the majority of support for this model is derived from pharmacological studies, it may be possible to observe quadratic effects for in vivo dopamine measures using an individual differences approach. Invoking the inverted-U model, Dreher et al. (2008) reported that aging shifts relationships between dopamine and reward-related neural activity, as the positive correlation between dopamine synthesis capacity and prefrontal activation observed in young people became negative in older adults. However, the behavioral relevance of this switch is unclear, as activation was collected independent of a performance measure.
Here, we examined relationships between age, dopamine synthesis capacity, fMRI, and behavioral measures of task switching. We examined whether DCA dopamine's influence on cognition and neural function was described by quadratic, inverted-U-shaped relationships across all participants. Next, we tested the hypothesis that the higher levels of dopamine synthesis capacity observed in older adults would transform relationships between DCA dopamine and task switching, revealing the effect of aging on dopaminergic neural mechanisms of cognitive flexibility.
Materials and Methods
Participants.
Fifteen older adults (70–83 years old; 9 female) and 21 young adults (20–31 years old; 10 female) participated in FMT PET and MRI sessions and are included in the present report (Table 1). Older adults were recruited from the Berkeley Aging Cohort Study, and young adults were recruited via fliers. The study was approved by the Institutional Review Boards at the University of California, Berkeley and Lawrence Berkeley National Laboratory. All participants provided written consent, and received monetary compensation for participating in the study.
Demographics and task switching performancea
Participants had normal or corrected to normal vision, did not take medication affecting cognition, and had no neurological or psychological conditions. Older adults were characterized as cognitively normal using an interview and test battery, including several subtests of the Wechsler Adult Intelligence Scale-Revised (Wechsler, 1987), and the Mini Mental State Exam (Folstein et al., 1975). Older adults scored at least 25 on the Mini Mental State Exam (range: 25–30; mean ± SD, 28.33 ± 1.59), and were screened for current depression using the Geriatric Depression Scale (range: 0–8; mean ± SD, 3.36 ± 2.76) (Yesavage et al., 1982). Data for an additional older adult and young adult were excluded from the present analysis due to sleeping during the fMRI task and structural abnormality.
Older adults previously underwent Pittsburg compound B ([11C]PIB) PET imaging under a different protocol. Distribution volume ratios were low (range 0.91–1.09; mean ± SD, 1.00 ± 0.04) with 14 of 15 PIB negative (Villeneuve et al., 2015).
Neuropsychological examination: cognitive factor score.
Neuropsychological factor scores for episodic memory, working memory, and executive function/processing speed were generated for each participant using a maximum likelihood analysis (Bentler and Kano, 1990) standardized by test results from a larger sample of 346 Berkeley Aging Cohort Study participants as described previously (Schöll et al., 2016) (Table 1). Briefly, 14 cognitive tests met criteria for factor score generation as they were normally distributed but not excessively collinear (r < 0.7), and allowed for minimal data loss (data present for >90% of total sample). The final solution (three factors) showed significant goodness-of-fit (χ2 = 67.85, df = 42, p < 0.01, RMSEA < 0.05) explaining 56% of total model variance. Factor 1 was interpreted as episodic memory and was comprised of the California Verbal Learning Test (Delis et al., 2000), WMS-III Logical Memory Story (A + B1) and Visual Reproduction (Wechsler, 1987), Listening Span (Daneman and Carpenter, 1980), and Category Fluency (Spreen and Benton, 1977). Factor 2 was interpreted as working memory and was comprised of WMS-II Digit Span (forward and backward), and Arithmetic tasks (Wechsler, 1987). Factor 3 was interpreted as processing speed/executive function and was composed of the Digit-Symbol test (Smith, 1982), Stroop (Stroop, 1938), WMS-III Mental Control test (Wechsler, 1987), Trails A, and taps per second on a finger tapping task (Reitan and Wolfson, 1985). Factor scores were generated by z-scoring each variable and, for each participant, multiplying by the MLE-derived weight specific to that variable for each factor. Factor scores for three young adults could not be generated due to incomplete neuropsychological battery data.
PET data acquisition.
The present study used the FMT PET tracer to measure dopamine synthesis capacity. Similar to 6-[18F]fluorodopa (FDOPA), FMT is a substrate for aromatic amino acid decarboxylase, an enzyme in the dopamine synthesis pathway. Although not the rate-limiting step in synthesis, its activity provides an estimate of dopamine synthesis capacity when provided with enough substrate (DeJesus, 2003). FDOPA's signal, unlike FMT's, is affected by catechol-O-methytransferase (DeJesus, 2003), decreasing signal-to-noise. FMT is oxidized to 6-[18F] fluorohydroxyphenylacetic acid, and its signal is trapped in the presynaptic terminal, unlike FDOPA, which is subject to postrelease processing.
FMT was synthesized at Lawrence Berkeley National Laboratory using methods previously described (VanBrocklin et al., 2004). Participants ingested 2.5 mg/kg of carbdopa ∼1 h before scanning to minimize the peripheral metabolism of FMT. Data were acquired using a Biograph Truepoint 6 PET/CT (Siemens Medical Systems). After a short CT scan, participants were injected with ∼2.5 mCi of FMT as a bolus in an antecubital vein. Dynamic acquisition frames were obtained over 90 min in 3D mode (25 frames total: 5 × 1 min, 3 × 2 min, 3 × 3 min, 14 × 5 min). Data were reconstructed using an ordered subset expectation maximization algorithm with weighted attenuation, corrected for scatter, and smoothed with a 4 mm FWHM kernel.
Structural MRI data acquisition and processing.
Images were acquired on a 1.5T Magnetom Avanto System (Siemens Medical Systems) with a 12-channel head coil. Each participant was scanned using a T1-weighted MPRAGE protocol (TR = 2110 ms; TE = 3.58 ms; FA = 15°; matrix = 256 × 256; FOV = 256; sagittal plane; voxel size = 1 × 1 × 1 mm; 160 slices). MPRAGE scans were segmented using FreeSurfer version 5.1 (http://surfer.nmr.mgh.harvard.edu/). Cerebellar gray matter was used as the reference region in FMT PET analyses because this region shows very little tracer uptake. The most anterior ¼ of cerebellar gray was removed from the reference region to limit contamination of signal from the substantia nigra and ventral tegmental area. PET analyses focused on dopamine synthesis capacity measured in regions of interest (ROIs) in striatum. DCA, dorsal putamen (DPUT), and ventral striatum (VST) ROIs were manually drawn on each participant's MPRAGE as previously described (Mawlawi et al., 2001) using Mango software (http://ric.uthscsa.edu/mango/) (see Fig. 1a). Intrarater reliability was high for all regions. For ROIs of 5 participants drawn by 2 raters, the Sorensen-Dice coefficient ranged from 0.80 to 0.86, and the intraclass correlation coefficient ranged from 0.95 to 0.99.
PET data analysis.
FMT PET data were preprocessed using SPM8 software (Friston et al., 2011). To correct for motion between frames, images were realigned to the middle (12th) frame. The first five images were summed before realignment. Structural images (including posterior cerebellar reference region and striatal ROIs) were coregistered to PET images using the mean image of frames corresponding to the first 20 min of acquisition as a target. Graphical analysis for irreversible tracer binding was performed using Patlak plotting (Patlak and Blasberg, 1985) with posterior cerebellar gray matter used as the reference region. Ki images were generated from PET frames corresponding to 25–90 min, which represent the amount of tracer accumulated in the brain relative to the cerebellum. These images are comparable with Ki images obtained using a blood input function but are scaled to the volume of tracer distribution in the reference region.
ROI-based partial volume correction was applied (Rousset et al., 1998) to generate the best approximation of FMT signal originating from DCA, DPUT, and VST, and to minimize differences in Ki driven by age-related changes in gray matter, white matter, and CSF volume. To apply the partial volume correction in native space, we used FreeSurfer-generated ROIs for gray matter cortical and subcortical regions, white matter, and CSF with hand-drawn striatal ROIs substituting for the automated striatal segmentation.
Ki and partial volume corrected Ki data were submitted to repeated-measures ANOVA in SPSS version 23 (IBM) with factors region (DCA, DPUT, VST) and age group (young, older). Greenhouse-Geisser sphericity correction was applied as needed for reporting p values, but degrees of freedom are reported as integers in the text for easier reading. Effect sizes for repeated-measures ANOVA are reported using ηG2 (Bakeman, 2005), which gives smaller values than the frequently used ηP2, but is preferable as it reduces error when comparing across studies (Fritz et al., 2012). Post hoc t tests were conducted with effect sizes computed using Cohen's d for between-subjects effects and dz for within-subjects effects.
fMRI behavioral task and analysis.
Participants performed a two-condition task switching experiment (see Fig. 2a). The task was generated using PsychoPy 1.80.06 (Peirce, 2007). On each trial, a cue (square or diamond) was presented for 300 ms after which a number between 1 and 9 was presented for 2000 ms. For square cues, participants determined whether the number was odd (right button press) or even (left button press). For diamond cues, participants determined whether the number was >5 (right button press) or <5 (left button press). The number 5 was never presented. Participants were given 2000 ms to respond and were instructed to respond as quickly as possible without sacrificing accuracy. Fixation periods of duration 2000, 4000, or 6000 ms separated each trial. Responses were limited to two buttons, rather than 4 buttons, to introduce response overlap, which is thought to increase cognitive control demands (Meiran et al., 2000). Response hand and cue-rule mapping were counterbalanced across participant.
Participants performed 6 functional runs comprised of 48 trials each. Four runs contained both conditions while two runs contained just a single condition. For the two-condition runs, stimuli were pseudo-randomized to ensure the same number stimulus was not presented two times in a row, the same cue was not presented more than four times in a row, consecutive right versus left hand responses did not repeat more than three times in a row, and there were no more than three consecutive rule repeat trials or rule switch trials. Participants performed out-of-scanner practice with and without trial-to-trial feedback, until they reached 60% accuracy. In the scanner, participants performed a brief practice to ensure they could see the stimuli.
The primary performance measure was response time (RT), although accuracy was also calculated. Analyses presented here focus on “switch cost,” the difference in RT for task switch trials (correct responses preceded by correct responses for the different condition) versus task repeat trials (correct responses preceded by correct responses for the same condition). High switch costs, in this framework, reflect cognitive inflexibility. Switch cost was the primary measure because it has been linked reliably to dopamine function (e.g., Samanez-Larkin et al., 2013). However, age effects for switch cost are minimal relative to the large and consistent age effects found for “mixing cost,” the difference in RT for task repeat trials in the two-condition runs versus RT for trials in the single condition run (Wasylyshyn et al., 2011).
RT and accuracy data were submitted to repeated-measures ANOVA with factors trial type (switch, repeat) and age group (young, older). Although we do not focus the present report on mixing cost, we include an additional repeated-measures ANOVA with factors trial type (repeat, single task) and age group (young, older) to confirm the current task switching paradigm replicates the broader literature.
Additionally, we evaluated the relationship between switch cost (% change in RT for switch trials relative to repeat trials) and DCA FMT Ki. Influential models posit dopaminergic modulation of cognitive function is best described as an inverted-U (Cools and Robbins, 2004; Cools and D'Esposito, 2011). To test the applicability of this model, we applied linear, quadratic, and cubic polynomial fits to data describing the relationship between FMT Ki and switch cost. Adjusted r2 values are reported. To preview our results, we found that data were best described by a quadratic, U-shaped fit. As a test of the robustness of the quadratic fit, a permutation procedure was applied to assess the likelihood that our results could have occurred by chance. Briefly, residuals from the model without the quadratic term were randomly shuffled without replacement using the mrg33k3a pseudorandom number generator (L'Ecuyer, 1999) for 1000 permutations. For each new sample, regression was performed using the quadratic term as the independent variable. p values indicate the likelihood that observed t values for the quadratic term were greater than random associations generated from shuffled data. p values from permutation testing are reported for all significant quadratic models.
To test whether the relationships between FMT Ki and switch cost were significantly different between age group, the age × FMT Ki × switch cost interaction was tested. Follow-up analyses tested Pearson correlations for individual groups. Regression, interaction, and correlation analyses were conducted using R version 3.2.2 (http://www.R-project.org). Permutation testing was conducted using Python software version 2.7 (http://www.python.org).
fMRI data acquisition and preprocessing.
Six functional task runs and resting state scans were collected in the same session as the T1 MPRAGE. The six functional task runs were acquired with at T2*-weighted EPI sequence (TR = 2200 ms; TE = 50; flip angle = 50°; matrix 64 × 64; FOV = 220 mm; voxel size = 3.44 × 3.44 × 3.91 mm). A total of 28 axial slices oriented to the AC-PC line were acquired in ascending order to give whole brain coverage. The resting state scan was acquired before task performance with a T2* EPI sequence (TR = 2200 ms; TE = 50; flip angle = 90°; matrix 64 × 64; FOV = 220 mm; voxel size = 3.44 × 3.44 × 3.91 mm). Twenty-eight axial slices were acquired in an interleaved order. The first 10 volumes acquired for each functional scan were discarded to allow for magnetization preparation and stabilization. In-plane turbo inversion recovery images (TR = 2000 ms, TE = 11, matrix = 256 × 256; FOV = 220; voxel size = 0.43 × 0.43 × 3.91 mm, 28 slices) were collected as an intermediate target for registration of MPRAGE and EPI images.
fMRI data were analyzed using SPM12 software. Functional images were corrected for differences in slice timing and were realigned to the first volume to correct for head movement. To spatially normalize functional images to the MNI template, the in-plane and MPRAGE images were used as intermediates. Following normalization, functional images were smoothed with an 8 mm FWHM isotropic Gaussian kernel.
fMRI GLM.
Data were analyzed using a multisession GLM implemented in SPM12 (Friston et al., 1995). Task switch (high-low rule, odd-even rule) and task repeat trials (high-low rule, odd-even rule) were modeled as separate predictors. All omissions, error responses, and trials following errors were modeled together as a single predictor and are not included in the present analysis. Predictors were time-locked to the onset of the number stimulus. To mitigate the effect of motion artifact, six motion regressors derived from individual subject realignment were included in the model. Additionally, a parametric regressor was created for individual trial RT to control for possible time-on-task effects on frontoparietal activation (Grinband et al., 2008). For each participant, we used a parametric regression approach to control for the effect of differences in mean RT on differences in brain activity between switch and repeat trials. Specifically, we created an index for each trial by subtracting the mean RT across all conditions and dividing by SD (Yarkoni et al., 2009).
A priori ROI analyses.
fMRI studies consistently implicate lateral prefrontal and posterior parietal cortex in task switching. Given that these frontoparietal regions are thought to be engaged transiently to support rapid updating of task set during switch trials, we selected ROIs from a previous mixed blocked/event-related design study of task switching (Braver et al., 2003). Regions showing transient increases in activation during switch trials included left lateral prefrontal cortex approximating BA 44/9 (MNI −46, 14, 24) and left posterior parietal cortex approximating BA 7 (MNI −28, −70, 45), and have been closely replicated (Jimura and Braver, 2010). We tested whether increased activation for switch relative to repeat trials in these ROIs (average of 8 mm radius spheres centered on peak coordinates from Braver et al., 2003 converted from Talairach space) showed linear, quadratic, or cubic relationships with DCA FMT Ki (averaged across left and right hemisphere), and conducted interaction analyses as described above. Control analyses were conducted using ROIs in right and left primary visual cortex (8 mm spheres surrounding ±28, −96, −6; Anderson et al., 2011; Kiviniemi et al., 2009), and left and right primary motor cortex (8 mm spheres surrounding ±37, −25, 63; Mayka et al., 2006) for which we hypothesized no relationship between DCA dopamine and activation, and no interaction with age group.
Exploratory whole-brain analyses.
First, we examined whether there were age-related differences in activation between young and older adults using a voxelwise independent t test on the switch > repeat contrast.
Second, we identified regions significantly activated for the switch > repeat contrast for all participants with age included as a covariate. We tested whether increased activation for switch relative to repeat trials in these regions showed a linear, quadratic, or cubic relationship with DCA FMT Ki, and conducted interaction analyses as described above. We first conducted analyses on all significantly activated regions, with complementary analyses limited to peaks in lateral prefrontal and posterior parietal cortex as these regions are strongly associated with task switching independent of stimulus type (Kimberg et al., 2000; Sohn et al., 2000; e.g., Dove et al., 2000).
For voxelwise significance testing, an initial cluster forming threshold of p < 0.001 was applied, which is more conservative than the commonly used p < 0.01 (Woo et al., 2014; Eklund et al., 2016). An additional minimum cluster extent threshold was applied using Monte Carlo simulation implemented in REST toolbox version 1.8 AlphaSim (www.restfmri.net/forum/; cluster-level threshold p < 0.05) (Song et al., 2011). This principled extent-thresholding method accounts for spatial correlations of voxels and does not make assumptions about data distribution and smoothness.
Resting state functional connectivity analysis.
Resting state scans were preprocessed as described above. Connectivity analyses were performed using the Functional Connectivity Toolbox version 14.p (CONN; www.nitrc.org/projects/conn) (Whitfield-Gabrieli and Nieto-Castanon, 2012) in SPM12. Anatomical images were segmented into gray matter, white matter and CSF to create masks for signal extractions. Regressors of no interest included motion, white matter, and CSF. Data underwent linear detrending and were bandpass filtered (0.008–0.09 Hz).
We created four functional connectivity maps for each participant using seeds identified in the univariate analysis in right and left inferior frontal gyrus (8 mm spheres surrounding MNI: 54, 6, 10; −54, 14, 18) and right and left precuneus (8 mm spheres surrounding MNI: 12, −68, 52; −12, −68, 48). The four seed-to-voxel maps were averaged to create a composite frontoparietal connectivity map for each subject. Connectivity values were extracted from bilateral DCA and tested for relationships with DCA FMT Ki as described above. Parallel control analyses were conducted for connectivity values extracted from right and left primary visual cortex (8 mm spheres surrounding MNI: ± 28, −96, −6) and primary motor cortex (8 mm spheres surrounding MNI: ±37, −25, 63).
Results
PET
Older adults had elevated FMT Ki relative to young adults (Fig. 1b). Repeated-measures ANOVA showed a main effect of age (F(1,34) = 18.94, p < 0.001, ηG2 = 0.24), a main effect of region (F(2,68) = 154.45, p < 0.001, ηG2 = 0.58), and a region × age interaction (F(2,68) = 3.96, p = 0.03, ηG2 = 0.04). Post hoc analyses determined the region × age interaction was driven by most elevated signal in DCA for older adults relative to young (t(34) = 5.07, p < 0.001, d = 1.68; independent t test), compared with DPUT (t(34) = 2.77, p = 0.009, d = 0.93) and VST (t(34) = 1.98, p = 0.056, d = 0.68).
FMT signal is elevated in older adults relative to young. a, ROIs in bilateral DCA, DPUT, and VST were drawn on each participant's structural scan and were coregistered with FMT Ki images. Above are one participant's ROIs overlaid on Ki image. The scale reflects Ki values. b, Dopamine synthesis capacity was elevated in older adults relative to young (F(1,34) = 18.94, p < 0.001). Partial volume corrected (PVC) Ki values are displayed. Data are mean ± SEM.
Although we use partial volume corrected FMT Ki as our primary measure, noncorrected data showed a similar pattern of results. There was a main effect of age (F(1,34) = 19.36, p < 0.001, ηG2 = 0.26), a main effect of region (F(2,68) = 224.75, p < 0.001, ηG2 = 0.62), but the region × age interaction was not significant (F(2,68) < 1).
fMRI behavior
Analyses focused on the comparison of task repeat and task switch trials in the two-condition runs. Task performance generally replicated previous findings (Table 1). Both groups showed robust and equivalent effects of task switching evidenced by main effect of trial type (F(1,34) = 40.13, p < 0.001, ηG2 = 0.05) and no trial × age interaction (F(1,34) < 1). Older adults were slower than young adults overall (F(1,34) = 14.06, p = 0.001, ηG2 = 0.28).
Accuracy data showed a similar pattern. Accuracies were lower for task switch trials relative to task repeat trials, at trend level (F(1,34) = 3.62, p = 0.07, ηG2 = 0.004). There was no trial type × age interaction (F(1,34) = 1.41, p = 0.24), although older adults were less accurate overall (F(1,34) = 7.07, p = 0.01, ηG2 = 0.17).
To confirm the present task was sensitive to effects reported in the broader literature, we conducted additional repeated-measures ANOVAs to detect age differences in global switch cost. Participants were slower and less accurate for task repeat trials relative to single task trials (main effect of task RT: F(1,34) = 330.30, p < 0.001, ηG2 = 0.54; accuracy: F(1,34) = 35.49, p < 0.001, ηG2 = 0.25). Older adults were slower but not less accurate than young adults overall (RT: F(1,34) = 12.77, p = 0.001, ηG2 = 0.23; accuracy: F(1,34) = 2.69, p = 0.11, ηG2 = 0.04), and there was an age × task interaction (RT: F(1,34) = 7.59, p < 0.01, ηG2 = 0.03; accuracy: F(1,34) = 6.50, p = 0.02, ηG2 = 0.06), indicating larger increases in RT and deceases in accuracy for repeat trials in older adults.
Association between FMT and switch cost
We submitted data to linear, quadratic, and cubic regression models testing the relationship between DCA FMT Ki and switch cost. Combining groups, the linear relationship was not significant (r2 = 0.03, F(1,34) = 2.10, p = 0.16), although quadratic (r2 = 0.36, F(2,33) = 10.89, p < 0.0005) and cubic were (r2 = 0.37, F(3,32) = 7.98, p < 0.0005). Comparative testing across models confirmed the quadratic fit outperformed the linear (F(1,33) = 18.59, p = 0.0001), and cubic fit did not explain more variance than quadratic (F(1,32) = 1.54, p = 0.22). Permutation testing confirmed the stability of the quadratic effect (p = 0.001; 1000 permutations).
Considering young and older adults separately, ANOVA revealed an interaction between age, DCA FMT Ki, and switch cost, indicating that these relationships differed as a function of age group (F(1,32) = 9.03, p = 0.005) (Fig. 2b). High synthesis appeared optimal for young adults (Pearson's r = −0.61, p = 0.003) as greater synthesis capacity was associated with increased flexibility. However, at the higher synthesis seen in older adults, this relationship was not evident (r = 0.24, p = 0.39).
Relationship between FMT and switch cost. a, Participants performed a two-condition task switching paradigm. Cues indicated whether the odd/even rule or the >5/<5 rule applied. Switch cost was calculated as the percentage increase in response time for task switch trials relative to task repeat trials. b, There was a U-shaped quadratic relationship between dorsal caudate FMT Ki and switch cost (r2 = 0.36, p = 0.0002; y = 231,948.60x2 − 13,754.936x + 208.85). Dotted line indicates quadratic fit. Gray represents 90% confidence interval. There was a significant interaction between age group, FMT Ki, and switch cost (F(1,32) = 9.03, p = 0.005). Solid line indicates young adult linear fit line. Dashed line indicates older adult linear fit line.
fMRI a priori ROI analyses
There was no age-related difference in frontoparietal activation in ROIs in left lateral prefrontal cortex (8 mm sphere surrounding MNI −46, 14, 24) and posterior parietal cortex (8 mm sphere surrounding MNI −28, −70, 45) defined by Braver et al. (2003) (average both ROIs, t(34) < 1). The effect of DCA FMT Ki on frontoparietal activation was tested in linear, quadratic, and cubic regression models. Across groups there was no linear relationship between DCA FMT Ki and averaged frontoparietal activation (F(1,34) < 1), although there was a significant quadratic (r2 = 0.26, F(2,33) = 7.29, p = 0.002) and cubic relationship (r2 = 0.24, F(3,32) = 4.73, p = 0.007). The quadratic fit outperformed the linear (F(1,33) = 13.27, p = 0.0009), with no significant improvement with the cubic fit (F(1,32) < 1). Permutation testing confirmed the stability of the quadratic effect (p = 0.001; 1000 permutations).
Considering young and older adults separately, ANOVA revealed an interaction between age, DCA FMT Ki, and activation, indicating these relationships differed as a function of age group (F(1,32) = 11.22, p = 0.002; Fig. 3a). This interaction remained significant after controlling for performance (F(1,31) = 8.50, p = 0.006). Closely mirroring dopamine effects observed for behavioral switch cost, young adults showed tight coupling between DCA FMT Ki and activation (r = −0.61, p = 0.003), with no significant relationship at the higher synthesis levels in older adults (r = 0.38, p = 0.16).
Relationship between FMT and frontoparietal activation. a, There was a quadratic relationship between dorsal caudate FMT Ki (average left and right) and frontoparietal activation in averaged a priori ROIs (r2 = 0.26, p = 0.002; y = 17477.89 x2 − 1027.98x + 15.11). a–c, Dotted line indicates quadratic fit. Gray represents 90% confidence interval. There was a significant interaction between age group, FMT Ki, and activation (F(1,32) = 11.22, p = 0.002). Solid line indicates young adult linear fit line. Dashed line indicates older adult linear fit line. ROIs were 8 mm spheres surrounding inferior frontal gyrus (−46, 14, 24), and superior parietal lobule (−28, −70, 45) defined by Braver et al. (2003). b, Similarly, there was a quadratic relationship between dorsal caudate FMT Ki and averaged frontoparietal activation in ROIs defined by the present study (r2 = 0.28, p = 0.002, y = 12,153.36x2 − 715.18x + 10.74). There was a significant interaction between age group, FMT Ki, and activation (F(1,32) = 8.53, p = 0.006). ROIs were 8 mm spheres surrounding bilateral inferior frontal gyrus (54, 6, 10; −54, 14, 18) and precuneus (12, −68, 52; −12, −68, 48). c, There was no quadratic relationship between dorsal caudate FMT Ki and primary visual cortex activation (r2 = 0.04, p = 0.20), and no interaction between age group, FMT Ki, and activation (F(1,32) < 1). ROIs were 8 mm spheres surrounding primary visual cortex (MNI ±28, −96, −6). All coordinates are in MNI space.
We performed a parallel analysis using ROIs in primary visual and motor cortex where we did not hypothesize a relationship between DCA FMT Ki and activation nor an interaction with age group. There was no age-related difference in primary visual activation (t(34) < 1). Linear, quadratic, and cubic fits were not significant (all F values <1.67, p > 0.20), and there was no age group × DCA FMT Ki × activation interaction (F(1,32) < 1; Fig. 3c). Similarly for motor cortex, there was no age-related difference in primary motor activation (t(34) < 1). Linear, quadratic, and cubic relationships were not significant (all F values <1), and there was no age group × DCA FMT Ki × activation interaction (F(1,32) < 1). As a further test of the anatomical specificity of the significant interaction effects reported for a priori frontoparietal ROIs, ANOVA probed the interaction between age group, DCA FMT Ki, activation, and region. The four-way interaction was significant for both primary visual (F(1,64) = 4.61, p = 0.04) and primary motor cortex (F(1,64) = 4.01, p = 0.05).
fMRI exploratory analyses
Voxelwise analyses, including all participants, identified a number of frontoparietal regions with significantly greater activation for task switch trials relative to task repeat trials (Fig. 4; Table 2), although no regions showed significantly greater activation for the opposite contrast. Likewise, there were no significant voxelwise age group differences in activation for either contrast as assessed by independent t tests.
Univariate activation during task switching. Activation is displayed for the contrast switch > repeat trials for young and older adults. T-maps are displayed on CARET semi-inflated templates with a combined height threshold of p < 0.001 and 50 voxel extent threshold (AlphaSim cluster-level threshold, p < 0.05). Scale reflects t values.
Results from the BOLD-fMRI analysis of the contrast switch > repeat
Analyses, including all activated regions, revealed no relationship between DCA FMT Ki and activation for the linear regression model (F(1,34) < 1), but significant quadratic (r2 = 0.18, F(2,33) = 4.81, p = 0.01) and cubic relationships (r2 = 0.15, F(3,32) = 3.11, p = 0.04). Comparative testing across models confirmed the quadratic fit outperformed the linear (F(1,33) = 8.51, p = 0.006), and cubic fit did not explain more variance than quadratic (F(1,32) < 1). Permutation testing confirmed the stability of the quadratic effect (p = 0.02; 1000 permutations). ANOVA specifically testing the interaction of age group, DCA FMT Ki, and activation indicated these relationships differed as a function of age group (F(1,32) = 4.66, p = 0.04), although the interaction did not remain significant after controlling for performance (F(1,31) = 1.97, p = 0.17). Young adults showed a correlation between DCA FMT Ki and activation (r = −0.47, p = 0.03), with no significant relationship at the higher synthesis levels in older adults (r = 0.22, p = 0.43).
Analyses restricted to activation in lateral prefrontal and posterior parietal cortex showed a similar pattern of results. Follow-up analyses focused on these regions as they have been consistently implicated in task switching (Dove et al., 2000; Kimberg et al., 2000; Sohn et al., 2000). Activation extracted from ROIs (8 mm spheres) surrounding peaks in bilateral inferior frontal gyrus (MNI: 54, 6, 10; −54, 14, 18) and precuneus (MNI: 12, −68, 52; −12, −68, 48) was tested for age group differences, and was submitted to linear, quadratic, and cubic regression. There were no age-related differences in activation (t(34) < 1). Additionally, there was not a significant linear relationship between activation and DCA FMT Ki (r2 = 0.001, F(1,34) = 1.04, p = 0.31), although the quadratic (r2 = 0.28, F(2,33) = 7.76, p = 0.002) and cubic relationships (r2 = 0.27, F(3,32) = 5.29, p = 0.004) were significant. Comparative testing across models confirmed the quadratic fit outperformed the linear (F(1,33) = 14.07, p < 0.0007), and cubic fit did not explain more variance than quadratic (F(1,32) < 1). Permutation testing confirmed the stability of the quadratic effect (p = 0.001; 1000 permutations). ANOVA specifically testing the interaction of age group, DCA FMT Ki, and activation indicated that these relationships differed as a function of age group (F(1,32) = 8.53, p = 0.006), which remained significant after controlling for performance (F(1,31) = 5.62, p = 0.02; Fig. 3b). Similar to effects reported above, young adults showed tight coupling between DCA FMT Ki and activation (r = −0.56, p = 0.008), with no significant relationship at the higher synthesis levels in older adults (r = 0.32, p = 0.24).
Dopamine synthesis capacity and functional connectivity
Task-relevant networks can be observed in the pattern of functional activity measured at rest (Fox et al., 2006). We investigated the influence of dopamine synthesis capacity on the strength of intrinsic functional connectivity within the corticostriatal networks thought to subserve cognitive flexibility. We generated frontoparietal functional connectivity maps using seeds defined from peak activation identified in the voxelwise analysis in lateral prefrontal and posterior parietal cortex (seeds from Fig. 3b, top, network in Fig. 5a), and measured the connectivity strength between this frontoparietal network and DCA.
Relationship between FMT and frontoparietal-dorsal caudate functional connectivity. a, Frontoparietal resting state functional connectivity maps were generated using seeds identified in the univariate analysis for the contrast switch > repeat. Seeds were 8 mm spheres surrounding bilateral inferior frontal gyrus (54, 6, 10; −54, 14, 18) and precuneus (12, −68, 52; −12, −68, 48). T-maps are displayed on CARET semi-inflated templates. Scale reflects t values. Coordinates are in MNI space. b, There was a quadratic relationship between dorsal caudate FMT Ki (average left and right) and frontoparietal-DCA functional connectivity (r2 = 0.16, p = 0.02, y = −1295.57x2 + 74.75 − 1.04). There was not a significant interaction between age group, FMT Ki, and connectivity (F(1,32) = 1.46, p = 0.24). b, c, Dotted line indicates quadratic fit. Gray represents 90% confidence interval. Solid line indicates young adult linear fit line. Dashed line indicates older adult linear fit line. c, There was no quadratic relationship between dorsal caudate FMT Ki and frontoparietal-primary visual cortex functional connectivity (r2 = 0.02, p = 0.27), and no interaction between age group, FMT Ki, and functional connectivity (F(1,32) < 1).
Our findings were consistent with the patterns of effects observed for neural and behavioral measures of task switching. We found no age-related connectivity differences (t(34) < 1). There was no linear relationship between DCA FMT Ki and connectivity (F(1,34) < 1), although there were significant quadratic (r2 = 0.16, F(2,33) = 4.44, p = 0.02) and cubic (r2 = 0.15, F(3,32) = 3.10, p = 0.04) relationships. Data were best described by the inverted-U-shaped quadratic function as tested against linear (F(1,33) = 8.81, p = 0.005) and cubic (F(1,32) < 1) fits (Fig. 5b). Permutation testing confirmed the stability of the quadratic effect (p = 0.004; 1000 permutations). Examining groups independently, there was not an interaction with age group (F(1,32) = 1.46, p = 0.24). Relationships between resting frontoparietal-DCA connectivity and dopamine were not significant (young r = 0.26, p = 0.26; older r = −0.16, p = 57) possibly reflecting less consistent coupling with task-independent measures relative to task-related measures. There was, however, a modest correlation between connectivity and activation across both groups (all participants: r = −0.38, p = 0.02), with stronger connectivity associated with smaller increases in frontoparietal activation during switch trials.
To confirm that our findings were specific to regions associated with cognitive flexibility, we measured the connectivity strength between the frontoparietal network and primary visual and motor cortex. We found a marginal age-related difference in connectivity with visual cortex (t(34) = 1.73, p = 0.09) with a trend for greater connectivity in older adults. Linear (r2 = 0.05, F(1,34) = 2.68, p = 0.11), quadratic (r2 = 0.02, F(2,33) = 1.36, p = 0.27), and cubic relationships (F(3,32) < 1) with DCA FMT Ki were not significant, and there was no age group × DCA FMT Ki × connectivity interaction (F(1,32) < 1) (Fig. 5c). Similarly, for motor cortex, we found no age-related difference in connectivity (t(34) < 1). Linear, quadratic, and cubic relationships were not significant (all F values ≤1), and there was no age group × DCA FMT Ki × connectivity interaction (F(1,32) < 1).
Discussion
We show that striatal dopamine synthesis capacity (PET FMT Ki) is elevated in healthy older adults relative to young. Although higher levels of dopamine synthesis capacity appeared to benefit cognitive flexibility in young adults, these relationships were altered in older adults as evidenced by interactions between age, dopamine, and task switching measures. Our findings contribute to a growing body of research associating age-related differences in endogenous dopamine markers with cognition. To date, the majority of studies have focused on the detrimental effect of decreases in postsynaptic measures of D1 (Wang et al., 1998) and D2 receptor density (Volkow et al., 1998; Bäckman et al., 2000) as well as presynaptic measures of the dopamine transporter (Mozley et al., 2001; Erixon-Lindroth et al., 2005) on cognition. Our findings offer a new perspective on how dopamine synthesis is related to cognitive flexibility, and how these mechanisms underlie variability in cognition.
Although there has been clear demonstration of decreases in receptor and transporter density (Seeman et al., 1987; Volkow et al., 1994; Frey et al., 1996; Ishibashi et al., 2009) underlying cognitive decline in aging (Bäckman et al., 2000; Erixon-Lindroth et al., 2005; Bäckman et al., 2011), there is accumulating evidence supporting upregulation of dopamine synthesis (Kish et al., 1995; Lee et al., 2000; Braskie et al., 2008; Nandhagopal et al., 2011), which could serve a compensatory function. Despite limited sample size, the present study replicates previous results demonstrating elevated FMT Ki in healthy older adults relative to young (Braskie et al., 2008). Importantly, a multimodal approach allowed us to assess age-related differences in dopamine function independent of changes involving amyloid deposition. While findings using FDOPA, another tracer targeting aromatic amino acid decarboxylase, have been mixed (Martin et al., 1989; Sawle et al., 1990; Bhatt et al., 1991; Wolf et al., 1991; Dreher et al., 2008), methodological considerations may contribute to inconsistencies. For example, FDOPA's affinity is 10-fold weaker than FMT's and is affected at longer scan times (Sossi et al., 2002), by vesicular uptake, synaptic release, and metabolism, decreasing signal-to-noise.
Our findings reveal strong links between DCA dopamine and cognitive flexibility and define the neural mechanisms of its influence in young adults. Specifically, higher synthesis capacity was associated with lower switch costs and smaller activation increases during switch trials, suggesting that frontoparietal efficiency underlies enhanced performance (Stelzel et al., 2010). These findings were consistent for a priori regions selected based on their transient response during switch trials (Braver et al., 2003), as well as for exploratory analyses. Neocortical activation has been previously linked to striatal measures of dopamine synthesis capacity (Landau et al., 2009; Dang et al., 2012; Deserno et al., 2015) and D1 (Bäckman et al., 2011) and D2 receptor binding (Nyberg et al., 2009). The influence of dopamine signaling in DCA on frontoparietal cortex is likely mediated by parallel cortico-striato-thalamo-cortical loops (Alexander et al., 1986). Transient activation of this circuitry may allow for a release in thalamic gating, enabling updating of representations held in cortex (Frank et al., 2001) in support of task switching.
Across groups, we found that relationships between DCA dopamine and cognitive flexibility measures were best described by quadratic, U-shaped regression models. These findings suggest that optimal cognitive flexibility is achieved only at midrange levels of dopamine synthesis. However, scrutiny of dopamine relationships within group revealed nonsignificant correlations between dopamine and cognitive flexibility measures for older adults. Although the interpretation of the within-group null effects is not entirely clear in view of the limited sample size, they temper the interpretation that elevated synthesis capacity in aging is detrimental. What is clear, however, is that age-related increases in dopamine synthesis capacity are associated with altered dopamine-cognition relationships. Supporting our findings, similar interactions with age and dopamine synthesis capacity have been demonstrated for prefrontal cortex activation during reward anticipation and outcome (Dreher et al., 2008).
Resting functional connectivity findings point to modulation of network tuning as a possible mechanism by which dopamine affects performance and frontoparietal activation. DCA dopamine's quadratic fit with frontoparietal-DCA connectivity suggested stronger connectivity strength for midrange synthesis values, although strong interpretation of these effects warrants caution as the correlation between synthesis capacity and connectivity strength was not significant for young adults. However, this proposed mechanism is in agreement with previous demonstrations in young adults that dopamine depletion regimens weaken functional connectivity in these networks (Nagano-Saito et al., 2008), likely making their successful recruitment more effortful during task performance. Further support is established in the present dataset by the correlation between connectivity strength and activation such that greater connectivity was associated with smaller activation increases for switch trials.
Future studies must try to resolve the mechanisms underlying the changes in dopamine-cognition relationships with age. Key questions remain regarding the relationship between elevated dopamine synthesis capacity in aging and the loss of receptors, transporters, and gray matter volume. PET studies of Parkinson's disease suggest there are compensatory increases in dopamine synthesis, but that increases are small and are quickly outpaced by losses incurred during disease progression (Lee et al., 2000; Nandhagopal et al., 2011). In the present study, it is possible that most elevated dopamine synthesis capacity tracks advanced dysfunction in the system. Any compensatory benefit of elevated dopamine synthesis for older adults may have provided a boost to performance and neural measures that precluded a pure reversal of effects observed in young adults. Pharmacological studies may explore this possibility by resolving whether dopamine depletion impairs flexible performance in older adults with highest synthesis values, or, alternatively, whether it improves performance. While the present study defined clear relationships between striatal dopamine and cognitive flexibility in young adults, the complexity of the relationships between age, presynaptic and postsynaptic dopamine measures, and brain activity will require further investigation using integrated neuroscience approaches to disentangle.
Footnotes
This work was supported by National Institutes of Health Grants AG044292 and AG047686. We thank the reviewers for valuable comments that greatly benefited this manuscript.
The authors declare no competing financial interests.
- Correspondence should be addressed to Dr. Anne S. Berry, Helen Wills Neuroscience Institute, University of California Berkeley, Berkeley, CA 94720. aberry{at}lbl.gov