# Effects of Simultaneously Performed Dual-Task Training with Aerobic Exercise and Working Memory Training on Cognitive Functions and Neural Systems in the Elderly

Authors: Hikaru Takeuchi; Daniele Magistro; Yuka Kotozaki; Kosuke Motoki; Keyvan Kashkouli Nejad; Rui Nouchi; Hyeonjeong Jeong; Chiho Sato; Salvatore Sessa; Ryoichi Nagatomi; Massimiliano Zecca; Atsuo Takanishi; Ryuta Kawashima
Journal: Neural Plasticity
 (2020)

---

## Abstract

Working memory (WM) training (WMT), aerobic exercise training (AET), and dual-task training improve cognitive functions and alter neural systems in older adults. In particular, the effects have been investigated of dual-task training that combines a walking or standing activity (balance exercise) simultaneously performed with cognitive training (which is ecologically difficult for the elderly). In this study, we investigated the effects of simultaneously performed dual-task training incorporating both AET and WMT (SDAEWMT), using a recumbent ergocycle bicycle and a WMT program that provided a portable console and made the training ecologically easy for the elderly. Older adults (65.9±13.7 years old) participated in 3 months of SDAEWMT, WMT, or AET after random allocation, and the effects of SDAEWMT were compared with those of WMT and AET. Prior to and after training, all the subjects underwent cognitive testing, magnetic resonance imaging (MRI) involving diffusion tensor imaging (DTI), and functional MRI during performance of an N-back WM task. SDAEWMT improved executive function (performance of a frontal assessment battery); however, there was no evidence of broader transfer effects or enhanced learning with WMT. SDAEWMT resulted in mean diffusivity changes in brain areas involving the dopaminergic system, suggesting that neural tissue changes occurred in these areas. SDAEWMT also resulted in an increase in brain activity during the 2-back working memory task in brain areas involved in attentional reorientation. These results suggest that SDAEWMT is effective for improving cognitive functions and inducing beneficial neural changes in older adults.


---

## Body

## 1. Introduction

The proportion of older individuals is increasing worldwide in parallel with improvements in the average life expectancy [1]. Despite an increase in life expectancy, various studies have highlighted that aging is associated with deterioration in both physical and cognitive functions [2]. With regard to the cognitive aspects of aging, previous studies have shown changes in brain structure and function with advancing age [3], particularly in the prefrontal areas involved in executive function and the attention system [4] as well as in areas involved in working memory [5]. Therefore, it is important to identify and implement effective strategies and activities that can maintain and even improve the cognitive abilities of older adults. Cognitive training and physical training have been primarily investigated for this purpose. The combination of these approaches is the topic of the present study.Working memory (WM) is the limited-capacity storage system involved in the maintenance and manipulation of information over a short time period [6]. Regions in the external attention system (EAS) that play a key role in external attention [7–9], such as the lateral prefrontal cortex (LPFC) and the posterior parietal cortex, are involved with WM and are activated during WM performance [6]. Dopaminergic transmission plays an important role in WM. Cortical dopamine release has been observed during WM [10], and dopaminergic neurotransmission and its dose are critical in WM performance and in the tuning of prefrontal activity during WM. A number of studies on WM training (WMT) have been conducted, and results are inconsistent; however, meta-analyses have shown that WMT improves the performance of untrained WM tasks as well as attention/inhibition [11] and nonverbal reasoning [12, 13], and the observed inconsistency could be due to methodological differences, as well as to statistical fluctuations and lack of statistical power in individual studies that were then combined with weak effects [12, 13]. Such training also changes the brain activity patterns observed during the performance of WM-related tasks and alters gray matter and white matter structures, particularly those within lateral frontal and lateral parietal areas, which are activated during WM tasks [14, 15]. Furthermore, changes in the dopamine D1 receptor density in the lateral frontal and lateral parietal areas are seen after WMT [16]. WMT that involves updating was shown to reduce transient changes in the binding potentials of D2 receptors in the striatum during updating, which in turn was suggested to reflect increases in striatal dopamine release during the updating [17].However, the extent of cognitive function plasticity subject to cognitive training is characterized by a low effect size [11, 12]. Furthermore, the potential magnitude of these changes in cognitive function and neural changes resulting from cognitive training deteriorates with increasing age [14].Aerobic exercise training (AET) is another approach for countering age-related cognitive decline. Numerous studies have been conducted on AET, and results are inconsistent; however, meta-analyses have shown that AET improves various cognitive functions in the elderly, including executive functions, processing speed, memory, and selective attention [18]. However, the effects of AET on WM performance were not clear in one meta-analysis [19], while another meta-analysis reported that intervention with various physical activities improved the WM performance using different WM outcomes [20]; thus, the findings were divided. Also, a recent study with the elderly investigated the effects of cognitive training involving several cognitive tasks, and physical training involving multiple physical activities including cardiovascular, and AET showed effects of cognitive training in one of the WM outcome measures, whereas physical training failed to show such effects [21]. Further, AET has been shown to increase prefrontal and temporal grey matter volume, anterior white matter volume, and neural activity of lateral frontal and parietal areas [22]. AET has also been reported to increase the efficiency of learning [22]. Part of this exercise effect may be mediated through increased secretion of brain-derived neurotrophic factor (BDNF), a key signal promoting activity-dependent neural plasticity [22].The term dual-task activity refers to the performance of two different tasks simultaneously or in close succession. Dual-task activity is one of the functions of the central executive function system [23]. Also, simultaneously performed dual-task training that involves exercise has been shown to enhance cognitive executive functions, including general cognitive functions [24]. Dual-task training that involves overlapping tasks, usually two-alternative forced choice tasks that have to be performed simultaneously or in close succession, alters brain activity during cognitive tasks [25]. In particular, the effects of simultaneously performed dual-task training that combines a walking or standing activity (balance exercise) with cognitive training have been investigated [24]. However, the methods of these studies, including outcome measures, have varied. some studies recruited healthy older adults, whereas others recruited cognitively impaired adults [24]. Many of these studies, but not all, have found that this type of simultaneously performed dual-task training that involves physical activity leads to improved performance of global function measures such as the frontal assessment battery (FAB) [24]. Another recent meta-analysis also concluded that simultaneously combined physical and cognitive training is more successful than the use of physical training alone or cognitive training alone, though the improvement was somewhat confined to trained cognitive functions, rather than generalizing to other cognitive or daily living skills [26].To summarize, we focused on SDAEWMT (combination of WMT, AET, and simultaneously performed dual-task training) because of the following reasons: (1) the development of an effective strategy to maintain cognitive abilities in older adults is important; (2) WMT, AET, and simultaneously performed dual-task training have all been shown to be effective for such purposes; and (3) AET leads to an increased secretion of BDNF, which facilitates neural plasticity. However, for older adults, this type of training (e.g., execution of walking and WM tasks simultaneously in real life) may be ecologically difficult because of the requirement of supervised exercise in combination with the implementation of effective computerized cognitive training tasks similar to previously developed systems [27]. In addition, previous studies have not investigated the effects of simultaneously performed dual-task training involving both AET and WMT (SDAEWMT) on neural systems.To address these unresolved issues, we first designed AET incorporating the use of a recumbent ergocycle bike (Figure1), which allows older subjects to sit in a stable position during training, thereby facilitating training that is both safe and ecologically applicable. Next, we designed a WMT program that used a portable console (3DS; Nintendo, Kyoto, Japan). Incorporating portable hardware (3DS) facilitated easy subject use during SDAEWMT. The combination of AET and WMT was adaptable to individual capacity and the physical condition of each older subject, enabling them to reach their optimal ecological level (for details, see Methods). We hypothesized that SDAEWMT is more effective than WMT or AET alone and/or (a) has effects on cognitive functions such as global cognitive functions, memory performance, attention, and executive functions, for which previous simultaneously performed dual-task training involving aerobic exercise or both WMT and AET have been effective, as described above; (b) alters brain structures and functional activities involving frontal and parietal areas where effects of WMT and AET have been observed, as described above, and that are critical for dual-task execution [25]; and (c) facilitates the improvement of WMT performance (greater improvement of training tasks’ performance of WMT). The aims of the present study were to test these hypotheses and to investigate the effects of SDAEWMT on older adult cognitive functions, neural systems, and learning processes during cognitive training.Figure 1
The SDAEWMT regimen conducted in the present study. Subjects were instructed to perform WMT implemented on a portable-type console (3DS) while riding a recumbent bike.


## 2. Methods

### 2.1. Subjects

For subject recruitment, study inclusion factors were the ability to walk without assistance, physician approval after assessment, right-handedness as determined using the Edinburgh Handedness Inventory (Oldfield, 1971), and age more than 65 years and less than 75 years. The right-handed subjects were recruited based on the custom of the MRI studies. Subjects were recruited from Sendai City. There were several exclusion criteria in this study and for details, see Supplemental Methods.The Ethics Committee of Tohoku University approved all study procedures.After an initial check of the inclusion and exclusion criteria, 102 right-handed older adults [65.9±13.7 (SD) years old] participated in a physician-administered health assessment involving an exercise electrocardiogram. After further exclusions and dropouts before and after random allocation based on the random number generated through the excel file, a total of 93 subjects successfully completed the intervention (see Figure 2 for a detailed subject flowchart), with 30 subjects in the WMT group, 30 in the SDAEWMT group, and 33 in the AET group. Random allocation was performed using the Microsoft Excel program. Group age and sex ratios are presented in Table 1. No participant had a history of neurological or psychiatric illness according to a routine questionnaire that was created and is used in this laboratory. Handedness was evaluated using the Edinburgh Handedness Inventory [28]. Written informed consent was obtained from each individual for the interventions in which he/she participated. The three groups did not differ significantly in age, sex ratio, MMSE score, or FAB score [P>0.1 by analysis of variance (ANOVA)]. Subjects who misunderstood the rules for certain psychological measures were excluded from analyses involving those measures. Misunderstanding of the rules was judged according to “answers based on wrong rules” or “answers of chance level in the tasks with multiple choices”.Figure 2
Study flowchart.Table 1
Characteristics of subjects who successfully completed the assigned training sessions.
WMTSDAEWMTAET(N=30)(N=30)(N=33)Age,mean±SD68.77±2.9468.03±3.2769.3±3.37Female,n191821MMSE-J,mean±SD28.27±1.3928.53±1.8628.00±1.81VO2max30.97±5.6430.53±6.7329.48±6.46Completed training sessions,mean±SD (range)33.33±1.56(31-37)33.67±1.60(31-37)33.48±2.18(27-37)Education length∗4.63±1.054.90±1.164.73±1.15∗For details of questions examining the education length of each participant, see our previous study’s description of the questions that established the data on length of education of the participants’ parents [1]. The same method was used.


### 2.2. Procedure

All the subjects were asked to engage in training three times each week in a laboratory located at the Institute of Development, Aging and Cancer, Tohoku University. If a subject was not available during a given week, he/she could make up for the lost sessions in the following weeks. Five subjects in the WMT group, 9 in the SDAEWMT group, and 8 in the AET group performed this procedure and attended four sessions in a week at least once. The first WMT and SDAEWMT group sessions were practice sessions to introduce the game software that would be utilized during training (see the next subsection for details). Each training session lasted for approximately 1 h, and sessions were held over a period of 12 weeks. All the subjects underwent MRI and completed psychological evaluations immediately before and after this 12-week training period.


