# Striatal and Hippocampal Involvement in Motor Sequence Chunking Depends on the Learning Strategy

Ovidiu Lungu, Oury Monchi, Geneviève Albouy, Thomas Jubault, Emanuelle Ballarin, Yves Burnod, Julien Doyon

► To cite this version:

Ovidiu Lungu, Oury Monchi, Geneviève Albouy, Thomas Jubault, Emanuelle Ballarin, et al. Striatal and Hippocampal Involvement in Motor Sequence Chunking Depends on the Learning Strategy. PLoS ONE, 2014, 9 (8), pp.e103885. 10.1371/journal.pone.0103885. hal-01366364

HAL Id: hal-01366364

https://hal.sorbonne-universite.fr/hal-01366364v1

Submitted on 14 Sep 2016

**HAL** is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire **HAL**, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Striatal and Hippocampal Involvement in Motor Sequence Chunking Depends on the Learning Strategy

Ovidiu Lungu¹,²,³,⁴*, Oury Monchi¹,²,⁸, Geneviève Albouy¹,²,⁵, Thomas Jubault¹,², Emanuelle Ballarin⁶, Yves Burnod⁷, Julien Doyon¹,²,⁵

¹Unité de Neuroimagerie Fonctionnelle (UNF), Montréal, Canada, ²Centre de recherche de l'Institut universitaire de gériatrie de Montréal, Montréal, Canada,

³Département de Psychiatrie, Université de Montréal, Montréal, Canada, ⁴Center for Research in Aging, Donald Berman Maimonides Geriatric Center, Montréal, Canada,

⁵Département de Psychologie, Université de Montréal, Montréal, Canada, ⁶Graduate Program in Neuropsychology, Université de Lyon 2, Lyon, France, ⁷INSERM U678, Université de Paris VI Jussieu, Paris, France, ⁸Département de Radiologie, Université de Montréal, Montréal, Canada

<figure><img src="image_2.png" /></figure>

CrossMark
click for updates

## Abstract

Motor sequences can be learned using an incremental approach by starting with a few elements and then adding more as training evolves (e.g., learning a piano piece); conversely, one can use a global approach and practice the whole sequence in every training session (e.g., shifting gears in an automobile). Yet, the neural correlates associated with such learning strategies in motor sequence learning remain largely unexplored to date. Here we used functional magnetic resonance imaging to measure the cerebral activity of individuals executing the same 8-element sequence after they completed a 4-days training regimen (2 sessions each day) following either a global or incremental strategy. A network comprised of striatal and fronto-parietal regions was engaged significantly regardless of the learning strategy, whereas the global training regimen led to additional cerebellar and temporal lobe recruitment. Analysis of chunking/grouping of sequence elements revealed a common prefrontal network in both conditions during the chunk initiation phase, whereas execution of chunk cores led to higher mediotemporal activity (involving the hippocampus) after global than incremental training. The novelty of our results relate to the recruitment of mediotemporal regions conditional of the learning strategy. Thus, the present findings may have clinical implications suggesting that the ability of patients with lesions to the medial temporal lobe to learn and consolidate new motor sequences may benefit from using an incremental strategy.

**Citation:** Lungu O, Monchi O, Albouy G, Jubault T, Ballarin E, et al. (2014) Striatal and Hippocampal Involvement in Motor Sequence Chunking Depends on the Learning Strategy. PLoS ONE 9(8): e103885. doi:10.1371/journal.pone.0103885

**Editor:** Edwin Robertson, Harvard Medical School, United States of America

**Received** September 20, 2013; **Accepted** July 8, 2014; **Published** August 22, 2014

**Copyright:** © 2014 Lungu et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** The work received the following support: Julien Doyon: Canadian Institutes of Health Research (MOP86463), Natural Sciences and Engineering Research Council of Canada (42909-2009); Ovidiu Lungu: Postdoctoral fellowship from the Ministère du Développement économique, de l'Innovation et de l'Exportation (MDEIE) du Québec, Québec, QC, Canada; Geneviève Albouy: Supported by a Postdoctoral fellowship from the Fondation FYSSEN, France and a C.I.H.R. postdoctoral grant; and Oury Monchi: Natural Sciences and Engineering Research Council of Canada (No 327518). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

* Email: o.lungu@umontreal.ca

## Introduction

Motor sequences are ubiquitous in everyday life, from simple behaviors such as preparing a cup of coffee to complex activities like speaking and dancing. As such, numerous studies in the past have investigated the neuronal correlates and mechanisms implicated in motor sequence acquisition [1,2]. As part of procedural memory processes [3], motor sequence learning is thought to take place in stages [4], whereby considerable improvements in motor performance are known to occur rapidly during the early stage of the acquisition process. These improvements are then followed by smaller gains spread out over longer periods of time across subsequent practice sessions in later stages. A large body of neuroimaging evidence has revealed that motor sequence learning is mediated by the cortico-striatal and cortico-cerebellar circuits in the early stage, regardless of whether the subject knows the sequence explicitly or not before practice begins, and by the cortico-striatal system during the later learning phases [2,5]. Yet changes in hippocampal activity has also been reported to be associated with both implicit and explicit motor sequence

learning [6–9], hence highlighting its ability to associate discontinuous but structured information.

Recently, there has been an increased interest in investigating the behavioral and neurophysiological determinants of long-term motor sequence memory consolidation. This process can be facilitated by sleep (both day and night sleep) and it was found to be dependent upon the cognitive functions recruited during the acquisition process [10–14]. Moreover, consolidation appears to be based on increased activity within the striatum [15,16] and mediotemporal lobe (MTL), especially the hippocampus [6,8,13,17–19]. Interestingly, the interaction of both striatum and hippocampus with frontal areas during initial training appears to facilitate the implementation of reproducible motor behavior [7]. While the role of the cortico-striatal circuit in long-term motor sequence acquisition seems to be in the grouping of sequence elements into single-action units (i.e., in creating chunks) [20–24], the functional contribution of the hippocampus could be in the detection and formation of higher-order sequential associations [9,25–27] and the stabilization of the motor performance [7]. Yet one important issue that has entirely been overlooked so far by the neuroimaging literature is the extent to which the involvement of
neural substrates found in later stages of motor sequence learning, such as the striatum and hippocampus, is modulated by the *type* of training regimens or learning strategies employed for acquiring a new motor skilled behavior.

In everyday life, when learning explicit sequences, such as steps of a dance or to play a new piece of music, we usually use an incremental approach, that is we start by practicing a part of the sequence of movements to be learned, and then expand to include more and more elements until we are able to perform it wholly. Yet in other occasions (e.g., learning to shift gears to manually drive an automobile), we may need to practice the whole sequence all at once (global approach) before we put it into use (e.g., driving in traffic). When using the incremental training regimen, the sequence representation changes and gradually builds up in complexity, whenever new elements are added to the sequence. This gradual build-up of sequence representation should be associated with a slowdown of motor performance during learning. For instance, in a series of behavioral experiments, Ganor-Stern and colleagues [28] showed that changes in earlier versus later sequence elements (i.e. the 3rd versus the 6th element of an 7-element sequence) were found to lead to a greater impairment in performance. These findings suggest, on the one hand, that participants built-up the sequence representation by chunking together the earlier elements first, then followed by the later elements, but also showed that motor performance is impaired whenever sequence representation is changed. By contrast, with the global learning strategy, stable and complex sequence representations are thought to be formed. For instance, it was shown that abstract, effector independent sequence representations can form rapidly and early on during the motor skill acquisition process, especially when the sequence is known explicitly [29,30], while motor, effector dependent representations are formed slowly after extended periods of training [29].

Given the scarcity of studies comparing directly these two learning strategies, some important knowledge gaps remain to be addressed. For instance, to what extent differences in sequence representation, as a result of incremental versus global training, are manifested behaviorally, and more specifically in the way sequence elements are chunked or grouped together? Also, neuroimaging studies on long-term explicit motor sequence learning employed exclusively a global learning strategy, whereby the same sequence was practiced over many training sessions, spanning days or weeks. In most cases, the contribution of striatum and hippocampus in this process has been assessed by comparing motor performance and brain activity between sessions (i.e. last vs. first). Thus, it is still unknown whether the long-term maintenance of the skill (i.e. producing the well-rehearsed and consolidated sequence at the end of multi-session training) is subserved by different neural substrates as a result of employing different learning strategies (i.e. incremental versus global). Here we directly compare the production of the same explicit motor sequence and we investigate, using fMRI, the contribution of the cortico-striatal and MTL regions in this process at the end of a 4-day incremental versus global behavioral training regimen of a sequence of movements.

# Materials and Methods

## Participants

Thirty-two subjects aged 19–36 years old were initially recruited for this study. About two thirds of the participants were female (23 out of 32) and the participants’ recruitment was done via advertisement on the university campus. None of the subjects had a history of neurological or psychiatric disorders as revealed

by anamnesis prior to the experiment. Out of the 32 participants, 4 did not complete the behavioral training regimen and consequently, they were excluded from the data analysis. Six others did terminate the behavioral training, but did not participate in the imaging session due to unanticipated technical problems with the scanner’s cooling system during the week they were trained; and thus their data were also excluded from the analysis. The data from the remaining 22 subjects is reported here.

## Ethics statement

All participants gave their voluntary consent and signed a form, which was approved by the joint research ethics committee of the Regroupement Neuroimagerie Québec at the “Centre de recherche, Institut universitaire de gériatrie de Montréal”. This committee, which approved the study and its experimental procedures, follows the guidelines of the Tri-Council Policy Statement of Canada, the civil code of Quebec, the Declaration of Helsinki and the code of Nuremberg.

## Experimental design and procedure

