Elisa Tatti1*
Serena Ricci1
Aaron B. Nelson1
Dave Mathew1
Henry Chen1
Angelo Quartarone2
Chiara Cirelli3
Giulio Tononi3
Maria Felice Ghilardi1*- 1CUNY School of Medicine, The City University of New York, New York, NY, United States
- 2Department of Biomedical, Dental, Morphological and Functional Imaging Sciences, University of Messina, Messina, Italy
- 3Department of Psychiatry, University of Wisconsin-Madison, Madison, WI, United States
Beta oscillations (13.5−25 Hz) over the sensorimotor areas are characterized by a power decrease during movement execution (event-related desynchronization, ERD) and a sharp rebound after the movement end (event-related synchronization, ERS). In previous studies, we demonstrated that movement-related beta modulation depth (peak ERS-ERD) during reaching increases within 1-h practice. This increase may represent plasticity processes within the sensorimotor network. If so, beta modulation during a reaching test should be affected by previous learning activity that engages the sensorimotor system but not by learning involving other systems. We thus recorded high-density EEG activity in a group of healthy subjects performing three 45-min blocks of motor adaptation task to a visually rotated display (ROT) and in another performing three blocks of visual sequence-learning (VSEQ). Each block of either ROT or VSEQ was followed by a simple reaching test (mov) without rotation. We found that beta modulation depth increased with practice across mov tests. However, such an increase was greater in the group performing ROT over both the left and frontal areas previously involved in ROT. Importantly, beta modulation values returned to baseline values after a 90-min of either nap or quiet wake. These results show that previous practice leaves a trace in movement-related beta modulation and therefore such increases are cumulative. Furthermore, as sleep is not necessary to bring beta modulation values to baseline, they could reflect local increases of neuronal activity and decrease of energy and supplies.
Introduction
Oscillations in the beta frequency range (13.5−25 Hz) are prominent in the entire sensorimotor network and show solid dynamics: desynchronization of beta power occurs with movement initiation (event-related desynchronization, ERD), while beta power rebounds (event-related synchronization, ERS) after the movement end (Lopes da Silva and Pfurtscheller, 1999; Neuper et al., 2006). Beta oscillations should function as a gatekeeper to sensorimotor information, with beta ERD representing a release of the sensorimotor network from the inhibitory power of beta activity (Lopes da Silva and Pfurtscheller, 1999; Toma, 2002; Kilavik et al., 2013). The amplitude of beta rebound is modulated by factors such as motor learning (Boonstra et al., 2007; Tan et al., 2016) and practice (Moisello et al., 2015; Nelson et al., 2017). In particular, recent studies from our laboratory have shown that extended practice in a reaching task is associated with a significant increase of the beta ERD-ERS peak-to-peak amplitude (i.e., beta modulation depth) over parietal and frontal regions. This measure is independent from mean power changes and is mostly due to amplitude increases of the peak ERS (Moisello et al., 2015; Nelson et al., 2017; Ricci et al., 2019a; Tatti et al., 2019). Interestingly, the beta modulation increase during the task is followed by a local beta power increase during the resting EEG (Moisello et al., 2015) and returns to baseline values when tested 24 h later (Nelson et al., 2017). In line with these observations, other works have demonstrated that beta ERS amplitude can be reduced by a sensorimotor perturbation (Boonstra et al., 2007) and enhanced after successful motor adaptation (Tan et al., 2016). The practice-related increase of beta modulation may thus represent use-dependent phenomena. If this were the case, first, previous learning activity that engages the sensorimotor system should be reflected in a progressive increase of beta modulation during successive performance in a simple motor test. Importantly, such an increase should be less evident after a learning task that does not involve the sensorimotor system. Second, if beta modulation enhancement was primarily caused by increased neuronal activity or temporary decrease of energy supply, a period of rest in quiet wake should restore beta modulation to baseline values. Conversely, if it were expression of long-term plasticity-related processes, a period of sleep would be needed for its renormalization. Therefore, we analyzed the changes of movement-related beta modulation depth during reaching tests (mov) recorded after three 1-h blocks of either a visuo-motor adaptation task (ROT) or a visual sequence learning task (VSEQ). After the third block of practice, one group of subjects took a 90-min nap and another rested quietly without sleeping for the same period of time. We found that, while beta modulation depth increased progressively across mov tests, the group performing ROT displayed greater values than the VSEQ group.
Materials and Methods
Subjects and Experimental Design
Two groups of right-handed healthy subjects with normal or corrected vision were enrolled for this study. During the morning hours, 28 participants (mean age ± SD: 24.4 ± 4.0 years, 16 women) performed three blocks with a visuo-motor adaptation task (ROT) and 23 subjects (mean age ± SD: 23.3 ± 4.6 years, 12 women) three blocks of visual sequence learning task (VSEQ).
All participants had no history of sleep or medical disorders. They reported an average of 7−8 h/night sleep for at least a week before the experiment, with consistent bed and rise times, as verified by their sleep diaries. Alcohol and caffeine-containing beverages were not allowed starting the night before and throughout each experiment. Briefly, around 8 am, subjects were fitted with a 256-channel HydroCel Geodesic Sensor Net (Electrical Geodesics Inc., Eugene, OR, United States). Participants were seated in a sound-shielded room in front of a computer display. As outlined in Figure 1A, in both ROT and VSEQ sessions, a baseline assessment (mov0) was run before both experimental conditions. Then, they performed three 45-min blocks of ROT, an implicit motor learning task (see description below), or VSEQ, a visual working memory task with a declarative learning component (see description below). After each 45-min block, participants completed a mov test for a total of four mov tests. Each block lasted about 1 h. Finally, after a brief lunch, all subjects were further divided in two groups: one group (15 subjects after ROT and 12 after VSEQ) was asked to take a nap and another (13 subjects after ROT and 11 after VSEQ) to rest quietly with eyes closed listening to an audio book. After such period of time, subjects were tested again with mov (mov4) (Figure 1A). The investigation was carried out in accordance with the latest version of the Declaration of Helsinki. The local Institutional Review Board approved the study and participants signed an IRB-approved consent form.
Figure 1. (A) Experimental design. Four blocks of a simple reaching task (mov) were interspaced by three 45-min blocks of practice of either ROT, an implicit motor learning task, or VSEQ, a visual working memory task with a declarative learning component. Another mov (mov4) was performed after lunch and a 90-min interval of either nap or quiet rest. (B) Mov test. One of 24 targets (three distances, eight directions) appeared in unpredictable order every 3 s. (C–E) Performance measures related to the movement.
mov Test
The characteristics of this test have been detailed in a recent paper (Tatti et al., 2019) and are illustrated in Figure 1B.
Briefly, subjects performed reaching movements with their right hand by moving a cursor on a digitizing tablet to targets appearing on a screen. The central starting point and the cursor were always visible. One out of 24 targets (three distances, eight directions) appeared on the screen every 3 s in randomized order. Instructions were to reach the target with out-and-back movements as fast as possible, with overlapping strokes, without stops in the target circle and without corrections. In addition, subjects were asked to minimize movement time and avoid anticipation. Each mov test entailed 96 movements.
As detailed in previous publications (Ghilardi et al., 2000, 2003; Perfetti et al., 2011a), we computed several measures for each movement. In this study, we focused on: reaction time, (i.e., the time from target appearance to movement onset); movement time, (i.e., the duration of the outgoing movement); total movement time, (i.e., the duration of the out and back movement); amplitude of peak velocity of the out-going segment (Figure 1C); normalized hand path area, a measure of interjoint coordination (i.e., the ratio between the area delimited by the out and back path and the squared path length, Figure 1D); directional error at peak velocity, (i.e., the difference between the direction of target and that of the trajectory at the time of peak velocity, Figure 1E). For each set, we also computed the percentage of correct movements, i.e., movements with values of reaction time, normalized hand path area and directional error within two standard deviations of the mov0 mean. Movements with any of such measures outside two standard deviations and those rejected from EEG preprocessing were excluded from EEG analyses.
ROT Task
Similarly to tasks used in previous studies (Huber et al., 2004; Perfetti et al., 2011a), in ROT we used an array of eight targets at a fixed distance in eight directions. One of targets blackened every 1.5 s in a random, unpredictable order. Instructions were the same as in mov. Participants performed three movement blocks of 21 sets of 56 movements each (30 s inter-set interval, 1176 movements per block). Differently from mov, every two sets, unbeknownst to the subjects, the direction of the cursor on the screen was rotated relative to the direction of the hand on the tablet to a maximum of 60°, in steps of 10°, 20°, or 30°. The small incremental rotation steps were implemented to minimize awareness and the use of cognitive strategies, thus triggering implicit learning processes. Importantly, the rotation steps were smaller in the first block and greater in the last one in order to keep a similar degree of adaptation across the entire session. Indeed, mean adaptation was similar in ROT1 (mean ± SD: 71.7% ± 4.2%), ROT2 (72.7% ± 3.7%), and ROT3 (71.9% ± 3.8%), suggesting that a similar learning rate occurred across ROT blocks. Importantly, each ROT block started and ended with two and three sets, respectively, without any imposed rotation; this was both to avoid and monitor possible interference of residual directional error on the subsequent mov test. Accordingly, we found that the mean directional error of the first and last set were similar in each block (mean ± SE: ROT1: 5.19 ± 0.16 vs. 5.49 ± 0.26; ROT2: 5.65 ± 0.23 vs. 5.36 ± 0.22; ROT3: 5.40 ± 0.17 vs. 5.64 ± 0.16; F(1,27) = 0.36, p = 0.55), without significant differences across blocks (F(2,54) = 0.72, p = 0.49) and Set × Block interaction (F(2,54) = 1.68, p = 0.20). This indicates that, at the end of each ROT block, performance returned on average to baseline levels and thus, interference of residual directional errors (due to ROT after-effects) on the subsequent mov tests should be considered as minimal.
VSEQ Task
In this visual sequence learning task participants were asked to memorize several 12-element spatial sequences (Ghilardi et al., 2003, 2009; Moisello et al., 2013; Steinemann et al., 2016). Sequences appeared on the screen with targets blackening every 1.5 s. Each sequence was presented three times per set (36 target presentation/set). At the end of each set, subjects reported the sequence order verbally. The same sequence was repeated until the subject correctly reported it, and then a new sequence was presented. The average number of sets required to learn a sequence was used to assess the subjects’ learning rate in each 45-min block.
EEG Recording and Analysis
High-density (HD) EEG data were acquired using a 257-channel HydroCel Geodesic Sensor Net (Electrical Geodesic Inc.) with a Net Amp 300 amplifier (250 Hz sampling rate, online reference electrode: Cz) and Net Station software (version 5.0). Sampling frequency was 250 Hz and channel impedances were maintained below 50 kΩ. All recorded data were preprocessed using EEGLAB v13.6.5b toolbox for MATLAB (v.2016b) (Delorme and Makeig, 2004; Makeig et al., 2004). The continuous signal was first filtered using a Finite Impulse Response Filter (FIR) between 1 and 80 Hz and Notch filtered at 60 Hz.
Recordings were then segmented in 4-s epochs centered on target onset and examined to remove sporadic artifacts and channels with poor signal quality. Additionally, Independent Component Analysis (ICA) with Principal Component Analysis (PCA)-based dimension reduction (max 108 components) was applied to identify stereotypical artifacts, such as eye movements and heartbeat. Electrodes with bad signal quality were reconstructed using spherical spline interpolation, whereas those located on the cheeks and neck were removed. Re-reference to overall signal average was finally applied on the resulting 180 channels.
mov Test EEG
After preprocessing, we discarded mov test epochs corresponding to “wrong” movements (see mov test description above). After trial rejection, the average number of trials per subject for the ROT and VSEQ sessions was 69.81 (±18.25 SD) and 78.75 (±5.65 SD), respectively. Data were then time-locked to movement onset (−1 to 2.5 s). Fieldtrip-based time-frequency representations within the beta frequency range (13.5−25 Hz) were computed using Complex Morlet Wavelets (0.5 Hz bins, 10 cycles). Data were normalized by the average beta power of the entire epoch. Afterward, the beta ERS-ERD peak-to-peak difference (beta modulation depth) topography was computed for each subject on mov0 to identify the electrode with the maximum beta modulation depth and the six neighbor ones. Specifically, the peak ERD and ERS amplitude was first determined over three broad regions corresponding to the frontal, left, and right located channels; peak ERD was defined as the minimum value of beta power within an interval between 100 ms before movement onset to 950 ms after; ERS amplitude was the maximum value in the interval from 700 to 2500 ms. Those values were finally used to find the electrode with the maximum beta modulation depth (peak ERS-peak ERD) and the six neighbors (see Supplementary Figures S1, S2 for a topological representation of the channels selection for each participants). Importantly, from now on, we are going to refer to this electrodes selection as Frontal, Left, and Right Regions of Interest (ROIs).
Time-frequency analyses were carried out on the selected ROIs (1:55 Hz, 0.5 Hz bins, 3:10 wavelet cycles) and normalized by the total power of the baseline test (mov0) according to this formula: (movn – mov0)/mov0. Peak beta ERS, ERD, modulation depth magnitude, as well as the ERS and ERD peak timing values were finally computed.
Tasks EEG
Time-frequency analyses on the task EEG signal were run using the MATLAB Toolbox Fieldtrip (Oostenveld et al., 2011). Time varying spectral components were estimated convolving the signal with complex Morlet Wavelets at linearly spaced frequencies (1−55 Hz, 0.5 Hz bins) and increasing number of wavelets cycles (3:10 cycles). To determine practice-related changes during ROT1, epochs recorded during the first (F) and last (L) sets of movements (both without imposed rotation) were normalized by the total power of the first block (all) according to the following formula: (TaskF/L−Task1all)/Task1all. The same approach was used for VSEQ, where the recordings of the first and the last sequence were used.
Analysis of the Nap and Quiet Rest Periods
Using standard guidelines (Berry et al., 2017), EEG recorded during the nap and the quiet rest periods was scored for sleep stages by trained experimenters with an open source, MATLAB-based, toolbox (Mensen et al., 2016). An experienced sleep scorer (AN) confirmed the scoring. Recordings were scored in 30-s epochs as: wakefulness (W), NREM sleep stage 1 (N1), NREM sleep stage 2 (N2), and NREM sleep stage 3 (N3). REM sleep was not present in either group. A mastoid reference was used and states were determined from classical derivations from the 10 to 20 montage (F4, F3, C4, C3, P3, P4, O1, and O2). The disappearance of posterior alpha oscillations and other rhythms associated with wakefulness as well as the occurrence of slow rolling eye movements were indicative of the transition to N1. Transition to N2 was marked by K complexes and sleep spindles, while transition and maintenance of N3 was determined by the occurrence of >75 uV slow waves for more than 20% of the epoch.
Statistical Analyses
Non-parametric permutation statistics were run to identify significant practice-related changes (Last-First) in beta oscillatory activity during the first block of ROT and VSEQ task. The reference distribution was created using the Monte Carlo method with 10000 permutations. The false-alarm rate was controlled by applying cluster-correction under each permutation distribution, with a threshold of three significant channels to form a cluster and a critical alpha of 0.05 (Maris and Oostenveld, 2007).
Analyses on mov peak ERD, ERS, modulation depth, latency timings and average beta amplitude were conducted on IBM SPSS statistics v.25.
We ascertained the effect of extensive practice on mov beta modulation depth, peak ERS and ERD magnitude, timings, and average beta amplitude using a mixed-model repeated measure ANOVAs with the three ROIs (left, frontal and right) and the four morning blocks (mov0, mov1, mov2, and mov3) as within-subjects factors and with Task (ROT and VSEQ) as between-subjects factor. The same approach was used to test the effects of nap and quiet wake with three blocks (mov0, mov3, and mov4) and the left and frontal regions as within-subject factors, and with task (VSEQ and ROT) and Nap and Quiet wake groups as between-subjects factors. For all these analyses, violation of sphericity was addressed with Greenhouse-Geisser correction. All post hoc pairwise comparisons that followed significant main effects were Bonferroni corrected.
Results
Learning in ROT and VSEQ Show Different EEG Correlates
All the subjects completed the three blocks of either VSEQ or ROT and the mov tests without difficulty. The performance of the subjects during both sessions showed signs of learning. In ROT, subjects successfully adapted their movements to the imposed rotation, learning more than 70% of the imposed rotation in each ROT block (see methods). In addition, the interjoint coordination improved from ROT1 to ROT3, as shown by a decrement of the hand path area (ROT1, mean ± SD: 0.057 ± 0.013; ROT3, 0.049 ± 0.009; t(27) = 3.60; p = 0.0006), suggesting the occurrence of meta-learning or a “learning how to learn” effect. Comparison of the recordings of the last and first set of movements in ROT1, both without imposed rotation, showed a power increase in the beta range in channels located over a broad region including the left temporo-parietal and frontal areas (22.4% ± 27.8%, cluster t: 138.45, p = 0.0003).
In VSEQ, the rate of learning improved from an average of 3.094 sets/sequence (SD: 0.667) in VSEQ1 to an average of 2.864 (SD: 0.716) in VSEQ3 (two-tailed paired t-test: t(22) = 1.90; p = 0.036), suggesting that subjects improved their ability to learn sequences. Cluster-based permutation analysis in VSEQ1 (last vs. first sequence) revealed a significant increase of beta power in a cluster of electrodes over the right temporo-parietal area (11.2% ± 17.9% cluster t: 65.01, p = 0.033) (Figure 2).
Figure 2. (Top) Topographic distribution of the difference in beta oscillatory activity between the last and first set/sequence VSEQ 1 (left) and ROT 1(right). (Bottom) Cluster t-values map. Dots indicate significant clusters of electrodes (p ≤ 0.05).
In summary, in both ROT and VSEQ beta power increased during learning; however, beta power increases in VSEQ and ROT was observed in different channels, with increases over regions associated to motor practice during ROT.
Previous Motor but Not Visual Learning Impairs Motor Performance in mov
Movements in mov were mostly straight and had bell-shaped velocity profiles both after VSEQ and ROT blocks. The percentage of correct movements decreased across blocks (F(3,147) = 15.82, p < 0.00001) and differently affected the two tasks (F(1,49) = 6.92, p = 0.011; block × task interaction: F(3,1) = 6.19, p = 0.001). In fact, correct movements significantly decreased from mov0 to mov1 (t(27) = 6.07; p < 0.0001), mov2 (t(27) = 4.12; p = 0.0018), and mov3 (t(27) = 6.88; p < 0.0001) only in ROT. In VSEQ, we found a significant difference only between mov0 and mov1 (t(22) = 3.044; p = 0.036). Reaction time, movement time, peak velocity, total movement duration, hand path area, and directional error of the corrected movements were similar in the two groups and without significant changes across blocks (Figure 3 and Table 1).
Figure 3. Mean and SE of performance measures for the correct movements in mov. No significant changes were observed across blocks and between practice groups for the kinematic measures.
Table 1. Results of mixed model ANOVAs for kinematic measures.
In summary, the small but significant decrease of correct movements after the ROT blocks but not after the VSEQ blocks supports the notion that error rate increased only when the practiced task shared some characteristics with mov. Nevertheless, the kinematic features of the mov correct movements were similar in both groups and were not significantly influenced by the preceding task.
Changes During mov Are Specific to the Previous Practice and the ROIs
As described in the methods, the selection of the ROIs to compute beta ERD, ERS and modulation depth were based on the topographical maps of mov0, recorded at baseline before either VSEQ or ROT learning task.
We first focus on the changes of ERD and ERS peak timings, and then on those of amplitudes of beta modulation, ERS and ERD. Importantly, analyses were performed on the data corresponding to the correct movements only in order to avoid possible contamination resulting from the inclusion of wrong movements. The results of the main mixed-model ANOVAs and post hoc comparisons are reported in Tables 2, 3. Briefly, the timings of ERS and ERD peaks differed in the three ROIs, independently of the practiced task. On average, peak ERD over the frontal ROI occurred 20 ms later than over the other regions, whereas peak ERS over the left ROI occurred 53 ms later than over the other two ROIs. The ERD timing decreased across blocks, peaking earlier in mov3 compared to mov2 (16 ± 0.5 ms, p = 0.006) and mov0 (20 ± 0.7 ms, p = 0.022) (Figure 4 and Table 3).
Table 2. Results of mixed model ANOVAs for the beta modulation depth, peak beta ERS, ERD, and average beta amplitude, and the timing of the peak ERD and ERS.
Table 3. Bonferroni-corrected post hoc comparisons for the magnitude of beta ERD and ERS and modulation depth, peak ERS and ERD timing, and average beta power.
Figure 4. Mean and SE of beta ERD and ERS peak timings in the Frontal, Left and Right ROIs.
ANOVAs on beta modulation depth revealed effects of ROI, block, and task. Similar results were found for ERS, but not for ERD magnitude (Tables 2, 3). As in a previous work (Moisello et al., 2015), beta modulation depth increased across blocks (Tables 2, 3) and was greater over the left ROI, followed by the frontal and the right ROIs (Figures 5, 6). In both ROT and VSEQ sessions, the right ROI reached significantly lower values than the other two ROIs (Table 3). The magnitude of beta modulation was different in the two sessions: greater values were found when mov was preceded by ROT than by VSEQ (Figures 5, 6 and Table 2).
Figure 5. Mean and SE of the magnitude of beta ERD, ERS, modulation depth, and average beta (dimensionless) in the Frontal, Left and Right ROIs.
Figure 6. Average mov VSEQ (top) and mov ROT (bottom) beta modulation depth topographies for the Left, Right, and Frontal ROIs.
Separate analyses for each ROI (Table 4) revealed a greater beta modulation depth after ROT than after VSEQ for the frontal and left ROIs, but not for the right ROI. Importantly, despite mean beta amplitude increased across the four blocks, we did not find differences between ROIs and tasks (Table 2).
Table 4. Results of mixed model ANOVAs for the magnitude of beta ERD and ERS, and modulation depth for each ROI.
In summary, these results confirm that beta modulation depth increases with practice. The increase is more evident over the left and the frontal ROIs following a motor learning task, ROT, which previously engaged such areas.
Both Nap and Quiet Rest Restore Beta Modulation Depth to Baseline Values
We then ascertained whether a period of either nap or quiet wake changed performance and beta modulation values. Therefore, after the three morning blocks of practice in VSEQ and ROT, a group of subjects took a nap while another group rested quietly but without sleeping for 90 min. The nap group after VSEQ (N = 12) slept for an average of 75% (SE: 5%) of the time not differently from the group that napped after ROT (N = 15, 76% ± 4%; two-tailed t-test for independent means: t25 = 0.13, p = 0.45, (μ1−μ2) = 1.2). NREM stage N2 was evident in all the subjects (average time in N2: after VSEQ: 35% ± 4%, after ROT: 37% ± 3%; t25 = 0.33, p = 0.37, (μ1−μ2) = 1.72), while N3 was present in 10 subjects after VSEQ (average time: 23% ± 7%) and in 12 after ROT (29% ± 5%; t20 = 0.75, p = 0.23, (μ1−μ2) = 6.07). This indicates that sleep was equally consolidated and deep after both VSEQ and ROT. Importantly, no differences between VSEQ and ROT conditions were observed in both delta (1−4 Hz) and theta (4.5−8 Hz) amplitude in the ROIs showing increased beta modulation (see Supplementary Figure S3). In the quiet wake group, five subjects out of 24 reached N1 stage (VSEQ, N = 3 out of 11, average time: 9% ± 4%; ROT, N = 2 out of 13, average time: 9% ± 2%) and two of them (one in VSEQ and the other in ROT) reached N2 stage but only for a short period of time (both subjects: 10%).
We first compared the number of correct movements in mov4, recorded after the 90-min interval, with mov0, recorded in the morning at baseline, and mov3, recorded at the end of the morning blocks, with a mixed model ANOVA (Block as repeated-measure factor, Task, and Group as between-subjects factors) followed by two one-way ANOVAs to isolate the effects for ROT and VSEQ (Table 5). Briefly, such analyses revealed that mov performance that was degraded in mov3 after ROT (Figure 3) improved in mov4 reaching baseline levels only after a nap but not after quiet wake (Figure 7, ROT awake vs. nap, mov0: mean difference = 2.01 ± 1.42, p = 0.165; mov3: mean difference = −6.18, p = 0.151; mov4: mean difference = −14.94 ± 4.18, p = 0.001). The increase of correct movements from mov3 to mov4 was positively correlated with delta power in the frontal ROI during the nap (N = 11, N2: r = 0.629, p = 0.038, 95% CI [0.047, 0.892]; N = 8, N3: r = 0.865, p = 0.005, 95% CI [0.411, 0.975]). Performance in mov4 after VSEQ (Figure 7) remained at the same high levels of the previous testing times (Figure 3) both after nap and quiet wake.
Table 5. (Left) Results of mixed model ANOVAs for the % of correct movements in mov0, mov3, and mov4 as Blocks, ROT vs. VSEQ as Task, and Nap vs. Awake as Group. (Right) Results of mixed model ANOVA for the ROT Task condition only.
Figure 7. (Left) Mean and SE of beta modulation depth (dimensionless) over the Left ROI after 90 min of either Nap or Quiet wake (Blocks: mov0, mov3, and mov4), in both ROT and VSEQ groups. (Right) Number of correct mov movements (%) during the baseline (mov0), the last mov of the morning (mov3), and after 90 min of either Nap or Quiet Wake (mov4).
We finally compared the magnitude of beta modulation in mov4, mov0 and mov3. The results of the ANOVA (Table 6) revealed only a significant effect of block and a borderline effect of task. Importantly, post hoc tests showed that, following the morning increase (mov3 vs. mov0: mean difference = 1.01 ± 0.22, p < 0.0001), beta modulation depth decreased after the 90-min interval (mov4 vs. mov3: mean difference = 0.91 ± 0.09, p < 0.00001), reaching baseline values of mov0 (mov4 vs. mov0: mean difference = 0.10 ± 0.19, p = 1, Figure 7). We obtained similar results for the magnitude of ERD, ERS and mean beta power (Table 6). Interestingly, the decrease of beta modulation, ERS and ERD values in the two subjects of the awake group that reached N2 was smaller than the average of the awake group. Indeed, if sleep contributed to decrease beta modulation values, we would have expected greater decreases in these two subjects.
Table 6. Results of mixed model ANOVAs for the changes in the magnitude of beta ERD and ERS and modulation depth after 90 min of either Nap or Quiet wake (Blocks: mov0, mov3, mov4; Task: ROT, VSEQ; Group: Nap, Quiet wake).
These results suggest that, while sleep is necessary to restore performance, a period of quiet rest without sleep is enough to reestablish beta modulation depth to baseline levels.
Discussion
The present study shows that increases of movement-related beta modulation depth are affected by former practice: in fact, they were greater after a visuo-motor learning task (ROT) than after a visual sequence-learning task (VSEQ). This task difference was observed over the left and frontal regions, areas that showed sustained increase of beta power during the ROT task. Crucially, we also found that beta modulation depth returned to the morning baseline values after a period of either quiet wake or sleep. Altogether, these findings suggest that former practice influences the magnitude of beta modulation increase, as a sort of performance “signature” that, however, is short lasting and does not require sleep to fade away. This implies that the practice-dependent increases may reflect increased neuronal activity or temporary decrease of energy supply, rather than full-fledged plasticity phenomena per se.
Regional Differences of Beta ERD/ERS Dynamics
The increases of beta modulation depth were different in the three ROIs. Beta oscillations in the motor, somatosensory, and frontal regions likely reflect different aspects of motor planning and execution that involve sensory, motor, and cognitive processes. Indeed, beta desynchronization during movement preparation and execution (ERD) and its rebound after movement completion (ERS) are predominant in the sensorimotor region contralateral to the moving effector, but are also observed in the ipsilateral and frontal regions (Rektor et al., 2006; Zaepffel et al., 2013; Meziane et al., 2015; Moisello et al., 2015; Ricci et al., 2019b) as also shown in the present work. However, our results further show that, while ERD magnitude was similar across ROIs, beta ERS magnitude was greater over the left and frontal areas and significantly lower over the right ROI, where practice-related increase was minimal. The discrepant results between ROIs suggest that beta rebound is a phenomenon that mostly involves the contralateral sensorimotor and frontal areas and less so the ipsilateral area. This finding is in agreement with other results showing that ipsilateral sensorimotor activity is linked mainly to movement selection and planning (Rao et al., 1993; Haaland et al., 2004) rather than to the processes occurring after movement cessation. Also, in line with past studies (Ohara et al., 2000; Szurhaj et al., 2003; Parkes et al., 2006; Wilson et al., 2010; Moisello et al., 2015; Heinrichs-Graham et al., 2017), the present results highlight the contribution of the frontal region not just to motor planning (Tombini et al., 2009; Perfetti et al., 2011b) but also to feedforward processes that occur after the movement; these are mostly reflected in the ERS magnitude, and are essential for updating internal models and learning. Indeed, beta ERS has been associated with many frontal functions, such as the maintenance of sensorimotor and cognitive sets (Gilbertson et al., 2005; Pogosyan et al., 2009; Engel and Fries, 2010; Tan et al., 2014a, b, 2016), the processing of sensory reafference (Cassim et al., 2001; Alegre et al., 2002), top-down executive control (Buschman and Miller, 2007; Spitzer and Haegens, 2017), as well as visuomotor attention (Classen et al., 1998; Kilavik et al., 2013).
Increases of Beta Modulation Depth Are Affected by Previous Practice
As in our previous work (Nelson et al., 2017; Ricci et al., 2019a, b), we found that beta modulation depth increased with practice over the three ROIs. Additionally, such increases were present, although to a lesser extent, after a visual learning task, a finding likely due to a carry-over effect from a mov test to the following one. This is in agreement with previous work showing that increased beta modulation can occur within a set with less than the 96 movements of each mov test (Nelson et al., 2017). The present results further demonstrate that such increases may be carried over to another test almost 1 h later, i.e., the time between two successive mov tests. The increases were mostly driven by beta ERS as discussed above, thus confirming that practice mainly affects the magnitude of the post-movement beta rebound.
A novel result is that beta modulation depth differently evolved after the two tasks. Starting from comparable baseline values, the magnitude beta modulation over the frontal and left parietal regions increased more when the test was preceded by the visuo-motor learning task (ROT). Indeed, while beta increases during the VSEQ task were confined to the right region, beta power during the ROT task increased locally over frontal and left parietal channels (Figure 3). The patterns of beta increase during the tasks indicate the occurrence of a cumulative effect of motor practice during the ROT task to beta modulation depth on successive mov tests. Therefore, the increases of beta modulation depth during mov may be the reflection of local neuronal use-dependent phenomena related to the development of the early phases of long-term potentiation. Support to this notion comes from studies in humans showing that intermittent theta burst stimulation produced an increase of cortical excitability (measured using MEPs) together with an increase of the amplitude of post-movement beta synchronization (Hsu et al., 2011). Moreover, since beta ERD likely reflects activation of the motor cortex and depression of the sensory cortex and beta ERS the reactivation of the sensory network (Lee and Schmit, 2018), the increase of beta modulation may express functional neural changes due to induction of long term potentiation by a repetitive pattern of activation and inactivation of the sensory and motor areas. Thus, the increases of beta modulation depth across blocks could represent the local and progressive saturation of the capacity for plasticity within the cortical areas involved in movement planning and execution (Moisello et al., 2015; Nelson et al., 2017).
A Period of Quiet Rest Returns Beta Modulation Depth to Baseline Values
Previous studies using a TMS-based repetitive pattern of sensorimotor activation have shown that the amplitude of beta modulation may express functional neural changes linked to the induction of long-term plasticity (LTP) processes (Hsu et al., 2011). These results suggest that our practice-dependent increases of beta modulation depth could be expression of plasticity-related phenomena. The early phases of LTP are short lasting (minutes to hours) and can be induced in brain slice preparations with a weak protocol (Huang, 1998). Conversely, the achievement of later phases needs a stronger induction protocol, the synthesis of new protein (Muller et al., 2002) and the decay of LTP early phases (Xiao et al., 1996; Villarreal et al., 2002) with the final result of increase of synaptic strengthening. All LTP phases require great energy availability and demand for delivery of cellular supplies. There is now strong evidence that only sleep can return the system to baseline levels in terms of net synaptic strength and cellular homeostasis while enhancing skill formation, thus overcoming the cellular consequences of full-fledged LTP (Tononi and Cirelli, 2014). However, sleep was not required to restore movement-related beta modulation depth to baseline values. It is therefore tempting to speculate that the increases of beta modulation depth may parallel some mechanisms involved in the early induction of plasticity-related phenomena, as previously indicated (Nelson et al., 2017). Alternatively but not exclusively, such increases might signal transient depletion of energy resources and cellular supplies needed to maintain LTP processes, a situation that can be restored to baseline levels by a period of rest and do not necessarily require sleep (Tononi and Cirelli, 2014).
Interestingly, as noted above, in the morning VSEQ practice we did not find reduction of beta modulation depth between 1-h apart mov tests. It is thus possible that one hour may not be sufficient to restore completely neural resources; alternatively, it is possible that the addition of activity, even when involving other brain areas, might have delayed the recovery process. In any event, while the decrease of beta modulation depth occurred even after a quiet wake period, sleep was necessary to improve the rate of correct movements, as shown by the comparison of the group that napped and the one that quietly rested after the ROT task. The discrepancy between performance and beta modulation suggests that the changes of beta modulation depth per se do not purely reflect learning-related changes. Therefore, as discussed above, the increase of movement-related beta modulation during the practice may not parallel the induction of full-scale long-term potentiation with all its benefits on performance but may signal greater neural activity and reduced energy and supplies.
Conclusion
In conclusion, these results confirmed our previous findings by showing that practice of reaching movements induces beta modulation depth increases. More importantly, they show that the magnitude of the increase is amplified by previous practice that involves similar brain areas and that a period of quiet wake without sleep can restore beta modulation to baseline values.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by the CUNY, Institutional Review Board (IRB). The patients/participants provided their written informed consent to participate in this study.
Author Contributions
MG, CC, and GT conceived and designed the study. ET, SR, and AN contributed to data acquisition. ET, MG, DM, HC, SR, and AN carried out the data analyses. ET, MG, AQ, CC, and GT interpreted the results and revised the manuscript. ET, MG, and AQ drafted the manuscript. All the authors agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.
Funding
This work was supported by NIH P01 NS083514 (GT, CC, and MG) and by DOD W81XWH-19-1-0810 (AQ and MG). Kinematic data were collected with custom-designed software, MotorTaskManager, produced by the E.T.T. s.r.l.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank Martina Bossini Baroggi for the implementation of the kinematic analysis program (Marky 1.0) that was used to analyze kinematic data. We also thank Shaina George, Cameron Chan, and Francesco Della Gatta for helping in analysis of kinematic data.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnsys.2020.00061/full#supplementary-material
FIGURE S1 | Topographies showing the selected channels for the Left, Frontal, and Right ROIs for each subject in the ROT group.
FIGURE S2 | Topographies showing the selected channels for the Left, Frontal, and Right ROIs for each subject in the VSEQ group.
FIGURE S3 | Delta and theta topographies during N2 and N3 stages in the ROT and VSEQ groups.
TABLE S1 | Results of mixed model ANOVAs for the delta and theta amplitudes in the selected ROIs during N2 and N3 stages in the ROT and VSEQ Nap groups.
References
Alegre, M., Labarga, A., Gurtubay, I. G., Iriarte, J., Malanda, A., and Artieda, J. (2002). Beta electroencephalograph changes during passive movements: sensory afferences contribute to beta event-related desynchronization in humans. Neurosci. Lett. 331, 29–32. doi: 10.1016/S0304-3940(02)00825-X
Berry, R. B., Brooks, R., Gamaldo, C., Harding, S. M., Lloyd, R. M., Quan, S. F., et al. (2017). AASM Scoring Manual Updates for 2017 (Version 2.4). J. Clin. Sleep Med. 13, 665–666. doi: 10.5664/jcsm.6576
Boonstra, T. W., Daffertshofer, A., Breakspear, M., and Beek, P. J. (2007). Multivariate time-frequency analysis of electromagnetic brain activity during bimanual motor learning. Neuroimage 36, 370–377. doi: 10.1016/j.neuroimage.2007.03.012
Buschman, T. J., and Miller, E. K. (2007). Top-down versus bottom-up control of attention in the prefrontal and posterior parietal cortices. Science 315, 1860–1864. doi: 10.1126/science.1138071
Cassim, F., Monaca, C., Szurhaj, W., Bourriez, J. L., Defebvre, L., Derambure, P., et al. (2001). Does post-movement beta synchronization reflect an idling motor cortex? Neuroreport 12, 3859–3863.
Classen, J., Gerloff, C., Honda, M., and Hallett, M. (1998). Integrative visuomotor behavior is associated with interregionally coherent oscillations in the human brain. J. Neurophysiol. 79, 1567–1573.
Delorme, A., and Makeig, S. (2004). EEGLAB: an open source toolbox for analysis of single-trial EEG dynamics including independent component analysis. J. Neurosci. Methods 134, 9–21.
Engel, A. K., and Fries, P. (2010). Beta-band oscillations-signalling the status quo? Curr. Opin. Neurobiol. 20, 156–165. doi: 10.1016/j.conb.2010.02.015
Ghilardi, M., Ghez, C., Dhawan, V., Moeller, J., Mentis, M., Nakamura, T., et al. (2000). Patterns of regional brain activation associated with different forms of motor learning. Brain Res. 871, 127–145.
Ghilardi, M. F., Eidelberg, D., Silvestri, G., and Ghez, C. (2003). The differential effect of PD and normal aging on early explicit sequence learning. Neurology 60, 1313–1319.
Ghilardi, M. F., Moisello, C., Silvestri, G., Ghez, C., and Krakauer, J. W. (2009). Learning of a sequential motor skill comprises explicit and implicit components that consolidate differently. J. Neurophysiol. 101, 2218–2229. doi: 10.1152/jn.01138.2007
Gilbertson, T., Lalo, E., Doyle, L., Di Lazzaro, V., Cioni, B., and Brown, P. (2005). Existing motor state is favored at the expense of new movement during 13-35 Hz oscillatory synchrony in the human corticospinal system. J. Neurosci. 25, 7771–7779. doi: 10.1523/JNEUROSCI.1762-05.2005
Haaland, K. Y., Elsinger, C. L., Mayer, A. R., Durgerian, S., and Rao, S. M. (2004). Motor sequence complexity and performing hand produce differential patterns of hemispheric lateralization. J. Cogn. Neurosci. 16, 621–636. doi: 10.1162/089892904323057344
Heinrichs-Graham, E., Kurz, M. J., Gehringer, J. E., and Wilson, T. W. (2017). The functional role of post-movement beta oscillations in motor termination. Brain Struct. Funct. 222, 3075–3086. doi: 10.1007/s00429-017-1387-1
Hsu, Y. F., Liao, K. K., Lee, P. L., Tsai, Y. A., Yeh, C. L., Lai, K. L., et al. (2011). Intermittent theta burst stimulation over primary motor cortex enhances movement-related beta synchronisation. Clin. Neurophysiol. 122, 2260–2267. doi: 10.1016/j.clinph.2011.03.027
Huang, E. P. (1998). Synaptic plasticity: going through phases with LTP. Curr. Biol. 8, 350–352. doi: 10.1016/s0960-9822(98)70219-2
Huber, R., Ghilardi, M. F., Massimini, M., and Tononi, G. (2004). Local sleep and learning. Nature 430, 78–81. doi: 10.1038/nature02663
Kilavik, B. E., Zaepffel, M., Brovelli, A., MacKay, W. A., and Riehle, A. (2013). The ups and downs of beta oscillations in sensorimotor cortex. Exp. Neurol. 245, 15–26. doi: 10.1016/j.expneurol.2012.09.014
Lee, J. J., and Schmit, B. D. (2018). Effect of sensory attenuation on cortical movement-related oscillations. J. Neurophysiol. 119, 971–978. doi: 10.1152/jn.00171.2017
Lopes da Silva, F. H., and Pfurtscheller, G. (1999). Event-related EEG/MEG synchronization and desynchronization: basic principles. Clin. Neurophysiol. 110, 1842–1857. doi: 10.1016/S1388-2457(99)00141-8
Makeig, S., Debener, S., Onton, J., and Delorme, A. (2004). Mining event-related brain dynamics. Trends Cogn. Sci. 8, 204–210.
Maris, E., and Oostenveld, R. (2007). Nonparametric statistical testing of EEG-and MEG-data. J. Neurosci. Methods 164, 177–190.
Mensen, A., Riedner, B., and Tononi, G. (2016). Optimizing detection and analysis of slow waves in sleep EEG. J. Neurosci. Methods 274, 1–12. doi: 10.1016/j.jneumeth.2016.09.006
Meziane, H. B., Moisello, C., Perfetti, B., Kvint, S., Isaias, I. U., Quartarone, A., et al. (2015). Movement preparation and bilateral modulation of beta activity in aging and Parkinson’s disease. PLoS One 10:e0114817. doi: 10.1371/journal.pone.0114817
Moisello, C., Blanco, D., Lin, J., Panday, P., Kelly, S. P., Quartarone, A., et al. (2015). Practice changes beta power at rest and its modulation during movement in healthy subjects but not in patients with Parkinson’s disease. Brain Behav. 5:e00374. doi: 10.1002/brb3.374
Moisello, C., Meziane, H. B., Kelly, S., Perfetti, B., Kvint, S., Voutsinas, N., et al. (2013). Neural activations during visual sequence learning leave a trace in post-training spontaneous EEG. PLoS One 8:e65882. doi: 10.1371/journal.pone.0065882
Muller, D., Nikonenko, I., Jourdain, P., and Alberi, S. (2002). LTP, memory and structural plasticity. Curr. Mol. Med. 2, 605–611. doi: 10.2174/1566524023362041
Nelson, A. B., Moisello, C., Lin, J., Panday, P., Ricci, S., Canessa, A., et al. (2017). Beta oscillatory changes and retention of motor skills during practice in healthy subjects and in patients with Parkinson’s disease. Front. Hum. Neurosci. 11:104. doi: 10.3389/fnhum.2017.00104
Neuper, C., Wortz, M., and Pfurtscheller, G. (2006). ERD/ERS patterns reflecting sensorimotor activation and deactivation. Prog. Brain Res. 159, 211–222. doi: 10.1016/S0079-6123(06)59014-4
Ohara, S., Ikeda, A., Kunieda, T., Yazawa, S., Baba, K., Nagamine, T., et al. (2000). Movement-related change of electrocorticographic activity in human supplementary motor area proper. Brain 123(Pt 6), 1203–1215. doi: 10.1093/brain/123.6.1203
Oostenveld, R., Fries, P., Maris, E., and Schoffelen, J.-M. (2011). FieldTrip: open source software for advanced analysis of MEG, EEG, and invasive electrophysiological data. Comput. Intell. Neurosci. 2011:156869. doi: 10.1155/2011/156869
Parkes, L. M., Bastiaansen, M. C. M., and Norris, D. G. (2006). Combining EEG and fMRI to investigate the post-movement beta rebound. Neuroimage 29, 685–696. doi: 10.1016/j.neuroimage.2005.08.018
Perfetti, B., Moisello, C., Landsness, E. C., Kvint, S., Lanzafame, S., Onofrj, M., et al. (2011a). Modulation of gamma and theta spectral amplitude and phase synchronization is associated with the development of visuo-motor learning. J. Neurosci. 31, 14810–14819. doi: 10.1523/JNEUROSCI.1319-11.2011
Perfetti, B., Moisello, C., Landsness, E. C., Kvint, S., Pruski, A., Onofrj, M., et al. (2011b). Temporal evolution of oscillatory activity predicts performance in a choice-reaction time reaching task. J. Neurophysiol. 105, 18–27. doi: 10.1152/jn.00778.2010
Pogosyan, A., Gaynor, L. D., Eusebio, A., and Brown, P. (2009). Boosting cortical activity at Beta-band frequencies slows movement in humans. Curr. Biol. 19, 1637–1641. doi: 10.1016/j.cub.2009.07.074
Rao, S. M., Binder, J. R., Bandettini, P. A., Hammeke, T. A., Yetkin, F. Z., Jesmanowicz, A., et al. (1993). Functional magnetic resonance imaging of complex human movements. Neurology 43, 2311–2311. doi: 10.1212/WNL.43.11.2311
Rektor, I., Sochurkova, D., and Bockova, M. (2006). Intracerebral ERD/ERS in voluntary movement and in cognitive visuomotor task. Prog. Brain Res. 159, 311–330. doi: 10.1016/S0079-6123(06)59021-1
Ricci, S., Mehraram, R., Tatti, E., Nelson, A. B., Bossini-Baroggi, M., Panday, P., et al. (2019a). Aging does not affect beta modulation during reaching movements. Neural Plast. 2019:1619290. doi: 10.1155/2019/1619290
Ricci, S., Tatti, E., Mehraram, R., Panday, P., and Ghilardi, M. F. (2019b). Beta band frequency differences between motor and frontal cortices in reaching movements. IEEE Int. Conf. Rehabil. Robot. 2019, 1254–1259. doi: 10.1109/ICORR.2019.8779373
Spitzer, B., and Haegens, S. (2017). Beyond the status quo: a role for beta oscillations in endogenous content (RE)activation. eNeuro 4:ENEURO.0170-17.2017. doi: 10.1523/ENEURO.0170-17.2017
Steinemann, N. A., Moisello, C., Ghilardi, M. F., and Kelly, S. P. (2016). Tracking neural correlates of successful learning over repeated sequence observations. Neuroimage 137, 152–164. doi: 10.1016/j.neuroimage.2016.05.001
Szurhaj, W., Derambure, P., Labyt, E., Cassim, F., Bourriez, J.-L., Isnard, J., et al. (2003). Basic mechanisms of central rhythms reactivity to preparation and execution of a voluntary movement: a stereoelectroencephalographic study. Clin. Neurophysiol. 114, 107–119. doi: 10.1016/s1388-2457(02)00333-4
Tan, H., Jenkinson, N., and Brown, P. (2014a). Dynamic neural correlates of motor error monitoring and adaptation during trial-to-trial learning. J. Neurosci. 34, 5678–5688. doi: 10.1523/JNEUROSCI.4739-13.2014
Tan, H., Wade, C., and Brown, P. (2016). Post-Movement Beta Activity in Sensorimotor Cortex Indexes Confidence in the Estimations from Internal Models. J. Neurosci. 36, 1516–1528. doi: 10.1523/JNEUROSCI.3204-15.2016
Tan, H., Zavala, B., Pogosyan, A., Ashkan, K., Zrinzo, L., Foltynie, T., et al. (2014b). Human subthalamic nucleus in movement error detection and its evaluation during visuomotor adaptation. J. Neurosci. 34, 16744–16754. doi: 10.1523/JNEUROSCI.3414-14.2014
Tatti, E., Ricci, S., Mehraram, R., Lin, N., George, S., Nelson, A. B., et al. (2019). Beta modulation depth is not linked to movement features. Front. Behav. Neurosci. 13:49. doi: 10.3389/fnbeh.2019.00049
Toma, K. (2002). Movement rate effect on activation and functional coupling of motor cortical areas. J. Neurophysiol. 88, 3377–3385. doi: 10.1152/jn.00281.2002
Tombini, M., Zappasodi, F., Zollo, L., Pellegrino, G., Cavallo, G., Tecchio, F., et al. (2009). Brain activity preceding a 2D manual catching task. Neuroimage 47, 1735–1746. doi: 10.1016/j.neuroimage.2009.04.046
Tononi, G., and Cirelli, C. (2014). Sleep and the price of plasticity: from synaptic and cellular homeostasis to memory consolidation and integration. Neuron 81, 12–34. doi: 10.1016/j.neuron.2013.12.025
Villarreal, D. M., Do, V., Haddad, E., and Derrick, B. E. (2002). NMDA receptor antagonists sustain LTP and spatial memory: active processes mediate LTP decay. Nat. Neurosci. 5, 48–52. doi: 10.1038/nn776
Wilson, T. W., Slason, E., Asherin, R., Kronberg, E., Reite, M. L., Teale, P. D., et al. (2010). An extended motor network generates beta and gamma oscillatory perturbations during development. Brain Cogn. 73, 75–84. doi: 10.1016/j.bandc.2010.03.001
Xiao, M. Y., Niu, Y. P., and Wigstrom, H. (1996). Activity-dependent decay of early LTP revealed by dual EPSP recording in hippocampal slices from young rats. Eur. J. Neurosci. 8, 1916–1923. doi: 10.1111/j.1460-9568.1996.tb01335.x
Keywords: EEG, beta, oscillations, ERD/ERS, motor control
Citation: Tatti E, Ricci S, Nelson AB, Mathew D, Chen H, Quartarone A, Cirelli C, Tononi G and Ghilardi MF (2020) Prior Practice Affects Movement-Related Beta Modulation and Quiet Wake Restores It to Baseline. Front. Syst. Neurosci. 14:61. doi: 10.3389/fnsys.2020.00061
Received: 11 June 2020; Accepted: 24 July 2020;
Published: 18 August 2020.
Edited by:
Bernat Kocsis, Harvard Medical School, United StatesReviewed by:
Vladyslav Vyazovskiy, University of Oxford, United KingdomJose L. Cantero, Universidad Pablo de Olavide, Spain
Copyright © 2020 Tatti, Ricci, Nelson, Mathew, Chen, Quartarone, Cirelli, Tononi and Ghilardi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Elisa Tatti, elisatatti@msn.com; Maria Felice Ghilardi, lice.mg79@gmail.com