### 2.3. Training Tasks

For WMT, we used the game suite “Brain Age: Concentration Training” published by Nintendo (Nintendo, Kyoto, Japan). This game includes eight WMT tasks, of which we used seven. Devilish reading, a reading span-like task, was excluded because it required the subject to read aloud, which may have disturbed other subjects playing in the same room. One session on an individual task ended when one trial ended after 5 min of training time (excluding interval times) had elapsed, after which the subjects proceeded to the next task. The order of the tasks was not specified, and the subjects could perform each task in an arbitrary order. For details of the training tasks, see Supplemental Methods. The schema of the tasks was presented in Figure3.Figure 3
Schemata of four complex WMT tasks.The exact time taken for the subjects to complete the training session on each day was not measured. However, it takes at least 5 min to complete each task, and there are intervals between each trial during which the subjects could rest. Therefore, it usually took more than 45 min to complete the training session.For AET, we used a recumbent ergocycle (DK-8718RP, Daikou, Japan) that electronically controls the workload. Using this ergocycle, the subjects could sit leaning back in a stable position during training. The training consisted of three sessions per week of 60 min each that included (i) 5 min of activity at low intensity so that all the subjects could achieve 35%–40% of their training intensity (see below for this calculation), (ii) 45 min of aerobic exercise on the ergocycle at moderate intensity based on the individual’s capabilities and clinical status, and (iii) 10 min for cooling down and resting. Each subject’s heart rate was monitored during all the training sessions by means of a wearable system which is attached to the left hand and measures the heart rate as presented in the index finger. The target was between 40% and 50% of training intensity by Karovonen’s formula [29] using the individual maximal heart rate estimated using Tanaka’s formula [30]. The range was 40% to 50% training intensities (TI) which were calculated as follows.
(1)
Estimation of one’s own maximal heart rate (MHR).MHR=208–.70×ageinyearsbpm [30](2)
Measurement of resting Heart Rate (RHR)(3)
Calculation of Heart Rate Reserve (HRR).HRR=MHR–RHR [29](1)
Calculation of Heart rate.40%TI=HRR×0.4+RHR. 50%TI=HPR×0.5+RHR [29]40% to 50% TI may seem weak. However, we used 40% to 50% TI in this study for the following reasons. First, prior to the experiment, we confirmed it was difficult to combine aerobic exercise of stronger TI with WMT in the setting of SDAEWMT, and this study is the study to investigate the effects of SDAEWMT. Second, as described above, previous simultaneously performed dual-task training paradigms involving exercise used walking (which is supposed to have weak TI) as an exercise to combine with cognitive training [24]. Third, a meta-analysis using data from 29 studies, involving 2049 participants, showed that there are no significant correlations between the TI of aerobic exercise on cognitive measures that showed the effects of AET (measures of attention and processing speed, executive function, and memory) and the effect size [19]; however, another meta-analysis that investigated the association between the status of aerobic exercise and cognitive functions in cross-sectional studies showed that in cross-sectional studies, either an intensity of aerobic exercise that is too low or an intensity that is too high is not effective for cognitive functions [31]. Another meta-analysis, using the data from randomized controlled trial or quasi-experimental studies with appropriate comparison groups, targeted participants aged 65 years or above and who had been prescribed with low-intensity exercise in at least one study arm and reported that low-intensity exercise can improve both cognitive functions and affective state [32]. These previous studies demonstrate that our results cannot be dismissed due to the fact that our study used weaker TI [33]. The group difference of the preintervention aerobic capacity cannot explain the present findings, because this is a randomized controlled study and such differences should not occur and preintervention VO2max level did not differ among 3 groups (ANOVA, P=0.630, Table 1).In SDAEWMT, the subjects were asked to perform the same cognitive training tasks as those described for WMT while performing the same aerobic exercises as those described for AET. If they did not finish WMT within 55 min of exercise, they were asked to keep cooling down while simultaneously completing WMT.The length of one session of AET and SDAEWMT was determined from how long it takes to complete all the WMT tasks (about 45 min) in the WMT group. Moreover, at least one meta-analysis of aerobic exercise showed that a training length of 30–45 min/session is more effective than training with longer sessions and training with shorter sessions [34].Furthermore, after the experiment, a questionnaire with a 9-point Likert scale was administered to some participants to subjectively evaluate each training task (WMT:N=21, SDAEWMT: N=23, AET: N=24). There was no significant difference among the three groups in motivation toward the training tasks [P=0.746, F=0.294, WMT: 7.381±0.326 (SEM), SDAEWMT: 7.304±0.312 (SEM), AET: 7.625±0.305 (SEM)], fatigue of the subjects after the training [P=0.618, F=0.486, WMT: 5.190±0.460 (SEM), SDAEWMT: 5.435±0.439 (SEM), AET: 4.833±0.430 (SEM)], satisfaction with the training tasks [P=0.792, F=0.234, WMT: 6.476±0.444 (SEM), SDAEWMT: 6.261±0.424 (SEM), AET: 6.667±0.415 (SEM)], or enjoyment during the training tasks [P=0.402, F=0.925, WMT: 7.143±0.348 (SEM), SDAEWMT: 6.870±0.332 (SEM), AET: 7.500±0.325 (SEM)]. Therefore, these group differences do not seem to be substantial.


### 2.4. Psychological Outcome Measures

The following neuropsychological measures assessing specific cognitive domains were administered as pre- and posttraining evaluations. They do not include computerized tasks. All measures were administered to test the primary hypotheses of the present study (superiority of SDAEWMT over WMT and AET).The details of tests (a), (b), (d), (f), (g), and (h) are provided in Supplemental Methods.(a) Nonverbal reasoning. Nonverbal reasoning was assessed using Raven’s advanced progressive matrix [35].(b) Verbal WM. Verbal WM was tested by a digit span task.(c) Immediate recall. Logical memory I (LM), which involves immediate recall after hearing a story, was evaluated using a subtest of the Wechsler Memory Scale-Revised (WMS-R) [36]. The test material consisted of two short paragraph-length stories (story A and story B). In the LM task, the subjects were required to memorize short story A, and soon after that, subjects were told to tell out loud the story they remembered. Then, the same procedure was repeated for story B. Performance was scored by the number of story units recalled as specified in the WMS-R scoring protocol. The primary measure for this task was the sum of the number of correct story units.(d) Attention. Attention was measured by the digit cancellation task (D-CAT) [37]. For this measure, we used the mean number of successfully cancelled items from the three trials.(e) Frontal lobe and executive function. FAB was employed for this evaluation [38]. The following description of FAB is a summary from our previous study [39]. FAB consists of six subtests, namely, those for similarities (conceptualization), lexical fluency (mental flexibility), motor series (programming), conflicting instructions (sensitivity to interference), go–no go (inhibitory control), and prehension behavior (environmental autonomy). FAB is scored from 0 to 18. Lower scores of FAB indicate greater degrees of executive dysfunction. The primary measure is the total score of this task (max=18). The full details of these subtests have been described extensively in a previous study [38].(f) Processing speed. Processing speed was evaluated by the Symbol Search [40] and Digit Symbol tasks [40].(g) Creativity measured by divergent thinking. The SA creativity test was employed to evaluate creativity [41–43]. This task involves three subtasks: (i) generate unique ways of using typical objects, (ii) imagine desirable functions for ordinary objects, and (iii) imagine the consequences of unimaginable things happening. The SA test scores the four dimensions of the creative process (fluency, originality, elaboration, and flexibility). In this study, the sum of the graded scores of the four dimensions was used for analysis.(h) General cognitive functions. MMSE [44], the screening instrument most widely used for detecting cognitive impairment in older adults, was used to evaluate general cognition. Because of an apparent ceiling effect, this task was not used for the assessment of training effects.(i) Inhibition. The Stroop task (Hakoda’s version) [45] was used to evaluate response inhibition and impulsivity. However, results were not analyzed because a substantial proportion of the subjects could not see yellow. Subjects were elderly, and with aging, sometimes the lenses of their eyes had become more yellow, making it difficult for such subjects to differentiate yellow from white (which is the background color of the test). In addition, we could not provide instructions to testers on how to address this problem; therefore, different testers provided different instructions to the subjects.Because of misunderstanding of the rules, one subject in the WMT group was excluded from the analysis of the LM task. Two subjects in the WMT group were excluded from the analysis of the Symbol Search task. One subject in the WMT group, one subject in the SDAEWMT group, and three subjects in the AET group were excluded from the analysis of D-CAT.Several questionnaires designed to assess the traits or states of the subjects were collected but are not described here. Physical tasks involving walking and dual tasks were also performed, but these results will be reported by other authors.The preassessment of psychological outcome measures was conducted in a double-blind manner. We used exactly the same versions of the test for the pre- and postmeasures, except for the S-A creativity test, which had different versions for pre and post [41].


### 2.5. Functional Magnetic Resonance Imaging (fMRI) Tasks

fMRI was used to map training-induced changes in brain activity related to WM. AnN-back task commonly used for fMRI studies was used to detect training-related changes in brain activity patterns related to these cognitive processes. In this version using digits as stimuli, there were two conditions: 0-back for testing simple cognitive processes and 2-back for testing WM. These are described in our previous publications [46, 47]. In brief, the subjects performed an N-back task using four numbers (i.e., 1–4). The subjects were asked to memorize a series of numbers and their temporal order, update the list of recent items, and select the responses that corresponded to the previously observed stimuli according to the N-back rule. The stimuli were shown in a random sequence on the screen. Each stimulus was presented for 0.5 s. Following this, a fixation symbol was presented for 2.5 s, and the next stimulus was then presented. The subjects had to respond either directly following the stimulus (0-back) or after a delay of two stimulus presentations (2-back). The schema of the task is presented in Figure 4. The task level of the memory load was displayed above the stimuli starting 2 s before task onset and remained visible during the task period. A single-condition block consisted of a 2-s presentation of the task instructions (e.g., “2-back”), followed by the 18-s actual task period. There were five blocks for each period. Because there are various strategies for performing N-back tasks, the subjects were requested to follow a specific strategy as instructed. The subjects were told to update memory N by N (e.g., in the 2-back task, they were asked to update their memory two by two) and not one by one. For instance, in the 2-back task, they first had to remember two numbers presented, and then, while they pushed the corresponding two buttons, they were asked to remember the two incoming numbers.Figure 4
Illustration of the N-back task, showing six trials of the 0-back level (left) and six trials of the 2-back level (right), together with the appropriate responses. At the 0-back level, subjects only have to push a button spatially corresponding to the location of a number on the screen. At the 2-back level, subjects must push a button spatially corresponding to a number that appeared two trials earlier.Visual stimuli were presented using Presentation Software (Neurobehavioral Systems, Inc., Albany, CA, USA). Stimuli were projected to the subjects via a screen positioned at the head-end of the board. The subjects viewed the screen through a mirror attached to the head coil. A fiber-optic light-sensitive key press button box was used to record the subjects’ behavior. In contrast to other variants of theN-back task, only accuracy was used as a measure of performance in this study. The reaction time did not reflect performance because in the 2-back task, prior to the presentation of the next stimuli, subjects already knew which button they should push (as determined by the preceding stimuli) after the presentation of the following stimuli and the timing of the presentation of the next stimuli before they were allowed to respond (see Figure 4).In this study, our focus was on the frontal and parietal areas, and activation in these regions occurs regardless whether the task is 2-back or 0-back in mostly similar areas, although the magnitude may differ. Furthermore, differences in brain activity between patients with schizophrenia and control subjects were mostly similar, regardless of whether the task was a 0-back task or 2-back task, especially in these areas (i.e., subtracting the activity during the 0-back task from the brain activity during the 2-back task substantially eliminates the group differences) [48, 49].