All participants followed a four-day behavioral training regimen with the aim of learning an explicit 8-element motor sequence by the end of it. In the fifth day, learning was assessed in a test session performed in the MRI scanner. Thus, the training routine always started on Monday and finished on Thursday, with two daily sessions (morning and afternoon) of behavioral training. The following day in the morning (always on Friday), subjects were tested on the motor sequence learning task while their brain activity was recorded. One group, ‘Global’ (n₁ = 11), always practiced the entire 8-element sequence in each of the training sessions, except the 4th (Global training). By contrast, the other group, ‘Incremental’ (n₂ = 11), started off with only two elements in the first training session, and was then asked to practice an increasingly more difficult sequence as one more element was added to the sequence on each subsequent session (except the 4th) (Incremental training) (see Figure 1A). Also, during behavioral training, each group was exposed to two conditions: “Sequence” where subjects practiced the 8-element target sequence, either globally or incrementally, and “Control” where participants practiced a simpler, 4-element sequence (for details, see task description). In all but the first and last training sessions of the Sequence condition subjects executed 16 blocks of 96 sequential movements each, interspersed with 10 seconds rest periods. In the first and last sessions subjects also performed separately 16 blocks with 96 trials each in the Sequence condition and 16 blocks with the same number of trials in the Control condition. In the 4th session, participants were exposed only to the Control condition (Figure 1A), again for 16 blocks with 96 trials. The next morning, after the 4-day behavioral training, subjects were tested in a single session with 5 imaging runs of 8 blocks each, interspersed with 30 seconds rest periods. In each scanning run, they were exposed to alternating pairs of blocks within the Sequence or Control condition. This particular distribution of the Sequence and Control training sessions over the 4-day training schedule was chosen to ensure that subjects in the Incremental group trained on a different sequence in each session. In addition, the decision to have morning and afternoon training sessions, rather than one per day or at bigger intervals, was made to minimize subject attrition rate because, otherwise, the training would have required subjects to commit to the study for more than a week.

## Task description

Subjects were tested using a version of the explicit motor sequence learning task similar to that developed by Karni and
<figure><img src="image_3.png" /><figcaption>B.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>C.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Experimental groups and sessions, the task and the automatic chunking procedure. **A.** Control (light gray) and target (dark gray) sequences practiced by the two groups of subjects across training sessions. The numbers 1 to 4 correspond to four fingers of the left hand, from the little (1) to the index (4) finger. **B.** The task interface presenting the executing instructions, the block information and the type of feedback after each trial. **C.** Reaction time (RT) in a block for a given subject and the classification of trials into initial element of a chunk (red squares) and the body or core of the chunk (blue triangles) by our automatic clustering algorithm.
doi:10.1371/journal.pone.0103885.g001</figcaption></figure>

colleagues [4,31,32]. In this task, participants were required to perform self-generated finger movements with their non-dominant (left) hand. They were shown four grey discs arranged in a row on a computer screen, either directly (during behavioral training) or through a mirror embedded within the head coil (in the scanner, at retest). Each disk corresponded to one of four buttons on a numerical keyboard (behavioral training) or on a MR compatible button box (in the scanner). Under the row of discs, subjects were also shown a small window comprising the following written instruction cues: “Pause” when subjects had to rest between

blocks, or “Go”, when they had to execute the finger movements (Figure 1B). Before the onset of each session, subjects were explicitly informed about the sequence of movements that had to be executed in each block. In the Sequence condition, participants in the Global group were asked to practice the following finger sequence: 2-4-1-3-4-2-3-1, where 1 represents the small finger and 4 represents the index finger of the non-dominant hand. By contrast, participants in the Incremental group were required to practice only parts of this sequence, with one additional element being added to the sequence on each subsequent session. For
<table><caption>**Table 1.** Error rates (percentages) across sessions, type of block, experimental condition and Group.</caption><thead><tr><th rowspan="3"></th><th colspan="4">GROUP</th></tr><tr><th colspan="2">Incremental</th><th colspan="2">Global</th></tr><tr><th>Mean</th><th>SEM</th><th>Mean</th><th>SEM</th></tr></thead><tbody><tr><td colspan="5"><strong>Behavioral training</strong></td></tr><tr><td colspan="5"><strong>Control condition</strong></td></tr><tr><td>Session 1</td><td>3.56</td><td>0.91</td><td>3.27</td><td>0.97</td></tr><tr><td>Session 4</td><td>3.37</td><td>0.65</td><td>1.88</td><td>0.47</td></tr><tr><td>Session 8</td><td>4.78</td><td>1.22</td><td>6.60</td><td>4.26</td></tr><tr><td colspan="5"><strong>Sequence condition</strong></td></tr><tr><td>Session 1</td><td>4.34</td><td>0.84</td><td>6.96</td><td>1.92</td></tr><tr><td>Session 2</td><td>4.02</td><td>1.10</td><td>5.00</td><td>1.44</td></tr><tr><td>Session 3</td><td>4.04</td><td>0.79</td><td>3.69</td><td>1.35</td></tr><tr><td>Session 5</td><td>3.71</td><td>0.59</td><td>2.37</td><td>0.46</td></tr><tr><td>Session 6</td><td>5.99</td><td>1.17</td><td>2.56</td><td>0.94</td></tr><tr><td>Session 7</td><td>6.01</td><td>1.59</td><td>2.23</td><td>0.57</td></tr><tr><td>Session 8</td><td>5.37</td><td>1.26</td><td>6.26</td><td>4.06</td></tr><tr><td colspan="5"><strong>Test (fMRI session)</strong></td></tr><tr><td>Control condition</td><td>3.60</td><td>0.74</td><td>1.70</td><td>0.30</td></tr><tr><td>Sequence condition</td><td>3.21</td><td>0.43</td><td>2.05</td><td>0.33</td></tr></tbody></table>

**Control condition during behavioral training.**

Session main effect: $F_{2,40} = 1.33; p = 0.27$.

Group main effect: $F_{1,20} = 0.95; p = 0.76$.

Session *Group interaction: $F_{2,40} = 0.59; p = 0.55$.

**Sequence condition during behavioral training.**

Session main effect: $F_{6,120} = 1.11; p = 0.35$.

Group main effect: $F_{1,20} = 0.22; p = 0.63$.

Session *Group interaction: $F_{6,120} = 1.63; p = 0.14$.

**Control and Sequence condition during test (fMRI session).**

Condition main effect: $F_{1,20} = 0.005; p = 0.94$.

Group main effect: $F_{1,20} = 3.63; p = 0.07$.

Condition *Group interaction: $F_{1,20} = 1.21; p = 0.28$.

doi:10.1371/journal.pone.0103885.t001

instance, in the first session they were asked to perform the first two elements of the sequence (2-4) only; in the second session – the first three elements (2-4-1), and so on, until the last training session where subjects had to perform all 8 elements of the sequence. This sequence was chosen for the following reasons: (1) all elements appear with the same frequency, (2) it is complex in structure because there are only second-order conditional transitions between elements (i.e. one needs to know a minimum of two immediately preceding elements to predict the current element), and (3) it needs some amount of practice in order to achieve asymptotic performance because there are no easy consecutive triplets (e.g. 1-2-3 or 4-3-2), which can be performed automatically. A simple sequence (1-2-3-4) was also employed as a control condition, which was administered during the 1st, 4th and 8th practice sessions (see Figure 1A). This sequence was chosen as control because it can automatize rapidly (i.e. asymptotic performance) due to the presence of consecutive elements and single-order conditional transitions. Subjects started each session with a 10 seconds rest period (30 seconds during the imaging session), indicated by the instruction cue. Five seconds prior to each block, the instruction cue informed the subject about the type of movements to be executed next: “Sequence” for the target sequence or “Control” for the control condition (Figure 1B). Each block of practice was initiated as soon as the cue changed to “Go”, hence instructing the subjects that they had to start executing the

sequence of movements at will. They were instructed to execute the movements as quickly as possible, while making as few errors as possible. Immediately after the subject pressed a button, the corresponding discs on the screen were switched to green (if subjects chose the correct button, a hit) or to red (if subjects chose a different button, an error). The location of the illuminated disc always indicated the correct location for that particular trial (Figure 1B). If an error was committed, subjects were instructed to continue the sequence in the proper order, with the next element. Therefore, the task design provided subjects in each trial with a visual feedback regarding their own performance.

**Behavioral data analysis**

**Reaction time.** Several dependent measures have been computed in past studies using serial reaction time task paradigms. For instance, in the classical serial reaction time task [3], learning-specific effects have been assessed by comparing the difference in reaction time between sequential and randomly presented elements. In our task, because motor responses were self-generated, it was impossible to have a random condition. Therefore, different learning metrics were employed. One was to compare the improvements in sequential motor performance time across sessions. Yet, by design, this was only possible in the Global group where the sequence was unchanged from one session to the next. Alternatively, the performance in the Control condition (the
sequence 1-2-3-4) can be taken as a benchmark for automaticity or asymptotic sequential performance; thus, the smaller the difference in performance between Sequence and Control conditions, the better the learning. By design, the two groups, Global and Incremental, did not practice the same number of elements from the target sequence during the behavioral training period (i.e. the Global group practiced always all 8 elements, whereas the Incremental group started with 2 elements in the first session and arrived to practice the whole 8-element sequence only in the last training session). Thus, the behavioral data could only be analyzed in situations where a match between the two groups was found in this regard. To this end, we analyzed the blocks from the Control condition during Sessions 1 and 4 and both the blocks from the Control and Sequence condition in the last behavioral training session (Session 8) and at test (fMRI session). For this analysis, we employed a general linear model (GLM) approach, where the Group (Global vs. Incremental), Session (Sessions 1st vs 4th vs 8th for Control condition only) or Group and Condition (Control vs. Sequence, during the 8th training or fMRI sessions) were used as fixed independent factors, and the reaction time (RT) of correct trials was the dependent measure. We used RTs per trial rather than the correct number of sequences per block or the average time to execute a sequence as a dependent variable, because we were also interested in obtaining a measure of the subject's ability to chunk sequence elements together based on differences in RTs from one trial to the next (see the "Trial chunking" section below). Subjects were always entered as a random factor in each analysis. Sidak tests were employed to test post-hoc differences between means in order to account for the number of multiple pairwise comparisons.

