Introduction {#Sec1}
============

The last 10--15 years have witnessed an explosion of studies examining whether working memory (WM) performance can be enhanced through various training procedures. Different approaches have been employed in this work, ranging from training of specific components of WM such as updating^[@CR1]--[@CR3]^, shifting^[@CR3],\ [@CR4]^, and inhibition^[@CR5]^ to composite programs that involve practicing several WM functions^[@CR6],\ [@CR7]^. Typically, the training runs for 4--5 weeks, with 3 sessions per week, each lasting for 45--60 min. Several recent reviews^[@CR8]--[@CR11]^ give rise to the following conclusions: (a) WM training results in gains on the trained tasks; (b) transfer effects to tasks tapping untrained cognitive domains (e.g., reasoning, intelligence, multitasking) are small or non-existent; and (c) although transfer to untrained WM tasks that share constituent processes with the trained WM tasks may be seen, generalizability to other WM tasks is also quite limited.

Several fMRI studies examining neural correlates of training-related WM gains report increased blood-oxygen-level dependent (BOLD) activity in striatum^[@CR1]--[@CR3],\ [@CR12]^. This is an interesting observation in view of the fact that striatal neurons are thought to serve a gating function in letting new information enter into WM^[@CR13]--[@CR15]^. The functional relevance of increased striatal BOLD activity following WM training was demonstrated by Brehmer and colleagues^[@CR16]^, who reported a strong relationship between the magnitude of training-related increases in striatal BOLD activity and the size of WM improvement post training.

Animal^[@CR17],\ [@CR18]^ and human^[@CR19],\ [@CR20]^ data reveal a link between BOLD activity and measures of dopamine (DA) release. These associations open up for the possibility that increased striatal BOLD activity after WM training is related to an increased release of DA. To infer DA release during task performance, typically quantified using positron emission tomography (PET) and specific radioligands, two conditions that vary in cognitive demands are contrasted. The idea is that binding of the radioligand to DA receptors should be reduced during the more challenging condition relative to the control condition. This is so because, in the former case, binding of the ligand to receptors competes with binding of endogenous DA to the same receptors to a greater extent than in the control condition. Thus, reduced binding is assumed to reflect an increased release of DA. This displacement principle was initially formulated in the context of pharmacological DA challenges^[@CR21]^, and evidence for displacement has been observed in both striatal and extrastriatal regions during verbal^[@CR22]^ and spatial^[@CR23],\ [@CR24]^ WM tasks.

Using the same updating training procedure as Dahlin *et al*.^[@CR1]^, Bäckman *et al*.^[@CR19]^ conducted a PET study that investigated whether 5 weeks of WM training was associated with increased striatal DA release. In that study, the radioligand raclopride was used to measure DAD2 receptor binding. There were two main PET findings: (a) increased DA release was observed during a trained updating task relative to a low-level control task before training, and (b) WM training resulted in a further increase of DA release. The locus of the DA effects was in left caudate, close by the training-related BOLD changes reported by Dahlin *et al*.^[@CR1]^ However, as this is the only study documenting increased DA release following WM training, replication of this finding is warranted. Replication of intriguing data is a rare happening in many fields, including cognitive neuroscience. When replication is attempted, a considerable number of findings reported in journals like *Science* and *Nature* cannot be reproduced^[@CR25],\ [@CR26]^. In the current work, a major goal was therefore to seek to replicate the Bäckman *et al*.^[@CR19]^ findings. The training procedure, the criterion task, as well as the PET assessment of DA release were identical to those employed by Bäckman *et al*.^[@CR19]^.

As noted, transfer effects from this type of WM training to untrained tasks are small or non-existent. However, we observed transfer from WM updating training to an untrained n-back task that shares a similar demand on updating operations^[@CR1],\ [@CR19]^. An interesting finding in Dahlin *et al*.^[@CR1]^ was that transfer effects were also observed in terms of neural activation patterns, as measured by BOLD fMRI. Specifically, increased BOLD activity was seen for both the letter-memory and n-back tasks in an overlapping part of the caudate. That said, the magnitude of the increases in the n-back task was considerably smaller compared to the letter-memory criterion task, both behaviorally and neurally. An additional objective of the present research was to investigate whether updating training would result in increased DA release also in an untrained n-back task, thus providing evidence for a neurochemical correlate of transfer of learning.

Results {#Sec2}
=======

Letter-memory criterion task {#Sec3}
----------------------------

Two participants in the updating training group and one participant in the control group were not included in the letter-memory analysis due to technical problems leading to loss of post-training behavioral data. The main analysis of the letter-memory criterion task revealed a significant group x time interaction (*F*(1, 23) = 24.579, *p* \< 0.001, η^2^partial = 0.52; *d* = 2.07), reflecting the fact that the training group showed larger performance gains after training than the control group (Fig. [1A](#Fig1){ref-type="fig"}). An identical analysis involving only those participants who were included in the PET analyses yielded the same significant interaction (one participant could not be included in this analysis due to a technical error). Thus, the present behavioral findings replicated exactly those from our two previous studies^[@CR1],\ [@CR19]^.Figure 1Behavioral and dopamine effects for the letter-memory criterion task. (**A**) Percent pre-post changes in letter memory. (**B**) Lower raclopride binding to striatal D2 receptors during letter memory compared to the control task before training in right striatum, reflecting greater DA release in response to the cognitive challenge (x,y,z = 14,−3,18; 30,−3,−2). (**C**,**D**) A training-induced decrease of raclopride binding to D2 receptors was found in left striatum (−18, 6, −12). (**E**) Training-related changes in raclopride binding to striatal D2 receptors in the peak region observed by Bäckman *et al*.^[@CR18]^. Error bars are standard errors.

In all analyses of the PET data, we only report effects that exceed 5 voxels to reduce the risk of spurious findings. We first examined a task effect on D2DR BP comparing the letter-memory and baseline tasks before training across all participants with a paired *t* test. We observed decreased D2DR binding during letter memory in bilateral striatum (x,y,z = 14,−3,18; 30,−3,−2, *p* \< 0.01; Fig. [1B](#Fig1){ref-type="fig"}; x,y,z = −17,21,0; −23,2,11, *p* \< 0.05). Critically, an ANOVA on the D2DR BP data yielded a significant group x time interaction in bilateral striatum (x,y,z = 26,−5,8; −18,6,−12, *p* \< 0.05). Figures [1C and D](#Fig1){ref-type="fig"} highlight this effect for left striatum. This interaction reflected the fact that the trained persons, but not the controls, showed decreased BP after training. In a peak-voxel analysis of a potential training-related change in the exact same region as originally observed (x,y,z = −17,−7,22;^[@CR18]^), a similar trend of greater reduction of BP in the trained group compared to the control group was observed, although this effect did not reach conventional significance (Fig. [1E](#Fig1){ref-type="fig"}; *p* \> 0.05).

Transfer tasks {#Sec4}
--------------

For several transfer tasks (in-scanner digit n-back and the offline tasks spatial n-back, letter-number sequencing, digit symbol, word recall, and number-letter RT), there were main effects of time (*p*s \< 0.05), reflecting general retest effects. However, for the in-scanner 3-back task (Fig. [2A](#Fig2){ref-type="fig"}) as well as for all off-line transfer tasks (Table [1](#Tab1){ref-type="table"}), there were no disproportionate time effects as a function of group (*p* \> 0.05; *d* = 0.00). Thus, no behavioral transfer effects were observed.Figure 2Behavioral and dopamine effects for the n-back transfer task. (**A**) Percent pre-post changes in 3-back. (**B**) Lower raclopride binding to striatal D2 receptors during 3-back compared to 1-back before training in bilateral striatum, reflecting greater DA release in response to the cognitive challenge (x,y,z = 27,0,0; 11,18,6; −23,0,8; −12,15,11). (**C**) Effects in bilateral striatum (x,y,z = 11,14,−9; −17,15,14) showing a training-induced decrease of raclopride binding to D2 receptors during 3-back. (**D**) Bar graph showing a selective training-related reduction in right striatum for trained subjects. Error bars are standard errors. Table 1Group means (*SD*) on the pre- and post-training cognitive measures.WM training groupControl groupPre-trainingPost-trainingPre-trainingPost-trainingLetter-memory accuracy %*n* = 12*n* = 1351.54 (11.63)83.33 (31.67)47.57 (12.87)56.28 (17.01)Digit n-back*n* = 10*n* = 12 1-back target accuracy %91.23 (7.78)91.58 (7.71)93.62 (3.05)89.14 (8.42) 3-back target accuracy %74.53 (18.24)86.12 (14.51)79.44 (20.66)91.10 (7.30)Letter-number sequencing*n* = 14*n* = 1413.50 (2.79)14.79 (2.61)14.00 (3.21)14.64 (3.84)Digit span forward*n* = 13*n* = 147.00 (1.23)7.38 (0.87)6.43 (1.22)6.36 (1.15)Digit span backward*n* = 14*n* = 146.15 (1.21)6.08 (1.50)6.36 (1.50)6.36 (1.28)Digit symbol*n* = 14*n* = 1471.00 (12.27)77.29 (13.98)68.07 (11.53)75.43 (11.65)Word recall accuracy %*n* = 14*n* = 1475.13 (10.54)80.69 (11.84)80.16 (15.70)82.14 (13.33)Simon task*n* = 14*n* = 13 Congruent RT413 (45)381 (49)430 (58)400 (42) Incongruent RT434 (38)398 (43)460 (48)420 (35) Congruent accuracy %99.00 (1.52)98.86 (1.51)98.62 (1.89)98.31 (1.97) Incongruent accuracy %96.71 (2.16)97.71 (2.33)96.92 (2.78)97.69 (2.29)Number-letter task*n* = 14*n* = 14 Mixed task: no-switch RT680 (104)599 (78)767 (116)642 (97) Mixed task: switch RT955 (288)816 (234)1115 (310)866 (193) Mixed task: no-switch accuracy %98.94 (2.00)99.70 (1.14)98.48 (2.27)98.78 (1.37) Mixed task: switch Accuracy %97.99 (3.15)97.99 (3.15)96.21 (3.51)96.21 (3.50)Spatial n-back*n* = 14*n* = 13 1-back accuracy in %96.80 (2.81)98.31 (2.09)96.96 (2.91)98.08 (1.67) 2-back accuracy in %94.84 (5.15)97.82 (2.97)93.27 (5.67)98.40 (2.53)

However, as with letter-memory, comparing the 3-back and 1-back tasks before training across all participants with a paired *t* test, we observed decreased D2DR binding during 3-back in bilateral striatum (Fig. [2B](#Fig2){ref-type="fig"}; x,y,z = 27,0,0; 11,18,6; −23,0,8; −12,15,11, *p* \< 0.01). Importantly, a group x time ANOVA on the D2DR BP n-back data yielded a significant interaction in bilateral striatum (Fig. [2C](#Fig2){ref-type="fig"}; x,y,z = 11,14, −9; 9,14,8; −17,15,14, *p* \< 0.05). The interaction in right striatum reflected the fact that the training group showed decreased BP after training, whereas the control group's level was stable (Fig. [2D](#Fig2){ref-type="fig"}).

Discussion {#Sec5}
==========

The main goal of this study was to replicate the Bäckman *et al*.^[@CR19]^ finding of increased DA release following WM updating training. First, replicating and extending prior work^[@CR19],\ [@CR22]--[@CR24]^, DAD2 BP was lower during both letter memory and n-back than in the control conditions before training. This pattern indicates increased DA release as a function of the cognitive challenge. A key observation was reduced striatal DAD2 BP during letter memory in the training group at post test, suggesting increased DA release after training. The intervention effect on neurotransmission was accompanied by pronounced training-related behavioral gains in the letter-memory criterion task^[@CR1],\ [@CR19]^. Although the DA effect did not overlap with the peak effect obtained in our previous study, both were located in left striatum. Also, a similar trend was seen in the same peak voxel as originally reported^[@CR19]^.

A unique feature of this study was the use of four PET examinations, allowing us to assess neurochemical effects during an untrained n-back task. Interestingly, our training was associated with increased striatal DA release at the post-training assessment also during the 3-back transfer task. To the best of our knowledge, this is the first demonstration of changes in neurotransmission in relation to transfer of learning. For n-back the peak response was in right striatum, whereas the corresponding effect for letter memory was in left striatum. The reason for the difference in laterality remains unclear, and was not expected from our previous observation of a training-related overlap in left caudate BOLD signal change for letter memory and n-back^[@CR1]^.

In agreement with several recent meta-analyses and qualitative reviews^[@CR8]--[@CR11]^, we found no evidence that updating training transferred to performance on the tasks included in the off-line battery tapping verbal and spatial working memory, motor speed, episodic memory, set shifting, and inhibitory control. Moreover, unlike the findings reported by Dahlin *et al*.^[@CR1]^ and Bäckman *et al*.^[@CR19]^, we did not observe behavioral generalization to the 3-back transfer task. The reason thereof may be that the present participants performed at a very high level for 3-back already at the pre-training assessment. The high level may reflect the length (more room for task practice) and the nature of the task (performed during scanning, which likely increases the perceived performance demands compared to an off-line testing situation). Note also that recent meta-analyses indicate that effect sizes for this form of near-transfer are quite small^[@CR9]--[@CR11]^. Thus, it is unsurprising that such transfer effects may or may not be observed in a specific study. This is especially true in PET research with relatively small sample sizes and limited statistical power. In view of this concern and even though the measurement scale might not have been sensitive enough to capture behavioral transfer effects, the PET data nevertheless demonstrated DA alterations during the untrained n-back task. The present pattern of results raises the possibility that increased DA release may be a necessary, but not sufficient, condition for behavioral transfer to occur.

The present research extends previous observations that the DA system is plastic. Such plasticity has been observed during both pharmacological and cognitive challenges^[@CR18],\ [@CR20],\ [@CR21],\ [@CR27]--[@CR31]^. The current replication of increased striatal DA release during letter memory after WM updating training, and the novel result of a corresponding increase during an untrained 3-back task, provide additional evidence that the DA system is malleable.

Methods {#Sec6}
=======

Participants {#Sec7}
------------

The effective sample included 28 right-handed healthy, non-smoking, and non-medicated Finnish male university students (19--26 years). The participants underwent structural MRI and medical screening. They were randomized into training and control groups (*n* = 14 for each). The two groups were comparable regarding years of education, age, and BDI-II scores, and on all pre-training neuropsychological measures, except for TMT B, *t*(25) = 2.35, *p* \< 0.05, where the control group showed significantly worse performance (Table [2](#Tab2){ref-type="table"}).Table 2Descriptive data on the study groups (means and *SD*).WM training group (*n* = 14)Control group (*n* = 14/13)Age22.21 (1.72)22.79 (1.48)Years of education14.21 (1.19)15.15 (1.45)BDI-II2.64 (3.13)1.79 (1.85)WAIS-III Vocabulary52.50 (6.97)49.86 (7.82)TMT A21.64 (5.23)22.93 (4.57)TMT B42.64 (6.79)50.69 (10.74)Pattern Comparison19.61 (3.70)20.04 (3.20)Number Copying51.46 (5.76)51.04 (6.25)Note. BDI = Beck Depression Inventory, WAIS = Wechsler Adult Intelligence Scale, TMT = Trail Making Test. Age in years, BDI-II scores, WAIS-III Vocabulary raw scores, TMT completion time in sec, Pattern Comparison and Number Copying raw scores (*SD*s in parentheses). One participant in the control group had missing data on years of education, one participant in the control group had unreliable performance data on the Pattern- Comparison task, and one participant in the control group was an extreme outlier in TMT B (hence, *n = *13 in the control group for these variables).

The study was approved by the Ethics Review Board of the Turku University Hospital District, the methods were carried out in accordance with relevant guidelines and regulations, and written informed consent was obtained from all participants.

Procedure {#Sec8}
---------

We employed a pre-training--intervention--post-training control group design. All participants took part in the pre- and post-training assessments, whereas only participants in the training group received training between these assessments. The training group practiced three times per week (45 min/session) for five weeks. The pre-training assessments included a structural MRI scan, neuropsychological testing, and two consecutive PET scans performed during the same day. The post-training assessments included neuropsychological testing and two PET scans, but no MRI. Time in-between the pre- and post-training assessments was 6--9 weeks.

The training program was a computerized in-house developed Visual Basic program and consisted of the letter memory criterion task and five other updating tasks^[@CR1],\ [@CR19]^. Four of the additional training tasks were similar to the letter memory task and involved updating of single items, but to foster generality different kinds of stimuli (i.e. numbers, letters, colors, and spatial locations) were used. In these training tasks, five lists of items were randomly presented and the task was to recall the four last presented items. Across the 5-week training period, list length varied to manipulate level of difficulty and thereby ensure that the training was sufficiently demanding (low level = 4--7 items; medium level = 6--11 items; high level = 5--15 items). Performance was monitored and at the end of each training week, level of difficulty was raised when the participant scored 80% or higher in the letter memory task. All subjects reached the most difficult level by the end of training. The final training task was a keep-track task. This task, too, taxes updating, but is structurally different from the other updating tasks in the training battery. The inclusion of the keep-track task should contribute toward strengthening of a general updating skill. In each trial, 15 words from different semantic categories were presented serially in random order (2.0 s per word) and participants were instructed to place the words into categories (animals, articles of clothing, countries, relatives, sports, professions), indicated by boxes at the bottom of the screen. They had to continuously update their working memory content and remember the last presented word in each category at the end of the presentation. Participants responded by typing the last presented word under each category box when the trial ended.

PET scanning {#Sec9}
------------

The pre- and post-training PET scans were identical: both involved two scans of which one was performed in the morning and the other in the afternoon. The letter-memory task was administered during the morning scan, whereas the digit n-back task was administered during the afternoon scan. The two assessments were separated by a 90-min lunch break.

During the first PET scan, we administered a computerized letter-memory task that taps verbal WM updating. Participants were shown 7--15 letters in a sequence; when a sequence suddenly ended, they were asked to report the last four letters in correct order by pressing buttons corresponding to A, B, C, and D. This task was preceded by a structurally equivalent control task not taxing updating (all letters in a sequence were identical and participants reported that letter). Prior to PET scanning, participants acquainted themselves with the task during a practice session. The PET session started with the control task (5--10 min prior to bolus injection continuing for 55 min post-injection), followed by the letter-memory task (25 min). The sequence began with instructions shown for 7000 ms. Stimulus duration was 2000 ms with a fixation cross in-between for 1000 ms. The participants were allowed 7000 ms to respond.

During the second PET scan, the participants performed a digit n-back task that measures verbal WM updating. In this task, they were to determine if a currently visible digit was the same as the previous digit (1-back) or the digit that was presented three steps back (3-back). The dependent measure was target accuracy on the 3-back vs. 1-back condition. After participants had read the task instructions and completed one 1-back and one 3-back practice task, the PET session was initiated. During scanning, they completed about 55 min of the 1-back (baseline) task, after which they immediately continued with a period of 25 min performing the 3-back task. Task instructions, informing the participant of which n-back version to perform, were shown at the beginning of each task version for 5000 ms, after which a digit was shown for 1500 ms. The digit was replaced by a fixation point that was visible for 450 ms, and after this the fixation point was once again replaced by a digit. The alternation of fixation points and digits continued until the end of the task. Targets encompassed 45% of the trials (i.e., the current digit and the digit n steps back matched), whereas 55% comprised non-targets. The participants responded by pressing a match button for targets and a non-match button for non-targets.

The study protocol and procedure was virtually identical to those employed by Bäckman *et al*.^[@CR19]^. The only differences between the studies concerned the duration of the letter- memory and digit n-back tasks performed during the PET scans. Unlike Bäckman *et al*.^[@CR19]^, the current letter-memory task included a 55-minute baseline phase and a 25-minute test phase, with no return to baseline. In other respects, the letter-memory tasks were identical. The digit n-back task was included in the neuropsychological off-line battery in Bäckman *et al*.^[@CR19]^, whereas here it was altered (i.e., much longer) and performed during PET scanning to mirror the design of the letter-memory task.

Imaging methods {#Sec10}
---------------

We used HRRT-PET (High Resolution Research Tomograph; Siemens Medical Solutions, Knoxville, TN, USA) and \[^11^C\]raclopride with bolus-plus-infusion^[@CR32]^ to measure striatal dopamine D2 receptor availability during WM vs. control task performance. Raclopride was prepared from \[^11^C\]methyl triflate following an established procedure^[@CR33]^. Its radiochemical purity and specific radioactivity were determined using high-performance liquid chromatography and ultraviolet detection at 214 nm. The aimed magnitude of the bolus was 50% of the total tracer volume^[@CR34]^, and constant infusion continued until 80 min after bolus injection (Kbol = 80 min). At scan start, the actual bolus injection and infused doses did not differ between groups or scans (Table [3](#Tab3){ref-type="table"}). Emission list-mode data were histogrammed into 3D sinograms in 20 time frames of variable length (8 × 2 min, 4 × 3 min, 2 × 4 min, 1 × 5 min, 1 × 6 min, 1 × 8 min, 3 × 8.3 min), taking declining radioactivity and task timing into account. Before each emission scan, a transmission scan was performed using a ^137^Cs point source. Tissue-attenuation maps were reconstructed using the maximum a posteriori transmission data (MAP-TR) algorithm with segmentation. Scattered events were estimated using the single-scatter simulation algorithm, and randoms were estimated from the block singles with a variance-reduction algorithm. All corrections were applied within the statistical image-reconstruction algorithm^[@CR35]^.Table 3Administered radioactivities (MBq; mean ± *SD*) in fast bolus and constant infusion formulations, and *p*-values in *t*-tests across groups and scans.BolusInfusionBolusInfusionBolusInfusionControl groupTraining group*t*-test (groups)Letter memoryPre Training262 ± 18 (*n* = 9)252 ± 14272 ± 31 (*n* = 13)255 ± 280.430.77Post training267 ± 19 (*n* = 10)253 ± 15280 ± 25 (*n* = 11)265 ± 250.520.22*t* test (scans)0.760.720.880.43n-backPre training267 ± 36 (*n* = 10)249 ± 32262 ± 29 (*n* = 13)264 ± 250.740.22Post training265 ± 24 (*n* = 10)260 ± 17269 ± 17 (*n* = 13)264 ± 160.440.51*t* test (scans)0.930.430.480.96

Special attention was paid to minimization of head-motion artifacts based on earlier observations of its detrimental effects on the PET signal and its interpretation^[@CR36]^. Specifically, head motion was minimized using an individually molded thermoplastic mask, and image reconstructions were made using an in-house version of the multiple-acquisition frame (MAF)^[@CR36]^ based motion-compensated image reconstruction algorithm. The algorithm employs external motion tracking (MT) as given by Vicra (Northern Digital) infrared detection, to define motion-free (amplitude less than 2.5 mm) framing of list-mode PET data. This procedure compensates for attenuation map misalignments using mutual information- based image registration and finally combines the motion-free subframes into original, desired framing. In a phantom study by Johansson *et al*. \[in preparation\], we showed good performance of the MAF-based algorithm, but a residual bias was seen in the presence of rapidly oscillating motion. Rapidly oscillating, high amplitude motion is not often observed in human PET scans but, if present, may introduce erroneous perturbations in the PET signal, which in turn could be incorrectly interpreted as activation effects^[@CR36],\ [@CR37]^.

Movements that are beyond the capabilities of the motion-compensation algorithm were sought through analysis of the number of subframes, which was identified as an important factor in our previous phantom experiment with regard to the algorithm's performance \[Johansson *et al*., in preparation\]. Analysis showed that in 84% (89/106) of the sessions, no additional acquisition frames were generated, indicating subthreshold (2.5 mm) amplitude of motion in all frames. The remaining cases were investigated further through detailed assessment of the time-motion graphs representing three translations and three rotations of the motion-recording target at the forehead at 1 Hz frequency. On the basis of visual inspection of the motion data, five letter-memory scans from four individuals, and three n-back scans from two individuals were deemed unrecoverable using the motion-compensation algorithm, and were thus excluded from analysis. In addition to those excluded due to extensive motion, the letter-memory PET variables of four subjects were missing due to other technical reasons. In one case, the ^11^C-raclopride infusion was started approximately 20 min overdue, in one case the PET data acquisition failed, and two participants were not PET scanned post training. Thus, a total of eight participants were missing from the PET letter-memory analysis, resulting in an *n* of 20. In turn, the n-back PET variables were missing from three subjects, two of whom were not scanned after training, and in one case the data acquisition failed. Hence, altogether five participants were missing from the n-back analysis, resulting in an *n* of 23.

Structural MR imaging was performed for excluding anatomical abnormalities and for anatomical reference of the PET data. T1-weighted MRI data were acquired using a 3T scanner (Philips Ingenuity TF PET/MR) with a 1 mm × 1 mm × 1 mm voxel size covering the whole brain. The T1-weighted MRI data were pre-processed using FreeSurfer (Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, US) to obtain automated region-of-interest (ROI) definitions^[@CR38]^. Subsequently, a PET sum image was registered to the skull-stripped MR image, using the mutual-information optimization algorithm in Statistical Parametric Mapping (version 8, SPM8; Wellcome Institute, London, UK). Spatial normalization parameters to the Montreal Neurological Institute (MNI) space coordinates were determined on the basis of MRI data using the unified-segmentation algorithm in SPM8^[@CR39]^ yielding a mapping for registered PET data normalization. The PET-image analysis was restricted to striatum, which serves as a central hub of dopaminergic activation. Striatal subvolumes were defined on the basis of FreeSurfer segmentations, combining segments of caudate, putamen, and nucleus accumbens into one composite striatal structure in each subject's individual space. Striatal masks were used in edge-preserving spatial smoothing of the PET data; outside-mask voxels were omitted in the Gaussian filter kernel (10 mm FWHM) calculation, attenuating the impact of partial-volume effects near the edges. Extensive spatial smoothing was deemed necessary for successful calculation of binding potential (BP) from the PET activation data.

Modeling of the dynamic PET data was based on an extension of the conventional simplified reference tissue model (SRTM)^[@CR40]^, including the effect of activation as an additional parameter. A similar approach as described by Alpert and colleagues^[@CR41]^ was adopted, where the activation effect on the ^11^C-raclopride signal was modeled by the time-dependent activation function *h*(*t*), and the relevant linear equations are:$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$[\begin{array}{c}{C}_{T}({t}_{1})\\ \vdots \\ {C}_{T}({t}_{m})\end{array}]=[\begin{array}{c}{C}_{R}({t}_{1})\,\\ \vdots \\ {C}_{R}({t}_{m})\end{array}\begin{array}{c}\,{\int }_{0}^{{t}_{1}}{C}_{R}(u)du\,\\ \vdots \\ \,{\int }_{0}^{{t}_{m}}{C}_{R}(u)du\,\end{array}\begin{array}{c}\,-{\int }_{0}^{{t}_{1}}{C}_{T}(u)du\,\\ \vdots \\ \,-{\int }_{0}^{{t}_{m}}{C}_{T}(u)du\end{array}\begin{array}{c}\,-{\int }_{0}^{{t}_{1}}{C}_{T}(u)h(u)du\\ \vdots \\ \,-{\int }_{0}^{{t}_{m}}{C}_{T}(u)h(u)du\end{array}]\,[\begin{array}{c}\begin{array}{c}{R}_{1}\\ {k}_{2}\end{array}\\ {k}_{2a}\\ \gamma \end{array}]$$\end{document}$$where the weighting factor $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\,\gamma $$\end{document}$ is depicted as the activation effect, $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${C}_{T}({t}_{i})$$\end{document}$ and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${C}_{R}({t}_{i})$$\end{document}$ are the instantaneously measured radioactivity concentrations in target and reference tissue, respectively, and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${R}_{1}$$\end{document}$, $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${k}_{2}$$\end{document}$, $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${k}_{2a}$$\end{document}$ are the conventional SRTM rate constants. Compared to our previous work^[@CR19]^, we now employed an explicit BP calculation algorithm for the activation phase, to quantitatively assess the magnitude of this effect relative to baseline. First, baseline BP (BP~0~) was determined from the fit parameters $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${k}_{2}$$\end{document}$, and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${k}_{2a}$$\end{document}$, which in turn were estimated using a linear least-squared solver in Matlab (version R2011b, Mathworks, US). Second, the activation BP (BP~1~) was determined from the fit parameters $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${k}_{2}$$\end{document}$, $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${k}_{2a}$$\end{document}$, and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\gamma $$\end{document}$, using an integral of $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\gamma h(t)$$\end{document}$ and averaged over the activation period. As the activation function, we employed a simplified gamma-variate function^[@CR42]^, with a fixed peak time of 8 min and 20 sec from active task initiation, λ = 3, and onset of activation at 55 min from the injection. Model calculation was restricted to the aforementioned striatal volume, and the resulting parametric images were normalized using the MRI-based deformation into MNI space coordinates. Voxel-wise statistical analysis was performed in SPM8 within a striatal search volume (2934 voxels, \~10 cm^3^), determined as the intersection of individual striatal subvolumes warped in MNI space. Because the images could not be assumed to be normally distributed, non-parametric analyses were conducted using SnPM (version SnPM12, <http://www2.warwick.ac.uk/fac/sci/statistics/staff/academic-research/nichols/software/snpm>.

A 2 (group) × 2 (time) ANOVA was conducted on BP = BP~1~−BP~0~. Variance smoothing was used with an FWHM of 10 mm. Five thousand permutations were used to determine a *p* \< 0.05 threshold within the striatal search space. Paired *t*-tests with a *p* \< 0.01 threshold were used to assess main effects of task (letter memory and n-back) on BP, combining both groups.

Off-line transfer tasks {#Sec11}
-----------------------

The pre- and post-training neuropsychological assessments included three WAIS-III subtests (digit span, letter-number sequencing, and digit symbol), a number-letter task, a Simon task, a visuospatial n-back task, and an episodic recall task. A pattern-comparison task, a number-copying task, Trail Making parts A and B, and the WAIS-III vocabulary subtest were only included before training. The pre-training assessment also included a background questionnaire, the Edinburgh Handedness Inventory, and the BDI-II.

### Digit span {#Sec12}

The digit span test was used to assess auditory attention as well as passive (forward span) and active (backward span) WM. Participants were asked to recall sequences of digits in the same or reversed order in which they were presented. Test administration, items, and scoring followed standard WAIS-III procedures^[@CR43]^. Maximum forward and backward spans were used as dependent measures.

### Letter-number sequencing {#Sec13}

The letter-number sequencing test was also used to assess auditory attention as well as passive and active working memory. Participants were instructed to recall and sort number-letter sequences. Test administration, items, and scoring followed standard WAIS-III procedures^[@CR43]^. The total score was used as the dependent measure.

### Visuospatial n-back {#Sec14}

The visuospatial n-back task was used to assess visuospatial WM updating^[@CR44]^. In this task, participants were asked to determine if a currently visible box was in the same location as the previous box (1-back) or the box that was presented two steps back (2-back). A matrix with eight stimulus locations was used (a 3 × 3 matrix without the middle square). Accuracy was used as the dependent measure.

### Digit-symbol coding {#Sec15}

The digit-symbol coding test was used to assess perceptual speed. Participants were asked to draw, as quickly as possible, different symbols in empty spaces on the basis of digits that were located above each empty space. A time limit of 90 sec was used here, but in all other respects test administration, items, and scoring followed standard WAIS-III procedures^[@CR43]^. The total score was used as the dependent measure.

### Number-letter task {#Sec16}

The number-letter task was used to assess set shifting. Here, participants were asked to categorize either the number or the letter in a number-letter pair depending on in which of two vertically aligned boxes the number-letter pair appeared. The task included both switch and no-switch trials, from which mixing and switching costs were calculated^[@CR44]^. The switching cost reflects the temporary cognitive load that is related to a task shift, whereas the mixing cost reflects the cost of maintaining attentional control in a situation where two task sets are active. The mixing and switching costs for RTs and accuracy were used as dependent measures.

### Simon task {#Sec17}

The Simon task was used to assess inhibitory control. In this task, participants were to determine the color of a square while disregarding where the square appeared^[@CR45]^. On congruent trials, the color and location of the square matched the positioning of the appropriate response key, whereas the stimulus and correct response key were crossed on incongruent trials. The Simon effect was calculated by comparing the performance on congruent and incongruent trials. Simon effects for accuracy and RT were used as dependent measures.

### Episodic word recall {#Sec18}

The word recall task was used to assess verbal episodic memory^[@CR46]^. Participants were instructed to read and remember 18 concrete nouns presented on a computer screen. They were instructed to recall as many words as possible, and the test leader recited any words that had been omitted. After hearing the omitted words, participants again attempted to recall as many words as possible, and once again the test leader recited any omissions. Finally, participants were again asked to recall as many words as possible. The total number of correctly recalled words was used as the dependent measure.

### Pattern comparison {#Sec19}

The pattern-comparison task was used to assess motor speed^[@CR47]^. In this task, participants were to determine, as quickly as possible, if two symbols were identical. The total number of correctly evaluated symbol-pairs was used as the dependent measure.

### Number copying {#Sec20}

The number-copying task was also used to assess motor speed^[@CR47]^. Here, participants were asked to copy, as quickly as possible, a set of numbers. The total number of correctly copied numbers was used as the dependent measure.

### Trail Making {#Sec21}

This task was used to assess visual search, hand-eye coordination, and set shifting^[@CR48]^. In part A, participants were to connect, as quickly as possible, a set of numbers in ascending order. In part B, they were to connect numbers and letters in both ascending and alphabetical order, constantly switching between the two types of stimuli (1-A-2-B, etc.). The dependent variable was the difference score obtained by subtracting the completion time of part A from part B.

### Vocabulary {#Sec22}

The vocabulary test was used to assess word comprehension. Participants were asked to explain the meaning of words of increasing complexity. Test administration, items, and scoring followed standard WAIS-III procedures^[@CR43]^. The total raw score was used as the dependent measure.

Statistical analyses {#Sec23}
--------------------

Those participants who were extreme outliers (i.e., showing a deviation of at least three times the interquartile range) before training were removed from the specific statistical analyses on which they deviated. Training effects were analyzed by mixed-model ANOVAs.

**Publisher\'s note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

The study was financially supported by grants from the Academy of Finland, the Åbo Akademi University endowment (the BrainTrain project), and the Turku University Hospital. LB was supported by an Alexander von Humboldt Research Award and a donation from the af Jochnick Foundation. LN was supported by a Scholar grant from the Knut and Alice Wallenberg Foundation.

A.S., K.A., and O.W. contributed to study design, acquisition of data, analyzed the neuropsychological data, interpretation of data, and revising the article. J.J. and M.A. analyzed and interpreted the PET data, and revised the article. L.B., M.L., J.O.R. and L.N. all contributed to study design, interpretation of data, and drafting of the article. All authors have given their approval of the final version.

Competing Interests {#FPar1}
===================

The authors declare that they have no competing interests.