### 2.6. Image Acquisition and Analysis

MRI data were acquired using a 3 T Philips Intera Achieva scanner. Diffusion-weighted data were acquired using a spin-echo EPI sequence [repetitiontimeTR=10293ms, echotimeTE=55ms, fieldofviewFOV=22.4cm, 2×2×2mm3 voxels, 60 slices, SENSEreductionfactor=2, numberofacquisitions=1]. The diffusion weighting was isotropically distributed along 32 directions (bvalue=1,000s/mm2). In addition, three images with no diffusion weighting (bvalue=0s/mm2) (b=0 images) were acquired using a spin-echo EPI sequence (TR=10293ms, TE=55ms, FOV=22.4cm, 2×2×2mm3 voxels, 60 slices). From the collected images, fractional anisotropy (FA) and mean diffusivity (MD) maps were calculated. Forty-two transaxial gradient-echo images (TE=30ms, flipangle=90°, slicethickness=3mm, FOV=192mm, matrix=64×64) covering the entire brain were acquired at TR of 2.5 s using an echo planar sequence. For the N-back session, 135 functional volumes were obtained.


### 2.7. Preprocessing of Diffusion Data

The MD and FA measures of diffusion tensor imaging (DTI) [50] can detect unique aspects of brain microstructural neural plasticity [15, 51], including those caused by WMT [15, 52, 53]. A reduction in MD reflects an increase in the tissue density of various cellular structures, including capillaries, synapses, spines, and macromolecular proteins, as well as changes in the properties of myelin, membranes, and axons; the shapes of neurons or glia; and tissue organization. However, MD is not specific to any one of these microlevel changes [50, 51]. On the other hand, an increase in FA is thought to have a relatively greater association with changes in microstructural properties related to brain connectivity, including myelination, axonal membrane thickness and diameter, and degree of axonal parallel organization [50, 52].Preprocessing of imaging data was performed using Parametric Mapping software (SPM8; Wellcome Department of Cognitive Neurology, London, UK) implemented in Matlab (Mathworks Inc., Natick, MA, USA). The descriptions in this subsection are a summary from our previous publication [15]. Pre- and posttraining MD and FA maps were independently segmented and normalized using the previously described diffeomorphic anatomical registration through an exponentiated Lie algebra (DARTEL)-based method (Takeuchi et al., 2013a). This method accounts not only for the outer edge but also for the signal distribution within the white matter during registration. The reason for the independent normalization has been described previously (Takeuchi et al., 2013a).From the pre- and postintervention normalized images of the normalized MD map, areas not strongly likely to be gray or white matter in our custom template (defined by “gray matter tissue probability + white matter tissue probability < 0.99”) were removed. Then, from the pre- and postintervention normalized images of normalized FA map, areas not strongly likely to be white matter in our custom template (defined by “white matter tissue probability < 0.99”) were removed. Subsequently, normalized MD and FA images were smoothed by convolution with an isotropic Gaussian kernel of 10- and 8-mm full-width at half maximum, respectively. These images were used for the group analysis described below.


### 2.8. Preprocessing and Data Analysis for Functional Activation Data

Preprocessing and data analysis were performed using SPM8 implemented in Matlab. Prior to analysis, BOLD images from the pretraining scan and the posttraining scan were realigned and resliced to the mean BOLD image of all pretraining scans. The images were then corrected for slice timing, coregistered to ab=0 diffusion-weighted image, and spatially normalized using the parameter of the abovementioned normalization of FA images to form images with 3×3×3mm3 voxels.A design matrix was fitted to each subject with one regressor for each task condition (0-back and 2-back in theN-back task) using the standard hemodynamic response function (HRF). The design matrix weighted each raw image according to its overall variability to reduce the impact of movement artifacts [54]. After estimation, beta images were smoothed (8-mm full-width at half maximum) and either taken to the second level or subjected to random-effect analysis. We removed low-frequency fluctuations using a high-pass filter with a cut-off value of 128 s. Individual-level statistical analyses were performed using a general linear model.In individual analyses, we focused on the activation change and compared activation relative to the conditions (0-/2-back versus rest) before and after an intervention period. The resulting maps for each subject represented changes between the pre- and postmeasures of brain activity during the corresponding conditions as well as the brain activity during the corresponding conditions in the premeasures, which were then included in a group analysis.


### 2.9. Preprocessing of Structural Data

Preprocessing of imaging data was performed using SPM12 implemented in Matlab. Using the Pairwise Longitudinal Registration toolbox [55], the midpoint average template images were generated for each participant as well as maps of differences between the Jacobian determinants and of pre- and postintervention T1-weighted structural images using default parameter settings. The midpoint average image of each subject was then segmented, and the Thorough Clean option was used to remove any odd voxels. Affine regularization was performed in accordance with the East Asian brain with a sampling distance (the approximate distance between the sampled points when estimating the model parameters) of 1 mm. Default parameter settings were used for other parameters in the segmentation process. Then, by multiplying maps of differences between the Jacobian determinants of pre- and postintervention images and gray matter segments of the midpoint average image of each subject, images of gray matter volume change were generated. Then, diffeomorphic anatomical registration was performed using DARTEL processes. Here, we used imported DARTEL images of gray matter and white matter for the processes created above. The DARTEL template was created from the entire set of subjects in this experiment, and the default parameter settings were used for the DARTEL process. Then, using the parameters generated for each participant, the images of gray matter change were normalized and smoothed by convolving each with an isotropic Gaussian kernel of 8 mm FWHM.


### 2.10. Statistical Group-Level Analysis of Imaging and Behavioral Data