**Error rates.** In most serial reaction time paradigms, error rates (i.e. proportion of errors over the total number of trials) are usually low and rarely constitute the main variable of interest. Nevertheless, we measured, reported and analyzed them in the current study. We computed the error rate for each block of 96 trials, in each condition, session and for each participant. Then, these block-related error rates were averaged by session, condition and group. We then used GLM models to assess differences in error rates as a function of session, condition and group both for the behavioral training data, as well as for the data acquired during the imaging session.

**Trial chunking.** A separate analysis was performed on the chunking of movements within a given block; the number and length of these chunks provide insights into the learning strategy that subjects used based upon the type of training (Global vs. Incremental) and allows for group comparison with respect to this aspect. As shown in previous studies [20,22,33-35] people tend to group the elements of a motor sequence in chunks or clusters, like dialing a long phone number by chunking its different parts such as the country code, area code, and number, for example. These chunks are characterized by a fast and similar reaction time for each of its elements, except the first. For example, previous investigations have shown that, during motor chunking, the subject's reaction time to execute the first sequence element is always slower due to transitions between chunks (i.e. the end of one and beginning of another), probably reflecting a starting cost or a higher memory load due to the fact that the previous cluster is discarded and another one is loaded [20,36,37]. In order to analyze the difference in chunking strategy between the two Global and Incremental groups, we devised an automatic algorithm that can categorize a given trial in a block as being the first element of a chunk [INITIAL], part of the rest of the elements in that chunk [BODY], or as being in neither one category or the other [OTHER]. This algorithm first computed the mean RT and the standard error of the mean for each block of

trials and for each subject. Then, all correct first trials in each block and all the correct trials immediately after an error were automatically classified as [INITIAL]. If the differences in reaction time between the ($i$)th correct trial and the two adjacent correct trials, ($i-1$)th and ($i+1$)th, were both higher than three standard errors of the mean for that particular block AND higher than two standard errors of the mean when compared with the ($i+2$)th correct trial, then the ($i$)th trial was categorized as INITIAL. Correct trials following an INITIAL trial were categorized as part of the BODY as long as the difference in reaction time between them remained equal or less than two standard errors of the mean. The errors and the rest of the correct trials that could not be classified as [INITIAL] or [BODY] were categorized as OTHER. Panel C in Figure 1 shows the reaction time in a block of trials for one participant, as well as the trial's categorization in either an [INITIAL] or [BODY] element. It is interesting to note that the first element of the sequence (trial #1, #9, #17, etc.) was overwhelmingly categorized as an [INITIAL] element by this automatic algorithm.

**Variability in performance.** Another way to measure the effects of our experimental manipulation on learning the motor sequence is to analyze the subjects' variability in movement execution. The more the sequence becomes consolidated and the more subjects are able to reproduce it in the same fashion, the lower the variability [38]. The evolution in performance variability is believed to reflect the implementation of a performance mode that would represent the whole sequence of movements in motor memory [38]. In order to perform this analysis for the behavioral training sessions, we computed the standard deviation for the RTs of correct trials only for each Sequence block, in each session and for each subject. We averaged the data across blocks for each subject in each session, hence obtaining a single value per subject and per session. These standard deviation measures were then used as dependent variable in a mixed repeated measures analysis of variance (ANOVA), with sessions as a within-subject factor and groups (Global vs. Incremental) as a between-subjects factor. A similar analysis approach was used for the behavioral data acquired during the fMRI session, except that instead of averaging the measures for the entire session, we averaged the standard deviations of RTs across the blocks of a given fMRI run. Therefore, in the mixed repeated-measures ANOVA, we used the run as a within-subject factor, while preserving the group as between-subjects factor.

## Imaging parameters

A 3 Tesla whole Body MR System (Magnetom TIM, Siemens Medical Systems, and Erlangen, Germany) was used for image acquisition. Prior to the functional run, 176 structural images were acquired in sagittal plane by using an MPRAGE imaging sequence (TR = 13 ms; TE = 4.92 ms; FA = 25°; FoV = 256 mm²; matrix size = 256×256; slice thickness = 1 mm, voxel size = 1 mm³). Then, whole brain fMRI was performed using an echo-planar imaging (EPI) sequence measuring blood oxygenation level dependent (BOLD) signal (TR = 2510 ms; TE = 30 ms; FA = 90°; FoV = 220 mm²; matrix size = 64×64; slice thickness = 3.4 mm, voxel size = 3.4 mm³; 41 slices). The functional slices covered the whole brain, were oriented in transverse plane and were angled to be parallel to the AC-PC line. An inline retrospective motion correction algorithm was employed while the EPI images were acquired. A total of 205 functional volumes were recorded in each functional run.
## Preprocessing of fMRI data

Brain Voyager QX (Brain Innovation B.V., Maastricht, the Netherlands) software was used for fMRI data preprocessing and analysis. The functional bi-dimensional images of every subject were preprocessed to correct for the difference in time slice acquisition (slice scan time correction). In addition to linear detrending, a high-pass filter of three cycles per time course (frequency domain) was applied to the corrected 2D slices. The functional series were then preprocessed to correct for possible motion artifacts in 3D space, and to ensure that movements in any plane did not exceed 3 mm. These functional images were subsequently used to reconstruct the 3D functional volume for every subject, and for every run. The 3D functional volumes were then aligned with the corresponding 3D anatomical volume, and both were normalized to standard Talairach space (Talairach and Tournoux, 1988. Spatial smoothing using a Gaussian kernel at 8 mm full width at half maximum (FWHM) was applied to the 3D functional data.

## Imaging data analysis

Both block-design and event-related analysis approaches were employed.

**Block design analysis.** In each group and within each run, two predictors describing the type of block or experimental conditions were defined (i.e. Sequence vs. Control). Then, three statistical contrasts were tested: one to assess group differences (Global vs. Incremental) regardless of experimental condition, one to evaluate the main effect of condition (Sequence vs. Control), regardless of the group membership and one contrast assessing the Group\*Condition interaction.

**Event-related design analysis.** For the event-related design, the following predictors were defined, based on the analysis of behavioral data at individual trial level (see the clustering of trials in the behavioral data analysis): (1) first element of a chunk in the Sequence blocks [*initial SEQ*], (2) the rest of the elements of a chunk in the Sequence blocks [*body SEQ*], (3) first element of the chunk in the Control blocks [*initial CTR*], (4) the rest of the elements within a chunk in the Control blocks [*body CTR*] and (5) all other elements that could not be categorized by the algorithm in one the previous four categories [*OTHER*]. The contrasts of interest in this analysis were *initial SEQ > body SEQ* and *initial SEQ < body SEQ*, in order to identify the brain regions in which the activity was higher during chunk initiation than during the execution of main chunk core or vice versa. These contrasts were tested (1) combining both groups together, (2) as a group difference (interaction contrast) and (3) separately within each group.

In both block and event-related analyses, the predictors were first entered as fixed factors in single subject GLM's, and then the parameters of this model were subsequently entered into a second level of analysis corresponding to a random-effect GLM model that was used for group analysis [39]. The statistical parameters of this latter model were estimated voxelwise for the entire brain and activation maps were computed for various contrasts between the predictors. While these contrasts were computed and will be reported for the whole brain, our discussion of the imaging results section will focus mostly on the cortico-striatal network and mediotemporal areas known to be involved in motor sequence learning. When displaying activation maps reflecting contrasts between predictors for all of the subjects, regardless of the group membership, the false discovery rate [$q(FDR)<0.05$] was employed as a criterion to correct for multiple comparisons with a minimum cluster size of 108 adjacent significant isovoxels (108 cubic millimeters in volume) that surpassed this threshold. When

displaying the activation maps for the contrasts reflecting group differences or performed separately for each group of subjects, we employed the same cluster size and a statistical threshold for each voxel in the cluster of at least $p<0.005$ (uncorrected) for group differences contrasts and of $p<0.001$ (uncorrected) for contrasts performed within each group.

Further GLM analyses were performed within the cortico-striatal and mediotemporal regions of interest that were activated significantly in previous contrasts. In these analyses the average percentage of BOLD signal change in the whole cluster was considered as dependent measure, with the group and block type as independent factors. The data for this analysis was extracted from Brain Voyager and the GLMs were analyzed by using SPSS software (SPSS Inc., Chicago, IL, USA).

## Results

### Behavioral results: Error rates and reaction time

Error rates were not significantly different as a function of any of the independent variables (session, type of block, groups), either alone or in interaction with each other (Table 1). Analysis of the mean reaction time data during correct trials in the 8th session (the last behavioral practice session, where both groups performed the same 8-element target sequence) as a function of experimental group (Global vs. Incremental) and the type of block (Control vs. Sequence) revealed a significant interaction effect [$F_{1,20} = 15.13; p<0.05$], as well as significant main effects of group [$F_{1,20} = 7.12; p<0.05$] and block type [$F_{1,20} = 16.48; p<0.05$]. Pairwise comparisons revealed that the interaction was driven by a significant difference between control and sequence blocks among subjects in the Incremental group (187.23 msec. vs. 290.53 msec.), whereas those in the Global group displayed similar RTs (136.27 msec. vs. 138.48 msec.) in both types of blocks (Figure 2A). These results show that, at the end of behavioral training, the Global condition group developed some level of automaticity regarding the target sequence as indicated by the fact that their reaction time was similar to that of the control sequence. Also, the main effect of group suggests that this level of automaticity in performance among participants in the Incremental group was lower than that among participants in the Global condition. However, our results suggest that the Incremental group also learned the target sequence. This is illustrated by the fact that their performance in the 8th session, that is when they were required to execute the full sequence for the first time, was significantly better than that of the Global group in the first training session [$F_{1,20} = 4.34; p<0.05$]. In fact, performance of the Incremental group executing the whole 8-element sequence in the 8th training session did not differ from that of the Global group in their 2nd [$F_{1,19} = 0.77; p = 0.78$] session.

RT analysis of the control sequence during the 1st, 4th and 8th training sessions revealed no significant group x session interaction [$F_{2,40} = 0.09; p = 0.91$], nor any significant main effect of group [$F_{1,20} = 2.33; p = 0.14$]. Yet the results showed a significant main effect of session [$F_{2,40} = 50.85; p<0.05$] (Figure 2B). Post-hoc Sidak tests also revealed that the two groups did not differ in their reaction time during the execution of control sequence in any of the three sessions. The latter results suggest that while the two groups improved over time for the control session, indicating a practice effect, the rate of improvement for the control sequence was similar in each group.

The RT analysis during the imaging session (9th) showed a significant interaction effect [$F_{1,20} = 6.35; p<0.05$] between the group and block type in the first run only (Figure 2C). For the remaining 4 runs of the sequence and control blocks, there was no
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 2. Behavioral results.** A. Mean reaction time for the target sequence in the training sessions, for each group. B. Mean reaction time for the control sequence in sessions 1, 4 and 8, for each group. C. Mean reaction time for target and control sequences and for each group during the fMRI (9th) session.
doi:10.1371/journal.pone.0103885.g002</figcaption></figure>
A.

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

B.

<figure><img src="image_11.png" /><figcaption>**Figure 3. Sequence chunking and variability of performance.** A. Mean number of clusters (first row) and the percentage of clusters with 3 or more elements (second row) in session 8th (left, behavioral) and 9th (right, imaging), for participants in the two groups. B. The mean variability in motor performance during the execution of the target sequence during 8 behavioral training sessions (left) and neuroimaging runs, during 9th session (right).</figcaption></figure>

doi:10.1371/journal.pone.0103885.g003

significant effect of group [$F_{1,20} = 0.56$; *p* = 0.81 for control blocks and $F_{1,20} = 0.99$; *p* = 0.32 for sequence blocks], nor any group * run interaction [$F_{3,60} = 0.91$; *p* = 0.44 for control blocks and $F_{3,60} = 0.84$; *p* = 0.47 for sequence blocks]. In addition, no significant interaction was found between run and type of block (sequence vs. control), when considering the data from the Incremental [$F_{3,30} = 1.34$; *p* = 0.27] and Global [$F_{3,30} = 1.93$; *p* = 0.14] group separately. The overall 3-way interaction between group, type of block and run was also not significant [$F_{3,60} = 0.31$; *p* = 0.81]. Altogether, these results suggest that performance of the control and target sequences of the two groups in the last 4 fMRI runs was similar, and we thus subsequently

decided to exclude the first run from the analysis of the imaging data.

## Behavioral results: clustering and variability

The clustering analysis of the 8th training session (Figure 3A, left panels) showed that the number of clusters within a block detected by our automatic algorithm was not significantly affected by the Global vs. Incremental group conditions [$F_{1,20} = 0.38$; *p* = 0.54], the type of block [$F_{1,20} = 0.76$; *p* = 0.39], and the group x type of block interaction was also not significant [$F_{1,20} = 0.0004$; *p* = 0.98]. However, the percentage of clusters with more than 3 or 4 elements in length were both significantly higher among subjects in the Incremental than in the Global group
<table><caption>**Table 2.** Brain regions activated by contrast [Sequence blocks>Control blocks].</caption><thead><tr><th>Contrasts and regions</th><th>Talairach coordinates</th><th>Volume (mm³)</th><th>Maximum t</th></tr></thead><tbody><tr><td colspan="4"><strong>Sequence>Control blocks (all subjects)</strong></td></tr><tr><td>Right precuneus (BA7)</td><td>21</td><td>-62</td><td>49</td><td>2698</td><td>6.12</td></tr><tr><td>Right SMA (BA6)</td><td>24</td><td>-14</td><td>53</td><td>730</td><td>5.29</td></tr><tr><td>Left precuneus (BA7)</td><td>-26</td><td>-55</td><td>49</td><td>5952</td><td>7.00</td></tr><tr><td>Left SMA (BA6)</td><td>-24</td><td>-12</td><td>54</td><td>2147</td><td>7.68</td></tr><tr><td>Lentiform nucleus (putamen)</td><td>-21</td><td>6</td><td>16</td><td>140</td><td>4.87</td></tr><tr><td colspan="4"><strong>Sequence>Control (Global>Incremental)</strong></td></tr><tr><td>Right STG (BA21)</td><td>54</td><td>-22</td><td>-2</td><td>830</td><td>3.84</td></tr><tr><td>Right culmen (cerebellum)</td><td>39</td><td>-37</td><td>-26</td><td>465</td><td>4.58</td></tr><tr><td>Right VLPFC (BA47)</td><td>39</td><td>29</td><td>1</td><td>224</td><td>4.45</td></tr><tr><td>Right declive (cerebellum)</td><td>33</td><td>-61</td><td>-11</td><td>665</td><td>4.33</td></tr><tr><td>Right parahippocampal gyrus (BA35)</td><td>27</td><td>-28</td><td>-20</td><td>347</td><td>3.90</td></tr><tr><td>Right PCC (BA31)</td><td>18</td><td>-55</td><td>22</td><td>342</td><td>4.05</td></tr><tr><td>Left culmen (cerebellum)</td><td>0</td><td>-52</td><td>-5</td><td>532</td><td>4.18</td></tr><tr><td>Left declive (cerebellum)</td><td>-3</td><td>-64</td><td>-11</td><td>153</td><td>3.62</td></tr></tbody></table>

SMA – supplementary motor area.

STG – superior temporal gyrus.

VLPFC – ventro-lateral prefrontal cortex.

PCC – posterior cingulate cortex.

doi:10.1371/journal.pone.0103885.t002

$[F_{1,20} = 9.47; p<0.05]$ for the percentage of clusters with more than 3 elements and $F_{1,20} = 7.37; p<0.05$ for the percentage of clusters with more than 4 elements in length. Such findings suggest that the type of learning (global versus incremental) led, at the end of training, to the development of longer clusters among participants in the Incremental than in the Global group condition, despite the fact that the total number of clusters was comparable among the two groups.

During the last 4 runs of the imaging session, an analysis measuring trial clustering (number of clusters per block, and the percentage of clusters with more than 3 or 4 elements) revealed no significant difference between the two groups, across runs and as a function of block type [$F_{3,60} = 0.31; p = 0.81$ for the total number of clusters; $F_{3,60} = 0.68; p = 0.56$ for the percentage of clusters with more than 3 elements and $F_{3,60} = 1.53; p = 0.21$ for the percentage of clusters with more than 4 elements]. These results (Figure 3A, right panels) indicate that the two groups chunked the sequence elements to the same degree and in similar sizes during the imaging session.

Subjects' average performance variability within blocks during the different behavioral training sessions revealed a crossover pattern (Figure 3B, left panel), where the variability decreased across sessions in the Sequence blocks for subjects in the Global group, but increased for those in the Incremental group. This pattern was expressed in a significant interaction effect between the type of group and session [$F_{6,120} = 20.54; p<0.05$]. Pairwise comparisons showed that the two groups were different in variability in all sessions (all Sidak tests $p<0.05$), except in the second one, where the crossover actually occurs. Moreover, after the third session, the variability in performance remained at the same level in both groups (all Sidak tests $p>0.05$). These results suggest that participants in the Global group became stable across sessions and less variable within session in their performance of the motor sequence during last 4 sessions, whereas those in the Incremental group were more variable within session, albeit pretty stable across sessions. In the fMRI session, the same analysis

revealed again a significant interaction between the type of group and run [$F_{4,80} = 2.98; p<0.05$]. However, Sidak tests revealed that this effect came from a significant difference between the first run and the rest of them for the subjects in the Incremental group only (Figure 3B, right panel), whereas individuals in the Global group showed no change in variability across runs.

Taken together, the behavioral results suggest that despite differences observed at the end of training, the participants in the two experimental groups displayed similar motor performance, both in terms of the speed of execution and clustering during the last four runs of the imaging session. The only difference between participants in the two groups was observed at the level of motor sequence reproducibility, as measured by the level of variability level, which was always higher in the Global than in the Incremental condition, a carryover effect of the nature of training during procedural learning.

## Imaging results: block contrasts analysis

As mentioned above, data from the first run of the imaging session were excluded from the analysis because behavioral performance of the two groups differed in terms of reaction time (Fig. 2C), and because the participants in the Incremental group had a higher variability relative to the subsequent runs (Fig. 3B, right panel). This ensured that differences in speed of movement execution would not 'contaminate' the observed differences in brain activation, and that only the carryover effect of different training regimens would be expressed in the neuroimaging data. Yet another reason to exclude the first run stemmed from the habituation effects to the new environment. During the behavioral training sessions, the participants were seated in front of a computer screen and were asked to use a keyboard to respond. By contrast, during the fMRI session, subjects were required to lie down in the scanner in a supine position and to use a MR compatible keypad. Thus, we believe that whatever changes in the behavioral performance during the first run of the fMRI session,
## A. Areas commonly activated in the two conditions
Contrast: Sequence>Control

<figure><img src="image_12.png" /></figure>

## B. Areas differentially activated in the two conditions
Contrast: Sequence>Control

<figure><img src="image_13.png" /></figure>
**Figure 4. Imaging results.** A. Areas in the cortico-striatal network activated more during Sequence than Control blocks, and that are commonly found in the two groups. B. Areas in ventro-orbital prefrontal cortex and parahippocampus were more activated in participants in the Global than Incremental condition, during Sequence than Control blocks.
doi:10.1371/journal.pone.0103885.g004

this could be, at least in part, attributed to changes in the body position and response device.

While contrasts were computed on the whole brain and the complete pattern of results is reported in Table 2, our description and discussion of the results will focus mostly on the cortico-striatal network and mediotemporal areas. The contrast [Global vs. Incremental], performed for both experimental conditions combined, did not reveal any activation cluster which surpassed the statistical threshold (minimum t-value $t_{(21)}$ = 4.30; FDR corrected) in either direction (i.e. Global>Incremental or Incremental>Global). The contrast [Sequence>Control] performed for both groups together identified 5 clusters of activation (Figure 4A, Table 2), one in the left putamen, and the other four, respectively, in the supplementary motor areas (SMA) and precuneus (Brodmann area BA7) bilaterally. Furthermore, in each of these clusters, the BOLD signal change was significantly higher during

the sequence than the control blocks, for each group (main effect of block; all *p*<0.05). In addition, in the left SMA cluster only, there was a significant interaction effect [*F*<sub>1,20</sub> = 6.04; *p*<0.05], indicating that the difference between the two types of blocks was higher in the Incremental than in the Global condition.

The same contrast analysis carried out to assess group differences revealed no region in which subjects in the Incremental group showed more activation than those in the Global group. In 8 clusters, however, the difference in BOLD signal between sequence and control blocks was significantly greater in the Global than the Incremental group (Figure 4B, Table 2). Half of these clusters were found within the cerebellum, with the remaining being located – in ventro-orbito prefrontal (BA47), posterior cingulate cortices (BA31), superior temporal (BA21) and parahippocampal gyri (BA35). Furthermore, the contrast [Sequence>Control] performed separately within each group revealed

<table><caption>**Table 3.** Brain regions activated in both groups by contrasts [Initial element>Core chunk] and [Core chunk>Initial element].</caption><thead><tr><th>Contrasts and regions</th><th>Talairach coordinates</th><th>Volume (mm³)</th><th>Maximum t</th></tr></thead><tbody><tr><td colspan="4"><strong>Initial element>Core chunk (all subjects)</strong></td></tr><tr><td>Right IPL (BA40)</td><td>41 -41 38</td><td>9173</td><td>8.48</td></tr><tr><td>Right VLPFC (BA44/BA45/BA47)/INS</td><td>40 13 14</td><td>6384</td><td>7.69</td></tr><tr><td>Right DLPFC (BA9)</td><td>37 36 33</td><td>744</td><td>5.72</td></tr><tr><td>Right SFG (BA9)</td><td>2 6 50</td><td>5619</td><td>6.72</td></tr><tr><td>Left IPL (BA40)</td><td>-31 -51 39</td><td>6930</td><td>7.18</td></tr><tr><td>Left SMA (BA6)</td><td>-28 -10 49</td><td>779</td><td>6.50</td></tr><tr><td>Left VLPFC (BA44/BA45/BA47)/INS</td><td>-33 17 12</td><td>3011</td><td>7.23</td></tr><tr><td>Left PMC (BA6)</td><td>-42 -2 29</td><td>1002</td><td>5.27</td></tr><tr><td>Left DLPFC (BA9)</td><td>-44 29 37</td><td>445</td><td>5.00</td></tr><tr><td colspan="4"><strong>Core chunk>Initial element (all subjects)</strong></td></tr><tr><td>Right MTG (BA21)</td><td>61 -11 -8</td><td>869</td><td>5.58</td></tr><tr><td>Right semi-lunar lobule (cerebellum)</td><td>31 -69 -35</td><td>1007</td><td>5.57</td></tr><tr><td>Right IFG (BA47)</td><td>35 19 -17</td><td>368</td><td>5.87</td></tr><tr><td>Right hippocampus</td><td>24 -17 -18</td><td>555</td><td>6.74</td></tr><tr><td>Right medial FG (BA10)</td><td>2 52 4</td><td>421</td><td>4.73</td></tr><tr><td>Left medial FG (BA8)</td><td>-6 50 40</td><td>2091</td><td>6.98</td></tr><tr><td>Left amygdala</td><td>-21 -8 -14</td><td>1221</td><td>6.20</td></tr><tr><td>Left MFG (BA10)</td><td>-42 46 -3</td><td>859</td><td>5.08</td></tr><tr><td>Left STG (BA38)</td><td>-45 16 -22</td><td>1791</td><td>7.76</td></tr><tr><td>Left angular gyrus (BA39)</td><td>-46 -68 35</td><td>1041</td><td>5.55</td></tr><tr><td>Left MTG (BA21)</td><td>-58 -5 -10</td><td>1944</td><td>6.15</td></tr></tbody></table>

IPL – inferior parietal lobule.

VLPFC – ventro-lateral prefrontal cortex.

INS – anterior insular cortex.

DLPFC – dorso-lateral prefrontal cortex.

SFG – superior frontal gyrus.

SMA – supplementary motor area.

PMC – premotor cortex.

MTG – middle temporal gyrus.

IFG – inferior frontal gyrus.

Medial FG – medial frontal gyrus.

STG – superior temporal gyrus.

MFG – middle frontal gyrus.

doi:10.1371/journal.pone.0103885.t003
<figure><img src="image_14.png" /><figcaption>**Figure 5.** Areas more activated during the core than initial chunk element (blue regions). The bar graphs show the percentage change in blood oxygenation level dependent (BOLD) signal in selected regions, for each group and type of trials. Blue stars indicate significant statistical differences (p<0.05 for one star and p<0.001 for two stars), whereas blue circles indicate no significant difference.
doi:10.1371/journal.pone.0103885.g005</figcaption></figure>

significant bilateral striatal activation in the Global, but not in the Incremental group.

Overall, the results of the block-based contrasts indicate that the participants in both conditions engaged the cortico-striatal network (putamen, SMA, parietal) significantly more during the sequence than the control blocks. In addition, participants exposed to the Global learning strategy during their training engaged significantly more the ventro-orbital prefrontal cortex, cerebellum, as well as the medial and lateral temporal areas. Also, these participants seemed to activate more the striatum bilaterally during sequence than during control blocks.

## Imaging results: event-related contrasts analysis

The contrast [*initial SEQ* > *body SEQ*], between brain activation during the initial element versus that during the core elements of a chunk of the target sequence, revealed that the two groups shared an extensive network of commonly activated regions both when initiating a sequential cluster, as well as when executing the elements forming the core of a chunk (Table 3). On one hand, the common network activated during chunk initiation (the orange-yellow areas in Figure 5, left column) included frontal and prefrontal areas almost exclusively (with the exception of inferior parietal lobule bilaterally). Of these, the activation in the

ventrolateral prefrontal clusters (VLPFC – encompassing areas BA44, BA45 and BA47) and anterior insular cortex (INS – BA13), bilaterally, was the most interesting as these brain regions are known to be implicated in a wide variety of language and short-term memory functions that are very much sequential in nature. A detailed analysis of the BOLD signal change in these regions also revealed a significant effect of the type of trial (initial vs. body) [i.e. $F_{1,20} = 44.74, p<0.05$ in left VLPFC/INS marked in Figure 5, left column], but no significant main effect of group [$F_{1,20} = 1.46, p>0.05$], and no significant group * trial type interaction [$F_{1,20} = 0.05, p>0.05$]. The trial type effect was significant in both groups [$F_{1,10} = 27.12, p<0.05$ for Incremental and $F_{1,10} = 19.44, p<0.05$ for Global]. In addition, we found a significant negative correlation between the difference in BOLD signal in this ROI and the difference in reaction time for the type of trial (initial vs. body), both when we considered the raw reaction time ($R^2 = 0.29, p<0.01$), as well as within-group z-score for reaction time ($R^2 = 0.47, p<0.001$). On the other hand, the common network involved in executing the chunk's core (the blue areas in Figure 5, left column) included mostly frontal and temporal areas, both medially and laterally. Among these regions, the results revealed a cluster in the left hippocampus (marked in Figure 5, left and middle columns), where we not only observed
the expected significant main effect of trial type [$F_{1,20} = 95.67, p<0.05$], but also a significant trial type and group interaction [$F_{1,20} = 30.98, p<0.05$]. The latter interaction resulted from a significantly higher difference between activation during the core elements than during the initial element of a chunk among people trained in the Global condition [$F_{1,10} = 130.37, p<0.05$] compared to those who trained in the Incremental group [$F_{1,10} = 8.10, p<0.05$]. Such effect was observed both when performing the interaction contrast for the entire brain (trial type by condition – Table 4), as well as when analyzing the imaging data separately for each group (Table 5). Importantly, this hippocampal cluster could still be observed in the Global group (Figure 5, middle column), but did not survive the statistical threshold in the Incremental group condition (Figure 5, right column). In fact, with the exception of one cluster located in semilunar lobule in cerebellum, all other clusters observed when pooling the data from both group together, could be found in the Global, but not the Incremental group when performing the same contrast separately for each group.

Previous experiments and the current models of motor sequence
posit that the striatum is involved in this process, but these claims
are usually based on block-related contrasts. In our experiment we
did not detect any activations in the striatum in response to the
event-related contrasts, performed either when comparing the two
groups or pooling their data together. However, when performing
the contrast [*initial SEQ>body SEQ*] for subjects in the
Incremental group only, we observed, as expected, a cluster in
posterior right putamen. Detailed analysis in this region revealed
that even though the interaction between experimental conditions
and trial type was not significant [$F_{1,20} = 0.30, p>0.05$], the
individuals in the Incremental group showed significantly more
activation in this region when putting sequence elements into well
integrated chunks than when executing the movements preceding
these chunks [$F_{1,10} = 17.18, p<0.05$]. Participants in the Global
condition, however, did not show such a significant difference
between trial types [$F_{1,10} = 0.37, p<0.05$]. Nevertheless, in
Global condition, the changes in BOLD signal for both types of
events in striatum (i.e. initial element and cluster core) were
positive and significantly above baseline, indicating that the
participants in this condition recruited the striatum to integrate
together the initial and the core elements of the cluster.

In conclusion, we observed that the two groups shared a
common neuronal network during chunk initiation, located almost
exclusively in frontal and prefrontal regions, and another one
during the execution of the chunk core, located mainly in parieto-
temporal lobes. Regarding the dissociation between the roles of
striatum and hippocampus in chunking, our results suggest that
the hippocampus contributes to the execution of chunks, but not to
their initiation, and that such an effect was observed to a greater
extent when people had extensive practice with the same sequence
(Global condition), then when they were exposed to it gradually
(Incremental condition). The activity in posterior putamen was
significantly higher than baseline for both the first and the
remaining elements of a chunk in the Global group (indicating the
integration of the first chunk element with the others), and only
during the execution of the chunk core for subjects in the
Incremental group.

Discussion

We manipulated the training regimen of two groups of
participants in order to test the roles of cortico- striatal circuits
and MTL in chunking of sequence elements during the late phase
of motor learning. Although both groups learned the same motor

sequence by the end of training, they differed in terms of their
reaction time and variability in motor responses. Individuals in the
Global condition, who developed a stable sequence representation,
were faster and less variable than their counterparts in the
Incremental condition. During the last four runs of the imaging
session, which took place the day after the 4 days of training were
complete, the two groups were similar in terms of motor
performance. The only observed difference between the two
groups being observed at the level of motor sequence reproduc-
ibility, which was, as expected, higher in the Global than in the
Incremental condition. When performing contrasts taking into
account the whole block of motor performance, we noted that the
two groups activated significantly more the cortico- striatal areas
during the execution of the target sequence than during the
control condition. However, performing the active sequence in the
Global group resulted in greater putaminal activity, bilaterally.

The event-related analysis allowed us to dissociate between the
cerebral activity specific to the initiation of sequence chunks
compared to that related to execution of the chunk cores. Both
groups engaged, to the same extent, a prefrontal network when
initiating sequential chunks, with a large volume of activation seen
in the ventrolateral prefrontal cortex, bilaterally. In contrast, the
neural correlates specific to the execution of well integrated chunks
were localized predominantly in the MTL, including hippocampus
and amygdala; in most of these regions, the activation was greater
in the Global than in the Incremental condition.

Behavioral results indicate that despite having different training
regimens and different ways of learning the same sequence, both
groups reached a similar level of performance during the imaging
session with some differences in the variability of performance.
Changes in reaction time during learning of the sequence and its
variability within a single training session were used in the past as
indicators of learning or movement execution strategies employed
by the subjects [7,38,40]. The scope of having different behavioral
training regimens was, in our study, to foster a stable, albeit
complex representation of the motor sequence in the Global group
and to gradually build up such a representation to develop in the
Incremental group. Both measures of reaction time and variability
were lower in the Global than in the Incremental group in the 8th
training session, hence suggesting that the latter participants had a
more stable representation of the sequence, as well as a stable
execution strategy. Interestingly, chunk's length was significantly
smaller in the Global than the Incremental group, suggesting that
a stable and complex representation of the sequence can rest on
small clusters, similar to knowing by heart a long phone number,
which is divided in smaller parts.

During the imaging session, the mean reaction times were
comparable in the two groups; the only difference being observed
only at the level of reaction time variability. This indicates that, at
this stage, subjects in the two groups were able to execute the
sequence at the same speed. Yet participants in the Global group
displayed less variability, hence a more stable performance,
suggesting that they developed a different strategy to perform
the task [38].

Another view that may explain the difference between the two
groups that arose from the specificity of behavioral training could
relate to the extent of consolidation of motor memory. Previous
studies have shown that periods of diurnal or nocturnal sleep can
contribute to the consolidation of motor sequences, expressed as
spontaneous gains between sessions in the absence of practice,
stability of the performance or resistance to interference
[6,11–14,41]. Given that subjects in the Global group were
exposed to the same sequence during their behavioral training, it is
likely that they benefited to a greater extent from memory
<table><caption>**Table 4.** Brain regions activated by contrast [Core chunk>Initial element] more in Global than in Incremental group.</caption><thead><tr><th>Contrasts and regions</th><th>Talairach coordinates</th><th>Volume (mm³)</th><th>Maximum t</th></tr></thead><tbody><tr><td>Core chunk>Init. elem. (Global>Increm.)</td><td></td><td></td><td>t(21)</td></tr><tr><td>Right STG (BA22)</td><td>57 -7 2</td><td>402</td><td>4.57</td></tr><tr><td>Right posterior insular cortex (BA13)</td><td>39 -17 16</td><td>702</td><td>4.09</td></tr><tr><td>Right IFG (BA47)</td><td>31 16 -17</td><td>201</td><td>5.17</td></tr><tr><td>Right PCC (BA31)</td><td>12 -56 21</td><td>3030</td><td>6.67</td></tr><tr><td>Right hippocampus</td><td>25 -20 -19</td><td>113</td><td>4.74</td></tr><tr><td>Right medial FG (BA10)</td><td>21 52 12</td><td>181</td><td>3.57</td></tr><tr><td>Left precuneus (BA31)</td><td>-12 -50 35</td><td>703</td><td>4.07</td></tr><tr><td>Left SMA (BA6)</td><td>-26 10 66</td><td>269</td><td>4.60</td></tr><tr><td>Left PMC (BA6)</td><td>-31 -12 68</td><td>716</td><td>4.78</td></tr><tr><td>Left posterior insular cortex (BA13)</td><td>-36 -20 13</td><td>740</td><td>5.08</td></tr><tr><td>Left MTG (BA39)</td><td>-48 -64 27</td><td>4179</td><td>5.20</td></tr><tr><td>Left postcentral gyrus (BA2)</td><td>-49 -27 56</td><td>710</td><td>4.91</td></tr><tr><td>Left postcentral gyrus (BA43)</td><td>-53 -7 14</td><td>133</td><td>3.75</td></tr><tr><td>Left STG (BA22)</td><td>-60 -23 3</td><td>154</td><td>3.69</td></tr><tr><td>Left MTG (BA21)</td><td>-62 -3 -15</td><td>438</td><td>4.89</td></tr></tbody></table>

STG – superior temporal gyrus.

IFG – inferior frontal gyrus.

PCC – posterior cingulate cortex.

Medial FG – medial frontal gyrus.

SMA – supplementary motor area.

PMC – premotor cortex.

MTG – middle temporal gyrus.

doi:10.1371/journal.pone.0103885.t004

consolidation induced by time and night sleeping periods; those in the Incremental group could only benefit from the consolidation due to the night of sleep between the last behavioral and the subsequent imaging sessions. Thus, differences in stability of the motor performance (i.e. variability) may also reflect differences in the motor skill consolidation during sleep and in the amount of practice on the same sequence.

The imaging results from the block-related analysis are in accord with previous findings showing that the putamen is involved in the self-generation of multiple novel actions especially when no clear distinction is available between competitive alternatives [42,43] (i.e. in our case, each individual response was performed with the same frequency throughout the task; thus each response had the same probability of occurring). Furthermore, dorsal striatum activity increases have been reported for category judgments in the context of greater category uncertainty [44]. In the context of the present study, the nature of the training regimen for the Incremental group makes it so that the earlier movements of the sequence are more heavily biased (practiced) than the latter ones, while they carry an equal value in the Global group, possibly explaining the greater recruitment of putamen in the latter group.

The greater recruitment of the striatum and its cortical projections in both group is concordant with the neurobiological model of motor sequence learning proposed by Doyon and collaborators [2,5,45], which predicts an increased cortico-striatal activity in the later phases of motor sequence learning. The present results show that the plastic changes in this specific network are relatively independent of the subject's learning strategy, given that they were present in both groups; this finding is novel in the motor sequence learning literature. In contrast, the Global, but not the Incremental learning strategy, led to

recruitment of additional areas in the ventro-lateral prefrontal cortex, cerebellum and temporal lobe. This additional ventro-lateral prefrontal activation seems to parallel findings in the literature on executive processes, which propose that the fronto-striatal loops are implicated in processing feedback information requiring a change in strategy [46], in planning and rule shifting when no clear clues are provided by the environment [43], and in switching from a planned response to a novel [46] or a more difficult one [47]. The greater involvement of temporal and cerebellar areas in Global than in Incremental learning group is consistent with recent findings showing that synchronization of motor performance (i.e. higher temporal reproducibility of the sequence) correlates positively with the level of cerebellar activity in later stages of motor learning [19] and that inter-regional connectivity tends to increase with practice after weeks of training, despite the fact that the regional cerebral activity may decrease over time [48]. Thus, the recruitment of additional areas outside the cortico-striatal network seems to support the development of a stable long-term representation and execution strategy of the motor sequence.

The novelty of our findings relies, however, in the event-related analysis, which permitted us to dissociate between the neural substrate mediating chunk initiation *versus* chunk execution. As stated before, the first element in a chunk has always a slower reaction time relative to the subsequent chunk elements due to the fact that transitions between clusters reflects a starting cost or higher memory load as a previous cluster is discarded and another one is loaded [20,36,37]. Here, we showed that despite differences in variability or in execution strategies between the two experimental conditions, subjects in both groups recruited the ventrolateral PFC/anterior insula, bilaterally, when executing the first element of a sequential chunk. In the neuroimaging literature
<table><caption>**Table 5.** Brain regions activated by contrast [Core chunk>Initial element] separately within each group.</caption><thead><tr><th>Contrasts and regions</th><th>Talairach coordinates</th><th>Volume (mm³)</th><th>Maximum t</th></tr></thead><tbody><tr><td colspan="4"><strong>Core chunk>Initial element (Incremental)</strong></td></tr><tr><td>Right MFG (BA11)</td><td>36 34 -11 177</td><td></td><td>t(10) 6.17</td></tr><tr><td>Right lentiform nucleus (putamen)</td><td>26 -12 7 119</td><td></td><td>8.13</td></tr><tr><td>Right pyramis (cerebellum)</td><td>25 -73 -34 122</td><td></td><td>5.69</td></tr><tr><td>Left STG (BA38)</td><td>-45 19 -24 108</td><td></td><td>6.92</td></tr><tr><td colspan="4"><strong>Core chunk>Initial element (Global)</strong></td></tr><tr><td>Right MTG (BA21)</td><td>59 -9 -6 1245</td><td></td><td>t(10) 9.27</td></tr><tr><td>Right posterior insular cortex (BA13)</td><td>37 -19 16 305</td><td></td><td>6.56</td></tr><tr><td>Right IFG (BA47)</td><td>32 17 -18 396</td><td></td><td>9.22</td></tr><tr><td>Right hippocampus</td><td>24 -20 -18 521</td><td></td><td>10.18</td></tr><tr><td>Right SMA (BA6)</td><td>21 -14 72 122</td><td></td><td>6.52</td></tr><tr><td>Right SPL (BA7)</td><td>11 -52 71 231</td><td></td><td>5.62</td></tr><tr><td>Left medial FG (BA8)</td><td>-8 48 43 1211</td><td></td><td>9.73</td></tr><tr><td>Left medial FG (BA9)</td><td>-4 59 26 280</td><td></td><td>6.78</td></tr><tr><td>Left precuneus (BA31)</td><td>-15 -46 33 321</td><td></td><td>6.63</td></tr><tr><td>Left amygdala</td><td>-19 -8 -18 382</td><td></td><td>6.44</td></tr><tr><td>Left SMA (BA6)</td><td>-29 -15 70 302</td><td></td><td>8.75</td></tr><tr><td>Left posterior insular cortex (BA13)</td><td>-38 -17 15 441</td><td></td><td>6.82</td></tr><tr><td>Left MTG (BA21)</td><td>-49 4 -20 3376</td><td></td><td>9.25</td></tr><tr><td>Left angular gyrus (BA39)</td><td>-47 -67 30 4327</td><td></td><td>12.88</td></tr><tr><td>Left IFG(BA10)</td><td>-45 44 -5 180</td><td></td><td>5.37</td></tr><tr><td>Left STG (BA38)</td><td>-44 15 -35 135</td><td></td><td>6.00</td></tr><tr><td>Left postcentral gyrus (BA2)</td><td>-47 -28 58 147</td><td></td><td>5.82</td></tr></tbody></table>

MFG – middle frontal gyrus.

STG – superior temporal gyrus.

MTG – middle temporal gyrus.

IFG – inferior frontal gyrus.

SMA – supplementary motor area.

SPL – superior parietal lobule.

Medial FG – medial frontal gyrus.

doi:10.1371/journal.pone.0103885.t005

on motor learning to date, there are only a few studies employing an event-related analysis on sequence chunking [34] and the findings support the recruitment of frontoparietal networks in chunk initiation or segmentation. In contrast, the literature on set/ task-shifting [46,49,50] and on rule/goal selection [51–53] has more such studies and they identified the VLPFC as one of the regions involved in rule shifting/selection processes. Although conjectural, it is possible that the latter activation in this area is related to an executive process related to chunk management, independently of the sequence representation. The fact that the areas identified in chunk initiation also included anterior insula is in line with recent findings which suggest that this part of the insular cortex has an important role in high-level cognitive and attentional processes [54]. In this context, it is hypothesized that anterior insula detects salient events and acts as a hub between large scale networks in order to facilitate proper allocation of attention and working memory resources. Consistent with this role, there was a significant negative correlation between the difference in BOLD signal and the difference in reaction time between the initial element and the body of the chunk, suggesting that the better or smoother the chunking process, the higher the difference in activity in this region.

A clear dissociation between the two experimental conditions was also observed in regard to the neural substrate associated with

the chunking itself, where differences were observed in: hippo-campus and other MTL areas in the Global group and in posterior putamen in Incremental condition. Similarly, the sensorimotor putamen was shown elsewhere to be significantly involved in chunk concatenation [34], consistent with findings showing that chunking is heavily dependent on dopaminergic circuits in PD patients [35]. Unlike Wymbs and colleagues [34], who did not manipulate learning strategy, we did not find chunking-related activation in putamen when we combined the two experimental conditions, but only in the Incremental condition, which is a clear indication that it depends on the training regimen. The absence of striatal finding for subjects in the Global group may seem at odds with the neurobiological model of motor sequence learning [2,5,45], which predicts increased striatal activity in the later stages of learning. However, following the argument described above regarding differential striatal activation in both groups, we propose that since all movements are equally trained in the global group, and that the chunks are indeed small, it is equally difficult to distinguish the chunk body from the element initiating the chunk, and therefore striatal requirement is similar for both conditions. On the other hand, it is possible that the first trial of a chunk is significantly differentiated from the body of the chunk in the incremental group, which leads to more putaminal activity when choosing the next movement of the sequence within the
chunk. Therefore, the fact that we did not observe a differential activation in striatum for the Global group, when comparing the chunk initiation with the chunk body itself; does not mean that the striatum is not implicated in later stages of learning. In fact, a detailed analysis of the activity in the posterior putamen, which was found to be more activated during the execution of core elements than during the initiation of a chunk in Incremental condition, showed that the activity level in the same region was above baseline (rest period) during both types of events for the participants in the Global group (Figure 5). Thus, this finding does support and even expand the role of striatum in the later stages of learning, by showing the involvement of this structure in clustering or chunking of sequence elements in both experimental conditions. In contrast, the pattern of activation in hippocampus was specific for the subjects in Global condition; only in this case, the activity was higher than baseline during chunking indicating that the development of a stable sequence representation is accompanied by a specific change in hippocampal activity. In the past, the recruitment of hippocampus during motor learning tasks was associated with the formation of higher-order associations [9], with the behavioral performance changes as learning progresses [6–8,19,55] and with overnight, but not over day gains in performance [6,17]. Our results go one step further: we report here that the activity in the hippocampus is specifically associated with chunking of sequence elements, but only when a stable sequence representation and execution strategy is in place. Of course, given the fact that our scanning session was administered at the end of 4 days of training, it is certain that the learning process in the Global group may have benefited from the consolidation due to four nights of sleep. However, this does not discount the fact that hippocampus was more active only during the execution of the elements that were already well clustered and was actually ‘deactivated’ during the initial chunk elements, thus

supporting our claim that its role is in chunking the elements once the sequence representation is stable.

Finally, our results have two major implications for the motor sequence learning literature. First, from a methodological viewpoint, our data suggest that employing an event-related design organized around chunking or clustering of sequence elements may provide insights on the functional roles of various cerebral regions that would not be otherwise revealed by classical block-based contrasts. Second, our findings could be relevant from a clinical point of view concerning the procedural learning capacity of amnesic patients with damage to the hippocampus. It has been shown in the literature that these individuals are capable of learning implicit motor skills [56,57], including motor sequences. Yet, people with amnesia have also been found to have a deficit in implicitly learning high-order associative information, similar to the second-order conditional probability which exists between elements of the sequence in the current study [58]. Our results indicate that in procedural learning of explicit sequences hippocampus is involved in the development of a stable sequence representation and execution strategy, which may be the basis for actually developing the motor expertise. The explicit memory demands in our task arguably pose a problem for people with anterograde amnesia. Still, there is past evidence that these individuals could perform well an explicit motor sequence task with no visual guidance, provided that the order of the elements in the sequence is determined by a well-known rule [59]. In this case, our data suggests that the use of an incremental learning strategy may be beneficial for these patients.

## Author Contributions

Conceived and designed the experiments: OL OM YB JD. Performed the experiments: OL TJ EB. Analyzed the data: OL OM GA TJ EB YB JD. Wrote the paper: OL OM GA JD.

## References

1. Dayan E, Cohen LG (2011) Neuroplasticity subserving motor skill learning. Neuron 72: 443–454.

2. Doyon J, Bellec P, Amsel R, Penhune V, Monchi O, et al. (2009) Contributions of the basal ganglia and functionally related brain structures to motor learning. Behav Brain Res 199: 61–75.

3. Willingham DB, Nissen MJ, Bullemer P (1989) On the development of procedural knowledge. J Exp Psychol Learn Mem Cogn 15: 1047–1060.

4. Karni A, Meyer G, Rey-Hipolito C, Jezzard P, Adams MM, et al. (1998) The acquisition of skilled motor performance: fast and slow experience-driven changes in primary motor cortex. Proc Natl Acad Sci U S A 95: 861–868.

5. Doyon J, Benali H (2005) Reorganization and plasticity in the adult brain during learning of motor skills. Curr Opin Neurobiol 15: 161–167.

6. Albouy G, Sterpenich V, Balteau E, Vandewalle G, Desseilles M, et al. (2008) Both the hippocampus and striatum are involved in consolidation of motor sequence memory. Neuron 58: 261–272.

7. Albouy G, Sterpenich V, Vandewalle G, Darsaud A, Gais S, et al. (2012) Neural correlates of performance variability during motor sequence acquisition. NeuroImage 60: 324–331.

8. Gheysen F, Van Opstal F, Roggeman C, Van Waelvelde H, Fias W (2010) Hippocampal contribution to early and later stages of implicit motor sequence learning. Exp Brain Res 202: 795–807.

9. Schendan HE, Searl MM, Melrose RJ, Stern CE (2003) An FMRI study of the role of the medial temporal lobe in implicit and explicit sequence learning. Neuron 37: 1013–1025.

10. Gomez Beldarrain M, Astorgano AG, Gonzalez AB, Garcia-Monco JC (2008) Sleep improves sequential motor learning and performance in patients with prefrontal lobe lesions. Clin Neurol Neurosurg 110: 245–252.

11. Korman M, Doyon J, Doljansky J, Carrier J, Dagan Y, et al. (2007) Daytime sleep condenses the time course of motor memory consolidation. Nat Neurosci 10: 1206–1213.

12. Doyon J, Korman M, Morin A, Dostie V, Hadj Tahar A, et al. (2009) Contribution of night and day sleep vs. simple passage of time to the consolidation of motor sequence and visuomotor adaptation learning. Exp Brain Res 195: 15–26.

13. Albouy G, Fogel S, Pottiez H, Nguyen VA, Ray L, et al. (2013) Daytime sleep enhances consolidation of the spatial but not motoric representation of motor sequence memory. PloS one 8: e52805.

14. Barakat M, Doyon J, Debas K, Vandewalle G, Morin A, et al. (2011) Fast and slow spindle involvement in the consolidation of a new motor sequence. Behavioural brain research 217: 117–121.

15. Debas K, Carrier J, Orban P, Barakat M, Lungu O, et al. (2010) Brain plasticity related to the consolidation of motor sequence learning and motor adaptation. Proceedings of the National Academy of Sciences of the United States of America 107: 17839–17844.

16. Barakat M, Carrier J, Debas K, Lungu O, Fogel S, et al. (2012) Sleep spindles predict neural and behavioral changes in motor sequence consolidation. Human brain mapping.

17. Walker MP, Stickgold R, Alsop D, Gaab N, Schlaug G (2005) Sleep-dependent motor memory plasticity in the human brain. Neuroscience 133: 911–917.

18. Albouy G, Vandewalle G, Sterpenich V, Rauchs G, Desseilles M, et al. (2012) Sleep Stabilizes Visuomotor Adaptation Memory: a functional Magnetic Resonance Imaging Study. Journal of sleep research.

19. Steele CJ, Penhune VB (2010) Specific increases within global decreases: a functional magnetic resonance imaging investigation of five days of motor sequence learning. J Neurosci 30: 8332–8341.

20. Sakai K, Kitaguchi K, Hikosaka O (2003) Chunking during human visuomotor sequence learning. Exp Brain Res 152: 229–242.

21. Boyd LA, Edwards JD, Siengsukon CS, Vidoni ED, Wessel BD, et al. (2009) Motor sequence chunking is impaired by basal ganglia stroke. Neurobiol Learn Mem 92: 35–44.

22. Graybiel AM (1998) The basal ganglia and chunking of action repertoires. Neurobiol Learn Mem 70: 119–136.

23. Levesque M, Bedard MA, Courtemanche R, Tremblay PL, Scherzer P, et al. (2007) Raclopride-induced motor consolidation impairment in primates: role of the dopamine type-2 receptor in movement chunking into integrated sequences. Exp Brain Res 182: 499–508.

24. Tremblay PL, Bedard MA, Levesque M, Chebli M, Parent M, et al. (2009) Motor sequence learning in primate: role of the D2 receptor in movement chunking during consolidation. Behav Brain Res 198: 231–239.

25. Ergorul C, Eichenbaum H (2006) Essential role of the hippocampal formation in rapid learning of higher-order sequential associations. J Neurosci 26: 4111–4117.
26. DeCoteau WE, Kesner RP (2000) A double dissociation between the rat hippocampus and medial caudoputamen in processing two forms of knowledge. Behav Neurosci 114: 1096-1108.

27. Devan BD, Goad EH, Petri HL (1996) Dissociation of hippocampal and striatal contributions to spatial navigation in the water maze. Neurobiol Learn Mem 66: 305-323.

28. Ganor-Stern D, Plonsker R, Perlman A, Tzelgov.J (2013) Are all changes equal? comparing early and late changes in sequence learning. Acta Psychol (Amst) 144: 180-189.

29. Bapi RS, Doya K, Harner AM (2000) Evidence for effector independent and dependent representations and their differential time course of acquisition during motor sequence learning. Exp Brain Res 132: 149-162.

30. Tubau E, Lopez-Moliner J (2004) Spatial interference and response control in sequence learning: the role of explicit knowledge. Psychol Res 68: 55-63.

31. Korman M, Raz N, Flash T, Karni A (2003) Multiple shifts in the representation of a motor sequence during the acquisition of skilled performance. Proc Natl Acad Sci U S A 100: 12492-12497.

32. Karni A, Meyer G, Jezzard P, Adams MM, Turner R, et al. (1995) Functional MRI evidence for adult motor cortex plasticity during motor skill learning. Nature 377: 155-158.

33. Hansen S, Tremblay L, Elliott D (2005) Part and whole practice: chunking and online control in the acquisition of a serial motor task. Res Q Exerc Sport 76: 60-66.

34. Wymbs NF, Bassett DS, Mucha PJ, Porter MA, Grafton ST (2012) Differential recruitment of the sensorimotor putamen and frontoparietal cortex during motor chunking in humans. Neuron 74: 936-946.

35. Tremblay PL, Bedard MA, Langlois D, Blanchet PJ, Lemay M, et al. (2010) Movement chunking during sequence learning is a dopamine-dependant process: a study conducted in Parkinson's disease. Experimental brain research Experimentelle Hirnforschung Experimentation cerebrale 205: 375-385.

36. Rosenbaum DA, Kenny SB, Derr MA (1983) Hierarchical control of rapid movement sequences. J Exp Psychol Hum Percept Perform 9: 86-102.

37. Bo J, Seidler RD (2009) Visuospatial working memory capacity predicts the organization of acquired explicit motor sequences. J Neurophysiol 101: 3116-3125.

38. Adi-Japha E, Karni A, Parnes A, Loewenschuss I, Vakil E (2008) A shift in task routines during the learning of a motor skill: group-averaged data may mask critical phases in the individuals' acquisition of skilled performance. J Exp Psychol Learn Mem Cogn 34: 1544-1551.

39. Penny WD, Holmes AJ (2003) Random-Effects Analysis. In: Frackowiak RSJ, Friston KJ, Frith C, Dolan R, Friston KJ, et al., editors. Human Brain Function. 2nd. ed. London: Academic Press.

40. Touron DR (2006) Are item-level strategy shifts abrupt and collective? Age differences in cognitive skill acquisition. Psychon Bull Rev 13: 781-786.

41. Fischer S, Hallschmid M, Elsner AL, Born J (2002) Sleep forms memory for finger skills. Proc Natl Acad Sci U S A 99: 11987-11991.

42. Francois-Brosseau FE, Martinu K, Strafella AP, Petrides M, Simard F, et al. (2009) Basal ganglia and frontal involvement in self-generated and externally-triggered finger movements in the dominant and non-dominant hand. Eur J Neurosci 29: 1277-1286.

43. Monchi O, Petrides M, Strafella AP, Worsley KJ, Doyon J (2006) Functional role of the basal ganglia in the planning and execution of actions. Ann Neurol 59: 257-264.

44. Daniel R, Wagner G, Koch K, Reichenbach JR, Sauer H, et al. (2011) Assessing the neural basis of uncertainty in perceptual category learning through varying levels of distortion. J Cogn Neurosci 23: 1781-1793.

45. Lehericy S, Benali H, Van de Moortele PF, Pelegrini-Issac M, Waechter T, et al. (2005) Distinct basal ganglia territories are engaged in early and advanced motor sequence learning. Proc Natl Acad Sci U S A 102: 12566-12571.

46. Monchi O, Petrides M, Petre V, Worsley K, Dagher A (2001) Wisconsin Card Sorting revisited: distinct neural circuits participating in different stages of the task identified by event-related functional magnetic resonance imaging. J Neurosci 21: 7733-7741.

47. Cameron IG, Coe BC, Watanabe M, Stroman PW, Munoz DP (2009) Role of the basal ganglia in switching a planned response. Eur J Neurosci 29: 2413-2425.

48. Ma L, Wang B, Narayana S, Hazeltine E, Chen X, et al. (2010) Changes in regional activity are accompanied with changes in inter-regional connectivity during 4 weeks motor learning. Brain Res 1318: 64-76.

49. Jamadar S, Hughes M, Fulham WR, Michie PT, Karayanidis F (2010) The spatial and temporal dynamics of anticipatory preparation and response inhibition in task-switching. Neuroimage 51: 432-449.

50. Monchi O, Petrides M, Mejia-Constain B, Strafella AP (2007) Cortical activity in Parkinson's disease during executive processing depends on striatal involvement. Brain 130: 233-244.

51. Souza MJ, Donohue SE, Bunge SA (2009) Controlled retrieval and selection of action-relevant knowledge mediated by partially overlapping regions in left ventrolateral prefrontal cortex. Neuroimage 46: 299-307.

52. Nee DE, Jonides J (2009) Common and distinct neural correlates of perceptual and memorial selection. Neuroimage 45: 963-975.

53. Donohue SE, Wendelken C, Crone EA, Bunge SA (2005) Retrieving rules for behavior from long-term memory. Neuroimage 26: 1140-1149.

54. Menon V, Uddin LQ (2010) Saliency, switching, attention and control: a network model of insula function. Brain Struct Funct 214: 655-667.

55. Fernandez-Seara MA, Aznarez-Sanado M, Mengual E, Loayza FR, Pastor MA (2009) Continuous performance of a novel motor sequence leads to highly correlated striatal and hippocampal perfusion increases. Neuroimage 47: 1797-1808.

56. Spiers HJ, Maguire EA, Burgess N (2001) Hippocampal amnesia. Neurocase 7: 357-382.

57. Hopkins RO, Waldram K, Kesner RP (2004) Sequences assessed by declarative and procedural tests of memory in amnesic patients with hippocampal damage. Neuropsychologia 42: 1877-1886.

58. Curran T (1997) Higher-order associative learning in amnesia: evidence from the serial reaction time task. Journal of cognitive neuroscience 9: 522-533.

59. Cavaco S, Anderson SW, Allen JS, Castro-Caldas A, Damasio H (2004) The scope of preserved procedural memory in amnesia. Brain: a journal of neurology 127: 1853-1867.