There were numerous potential contrasts in this study, because there were three groups with two factors (aerobic exercise and WMT). However, we only tested two-step analyses to confirm our hypothesis. The first step was the superiority of SDAEWMT over single training (i.e., comparison of SDAEWMT with the combined groups of WMT and AET) which was followed by the second step that confirmed the superiority of SDAEWMT over AET as well as the superiority of SDAEWMT over WMT (the comparison between SDAEWMT and AET and the comparison between SDAEWMT and WMT). Behavioral data were analyzed using SPSS 22.0 (SPSS Inc., Chicago, IL). In behavioral analysis, because the superiority of SDAEWMT was the primary study hypothesis (and we are not aware of any study indicating that dual-task training is inferior to single-task training in terms of cognitive transfer), test–retest changes within the SDAEWMT group were compared with those in the WMT or AET group using one-tailed, one-way analysis of covariance (ANCOVA) with the posttest measures as dependent variables and pretest scores as independent variables (P<0.05). This ANCOVA design was selected to increase the sensitivity of analysis. As described in our previous study [56], we employed ANCOVAs instead of repeated measure ANOVAs to control the effects of pretest scores. Statistical experts strongly recommend to use ANCOVA instead of repeated measure ANOVA in this type of study design [57]. With randomized designs, the purpose of ANCOVA is to reduce error variance, whereas with nonrandomized designs (or with analyses involving substantial preexisting group differences), ANCOVA is used to adjust the posttest means for pretest differences among groups [57]. One might recommend that posttest scores be used, instead of using the differences between pre- and posttest measures. However, in practice, when the pretest scores are included as covariates, the two analyses return the same P value.The use of a one-tailed test was consistent with previous studies [27, 47, 58, 59]. In tests of the effects of SDAEWMT group compared with those of others, for post hoc comparisons between SDAEWMT and AET groups as well as comparisons between SDAEWMT and WMT groups, when ANCOVAs showed the same pattern (P<0.1) of superiority for the effects of SDAEWMT, the results were regarded as showing the superiority of SDAEWMT over single training (here the main contrast is SDAEWMT vs AET + WMT; however, the existence of similar trends in the contrasts of (SDAEWMT vs AET) and (SDAEWMT vs WMT) was confirmed with a more liberal threshold). In all nonwhole-brain analyses (i.e., behavioral analyses) performed in this study, results with a threshold of P<0.05 corrected for false discovery rate (FDR) using the graphically sharpened method [60] were considered statistically significant. The correction for multiple comparisons using this method was applied to the results of 10 ANCOVAs that tested the superiority of SDAEWMT over single training. ANOVAs were utilized to evaluate group differences in cognitive performance before the intervention. In group-level imaging analyses, we tested for group-wise differences in the change of each measure across the entire brain.For determining the MD, FA, and brain activity, analyses were performed using Biological Parametrical Mapping [61], the extension of SPM, which allow the inclusion of signal of other images at every voxel, to perform voxel-by-voxel ANCOVAs. In analyses of MD and FA, we performed voxel-wise ANCOVAs with the values of postscans at each voxel as dependent variables and the values for the prescans at each voxel as independent variables. In analysis of the N-back task, we performed voxel-wise ANCOVAs with the contrast of activity (0-back or 2-back) changes between pre- and postscans at each voxel as dependent variables and the corresponding value for the prescan at each voxel, the accuracy of the corresponding tasks in the prescan, and its pre- to postchanges (post–pre) as independent variables. For the N-back task, only the subjects who achieved an accuracy of 80% or higher during both pre- and postsessions were included in the analysis, because when the task difficulty was too high for the subjects, brain activity showed nonlinear patterns and such nonlinear changes could not be assessed by the statistical tests applied [62]. In other words, areas that showed an increase in activity as the load increased started to exhibit a decrease in activity when the task difficulty became too high. By this inclusion, we could focus only on linear changes and could ignore the effects of “not doing the tasks properly.” Consistently, Jansma et al. [48] showed that 75% accuracy is too low for this type of N-back task among healthy participants to detect linear changes (brain activity in relevant areas increased from 0-back to 2-back task in which accuracy was about 82%, but dropped substantially in the 3-back condition in which the accuracy was about 72.5% in the case of healthy participants). Since the random response can achieve 25% accuracy rate in this task, 75% accuracy is actually not high, and it means they are substantially failing. Consequently, 13 WMT group subjects (mean accuracy, pre: 90.0±6.8%, post: 91.9±7.2%), 12 SDWMAET group subjects (mean accuracy, pre: 90.8±8.1%, post: 94.2±5.3), and 13 AET group subjects (mean accuracy, pre: 89.6±6.6%, post: 91.5±6.0%) were included in analyses of the 2-back task. Without surprisingly, there were no significant differences of pre- to postchanges in accuracy among these subjects after correcting preintervention accuracy. 28 WMT group subjects, 26 SDWMAET group subjects, and 27 AET group subjects were included in analyses of the 0-back task. In behavioral analyses of the N-back task, only the subjects who performed the task without misunderstanding of the rules were included from analyses. In these neuroimaging analyses, the effects of SDAEWMT compared with those of single training were tested using the contrasts of −1×WMT2×SDAEWMT−1×AET or 1×WMT−2×SDAEWMT1×AET with the uncorrected masks of −1×WMT1×SDAEWMT or 1×WMT−1×SDAEWMT and 1×SDAEWMT−1×AET or −1×SDAEWMT1×AET (P<0.05, uncorrected). Regions with significance were inferred using cluster-level statistics [63]. Only clusters with P<.05, after correction for multiple comparisons at cluster size, with a voxel-level cluster-determining threshold of P<.001 (uncorrected) were considered statistically significant.In the group-level imaging analysis of rGMV, we tested for group-wise differences in the change in rGMV across the whole brain with the same contrasts. We used the one-way ANOVA option in SPM8. In the imaging analysis, the effects of training with TCSSL were estimated by comparing differences between the groups of rGMV, as calculated above at each voxel. This analysis was limited to areas where the average rGMV value for the segmented and normalized mean images of pre- and postscans of all participants was greater than 0.1. A multiple comparison correction was performed using threshold-free cluster enhancement (TFCE) [64] with randomized (5,000 permutations) nonparametric testing using the TFCE toolbox (http://dbm.neuro.uni-jena.de/tfce/). We applied a threshold of family-wise error corrected at P<0.05. The standard cluster size test is inappropriate in voxel-based morphometry [65], whereas a permutation-based test is assumed to be correct. However, the abovementioned analysis using BPM cannot be applied to rGMV analysis because the TFCE toolbox cannot be used in BPM, which is why we chose this statistical design for rGMV analysis.


## 2.1. Subjects

For subject recruitment, study inclusion factors were the ability to walk without assistance, physician approval after assessment, right-handedness as determined using the Edinburgh Handedness Inventory (Oldfield, 1971), and age more than 65 years and less than 75 years. The right-handed subjects were recruited based on the custom of the MRI studies. Subjects were recruited from Sendai City. There were several exclusion criteria in this study and for details, see Supplemental Methods.The Ethics Committee of Tohoku University approved all study procedures.After an initial check of the inclusion and exclusion criteria, 102 right-handed older adults [65.9±13.7 (SD) years old] participated in a physician-administered health assessment involving an exercise electrocardiogram. After further exclusions and dropouts before and after random allocation based on the random number generated through the excel file, a total of 93 subjects successfully completed the intervention (see Figure 2 for a detailed subject flowchart), with 30 subjects in the WMT group, 30 in the SDAEWMT group, and 33 in the AET group. Random allocation was performed using the Microsoft Excel program. Group age and sex ratios are presented in Table 1. No participant had a history of neurological or psychiatric illness according to a routine questionnaire that was created and is used in this laboratory. Handedness was evaluated using the Edinburgh Handedness Inventory [28]. Written informed consent was obtained from each individual for the interventions in which he/she participated. The three groups did not differ significantly in age, sex ratio, MMSE score, or FAB score [P>0.1 by analysis of variance (ANOVA)]. Subjects who misunderstood the rules for certain psychological measures were excluded from analyses involving those measures. Misunderstanding of the rules was judged according to “answers based on wrong rules” or “answers of chance level in the tasks with multiple choices”.Figure 2
Study flowchart.Table 1
Characteristics of subjects who successfully completed the assigned training sessions.
WMTSDAEWMTAET(N=30)(N=30)(N=33)Age,mean±SD68.77±2.9468.03±3.2769.3±3.37Female,n191821MMSE-J,mean±SD28.27±1.3928.53±1.8628.00±1.81VO2max30.97±5.6430.53±6.7329.48±6.46Completed training sessions,mean±SD (range)33.33±1.56(31-37)33.67±1.60(31-37)33.48±2.18(27-37)Education length∗4.63±1.054.90±1.164.73±1.15∗For details of questions examining the education length of each participant, see our previous study’s description of the questions that established the data on length of education of the participants’ parents [1]. The same method was used.


## 2.2. Procedure

All the subjects were asked to engage in training three times each week in a laboratory located at the Institute of Development, Aging and Cancer, Tohoku University. If a subject was not available during a given week, he/she could make up for the lost sessions in the following weeks. Five subjects in the WMT group, 9 in the SDAEWMT group, and 8 in the AET group performed this procedure and attended four sessions in a week at least once. The first WMT and SDAEWMT group sessions were practice sessions to introduce the game software that would be utilized during training (see the next subsection for details). Each training session lasted for approximately 1 h, and sessions were held over a period of 12 weeks. All the subjects underwent MRI and completed psychological evaluations immediately before and after this 12-week training period.


## 2.3. Training Tasks

For WMT, we used the game suite “Brain Age: Concentration Training” published by Nintendo (Nintendo, Kyoto, Japan). This game includes eight WMT tasks, of which we used seven. Devilish reading, a reading span-like task, was excluded because it required the subject to read aloud, which may have disturbed other subjects playing in the same room. One session on an individual task ended when one trial ended after 5 min of training time (excluding interval times) had elapsed, after which the subjects proceeded to the next task. The order of the tasks was not specified, and the subjects could perform each task in an arbitrary order. For details of the training tasks, see Supplemental Methods. The schema of the tasks was presented in Figure3.Figure 3
Schemata of four complex WMT tasks.The exact time taken for the subjects to complete the training session on each day was not measured. However, it takes at least 5 min to complete each task, and there are intervals between each trial during which the subjects could rest. Therefore, it usually took more than 45 min to complete the training session.For AET, we used a recumbent ergocycle (DK-8718RP, Daikou, Japan) that electronically controls the workload. Using this ergocycle, the subjects could sit leaning back in a stable position during training. The training consisted of three sessions per week of 60 min each that included (i) 5 min of activity at low intensity so that all the subjects could achieve 35%–40% of their training intensity (see below for this calculation), (ii) 45 min of aerobic exercise on the ergocycle at moderate intensity based on the individual’s capabilities and clinical status, and (iii) 10 min for cooling down and resting. Each subject’s heart rate was monitored during all the training sessions by means of a wearable system which is attached to the left hand and measures the heart rate as presented in the index finger. The target was between 40% and 50% of training intensity by Karovonen’s formula [29] using the individual maximal heart rate estimated using Tanaka’s formula [30]. The range was 40% to 50% training intensities (TI) which were calculated as follows.
(1)
Estimation of one’s own maximal heart rate (MHR).MHR=208–.70×ageinyearsbpm [30](2)
Measurement of resting Heart Rate (RHR)(3)
Calculation of Heart Rate Reserve (HRR).HRR=MHR–RHR [29](1)
Calculation of Heart rate.40%TI=HRR×0.4+RHR. 50%TI=HPR×0.5+RHR [29]40% to 50% TI may seem weak. However, we used 40% to 50% TI in this study for the following reasons. First, prior to the experiment, we confirmed it was difficult to combine aerobic exercise of stronger TI with WMT in the setting of SDAEWMT, and this study is the study to investigate the effects of SDAEWMT. Second, as described above, previous simultaneously performed dual-task training paradigms involving exercise used walking (which is supposed to have weak TI) as an exercise to combine with cognitive training [24]. Third, a meta-analysis using data from 29 studies, involving 2049 participants, showed that there are no significant correlations between the TI of aerobic exercise on cognitive measures that showed the effects of AET (measures of attention and processing speed, executive function, and memory) and the effect size [19]; however, another meta-analysis that investigated the association between the status of aerobic exercise and cognitive functions in cross-sectional studies showed that in cross-sectional studies, either an intensity of aerobic exercise that is too low or an intensity that is too high is not effective for cognitive functions [31]. Another meta-analysis, using the data from randomized controlled trial or quasi-experimental studies with appropriate comparison groups, targeted participants aged 65 years or above and who had been prescribed with low-intensity exercise in at least one study arm and reported that low-intensity exercise can improve both cognitive functions and affective state [32]. These previous studies demonstrate that our results cannot be dismissed due to the fact that our study used weaker TI [33]. The group difference of the preintervention aerobic capacity cannot explain the present findings, because this is a randomized controlled study and such differences should not occur and preintervention VO2max level did not differ among 3 groups (ANOVA, P=0.630, Table 1).In SDAEWMT, the subjects were asked to perform the same cognitive training tasks as those described for WMT while performing the same aerobic exercises as those described for AET. If they did not finish WMT within 55 min of exercise, they were asked to keep cooling down while simultaneously completing WMT.The length of one session of AET and SDAEWMT was determined from how long it takes to complete all the WMT tasks (about 45 min) in the WMT group. Moreover, at least one meta-analysis of aerobic exercise showed that a training length of 30–45 min/session is more effective than training with longer sessions and training with shorter sessions [34].Furthermore, after the experiment, a questionnaire with a 9-point Likert scale was administered to some participants to subjectively evaluate each training task (WMT:N=21, SDAEWMT: N=23, AET: N=24). There was no significant difference among the three groups in motivation toward the training tasks [P=0.746, F=0.294, WMT: 7.381±0.326 (SEM), SDAEWMT: 7.304±0.312 (SEM), AET: 7.625±0.305 (SEM)], fatigue of the subjects after the training [P=0.618, F=0.486, WMT: 5.190±0.460 (SEM), SDAEWMT: 5.435±0.439 (SEM), AET: 4.833±0.430 (SEM)], satisfaction with the training tasks [P=0.792, F=0.234, WMT: 6.476±0.444 (SEM), SDAEWMT: 6.261±0.424 (SEM), AET: 6.667±0.415 (SEM)], or enjoyment during the training tasks [P=0.402, F=0.925, WMT: 7.143±0.348 (SEM), SDAEWMT: 6.870±0.332 (SEM), AET: 7.500±0.325 (SEM)]. Therefore, these group differences do not seem to be substantial.


## 2.4. Psychological Outcome Measures

The following neuropsychological measures assessing specific cognitive domains were administered as pre- and posttraining evaluations. They do not include computerized tasks. All measures were administered to test the primary hypotheses of the present study (superiority of SDAEWMT over WMT and AET).The details of tests (a), (b), (d), (f), (g), and (h) are provided in Supplemental Methods.(a) Nonverbal reasoning. Nonverbal reasoning was assessed using Raven’s advanced progressive matrix [35].(b) Verbal WM. Verbal WM was tested by a digit span task.(c) Immediate recall. Logical memory I (LM), which involves immediate recall after hearing a story, was evaluated using a subtest of the Wechsler Memory Scale-Revised (WMS-R) [36]. The test material consisted of two short paragraph-length stories (story A and story B). In the LM task, the subjects were required to memorize short story A, and soon after that, subjects were told to tell out loud the story they remembered. Then, the same procedure was repeated for story B. Performance was scored by the number of story units recalled as specified in the WMS-R scoring protocol. The primary measure for this task was the sum of the number of correct story units.(d) Attention. Attention was measured by the digit cancellation task (D-CAT) [37]. For this measure, we used the mean number of successfully cancelled items from the three trials.(e) Frontal lobe and executive function. FAB was employed for this evaluation [38]. The following description of FAB is a summary from our previous study [39]. FAB consists of six subtests, namely, those for similarities (conceptualization), lexical fluency (mental flexibility), motor series (programming), conflicting instructions (sensitivity to interference), go–no go (inhibitory control), and prehension behavior (environmental autonomy). FAB is scored from 0 to 18. Lower scores of FAB indicate greater degrees of executive dysfunction. The primary measure is the total score of this task (max=18). The full details of these subtests have been described extensively in a previous study [38].(f) Processing speed. Processing speed was evaluated by the Symbol Search [40] and Digit Symbol tasks [40].(g) Creativity measured by divergent thinking. The SA creativity test was employed to evaluate creativity [41–43]. This task involves three subtasks: (i) generate unique ways of using typical objects, (ii) imagine desirable functions for ordinary objects, and (iii) imagine the consequences of unimaginable things happening. The SA test scores the four dimensions of the creative process (fluency, originality, elaboration, and flexibility). In this study, the sum of the graded scores of the four dimensions was used for analysis.(h) General cognitive functions. MMSE [44], the screening instrument most widely used for detecting cognitive impairment in older adults, was used to evaluate general cognition. Because of an apparent ceiling effect, this task was not used for the assessment of training effects.(i) Inhibition. The Stroop task (Hakoda’s version) [45] was used to evaluate response inhibition and impulsivity. However, results were not analyzed because a substantial proportion of the subjects could not see yellow. Subjects were elderly, and with aging, sometimes the lenses of their eyes had become more yellow, making it difficult for such subjects to differentiate yellow from white (which is the background color of the test). In addition, we could not provide instructions to testers on how to address this problem; therefore, different testers provided different instructions to the subjects.Because of misunderstanding of the rules, one subject in the WMT group was excluded from the analysis of the LM task. Two subjects in the WMT group were excluded from the analysis of the Symbol Search task. One subject in the WMT group, one subject in the SDAEWMT group, and three subjects in the AET group were excluded from the analysis of D-CAT.Several questionnaires designed to assess the traits or states of the subjects were collected but are not described here. Physical tasks involving walking and dual tasks were also performed, but these results will be reported by other authors.The preassessment of psychological outcome measures was conducted in a double-blind manner. We used exactly the same versions of the test for the pre- and postmeasures, except for the S-A creativity test, which had different versions for pre and post [41].


## 2.5. Functional Magnetic Resonance Imaging (fMRI) Tasks

fMRI was used to map training-induced changes in brain activity related to WM. AnN-back task commonly used for fMRI studies was used to detect training-related changes in brain activity patterns related to these cognitive processes. In this version using digits as stimuli, there were two conditions: 0-back for testing simple cognitive processes and 2-back for testing WM. These are described in our previous publications [46, 47]. In brief, the subjects performed an N-back task using four numbers (i.e., 1–4). The subjects were asked to memorize a series of numbers and their temporal order, update the list of recent items, and select the responses that corresponded to the previously observed stimuli according to the N-back rule. The stimuli were shown in a random sequence on the screen. Each stimulus was presented for 0.5 s. Following this, a fixation symbol was presented for 2.5 s, and the next stimulus was then presented. The subjects had to respond either directly following the stimulus (0-back) or after a delay of two stimulus presentations (2-back). The schema of the task is presented in Figure 4. The task level of the memory load was displayed above the stimuli starting 2 s before task onset and remained visible during the task period. A single-condition block consisted of a 2-s presentation of the task instructions (e.g., “2-back”), followed by the 18-s actual task period. There were five blocks for each period. Because there are various strategies for performing N-back tasks, the subjects were requested to follow a specific strategy as instructed. The subjects were told to update memory N by N (e.g., in the 2-back task, they were asked to update their memory two by two) and not one by one. For instance, in the 2-back task, they first had to remember two numbers presented, and then, while they pushed the corresponding two buttons, they were asked to remember the two incoming numbers.Figure 4
Illustration of the N-back task, showing six trials of the 0-back level (left) and six trials of the 2-back level (right), together with the appropriate responses. At the 0-back level, subjects only have to push a button spatially corresponding to the location of a number on the screen. At the 2-back level, subjects must push a button spatially corresponding to a number that appeared two trials earlier.Visual stimuli were presented using Presentation Software (Neurobehavioral Systems, Inc., Albany, CA, USA). Stimuli were projected to the subjects via a screen positioned at the head-end of the board. The subjects viewed the screen through a mirror attached to the head coil. A fiber-optic light-sensitive key press button box was used to record the subjects’ behavior. In contrast to other variants of theN-back task, only accuracy was used as a measure of performance in this study. The reaction time did not reflect performance because in the 2-back task, prior to the presentation of the next stimuli, subjects already knew which button they should push (as determined by the preceding stimuli) after the presentation of the following stimuli and the timing of the presentation of the next stimuli before they were allowed to respond (see Figure 4).In this study, our focus was on the frontal and parietal areas, and activation in these regions occurs regardless whether the task is 2-back or 0-back in mostly similar areas, although the magnitude may differ. Furthermore, differences in brain activity between patients with schizophrenia and control subjects were mostly similar, regardless of whether the task was a 0-back task or 2-back task, especially in these areas (i.e., subtracting the activity during the 0-back task from the brain activity during the 2-back task substantially eliminates the group differences) [48, 49].


## 2.6. Image Acquisition and Analysis

MRI data were acquired using a 3 T Philips Intera Achieva scanner. Diffusion-weighted data were acquired using a spin-echo EPI sequence [repetitiontimeTR=10293ms, echotimeTE=55ms, fieldofviewFOV=22.4cm, 2×2×2mm3 voxels, 60 slices, SENSEreductionfactor=2, numberofacquisitions=1]. The diffusion weighting was isotropically distributed along 32 directions (bvalue=1,000s/mm2). In addition, three images with no diffusion weighting (bvalue=0s/mm2) (b=0 images) were acquired using a spin-echo EPI sequence (TR=10293ms, TE=55ms, FOV=22.4cm, 2×2×2mm3 voxels, 60 slices). From the collected images, fractional anisotropy (FA) and mean diffusivity (MD) maps were calculated. Forty-two transaxial gradient-echo images (TE=30ms, flipangle=90°, slicethickness=3mm, FOV=192mm, matrix=64×64) covering the entire brain were acquired at TR of 2.5 s using an echo planar sequence. For the N-back session, 135 functional volumes were obtained.


## 2.7. Preprocessing of Diffusion Data

The MD and FA measures of diffusion tensor imaging (DTI) [50] can detect unique aspects of brain microstructural neural plasticity [15, 51], including those caused by WMT [15, 52, 53]. A reduction in MD reflects an increase in the tissue density of various cellular structures, including capillaries, synapses, spines, and macromolecular proteins, as well as changes in the properties of myelin, membranes, and axons; the shapes of neurons or glia; and tissue organization. However, MD is not specific to any one of these microlevel changes [50, 51]. On the other hand, an increase in FA is thought to have a relatively greater association with changes in microstructural properties related to brain connectivity, including myelination, axonal membrane thickness and diameter, and degree of axonal parallel organization [50, 52].Preprocessing of imaging data was performed using Parametric Mapping software (SPM8; Wellcome Department of Cognitive Neurology, London, UK) implemented in Matlab (Mathworks Inc., Natick, MA, USA). The descriptions in this subsection are a summary from our previous publication [15]. Pre- and posttraining MD and FA maps were independently segmented and normalized using the previously described diffeomorphic anatomical registration through an exponentiated Lie algebra (DARTEL)-based method (Takeuchi et al., 2013a). This method accounts not only for the outer edge but also for the signal distribution within the white matter during registration. The reason for the independent normalization has been described previously (Takeuchi et al., 2013a).From the pre- and postintervention normalized images of the normalized MD map, areas not strongly likely to be gray or white matter in our custom template (defined by “gray matter tissue probability + white matter tissue probability < 0.99”) were removed. Then, from the pre- and postintervention normalized images of normalized FA map, areas not strongly likely to be white matter in our custom template (defined by “white matter tissue probability < 0.99”) were removed. Subsequently, normalized MD and FA images were smoothed by convolution with an isotropic Gaussian kernel of 10- and 8-mm full-width at half maximum, respectively. These images were used for the group analysis described below.


## 2.8. Preprocessing and Data Analysis for Functional Activation Data

Preprocessing and data analysis were performed using SPM8 implemented in Matlab. Prior to analysis, BOLD images from the pretraining scan and the posttraining scan were realigned and resliced to the mean BOLD image of all pretraining scans. The images were then corrected for slice timing, coregistered to ab=0 diffusion-weighted image, and spatially normalized using the parameter of the abovementioned normalization of FA images to form images with 3×3×3mm3 voxels.A design matrix was fitted to each subject with one regressor for each task condition (0-back and 2-back in theN-back task) using the standard hemodynamic response function (HRF). The design matrix weighted each raw image according to its overall variability to reduce the impact of movement artifacts [54]. After estimation, beta images were smoothed (8-mm full-width at half maximum) and either taken to the second level or subjected to random-effect analysis. We removed low-frequency fluctuations using a high-pass filter with a cut-off value of 128 s. Individual-level statistical analyses were performed using a general linear model.In individual analyses, we focused on the activation change and compared activation relative to the conditions (0-/2-back versus rest) before and after an intervention period. The resulting maps for each subject represented changes between the pre- and postmeasures of brain activity during the corresponding conditions as well as the brain activity during the corresponding conditions in the premeasures, which were then included in a group analysis.


## 2.9. Preprocessing of Structural Data

Preprocessing of imaging data was performed using SPM12 implemented in Matlab. Using the Pairwise Longitudinal Registration toolbox [55], the midpoint average template images were generated for each participant as well as maps of differences between the Jacobian determinants and of pre- and postintervention T1-weighted structural images using default parameter settings. The midpoint average image of each subject was then segmented, and the Thorough Clean option was used to remove any odd voxels. Affine regularization was performed in accordance with the East Asian brain with a sampling distance (the approximate distance between the sampled points when estimating the model parameters) of 1 mm. Default parameter settings were used for other parameters in the segmentation process. Then, by multiplying maps of differences between the Jacobian determinants of pre- and postintervention images and gray matter segments of the midpoint average image of each subject, images of gray matter volume change were generated. Then, diffeomorphic anatomical registration was performed using DARTEL processes. Here, we used imported DARTEL images of gray matter and white matter for the processes created above. The DARTEL template was created from the entire set of subjects in this experiment, and the default parameter settings were used for the DARTEL process. Then, using the parameters generated for each participant, the images of gray matter change were normalized and smoothed by convolving each with an isotropic Gaussian kernel of 8 mm FWHM.


## 2.10. Statistical Group-Level Analysis of Imaging and Behavioral Data

There were numerous potential contrasts in this study, because there were three groups with two factors (aerobic exercise and WMT). However, we only tested two-step analyses to confirm our hypothesis. The first step was the superiority of SDAEWMT over single training (i.e., comparison of SDAEWMT with the combined groups of WMT and AET) which was followed by the second step that confirmed the superiority of SDAEWMT over AET as well as the superiority of SDAEWMT over WMT (the comparison between SDAEWMT and AET and the comparison between SDAEWMT and WMT). Behavioral data were analyzed using SPSS 22.0 (SPSS Inc., Chicago, IL). In behavioral analysis, because the superiority of SDAEWMT was the primary study hypothesis (and we are not aware of any study indicating that dual-task training is inferior to single-task training in terms of cognitive transfer), test–retest changes within the SDAEWMT group were compared with those in the WMT or AET group using one-tailed, one-way analysis of covariance (ANCOVA) with the posttest measures as dependent variables and pretest scores as independent variables (P<0.05). This ANCOVA design was selected to increase the sensitivity of analysis. As described in our previous study [56], we employed ANCOVAs instead of repeated measure ANOVAs to control the effects of pretest scores. Statistical experts strongly recommend to use ANCOVA instead of repeated measure ANOVA in this type of study design [57]. With randomized designs, the purpose of ANCOVA is to reduce error variance, whereas with nonrandomized designs (or with analyses involving substantial preexisting group differences), ANCOVA is used to adjust the posttest means for pretest differences among groups [57]. One might recommend that posttest scores be used, instead of using the differences between pre- and posttest measures. However, in practice, when the pretest scores are included as covariates, the two analyses return the same P value.The use of a one-tailed test was consistent with previous studies [27, 47, 58, 59]. In tests of the effects of SDAEWMT group compared with those of others, for post hoc comparisons between SDAEWMT and AET groups as well as comparisons between SDAEWMT and WMT groups, when ANCOVAs showed the same pattern (P<0.1) of superiority for the effects of SDAEWMT, the results were regarded as showing the superiority of SDAEWMT over single training (here the main contrast is SDAEWMT vs AET + WMT; however, the existence of similar trends in the contrasts of (SDAEWMT vs AET) and (SDAEWMT vs WMT) was confirmed with a more liberal threshold). In all nonwhole-brain analyses (i.e., behavioral analyses) performed in this study, results with a threshold of P<0.05 corrected for false discovery rate (FDR) using the graphically sharpened method [60] were considered statistically significant. The correction for multiple comparisons using this method was applied to the results of 10 ANCOVAs that tested the superiority of SDAEWMT over single training. ANOVAs were utilized to evaluate group differences in cognitive performance before the intervention. In group-level imaging analyses, we tested for group-wise differences in the change of each measure across the entire brain.For determining the MD, FA, and brain activity, analyses were performed using Biological Parametrical Mapping [61], the extension of SPM, which allow the inclusion of signal of other images at every voxel, to perform voxel-by-voxel ANCOVAs. In analyses of MD and FA, we performed voxel-wise ANCOVAs with the values of postscans at each voxel as dependent variables and the values for the prescans at each voxel as independent variables. In analysis of the N-back task, we performed voxel-wise ANCOVAs with the contrast of activity (0-back or 2-back) changes between pre- and postscans at each voxel as dependent variables and the corresponding value for the prescan at each voxel, the accuracy of the corresponding tasks in the prescan, and its pre- to postchanges (post–pre) as independent variables. For the N-back task, only the subjects who achieved an accuracy of 80% or higher during both pre- and postsessions were included in the analysis, because when the task difficulty was too high for the subjects, brain activity showed nonlinear patterns and such nonlinear changes could not be assessed by the statistical tests applied [62]. In other words, areas that showed an increase in activity as the load increased started to exhibit a decrease in activity when the task difficulty became too high. By this inclusion, we could focus only on linear changes and could ignore the effects of “not doing the tasks properly.” Consistently, Jansma et al. [48] showed that 75% accuracy is too low for this type of N-back task among healthy participants to detect linear changes (brain activity in relevant areas increased from 0-back to 2-back task in which accuracy was about 82%, but dropped substantially in the 3-back condition in which the accuracy was about 72.5% in the case of healthy participants). Since the random response can achieve 25% accuracy rate in this task, 75% accuracy is actually not high, and it means they are substantially failing. Consequently, 13 WMT group subjects (mean accuracy, pre: 90.0±6.8%, post: 91.9±7.2%), 12 SDWMAET group subjects (mean accuracy, pre: 90.8±8.1%, post: 94.2±5.3), and 13 AET group subjects (mean accuracy, pre: 89.6±6.6%, post: 91.5±6.0%) were included in analyses of the 2-back task. Without surprisingly, there were no significant differences of pre- to postchanges in accuracy among these subjects after correcting preintervention accuracy. 28 WMT group subjects, 26 SDWMAET group subjects, and 27 AET group subjects were included in analyses of the 0-back task. In behavioral analyses of the N-back task, only the subjects who performed the task without misunderstanding of the rules were included from analyses. In these neuroimaging analyses, the effects of SDAEWMT compared with those of single training were tested using the contrasts of −1×WMT2×SDAEWMT−1×AET or 1×WMT−2×SDAEWMT1×AET with the uncorrected masks of −1×WMT1×SDAEWMT or 1×WMT−1×SDAEWMT and 1×SDAEWMT−1×AET or −1×SDAEWMT1×AET (P<0.05, uncorrected). Regions with significance were inferred using cluster-level statistics [63]. Only clusters with P<.05, after correction for multiple comparisons at cluster size, with a voxel-level cluster-determining threshold of P<.001 (uncorrected) were considered statistically significant.In the group-level imaging analysis of rGMV, we tested for group-wise differences in the change in rGMV across the whole brain with the same contrasts. We used the one-way ANOVA option in SPM8. In the imaging analysis, the effects of training with TCSSL were estimated by comparing differences between the groups of rGMV, as calculated above at each voxel. This analysis was limited to areas where the average rGMV value for the segmented and normalized mean images of pre- and postscans of all participants was greater than 0.1. A multiple comparison correction was performed using threshold-free cluster enhancement (TFCE) [64] with randomized (5,000 permutations) nonparametric testing using the TFCE toolbox (http://dbm.neuro.uni-jena.de/tfce/). We applied a threshold of family-wise error corrected at P<0.05. The standard cluster size test is inappropriate in voxel-based morphometry [65], whereas a permutation-based test is assumed to be correct. However, the abovementioned analysis using BPM cannot be applied to rGMV analysis because the TFCE toolbox cannot be used in BPM, which is why we chose this statistical design for rGMV analysis.


## 3. Results

### 3.1. Basic Data and Results of Training Data

The numbers of training sessions for subjects who completed the intervention in each training group are presented in Table1. There were no significant preexisting group differences in the general cognitive function (MMSE score, FAB), VO2Max, education length, working memory, or short-term memory performance (digit span performance, 2-back task performance), as shown in Table 1 (ANOVA, P>0.05).Training data for WMT and SDAEWMT are presented in Table2. The increase in the average of the highest performance scores (N) for trained WM tasks from the first three training sessions (excluding the first session, which was completed as a practice and explanation session) to the last three sessions was not significantly greater with SDAEWMT than with WMT after correcting for the highest performances of trained WM tasks during the first three training sessions, excluding the first session (P=0.833 by one-tailed ANCOVA). The best performance among the three sessions was used because the brevity of each session for each training task would lead to instability of results if we used the data from just one training day. This same group procedure was applied previously [52, 56, 66].Table 2
Highest performance scores (mean±SD) of WMT and SDAEWMT groups on trained WM tasks among the second to fourth training sessions and final three training sessionsa.
WMTSDAEWMTFirst three sessionsLast three sessionsFirst three sessionsLast three sessionsDevilish calculations2.08±0.233.05±0.392.02±0.242.98±0.42Devilish pairs (oni-mekuri)5.93±0.937.63±1.256.13±0.927.57±1.45Devilish mice2.47±0.733.15±0.792.38±0.403.18±0.65Devilish shapes2.02±0.092.15±0.392±02.13±0.29Devilish blocks4.43±1.055.93±1.314.43±1.125.70±1.16Devilish cups4.68±0.565.62±0.804.62±0.565.62±0.8Devilish listening1.75±0.342.30±0.401.57±0.402.27±0.7aThe data from the first training session were excluded because the first training session mainly consisted of practice and explanation of the training tasks.Overall, the expected improvement of cognitive training task performance was not observed. On the other hand, most subjects successfully completed DAEMT training as the performance of each specific training task was comparable with that of WMT group subjects.We also recorded subjects’ heart rates during the resting state prior to exercise (HRRSBE) in each training session and used as one measure of cardiovascular function. Themean±SD lowest HRRSBE from the first three training sessions (excluding the first session, which was completed as a practice and explanation session) and the last three training sessions were 75.61±8.37 beats per minute (bpm) and 70.09±8.48bpm, respectively, in the AET group, and 76.63±10.46bpm and 74.07±6.59bpm, respectively, in the SDAEWMT group (we focused on the best performance among the three training days to prevent instability of the results based on a single trial. This method is consistent with those of our previous studies [15, 52, 56]). The decrease in HRRSBE in the AET group tended to be greater than that in the SDAEWMT group after correcting for the effect of the lowest HPRSBE value recorded during the first three training sessions, excluding the first session (P=0.052, two-tailed ANCOVA). This finding suggests that although all the subjects in AET and SDAEWMT groups received the same instructions to keep their heart rate within a fixed range, the subjects in the SDAEWMT group may have been unable to follow this instruction all the time when they also had to concentrate on cognitive training. Therefore, the effects of aerobic exercise on HRRSBE may have been reduced during SDAEWMT.


### 3.2. Analyses of Cognitive Test Data and fMRI Task Behavioral Data

Compared with the WMT or AET group, the SDAEWMT group showed a significantly greater increase in 2-back task accuracy and FAB score (for statistical values including rawP values, see Table 3). This difference remained significant after correcting for multiple comparisons (both P=0.040, corrected for FDR). However, the SDAEWMT group did not show a significantly larger improvement in 2-back task accuracy than the WMT group. The greater increase in FAB performance in the SDAEWMT group than in the WMT or AET group was moderated by the greater increase in the SDAEWMT group than in the WMT group as well as the significantly greater increase in the SDAEWMT group than in the AET group (Table 3).Table 3
Pre- and posttest scores for psychological measures (mean±SEM).
WMTSDAEWMTAETRawP valueaP valuebPrePostPrePostPrePostSDAEWMT vs. othersPretest differences0-back accuracy95.06±7.0996.89±4.3893.10±7.8896.21±8.0194.22±8.8293.87±9.770.2600.6562-back accuracy73.52±20.8188.52±11.4573.52±20.990.00±9.8174.68±20.1676.45±20.530.008c0.970RAPM12.97±4.4513.97±5.0214.67±5.4115.50±5.0713.67±4.5013.39±5.670.1790.401Digit span14.20±3.0515.00±2.7613.90±2.7014.87±3.0312.76±2.8714.52±3.190.6640.121Logical memory22.41±6.1226.69±6.6125.37±6.8526.07±6.521.39±6.5823.67±7.870.9380.056D-CAT41.63±4.9841.95±6.6842.33±6.9342.83±6.5341.59±6.8542.21±7.050.4490.884FAB15.57±1.6515.07±1.915.17±2.115.6±1.7015.03±1.8214.21±2.040.008d0.514Symbol search34.00±5.2135.46±5.5534.07±5.1434.1±5.4432.15±5.0034.36±6.210.9750.257Digit symbol66.83±14.3967.73±12.365.90±11.4669.63±12.5163.45±13.2669.48±14.230.4300.581SACT17.67±5.1518.80±6.0815.83±4.7916.83±4.8215.30±5.5716.83±4.820.8370.187aOne-way analysis of covariance (ANCOVA) tests comparing the SDAEWMT group with the other two groups for test–retest differences, with psychological measures as dependent variables and pretest scores for the psychological measures as covariates. For all tests, one-tailed tests were performed to assess the superiority of SDAEWMT. bOne-way ANOVAs comparing the pretest scores of the three groups. cOne-tailed, one-way ANCOVAs revealed that the increase in the SDAEWMT group was not significantly greater than that in the WMT group (P=0.302) but that there was a significantly greater increase in the SDAEWMT group than in the AET group (P=0.0003). dOne-tailed, one-way ANCOVAs revealed a trend toward a greater increase in the SDAEWMT group than in the WMT group (P=0.096) and a significantly greater increase in the SDAEWMT group than in the AET group (P=0.003).


### 3.3. MD Analyses

We also compared MD changes in the SDAEWMT group with those in the WMT and AET groups. This analysis revealed a significant SDAEWMT-related decrease from pre- to postmeasurement in MD regions, including areas adjacent to the right dorsolateral prefrontal cortex (DLPFC) and dorsal anterior cingulate cortex, areas in the right orbitofrontal cortex, areas in the left hippocampus, midbrain areas involving the bilateral substantia nigra, ventral tegmental area, periaqueductal gray and reticular formation, right ventral anterior cingulate cortex, and part of the left putamen and left pallidum (Figure5(a), Table 4). There were no significant results of SDAEWMT-related increase of MD from before measurement to after measurement.Figure 5
Change in diffusion imaging measures related to SDAEWMT. Decrease in mean diffusivity (MD) in the SDAEWMT group compared with that in the WMT and AET groups. The results are shown with a threshold ofP<0.05, corrected for multiple comparisons at the cluster level, with an underlying voxel-level P<0.001, uncorrected. Compared with the WMT and AET groups, SDAEWMT training resulted in decreased MD in areas adjacent to the right dorsolateral prefrontal cortex (DLPFC) and dorsal anterior cingulate cortex (dACC), areas in the right orbitofrontal cortex, areas in the left hippocampus, midbrain areas involving the bilateral substantia nigra, ventral tegmental area, periaqueductal gray and reticular formation, ventral anterior cingulate cortex (vACC), and parts of the left putamen and left pallidum. The color intensity represents the T value. The color intensity represents the T value.Table 4
SDAEWMT-related decrease in MD compared with the control group.
AreaMNI coordinatesT scoreCorrectedP value (cluster)Cluster sizexyzMidbrain (B)/hippocampus (L)/putamen (L)/pallidum (L)-13.5-16.5-13.55.18<0.0013071Middle frontal gyrus (L)-34.5946.54.180.014901Superior frontal gyrus (R)/dorsal anterior cingulate cortex (R)13.540.5333.960.0021337Orbital frontal gyrus (R)/ventral anterior cingulate cortex (R)1230-93.850.0031185MNI: Montreal Neurological Institute.


### 3.4. FA Analyses

When FA for the SDAEWMT group was compared with that for the WMT and AET groups, there were no significant results of SDAEWMT-related changes of FA from before measurement to after measurement.


### 3.5. Brain Activity Analyses

Brain activity was assessed during the 2-back task, and the difference between the pre- and posttraining measurements was compared among the groups. The SDAEWMT group showed a significant increase in brain activity (meaning that the brain activity corresponding to the contrast of post versus pre is greater than in other groups) in the left temporoparietal junction [MNI coordinates:x,y,z=−63,−42,15; t=4.68; P=0.002 (3159 mm3), corrected for multiple comparisons with the cluster determining threshold of P<0.001] and in anatomical clusters that included the right temporoparietal junction [MNI coordinates: x,y,z=69,−36,27; t=4.68; P=0.007 (2430 mm3), corrected for multiple comparisons with the cluster determining threshold of P<0.001] and the right superior temporal gyrus compared with the WMT and AET groups (Figure 6(a)). There were no significant results of SDAEWMT-related decrease of brain activity from before measurement to after measurement.Figure 6
Brain activity changes related to training. Increase in brain activity during the 2-back WM task in the SDAEWMT group compared with that in the WMT and AET groups. The results are shown with a threshold ofP<0.05, corrected for multiple comparisons at the cluster level, with an underlying voxel-level P<0.001, uncorrected. Compared with AET, WMT training resulted in increased brain activity during the 2-back WM task in bilateral areas of the temporoparietal junction. The color intensity represents the T value.


### 3.6. rGMV Analyses

The SDAEWMT-related changes in the rGMV from before measurement to after measurement were not statistically significant.


## 3.1. Basic Data and Results of Training Data

The numbers of training sessions for subjects who completed the intervention in each training group are presented in Table1. There were no significant preexisting group differences in the general cognitive function (MMSE score, FAB), VO2Max, education length, working memory, or short-term memory performance (digit span performance, 2-back task performance), as shown in Table 1 (ANOVA, P>0.05).Training data for WMT and SDAEWMT are presented in Table2. The increase in the average of the highest performance scores (N) for trained WM tasks from the first three training sessions (excluding the first session, which was completed as a practice and explanation session) to the last three sessions was not significantly greater with SDAEWMT than with WMT after correcting for the highest performances of trained WM tasks during the first three training sessions, excluding the first session (P=0.833 by one-tailed ANCOVA). The best performance among the three sessions was used because the brevity of each session for each training task would lead to instability of results if we used the data from just one training day. This same group procedure was applied previously [52, 56, 66].Table 2
Highest performance scores (mean±SD) of WMT and SDAEWMT groups on trained WM tasks among the second to fourth training sessions and final three training sessionsa.
WMTSDAEWMTFirst three sessionsLast three sessionsFirst three sessionsLast three sessionsDevilish calculations2.08±0.233.05±0.392.02±0.242.98±0.42Devilish pairs (oni-mekuri)5.93±0.937.63±1.256.13±0.927.57±1.45Devilish mice2.47±0.733.15±0.792.38±0.403.18±0.65Devilish shapes2.02±0.092.15±0.392±02.13±0.29Devilish blocks4.43±1.055.93±1.314.43±1.125.70±1.16Devilish cups4.68±0.565.62±0.804.62±0.565.62±0.8Devilish listening1.75±0.342.30±0.401.57±0.402.27±0.7aThe data from the first training session were excluded because the first training session mainly consisted of practice and explanation of the training tasks.Overall, the expected improvement of cognitive training task performance was not observed. On the other hand, most subjects successfully completed DAEMT training as the performance of each specific training task was comparable with that of WMT group subjects.We also recorded subjects’ heart rates during the resting state prior to exercise (HRRSBE) in each training session and used as one measure of cardiovascular function. Themean±SD lowest HRRSBE from the first three training sessions (excluding the first session, which was completed as a practice and explanation session) and the last three training sessions were 75.61±8.37 beats per minute (bpm) and 70.09±8.48bpm, respectively, in the AET group, and 76.63±10.46bpm and 74.07±6.59bpm, respectively, in the SDAEWMT group (we focused on the best performance among the three training days to prevent instability of the results based on a single trial. This method is consistent with those of our previous studies [15, 52, 56]). The decrease in HRRSBE in the AET group tended to be greater than that in the SDAEWMT group after correcting for the effect of the lowest HPRSBE value recorded during the first three training sessions, excluding the first session (P=0.052, two-tailed ANCOVA). This finding suggests that although all the subjects in AET and SDAEWMT groups received the same instructions to keep their heart rate within a fixed range, the subjects in the SDAEWMT group may have been unable to follow this instruction all the time when they also had to concentrate on cognitive training. Therefore, the effects of aerobic exercise on HRRSBE may have been reduced during SDAEWMT.


## 3.2. Analyses of Cognitive Test Data and fMRI Task Behavioral Data

Compared with the WMT or AET group, the SDAEWMT group showed a significantly greater increase in 2-back task accuracy and FAB score (for statistical values including rawP values, see Table 3). This difference remained significant after correcting for multiple comparisons (both P=0.040, corrected for FDR). However, the SDAEWMT group did not show a significantly larger improvement in 2-back task accuracy than the WMT group. The greater increase in FAB performance in the SDAEWMT group than in the WMT or AET group was moderated by the greater increase in the SDAEWMT group than in the WMT group as well as the significantly greater increase in the SDAEWMT group than in the AET group (Table 3).Table 3
Pre- and posttest scores for psychological measures (mean±SEM).
WMTSDAEWMTAETRawP valueaP valuebPrePostPrePostPrePostSDAEWMT vs. othersPretest differences0-back accuracy95.06±7.0996.89±4.3893.10±7.8896.21±8.0194.22±8.8293.87±9.770.2600.6562-back accuracy73.52±20.8188.52±11.4573.52±20.990.00±9.8174.68±20.1676.45±20.530.008c0.970RAPM12.97±4.4513.97±5.0214.67±5.4115.50±5.0713.67±4.5013.39±5.670.1790.401Digit span14.20±3.0515.00±2.7613.90±2.7014.87±3.0312.76±2.8714.52±3.190.6640.121Logical memory22.41±6.1226.69±6.6125.37±6.8526.07±6.521.39±6.5823.67±7.870.9380.056D-CAT41.63±4.9841.95±6.6842.33±6.9342.83±6.5341.59±6.8542.21±7.050.4490.884FAB15.57±1.6515.07±1.915.17±2.115.6±1.7015.03±1.8214.21±2.040.008d0.514Symbol search34.00±5.2135.46±5.5534.07±5.1434.1±5.4432.15±5.0034.36±6.210.9750.257Digit symbol66.83±14.3967.73±12.365.90±11.4669.63±12.5163.45±13.2669.48±14.230.4300.581SACT17.67±5.1518.80±6.0815.83±4.7916.83±4.8215.30±5.5716.83±4.820.8370.187aOne-way analysis of covariance (ANCOVA) tests comparing the SDAEWMT group with the other two groups for test–retest differences, with psychological measures as dependent variables and pretest scores for the psychological measures as covariates. For all tests, one-tailed tests were performed to assess the superiority of SDAEWMT. bOne-way ANOVAs comparing the pretest scores of the three groups. cOne-tailed, one-way ANCOVAs revealed that the increase in the SDAEWMT group was not significantly greater than that in the WMT group (P=0.302) but that there was a significantly greater increase in the SDAEWMT group than in the AET group (P=0.0003). dOne-tailed, one-way ANCOVAs revealed a trend toward a greater increase in the SDAEWMT group than in the WMT group (P=0.096) and a significantly greater increase in the SDAEWMT group than in the AET group (P=0.003).


## 3.3. MD Analyses

We also compared MD changes in the SDAEWMT group with those in the WMT and AET groups. This analysis revealed a significant SDAEWMT-related decrease from pre- to postmeasurement in MD regions, including areas adjacent to the right dorsolateral prefrontal cortex (DLPFC) and dorsal anterior cingulate cortex, areas in the right orbitofrontal cortex, areas in the left hippocampus, midbrain areas involving the bilateral substantia nigra, ventral tegmental area, periaqueductal gray and reticular formation, right ventral anterior cingulate cortex, and part of the left putamen and left pallidum (Figure5(a), Table 4). There were no significant results of SDAEWMT-related increase of MD from before measurement to after measurement.Figure 5
Change in diffusion imaging measures related to SDAEWMT. Decrease in mean diffusivity (MD) in the SDAEWMT group compared with that in the WMT and AET groups. The results are shown with a threshold ofP<0.05, corrected for multiple comparisons at the cluster level, with an underlying voxel-level P<0.001, uncorrected. Compared with the WMT and AET groups, SDAEWMT training resulted in decreased MD in areas adjacent to the right dorsolateral prefrontal cortex (DLPFC) and dorsal anterior cingulate cortex (dACC), areas in the right orbitofrontal cortex, areas in the left hippocampus, midbrain areas involving the bilateral substantia nigra, ventral tegmental area, periaqueductal gray and reticular formation, ventral anterior cingulate cortex (vACC), and parts of the left putamen and left pallidum. The color intensity represents the T value. The color intensity represents the T value.Table 4
SDAEWMT-related decrease in MD compared with the control group.
AreaMNI coordinatesT scoreCorrectedP value (cluster)Cluster sizexyzMidbrain (B)/hippocampus (L)/putamen (L)/pallidum (L)-13.5-16.5-13.55.18<0.0013071Middle frontal gyrus (L)-34.5946.54.180.014901Superior frontal gyrus (R)/dorsal anterior cingulate cortex (R)13.540.5333.960.0021337Orbital frontal gyrus (R)/ventral anterior cingulate cortex (R)1230-93.850.0031185MNI: Montreal Neurological Institute.


## 3.4. FA Analyses

When FA for the SDAEWMT group was compared with that for the WMT and AET groups, there were no significant results of SDAEWMT-related changes of FA from before measurement to after measurement.


## 3.5. Brain Activity Analyses

Brain activity was assessed during the 2-back task, and the difference between the pre- and posttraining measurements was compared among the groups. The SDAEWMT group showed a significant increase in brain activity (meaning that the brain activity corresponding to the contrast of post versus pre is greater than in other groups) in the left temporoparietal junction [MNI coordinates:x,y,z=−63,−42,15; t=4.68; P=0.002 (3159 mm3), corrected for multiple comparisons with the cluster determining threshold of P<0.001] and in anatomical clusters that included the right temporoparietal junction [MNI coordinates: x,y,z=69,−36,27; t=4.68; P=0.007 (2430 mm3), corrected for multiple comparisons with the cluster determining threshold of P<0.001] and the right superior temporal gyrus compared with the WMT and AET groups (Figure 6(a)). There were no significant results of SDAEWMT-related decrease of brain activity from before measurement to after measurement.Figure 6
Brain activity changes related to training. Increase in brain activity during the 2-back WM task in the SDAEWMT group compared with that in the WMT and AET groups. The results are shown with a threshold ofP<0.05, corrected for multiple comparisons at the cluster level, with an underlying voxel-level P<0.001, uncorrected. Compared with AET, WMT training resulted in increased brain activity during the 2-back WM task in bilateral areas of the temporoparietal junction. The color intensity represents the T value.


## 3.6. rGMV Analyses

The SDAEWMT-related changes in the rGMV from before measurement to after measurement were not statistically significant.


## 4. Discussion

In the present study, we developed and successfully implemented an ecologically applicable SDAEWMT (simultaneously performed dual-task training involving both aerobic exercise training and working memory training) method and examined its effects on cognitive function (subjects could successfully perform cognitive training while doing the exercise in this setting), brain structure, and brain activity patterns in comparison with those of WMT and AET methods. We showed that (a) SDAEWMT significantly improved the performance of the FAB, but there were no significant results in other broader transfer measures. These results were only partly consistent with our hypotheses that (a) SDAEWMT would be more effective than WMT or AET alone and/or that SDAEWMT would have effects on cognitive functions such as global cognitive functions, memory performance, attention, and executive functions, and (b) SDAEWMT resulted in decreases in the MD in the frontal and subcortical brain areas involving the dopaminergic system. And SDAEWMT produced an increase in WM-related brain activity in the bilateral temporoparietal junction. These results were partly overlapping and partly distinctive compared to our hypothesis (b), which SDAEWMT would alter brain structures and functional activities involving frontal and parietal areas. (c) SDAEWMT did not result in a greater improvement of the performance of WMT tasks compared with WMT alone, which is inconsistent with our hypothesis that SDAEWMT facilitates the improvement of WMT performance.The present behavioral results suggest that SDAEWMT has a specific effect on executive function; however, there was no evidence of broader transfer effects or enhanced learning ability. This result was based on results showing that SDAEWMT significantly improved the performance of FAB, which consisted of the battery for global executive function, compared with the other two training groups. This result is consistent with a previous review suggesting that simultaneously performed dual-task training incorporating physical training is effective in improving subject performance in this type of general cognitive battery [24]. However, there were no other cognitive functions that showed such patterns. Furthermore, SDAEWMT did not appear to facilitate learning processes for WMT performance, as indicated by the lack of significantly better performance of training tasks compared with the WMT group. We obtained this result despite previous studies reporting that aerobic exercise facilitates learning processes in mice [67]. The reason for this lack of facilitation and for the lack of broad transfer to a wide range of cognitive functions remains unclear. However, as described below, at least some of the WMT tasks employed in this study may have been too difficult for the older adults in our subject cohort. Furthermore, the simultaneously performed dual-task setting may have made it difficult for the exercising subjects to simultaneously focus on the cognitive training tasks. However, these are only speculations, and future studies are required to establish whether SDAEWMT can enhance cognitive training learning processes under different conditions.The observed MD changes in the SDAEWMT group were found in parts of three major dopaminergic systems. ACC and DLPFC are both involved in top-down attentional control and dual-task execution [7, 68], while the hippocampus is involved in long-term memory [69]. On the other hand, the majority of areas where the effects of SDAEWMT on MD were observed are included in the three main dopaminergic systems of the brain. One of these systems is the nigrostriatal system, which projects from the substantia nigra to the striatum and is involved in motor control as well as some higher cognitive functions [70]. The second system is the mesolimbic system, which projects from the ventral tegmental area to regions such as the hippocampus and nucleus accumbens, and is involved in the reward system [70]. The third system is the mesocortical system, which projects from the ventral tegmental area to the prefrontal cortex and is involved in WM and other higher-order cognitive functions [70]. The notion that the dopaminergic system plays at least a partial role in simultaneously performed dual tasks involving physical movement has been suggested by research on Parkinson’s disease, which is characterized by degeneration of the dopaminergic system and deficits in such dual-task performance [71]. In addition, while coordination of multiple cognitive processes is assumed to be important for WM [72], as already described, different types of WMT have been shown to affect the dopaminergic system in the frontal-parietal areas and the striatum [73, 74]. Lower MD is known to reflect higher tissue densities resulting from various conditions [50], and this measure is reported to be sensitive to use-dependent neural plasticity [51]. However, recent studies have shown that MD in the dopaminergic system is uniquely related to pathological, pharmacological, and cognitive differences or changes directly related to dopamine [15, 75–77]. The changes in MD observed in the present study may reflect alterations in the dopaminergic functions of the three dopaminergic systems. Through these changes, effective dual-task performance and improved frontal lobe functions may be possible. However, future physiological or imaging studies using positron emission topography are required for confirmation.Functional imaging analysis of WM-related brain activity revealed that SDAEWMT increased WM-related brain activity in areas mainly located around the bilateral temporoparietal junctions, which function in attentional mechanisms. These regions are critical areas of the ventral attention system [78] and are involved in reorienting attention to an outside event or switching attention between different matters [78]. This region is typically deactivated during WM task performance, potentially to prevent attention from reorienting to distracting events [78]. The elderly can learn to control attention and perform dual task efficiently as they are prepared for simultaneously performed dual-task training involving exercise [79]. Also, a dual-task process typically requires task switching [80]. We believe that the same process partially underlies how the elderly exhibited improved performance in both aerobic capacity and training tasks. Thus, the fact that activity in this region increased following SDAEWMT suggests that through this training, the subjects were able to focus better on outside events or to switch between different tasks. Through these changes, the effective execution of simultaneously performed dual tasks during SDAEWMT may become possible. However, this is speculation based on reverse inference. Future studies are required to elucidate the changes in attentional mechanism induced by SDAEWMT.The successful implementation of simultaneously performed dual-task training involving exercise in the present study may have implications in daily life. As described in theIntroduction, both AET and simultaneously performed dual-task training involving exercise have been shown to be effective in the elderly. However, there is generally a fall risk in the elderly when performing tasks that involve locomotion [81]. In addition, it is difficult to perform these WMTs in daily life while walking or running. Thus, the use of a recumbent bike may be applicable to a wide range of cognitive activities for the elderly in addition to the present dual-training regimen and may make a wide range of ecologically valid and safe training regimens possible.The reasons for the lack of significant results in rGMV and FA analyses despite significant effects of WMT on those measures in young adults [14] are not clear. However, previous studies generally suggest the reduction of neural plasticity in the elderly [82], and also, a previous study showed that 1-year AET led to greater rGMV increases than 6-month AET in some of the brain regions [83]. The length of training might not be long enough to obtain significant results among the different intervention methods in the rGMV analyses.There are several limitations to this study. First, we could not include passive or active control groups that would have revealed null effects in any of the measures because of our resource limitations. Accordingly, we could not detect potential common effects among the three groups, although this design can test the hypothesis that simultaneously combined physical and cognitive training would be more effective than single physical or single cognitive training, as was done in an earlier review [26]. In particular, we could not identify any significant effects of training on digit span or divergent thinking performance This lack of significance may be unexpected given previous studies reporting consistent WMT effects on untrained WM performance, including effects of N-back task training on digit span performance [11, 84] and the consistently reported effects of AET on divergent thinking performance [85]. However, a previous three-month intervention study involving older adults also did not show any increase in digit span performance in two control groups (0.3 and 0 point decreases from pre- to posttest, respectively) [86]. Furthermore, in our studies on young adults, including two separate one-month intervention studies, S-A creativity test performance did not show a mean increase in the passive control group [66]. On the other hand, the 93 subjects who completed the present intervention showed an increase in digit span performance (1.194 point increase, P=5.473×10−15, paired t-test) and S-A creativity test performance (1.527 point increase, P=1.931×10−18, paired t-test). In the present study, the game used for WMT was apparently too difficult for the older subjects to complete. In particular, in each of the three N-back WMT tasks, the final performance level was only slightly different from the lowest level for the WM condition (2-back task). The initial performance levels for these tasks were also low, and based on nonquantitative observation of training, the subjects appeared to have difficulty in adopting an effective strategy for completing the task, consistent with previous reports [87]. These observations reveal a less effective adaptive difficulty control (i.e., when the span increases from 3 to 4, the level increases by 33%, which is twice as large as when the span increases from 6 to 7), resulting in only modest improvement in tasks such as the digit span. Future studies evaluating young adult subjects or incorporating easier WMT tasks for older adults may be warranted to address this problem. Further, in this study, measurements of cardiovascular fitness were not performed postexperiment as they were not central to the purpose of the study and as it takes resources to perform the procedures. HRRSBE was gathered, but only from participants in the AET and SDAEWMT groups. Thus, the effects of aerobic exercise might not have been sufficiently evaluated in this study. Also, in the AET and SDAEWMT groups, participants could select the order of tasks that they conducted during the training sessions. Although this procedure might mitigate the toughness of training, it might increase the unwanted variables among training groups and reduce the sensitivity of the analysis. Further, there is a lack of follow-up experimentation after the training ceased, to evaluate how long the training effects lasted; this lack is due to the practical limitation of the various kinds of resources. Intervention studies in the field often take place without an opportunity to evaluate the duration of the training effects, and future studies should investigate this aspect. Finally, all the subjects were asked to engage in training three times each week in a laboratory, but if a subject was not available during a given week, he/she was allowed to make up for the lost sessions in the following weeks. This was necessary as participants have daily lives that do not always fit the laboratory schedule. This kind of imperfect control of the training protocol might affect the findings.From the present results, we showed the effects of SDAEWMT and made the following conclusions: (a) SDAEWMT improved executive function in the frontal lobe. (b) There was no evidence of broader transfer effects such as in memory performance and attention measures (in contrast to our hypothesis) or in processing speed, divergent thinking, or other effects that enhanced learning. (c) SDAEWMT resulted in neural tissue changes in brain areas involving the dopaminergic system. (d) SDAEWMT produced an increase in WM-related brain activity in brain areas involved in attentional reorienting.



---

*Source: 3859824-2020-06-05.xml*