Articles in PresS. J Neurophysiol (November 4, 2015). doi:10.1152/jn.00497.2015
1
Changes in cortical activity measured with EEG during a high intensity cycling exercise.
2
Hendrik Enders1, Filomeno Cortese2, Christian Maurer3, Jennifer Baltich1, Andrea B.
3
Protzner2,4, Benno M. Nigg1
4 5 6 7 8 9
1
Human Performance Laboratory, Faculty of Kinesiology, University of Calgary, Calgary, AB,
Canada 2
Hotchkiss Brain Institute, University of Calgary, Calgary, AB, Canada
3
Red Bull Diagnostic and Training Center, Thalgau, Austria
4
Department of Psychology, University of Calgary, Calgary, AB, Canada
Running Title: Cortical activity during cycling
10
Corresponding author and address:
11
Hendrik Enders, e-mail:
[email protected]
12
Human Performance Laboratory, Faculty of Kinesiology, University of Calgary
13
2500 University Drive NW, Calgary, Alberta, Canada, T2N 1N4
14
Phone: (403) 220-2413
15
Fax: (403) 284-3553
16
Research article submitted to Journal of Neurophysiology
17
Copyright © 2015 by the American Physiological Society.
18 19
Abstract
20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36
This study investigated the effects of a high intensity cycling exercise on changes in spectral and temporal aspects of electroencephalography (EEG) measured from ten experienced cyclists. Cyclists performed a maximum aerobic power test on the first testing day followed by a time-toexhaustion trial at 85% of their maximum power output on two subsequent days that were separated by approximately 48 hours. EEG was recorded using a 64 channel system at 500 Hz. Independent component analysis parsed the EEG scalp data into maximally independent components (IC). An equivalent current dipole model was calculated for each IC and results were clustered across subjects. A time-frequency analysis of the identified electrocortical clusters was performed to investigate the magnitude and timing of event-related spectral perturbations. Significant changes (p < 0.05) in electrocortical activity were found in frontal, supplementary motor and parietal areas of the cortex. Overall, there was a significant increase in EEG power as fatigue developed throughout the exercise. The strongest increase was found in the frontal area of the cortex. The timing of event-related desynchronization (ERD) within the supplementary motor area corresponds with the onset of force production and the transition from flexion to extension in the pedalling cycle. The results indicate an involvement of the cerebral cortex during the pedalling task that most likely involves executive control function as well as motor planning and execution.
37
Keywords: Electroencephalography, Locomotion, Motor control, Fatigue
38 39
40 41
Introduction Coordinated motor skills require the integration of information from peripheral sensors,
42
spinal locomotor networks as well as supraspinal commands (Grillner et al., 2008). However,
43
only in the last decade, researchers have begun to investigate the function of the human brain
44
with respect to voluntary, mostly rhythmic locomotor tasks in real time. Owing the dynamic
45
nature of human locomotion, a primary requirement to quantify cortical dynamics during
46
locomotion is a time resolution that allows the recording of modulations in cortical activity
47
within one cycle of a repeating movement sequence. Classical neuroimaging tools such as
48
functional magnetic resonance imaging (fMRI) or positron emission tomography (PET) are not
49
well suited to resolve such rapid modulations of activity. Contrary, electroencephalography
50
(EEG) is a mobile brain imaging technique with more-than-sufficient time resolution to capture
51
cortical dynamics during human locomotion. While historically, EEG has been considered too
52
prone to noise and artifacts for dynamic recordings, recent advancements in hardware and
53
software have opened the possibility for mobile brain imaging using EEG (Makeig et al., 2009;
54
Gramann et al., 2011).
55
In the last decade, researchers have begun to study time-dependent cortical dynamics in
56
locomotor and exercise-related situations with the majority of studies focusing on upright human
57
walking and slow jogging. It was proposed that the low frequency component (0.1 – 2 Hz) of
58
EEG recordings can be used to infer the overall time-dependent, linear and angular kinematics of
59
the ankle, knee and hip joints during walking (Presacco et al., 2011). Another study has shown
60
that specific brain regions such as the primary motor, the parietal, and frontal cortices exhibit
61
distinct time-frequency patterns that are closely linked to the overall movement kinematics of the
62
gait cycle (Gwin et al., 2011). Similarly, it was shown that task-dependent EEG can be extracted
63
during slow to moderate walking speeds (De Sanctis et al., 2012) and that cortical modulations
64
seen in EEG differ between active and passive walking conditions (Wagner et al., 2012). More
65
recently, Wagner and colleagues showed that providing interactive feedback during walking
66
enhances premotor and parietal areas in the brain that are thought to play an important role for
67
motor planning and intention (Wagner et al., 2014). These studies demonstrate that cortical
68
dynamics respond to changes in task parameters indicating that EEG can be used to study neural
69
control during human locomotor tasks. The typical frequency bands that are investigated in
70
movement related tasks are theta (θ: 4 – 8 Hz), alpha (α: 8 – 12 Hz), beta (β:12 – 30 Hz) and low
71
gamma (γ: 30 – 50 Hz) frequency bands as they are believed to reflect different aspects with
72
respect to motor planning, execution and control (Pfurtscheller et al., 2000; Taniguchi et al.,
73
2000; Jurkiewicz et al., 2006; Waldert et al., 2008; Gwin et al., 2011). This was supported by a
74
recent study that found sustained suppression of cortical oscillations in the α and β band while
75
cortical activity in the low γ band seemed to be modulated during the gait cycle phases (Seeber et
76
al., 2014) suggesting that these responses characterize different elements of human walking. A
77
follow up study found that low and high γ amplitudes are modulated throughout the gait cycle,
78
however, their amplitude envelopes were negatively correlated throughout the movement
79
(Seeber et al., 2015). The majority of the mentioned studies have been conducted during human
80
walking at low intensities. The aim of this study was to apply similar analysis techniques to high
81
effort locomotor tasks. Cycling avoids the high impact forces that are typically observed during
82
high effort running and has often been used in studies investigating neuromuscular strategies due
83
to the well-controlled movement and repeatability of the task (Enders et al., 2013). In fact, it has
84
previously been used for EEG measurements during high intensity cycling and evidence was
85
provided for increased communication between mid/anterior insular and motor cortex (Hilty et
86
al., 2011). This increased communication between cortical areas was found at the end of the
87
cycling exercise and returned to baseline after a rest period suggesting that (muscle) fatigue
88
alters interaction between cortical structures. Expanding upon these findings, the purpose of this
89
study was to localize electrical brain activity using EEG during cycling at a constant workload
90
and quantify changes in temporal brain activation throughout an exhaustive bout of exercise.
91
Based on previous observations in cycling (Schneider et al., 2009; Brümmer et al., 2011), we
92
hypothesized that electrocortical sources will be located in the motor areas, the frontal cortex and
93
the parietal cortex. Furthermore, it was hypothesized that EEG activity within motor areas will
94
demonstrate an oscillating pattern of increased and decreased activation throughout the pedalling
95
cycle. We further hypothesized that overall EEG activity would increase as fatigue develops
96
throughout the exercise duration.
97
Materials and methods
98
Subjects and ethics approval
99
Ten healthy, experienced male cyclists (Table 1) with no history of lower limb injury
100
within the past two years and no known neurological or muscular deficits volunteered for this
101
study. All subjects provided written informed consent prior to testing. All procedures complied
102
with the standards of the Declaration of Helsinki and were reviewed and approved by the
103
University of Calgary research ethics review board. Insert Table 1 here
104 105
Experimental protocol
106
Each subject performed testing on three different days. All experimental protocols were
107
supervised by a Canadian Society of Exercise Physiology Certified Exercise Physiologist
108
(CSEP-CEP) and participants were cleared for exercise using a Physical Activity Readiness
109
Questionnaire (Par-Q+). Resting blood pressure and heart rate were measured prior to daily
110
testing and standardized ceiling levels according to the Canadian Society of Exercise Physiology
111
were adhered to. All testing was completed using a Velotron Dynafit Pro cycle ergometer
112
(RacerMate Inc., Seattle, WA) that was calibrated for each testing session (Figure 1A). All
113
participants used their own pedals (clipless) and corresponding footwear that locked into the
114
pedals to accommodate their usual bike racing equipment. The first testing day was used to
115
establish the maximum aerobic power (MAP) capacity of each subject based on an incremental
116
test. Testing on the second day consisted of a time-to-exhaustion (TTE) test at 85% of the
117
individual’s MAP and was used for familiarization with the protocol and equipment. The third
118
day was a repeat of the second day and was used for the data analysis. The three testing sessions
119
for each subject were separated by 48 ± 1 hours. To increase repeatability of the results, caffeine
120
intake was restricted to two hours prior to testing and the subject was asked to maintain a normal
121
diet and refrain from any resistance training for the duration of the study.
122
Maximum aerobic power (MAP)
123
On the first day of testing an incremental cycling exercise test was used to establish blood
124
lactate thresholds LT1 and LT2 for each subject according to the double breakaway model
125
(Kindermann et al., 1979; Anderson and Rhodes, 1989). Subjects were given twenty minutes of
126
rest after the incremental test prior to starting the MAP protocol. The MAP test was started at a
127
workload slightly below the LT1 and the workload increased every minute by 25W. The test was
128
terminated when the cadence dropped by more than 15 revolutions per minute (RPM) or below
129
70 RPM (Figure 1B). Subjects were included in the study if they reached an MAP of at least
130
4.5W/kg. Insert Figure 1 here
131 132 133
Time-to-exhaustion (TTE) trial Prior to the TTE trial, each subject performed a standardized warm-up protocol based on
134
three five-minute stages corresponding to resistance values below, at, and above the measured
135
LT1. After the warm-up, the TTE trial was performed at a constant workload of 85% of the MAP
136
with the same seat and handlebar settings of the ergometer that were used to perform the MAP
137
test. On both days subjects were unaware of the workload and the elapsed time of the TTE trial
138
to ensure consistency between days. The preferred cadence was freely chosen by the cyclists
139
(typically between 90 and 100 RPM) and the test was terminated if the cadence dropped by 15
140
RPM or below 70 RPM. Verbal encouragement was provided throughout the test on both days to
141
ensure maximum performance.
142
EEG Data collection
143
Subjects cycled at constant intensity using the ergometer described above. A magnetic
144
switch was used to synchronize the data streams and identify the pedalling cycle. The switch sent
145
a square wave pulse that was recorded in synchrony with the EEG data in order to identify
146
individual pedalling revolutions. EEG was recorded using an active electrode 64-channel (10-20
147
positioning) BrainVision actiCHamp system (Brain Products GmbH, Germany). Subjects wore a
148
tight fitting electrode cap corresponding to the measured head circumference using a chinstrap to
149
reduce movement of the cap (Figure 1A). For each subject nasion, inion and preauricular points
150
were used as anatomical landmarks to position the electrode cap. Prior to data collection,
151
SuperVisc electrode gel (EasyCap GmbH, Germany) was used to ensure that the impedance was
152
less than 20 kΩ for each channel. During acquisition, all EEG channels were referenced to the
153
parietal midline (Pz) electrode and a frontal midline (AFz) electrode was used as ground. EEG
154
signals were recorded with a sampling frequency of 500 Hz per channel and stored on a
155
computer for offline processing. All data processing was performed in Matlab (The Mathworks,
156
Natick, MA) using custom written software in combination with scripts based on EEGLab
157
13.2.2b, an open source toolbox for the analysis of electrophysiological data sets (Delorme and
158
Makeig, 2004).
159
EEG Data analysis
160
Analysis of continuous EEG signals was carried out according to previously used
161
methodologies (Gwin et al., 2011; Wagner et al., 2012, 2014; Kline et al., 2014). Briefly, EEG
162
signals were band pass filtered (1-50 Hz) in order to remove signal drifts as well as line noise
163
and were re-referenced to a common average reference. Channels were removed if one of the
164
following criteria were true: (1) channels with standard deviation (SD) larger than 1000 μV, (2)
165
channels whose kurtosis was more than 5 SD from the mean and (3) channels that were
166
uncorrelated with neighbouring channels (r < 0.4) for more than 0.1% of the time. On average
167
51.6 channels were used for further analysis (SD: 1.96, range 49 – 54). Independent component
168
analysis (ICA) was used to decompose the signal into maximally independent components (IC)
169
(Makeig et al., 1996). ICA was applied to the data set of individual subjects. EEG data
170
corresponding to the entire time-to-exhaustion trial (approximately seven minutes) was used for
171
the ICA decomposition. ICs representing line noise, movement artifacts and eye activity were
172
rejected prior to further analysis (Jung et al., 2000). In order to locate the neural sources in brain
173
space, an equivalent current dipole model was computed for individual component scalp maps
174
using the DIPFIT plugin within EEGLab (Oostenveld and Oostendorp, 2002). Each IC was
175
modeled as a single dipole. We used a boundary element head model based on the Montreal
176
Neurological Institute (MNI) template that utilizes the average of 152 MRI scans from healthy,
177
young adult individuals. No specific spatial constraints were enforced as spatial accuracy is
178
limited without the individual anatomy of subjects. ICs were excluded from further analysis if
179
the projection of the dipole to the scalp accounted for less than 80% of the scalp map variance.
180
All remaining ICs of all individual subjects (approximately 100 in total) were clustered across all
181
subjects using k-means clustering algorithms implemented in EEGLab. The algorithm used the
182
dipole location in brain space and the power spectra of each individual IC for clustering purposes
183
and was reduced to the first ten principal dimensions (as in Gwin et al., 2011). Clusters and
184
dipoles indicating eye movement artifacts were not considered for further analysis. Additionally,
185
any components that were located near the lower back part of the head were not considered for
186
further analysis as they typically represent EMG activity from neck muscle that is characterized
187
by an increase in spectral power above 30 Hz. Lastly, any clusters that were located in deeper
188
brain structures near the cerebellum were not considered for analysis as EEG is not able to
189
reliable localize neural activity in these areas.
190
Spectral analysis of EEG data
191
For all individual components in the identified electrocortical clusters we calculated
192
power spectra using Welch’s method. The average power spectra for each cluster were compared
193
for different stages of the exercise corresponding to the first, middle and last 20% of all pedal
194
revolutions corresponding to the fresh, middle and fatigue condition during the time-to-
195
exhaustion trial. An analysis of variance (ANOVA) with 3 conditions (fresh, middle, fatigue)
196
was computed within the EEGLab toolbox in order to compare power spectra between
197
conditions controlling for false discovery rate (Benjamini and Hochberg, 1995) in case of
198
multiple comparisons (p < 0.05).
199
Time-frequency analysis
200
Specifically for electrocortical clusters in the motor areas, we calculated spectrograms for
201
each pedal revolution as a function of crank angle. Spectrograms were calculated based on a set
202
of sinusoidal wavelets (Delorme and Makeig, 2004). The number of cycles increased slowly with
203
increasing center frequency which allows for a better frequency resolution at higher frequencies
204
compared to a traditional wavelet approach using constant cycle length. The time-frequency
205
spectrograms were time-locked to the start of the pedal revolution and then linearly time warped
206
so that the onset and offset of each pedal revolution occurred at the same latencies (Gwin et al.,
207
2011; Gwin and Ferris, 2012; Castermans et al., 2014; Seeber et al., 2015). The time warping
208
procedure uses event latencies of each individual pedal revolution onset and offset and warps
209
these latencies to the average onset and offset latency, thereby, time normalizing each pedal
210
revolution. In order to display electrocortical fluctuations within a pedal revolution we expressed
211
the time-frequency transformed data relative to the average spectrogram (averaged across the
212
time domain of the pedal revolution) similar to previous studies (Gwin et al., 2011; Wagner et
213
al., 2012, 2014; Kline et al., 2014). This was done based on the following equation:
= 10 ∙ log
214
Where Yt is the power at time t and Y0 is the average power over the pedalling cycle.
215
These fluctuations from baseline are typically referred to as event-related spectral
216
perturbations (ERSP). We calculated ERSPs for each component within the motor area cluster. A
217
grand average ERSP was obtained by averaging the individual ERSPs. Significant pedalling
218
ERSPs were identified using bootstrapping methods available within EEGLab (200 iterations).
219
Bootstrapping is based on surrogate data distributions by randomly resampling the observed
220
data, re-computing the outcome variable and comparing the experimental data to the generated
221
surrogate data distribution. Multiple comparisons were controlled using the false discovery rate
222
with a significance level set at 0.05. All non-significant fluctuations from baseline in the ERSP
223
plots were masked for visualization in this manuscript. In order to reveal the relative timing of
224
brain activity we calculated waveforms representing spectral power fluctuations in the alpha (8-
225
12 Hz), beta (12-30 Hz) and low gamma (30-50 Hz) band. ERSPs and waveforms were
226
visualized for the first 20% of the exercise (fresh) and the last 20% of the exercise (fatigue).
227
These timeframes correspond to approximately 120 – 135 pedal revolutions. Comparisons
228
between the ERSP conditions was made within EEGLab using the bootstrapping procedure
229
outlined above using the false discovery rate to correct for multiple comparisons.
230
EMG data collection and analysis
231
While the analysis of muscle activity was not the primary focus of this manuscript,
232
surface electromyography (EMG) data were collected to ensure if potential cortical changes were
233
associated with altered muscular activation patterns. Details of the data collection and analysis
234
techniques is provided in a previously published manuscript that utilized the same techniques
235
(Enders et al., 2015). Briefly, muscle activity was collected from muscles spanning the ankle,
236
knee and hip joints with a sampling rate of 2000 Hz. Raw EMG signals were separated into
237
individual pedal revolutions according to a magnetic pedal switch. In order to obtain temporal
238
waveforms, the EMG activity was wavelet transformed (von Tscharner, 2000). The output of the
239
wavelet transform is an intensity pattern indicating the EMG activity at each time point for
240
specific center frequencies. Summing up the activity across all frequency bands results in the
241
total intensity for each time point. These temporal waveforms were normalized to the pedalling
242
cycle with a time resolution corresponding to 0.5° of the crank cycle. Grand average waveforms
243
were obtained by averaging the trials of each subject corresponding to the first and last 20% of
244
the time-to-exhaustion trial.
245
Results
246
Time-to-exhaustion
247
On average, participants were able to produce 85% of their maximum aerobic power for
248
7:04 minutes (range: 6:01 – 8:58 minutes). Participants started the trial at a cadence of
249
approximately 97 RPM. The cadence significantly dropped to 90 RPM during the last 20% of the
250
TTE trial. As per our definition of task failure, a drop in cadence of 15 RPM or below an
251
absolute number of 70 RPM, the average cadence for the last 5% of the pedal revolution
252
decreased to 86 RPM (Figure 1C).
253
Location of clusters in brain space
254
Across all subjects, we identified four electrocortical clusters that contained components
255
from more than half of our subjects. The corresponding Talairach coordinates (Lancaster et al.,
256
2000) and Brodmann areas (BA) are illustrated in Table 2. Two clusters (left and right
257
hemisphere) were identified in the parietal cortex with corresponding BA 39 and 7. One cluster
258
was identified in the superior frontal gyrus of the frontal cortex corresponding to BA 8. Lastly,
259
we identified one cluster in the supplementary motor area (SMA) of the brain in the precentral
260
gyrus of the frontal cortex corresponding to BA 6. Additionally, two more clusters were
261
identified corresponding to the right frontal and right premotor areas, however, only four and
262
five subjects, respectively, contributed to these clusters. According to our inclusion criteria if a
263
cluster did not contain dipoles of at least more than half of the tested subjects, it was not included
264
in any further analysis. Insert Table 2 here
265 266 267
Changes with exercise duration All clusters of cortical activity showed a significant increase in EEG power (p < 0.05).
268
No differences were found between the fresh and middle portion of the TTE trial and, therefore,
269
the results are focused on the fresh and fatigued states. The left frontal cortex showed an increase
270
in α-, β- and γ-band while the clusters in the supplementary motor area and left parietal cortex
271
only showed a significant increase in the α- and β-band and the right parietal cortex only for the
272
α-band (Figure 2). The strongest increase in EEG power was observed in BA 8, followed by
273
increases of smaller magnitude in the SMA (BA 6) and the parietal cortex.
274 275 276
Insert Figure 2 here Spectral fluctuations within motor area The ERSP pattern for BA 6 (SMA) is visualized with warm colors, indicating increased
277
brain activity, and cool colors indicating decreased brain activity, with respect to the average
278
spectrogram that was considered as baseline (Figure 3). All non-significant modulations of
279
spectral power were masked (green) for better visualization. All frequency bands showed an
280
oscillating pattern of increasing and decreasing activity throughout the pedalling cycle in both
281
conditions. In the fresh condition, significant modulations in α-, β- and γ-band power are
282
characterized by two local minima representing event-related desynchronization (ERD), which
283
occur before the onset of right and left down stroke, respectively. Modulations of γ-band power
284
remain similar in the fatigue condition, with the second negative peak occurring slightly later
285
compared to the fresh condition. On the contrary, modulations of α- and β-band power show an
286
altered pattern in the fatigue condition with a pronounced ERD occurring prior to the onset of the
287
right down stroke, followed by a phase of event-related synchronization (ERS) prior to the left
288
down stroke. Overall, there was a larger amplitude modulation in the α- and β-band in the fatigue
289
condition. Insert Figure 3 here
290 291 292
Muscle activation patterns Grand average temporal waveforms of EMG activity were visualized for both conditions
293
(Figure 4). The waveforms are in general agreement with previous results investigating muscle
294
activity during cycling (Ryan and Gregor, 1992; Baum and Li, 2003; Enders et al., 2015). Most
295
muscles showed a trend towards increased mucsle activity, however, across all subjects the
296
differences were not statistically different. The overall temporal pattern of muscle recruitment,
297
however, seemed very robust despite the development of significant fatigue. Insert Figure 4 here
298 299
Discussion
300
Modulation of cortical activity during cycling
301
This is the first study to quantify changes in the temporal pattern of electrocortical
302
activity during high intensity cycling in a time-to-exhaustion trial. This expands on previous
303
work investigating intra-cortical communication during a cycling time-to-exhaustion trial where
304
communication between mid/anterior insular and motor cortex significantly increased with the
305
development of fatigue (Hilty et al., 2011). Similar to previous studies on gait, the results of our
306
study highlight temporal patterns of human brain activity that are spatially localized to the motor
307
area within the frontal lobe (BA 6) (Gwin et al., 2011; Wagner et al., 2012). Significant
308
modulations from baseline were observed in EEG power for an electrocortical cluster that was
309
localized in BA 6 representing the SMA. Significant ERD in EEG power were observed across
310
the α, β and γ frequency range for the SMA during the transition from the recovery phase to the
311
power production phase for both legs indicating activation and increased neuron excitability for
312
these brain areas at these specific latencies (Pfurtscheller and Lopes da Silva, 1999; Neuper and
313
Pfurtscheller, 2001). Thus, neural activation within the SMA is greatest when muscles are
314
recruited that are critical in transitioning the lower limb from a flexed to an extended position.
315
This result is in agreement with previously described cortical activation patterns during cycling
316
(Jain et al., 2013). Interestingly, a study quantifying cortical dynamics during walking supports
317
the notion that modulation of activity in the BA 6 is related to movement with specific phases of
318
repetitive locomotor activity (Wagner et al., 2012). They also found a trend for increased neural
319
activation in the gamma frequency band in active compared to passive walking. Our results
320
during active cycling also show modulation of gamma activity in the SMA supporting the
321
suggestion by Wagner and colleagues (2012), that it might be an important contributor to motor
322
planning and sensorimotor processing.
323
The role of the motor areas in the cycling task used in this experimental protocol can be
324
explored in more detail when interpreting the alternating pattern of desynchronization and
325
synchronization of EEG power (Figure 3). Desynchronization refers to clusters of neurons that
326
are not firing in synchrony and reflect activity distributed across various brain areas. Neural
327
desynchronization is often interpreted as active information processing (Pfurtscheller and
328
Berghold, 1989; Jurkiewicz et al., 2006). An increase in desynchronization can be indicative of
329
an increase in task complexity with a higher demand for information processing (Pfurtscheller
330
and Lopes da Silva, 1999). The opposite phenomenon, an increase in EEG power at a given
331
frequency, is typically observed due to a synchronized firing of neurons. Although it is generally
332
believed that ERS in the alpha and beta band reflects coherent activity and a deactivated
333
communication between neural networks (Pfurtscheller and Lopes da Silva, 1999), in the context
334
of movement tasks it is often interpreted as movement-related sensory processing (Cassim et al.,
335
2001). In summary, ERD reflects active information processing with activity distributed across
336
neural networks while ERS reflects a more focal activation pattern indicating a neural cell
337
assembly responding in the same way. The results in this study show different responses for the
338
alpha and beta compared to the gamma frequency band and will, therefore, be discussed
339
separately.
340
In the alpha and beta band, during the fresh condition, we see a small ERD that occurs
341
during the phases of the right and left down stroke followed by ERS towards the end of the crank
342
cycle that is stronger in the beta compared to the alpha band. In both frequency bands this pattern
343
is changed during the fatigue condition and pronounced in amplitude. A clear phase of ERD
344
during the down stroke of the right leg is followed by a phase of ERS in the second half of the
345
cycle. The magnitude of ERD and ERS are increased in both alpha and beta frequency band.
346
Specifically, during the fresh condition the phases of ERD, representing an active information
347
processing behavior, correspond to the phases when the leg prepares for force production and
348
transitions from a flexed to an extended position (see Figure 3).
349
The gamma frequency band shows similar behavior to the alpha and beta frequency
350
bands during the fresh condition with oscillations of ERD and ERS throughout the pedalling
351
cycle. Again, phases of ERD, interpreted as active information processing, correlate well with
352
the timing of leg force production. While the temporal pattern seemed to be altered in the alpha
353
and beta band, the oscillations in the gamma band remained similar. In both conditions, the
354
gamma band shows a pattern of two local minima corresponding to the initiation of power
355
production and preparation for the extension phase of the right and left leg.
356
The oscillating pattern of ERD and ERS could be interpreted as an alternating behavior in
357
the motor areas reflecting a distributed activity and active information processing in preparation
358
for leg extension (ERD) followed by a focal activation behavior (ERS) during the main phase of
359
force production. This pattern was most pronounced in the gamma frequency band. This
360
observation was made in a cluster slightly localized in the left SMA. One may expect a mirror
361
effect in the other hemisphere or a cluster that is located over the central midline. Our identified
362
cluster contained dipoles located in both the left hemisphere as well as in the midline of BA 6.
363
While we did find a cluster in the right part of the premotor cortex, this cluster did not contain
364
ICs from more than half of our subjects.
365
Overall, we speculate that differences in EEG power modulation between the fresh and
366
fatigued state are due to an increased effort and task difficulty that participants need to cope
367
with. An alternative option would be a substantial change in neuromuscular strategies of the leg
368
muscles to keep up with the cycling task. However, preliminary analysis of EMG data that was
369
collected in synchrony with the EEG data does not support the hypothesis of altered lower limb
370
neuromuscular strategies as the overall EMG pattern remained similar (Figure 4).
371
Spectral power increase in EEG with fatigue
372
Our study observed a broadband increase in EEG power, even though for most identified
373
clusters this was only true in the α- and β-band. This increase across a broader frequency band is
374
in contrast to studies that have investigated different walking conditions and found changes
375
limited to the α- and β-band (Wagner et al., 2012, 2014). It is unclear, however, if our fatigue
376
protocol directly relates to comparisons of active and passive walking (Wagner et al., 2012) or
377
different visual feedback conditions (Wagner et al., 2014). Typically, active movement is
378
associated with desynchronization in the α- and β-band and, therefore, it is not surprising to find
379
such differences between active and passive movement. Instead our comparison was between
380
two conditions that were both based on intense movements. Interestingly, a study that has
381
investigated fatigue on a recumbent bicycle ergometer shows results that are very similar to this
382
study. They showed that as participants fatigued during the cycling exercise there was a
383
broadband increase in spectral power that was significant in the θ (4.5 – 7.99 Hz), low α (8 –
384
10.49 Hz), high α (10.5 – 12.99 Hz), low β 1 (13 – 17.99 Hz) and high β (18 – 30 Hz) frequency
385
bands (Bailey et al., 2008). Another possibility is noise in the acquired EEG that was not
386
removed by our data analysis. While this may partially contribute to our observation, the changes
387
in modulation of ERD and ERS also suggest a response in cortical involvement during the
388
exercise. Limitations with respect to noise and motion artifacts are further discussed at the end of
389
this section.
390
Interaction of electrocortical clusters
391 392
The four identified clusters of electrical activity in specific cortical domains (prefrontal, premotor and posterior parietal cortices) were consistently active in the majority of subjects in
393
this study. Common activation between these areas indicate a sequence of cognitive processing
394
that includes executive control, motor planning and execution as well as sensorimotor
395
integration. Indeed, previous research suggests that activation in the prefrontal cortex area relates
396
to the executive behavioural control function (Tanji and Hoshi, 2008) and motor planning (Wise
397
et al., 1997; Rizzolatti and Luppino, 2001; Tanji, 2001). Information about motor planning may
398
be passed on to the SMA and premotor cortex as it is well documented that these areas have
399
strong connections to the prefrontal cortex (Rizzolatti and Luppino, 2001; Chouinard and Paus,
400
2006), specifically the dorsolateral prefrontal cortex (Lu et al., 1994). Thus, the output of the
401
prefrontal cortex likely targets specific premotor areas. It has been proposed that the prefrontal
402
dependent motor areas (anterior parts of the premotor cortex and SMA) receive higher order
403
cognitive information regarding motor plans and motivation (Rizzolatti and Luppino, 2001)
404
which supports the executive control function of the frontal areas. Interestingly, the largest
405
response in neural activity accompanying fatigue is seen in the frontal areas. We speculate that
406
this is consistent with an increase in executive control to keep up with the challenging task
407
demand. However, this speculation needs careful evaluation in a future study.
408
In our analysis no cluster was directly associated to activity in the primary motor cortex.
409
Several technical as well as functional explanations may describe this observation. It has to be
410
noted that EEG has a limited spatial resolution limiting the exact location of activity within the
411
motor areas. Further, individual anatomy varies across participants and this study did not utilize
412
individual MRI head models which would improve source localization. However, it is well
413
known that strong interconnections between premotor and primary motor areas exist. It suggests
414
that information that is processed and passed on to the premotor cortex is directed further to the
415
primary motor cortex (Künzle, 1978; Chouinard and Paus, 2006) and the spinal cord (Dum and
416
Strick, 1991; Münchau et al., 2002). In fact, the majority of cortico-cortical connections of the
417
primary motor cortex are with BA 6 (Leichnetz, 1986; Stepniewska et al., 1993) implying that
418
reciprocal information processing between premotor and primary areas plays a crucial role in
419
motor execution. It is further known that both premotor and primary motor areas have strong
420
connections to the spinal cord through corticospinal neurons that are activated during movement
421
(Dum and Strick, 1991; Quallo et al., 2012).
422
Lastly, the connectivity between parietal areas and motor areas (specifically the SMA and
423
premotor cortex) has been demonstrated in many different scenarios. For a detailed review of
424
this network the reader is referred to more in depth literature (Desmurget and Sirigu, 2009).
425
Studies link activity in the parietal cortex to motor awareness and movement intention prior to
426
movement (Desmurget et al., 2009), coordinative aspects of movement, and functioning as an
427
interface between sensorimotor and visually guided aspects of movements (Taira et al., 1990;
428
Buneo and Andersen, 2006). The fact that the posterior parietal cortex plays a role in sensory and
429
motor function led to the proposal that it is heavily involved in sensory-motor interaction.
430
Interestingly, the parietal cortex has typically been identified in studies focusing on real-time
431
EEG recordings during repetitive locomotor tasks (Gwin et al., 2011; Wagner et al., 2014).
432
Specifically, a premotor-parietal network has been identified during rhythmic walking which is
433
believed to represent motor planning and intention and has been shown to increase activity in the
434
presence of task feedback (Wagner et al., 2014). In summary, connectivity between parietal and
435
premotor areas seems to be crucial for motor planning as well as sensory motor integration.
436
Increased EEG power was observed in this area which might suggest the importance of motor
437
planning and the integration of sensory input towards the end of an intense exercise.
438
Limitations: Muscle drive vs. sensory afferent integration
439
The spectral analysis performed in this study adds to the emerging understanding that the
440
cortex may play a more important role during rhythmic movement activity than previously
441
thought. It confirms similar findings of alternating latencies of event-related desynchronization
442
and synchronization within a movement cycle (Gwin et al., 2011; Wagner et al., 2012, 2014; Jain
443
et al., 2013). However, solely based on the spectral analysis it remains unknown if the observed
444
cortical activation patterns are directly involved in the transmission of motor commands to
445
facilitate muscle activation. Another possibility, a hypothesis that was previously discussed
446
(Gwin et al., 2011), is the integration of sensory afferent signals that are used to modulate
447
efferent corticospinal transmission. Most likely, the observed activity reflects involvement in
448
both descending efferent pathways as well as sensory information processing. In order to
449
investigate this proposal in more detail, a coherence analysis should be used in order to address
450
the idea of an anatomical coupling between cortical and muscular signals during movement
451
(Petersen et al., 2012).
452
Limitations: Movement-related artifacts
453
Movement-related and specifically gait-related artifacts have been a central issue to the
454
debate about the interpretation of EEG data recorded during locomotor activities. Following the
455
recommendation of previous reports, we excluded any components that were clearly indicative
456
of electromyographic artifacts (neck, jaw, facial). Additionally, we recorded the data in an
457
electrically shielded chamber to minimize electromagnetic interference from power lines and
458
devices. Due to the nature of the experimental protocol, we had to address the issue of
459
perspiration. Sweating typically results in very slow drifts of the EEG signal. We found these
460
drifts to be well below 1 Hz and, thus, we were able to remove these drifts using the applied
461
filters. Additionally, sweating was relatively constant throughout the exercise as subjects started
462
to sweat during the warm up period. Therefore, electrode impedance did not change significantly
463
throughout the TTE trial. Participants were not allowed to stand up from the saddle to minimize
464
head and upper body movement as much as possible. However, some movement still occurs
465
during an intense cycling trial and while we have addressed common noise issues from an
466
experimental and computational side, there is a remaining issue as artifact sources were shown to
467
greatly vary across subjects and efforts (Kline et al., 2015) making it difficult to reject artifacts
468
based on templates and models. Therefore, most likely, many reports of EEG during locomotion
469
contain a mixture of true electrocortical and artifact signals which should be kept in mind by the
470
reader.
471
Acknowledgements We would like to thank Jessi Temple for the data collection and supervision of the
472 473
exercise physiology portion of the study protocol as well as Elias Tomaras for technical
474
assistance. We would also like to acknowledge the time commitment and effort of all
475
participants who contributed to this study.
476
Grants
477
H.E. is supported by Vanier Canada (Natural Sciences and Engineering Research
478
Council), the Killam foundation as well as the CREATE program of the Natural Sciences and
479
Engineering Research council. J.B. is supported by Vanier Canada (Canadian Institute of Health
480
Research). A.B.P. is supported by Natural Sciences and Engineering Research Council of
481
Canada (NSERC), Canadian Foundation for Innovation Leaders Opportunity Fund, and Alberta
482
Enterprise and Advanced Education Research Capacity Program, Alberta Alignment Grants.
483
Disclosures
484 485 486 487 488
There are no professional relationships with companies or manufacturers to disclose for all authors.
489
References
490 491
Anderson GS, Rhodes EC. A review of blood lactate and ventilatory methods of detecting transition thresholds. Sports Med 8: 43–55, 1989.
492 493
Bailey SP, Hall EE, Folger SE, Miller PC. Changes in EEG during graded exercise on a recumbent cycle ergometer. J Sport Sci Med 7: 505–511, 2008.
494 495
Baum BS, Li L. Lower extremity muscle activities during cycling are influenced by load and frequency. J Electromyogr Kinesiol 13: 181–190, 2003.
496 497
Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J R Stat Soc Ser B 57: 289–300, 1995.
498 499
Brümmer V, Schneider S, Strüder H, Askew C. Primary motor cortex activity is elevated with incremental exercise intensity. Neuroscience 181: 150–62, 2011.
500 501
Buneo CA, Andersen RA. The posterior parietal cortex: sensorimotor interface for the planning and online control of visually guided movements. Neuropsychologia 44: 2594–606, 2006.
502 503 504
Cassim F, Monaca C, Szurhaj W, Bourriez JL, Defebvre L, Derambure P, Guieu JD. Does post-movement beta synchronization reflect an idling motor cortex? Neuroreport 12: 3859–63, 2001.
505 506 507
Castermans T, Duvinage M, Cheron G, Dutoit T. About the cortical origin of the low-delta and high-gamma rhythms observed in EEG signals during treadmill walking. Neurosci Lett 561: 166–70, 2014.
508 509
Chouinard PA, Paus T. The primary motor and premotor areas of the human cerebral cortex. Neuroscientist 12: 143–52, 2006.
510 511
Delorme A, Makeig S. EEGLAB: an open source toolbox for analysis of single-trial EEG dynamics including independent component analysis. J Neurosci Methods 134: 9–21, 2004.
512 513
Desmurget M, Reilly KT, Richard N, Szathmari A, Mottolese C, Sirigu A. Movement intention after parietal cortex stimulation in humans. Science 324: 811–3, 2009.
514 515
Desmurget M, Sirigu A. A parietal-premotor network for movement intention and motor awareness. Trends Cogn Sci 13: 411–9, 2009.
516 517
Dum RP, Strick PL. The origin of corticospinal projections from the premotor areas in the frontal lobe. J Neurosci 11: 667–89, 1991.
518 519
Enders H, Maurer C, Baltich J, Nigg BM. Task-Oriented Control of Muscle Coordination during Cycling. Med Sci Sports Exerc 45: 2298–2305, 2013.
520 521
Enders H, von Tscharner V, Nigg BM. Neuromuscular Strategies during Cycling at Different Muscular Demands. Med Sci Sport Exerc 47: 1450–1459, 2015.
522 523
Gramann K, Gwin JT, Ferris D, Oie K, Jung T-P, Lin C-T, Liao L-D, Makeig S. Cognition in action: imaging brain/body dynamics in mobile humans. Rev Neurosci 22: 593–608, 2011.
524 525
Grillner S, Wallén P, Saitoh K, Kozlov A, Robertson B. Neural bases of goal-directed locomotion in vertebrates--an overview. Brain Res Rev 57: 2–12, 2008.
526 527
Gwin J, Gramann K, Makeig S, Ferris D. Electrocortical activity is coupled to gait cycle phase during treadmill walking. Neuroimage 54: 1289–96, 2011.
528 529
Gwin JT, Ferris D. An EEG-based study of discrete isometric and isotonic human lower limb muscle contractions. J Neuroeng Rehabil 9: 1–11, 2012.
530 531 532
Hilty L, Langer N, Pascual-Marqui R, Boutellier U, Lutz K. Fatigue-induced increase in intracortical communication between mid/anterior insular and motor cortex during cycling exercise. Eur J Neurosci 34: 2035–42, 2011.
533 534
Jain S, Gourab K, Schindler-Ivens S, Schmit B. EEG during pedaling: evidence for cortical control of locomotor tasks. Clin Neurophysiol 124: 379–90, 2013.
535 536 537
Jung TP, Makeig S, Humphries C, Lee TW, McKeown MJ, Iragui V, Sejnowski TJ. Removing electroencephalographic artifacts by blind source separation. Psychophysiology 37: 163–78, 2000.
538 539
Jurkiewicz MT, Gaetz WC, Bostan AC, Cheyne D. Post-movement beta rebound is generated in motor cortex: evidence from neuromagnetic recordings. Neuroimage 32: 1281–9, 2006.
540 541 542
Kindermann W, Simon G, Keul J. The significance of the aerobic-anaerobic transition for the determination of work load intensities during endurance training. Eur J Appl Physiol Occup Physiol 42: 25–34, 1979.
543 544
Kline JE, Huang HJ, Snyder KL, Ferris DP. Isolating gait-related movement artifacts in electroencephalography during human walking. J Neural Eng 12: 046022, 2015.
545 546
Kline JE, Poggensee K, Ferris DP. Your brain on speed: cognitive performance of a spatial working memory task is not affected by walking speed. Front Hum Neurosci 8: 288, 2014.
547 548
Künzle H. An autoradiographic analysis of the efferent connections from premotor and adjacent prefrontal regions (areas 6 and 9) in macaca fascicularis. Brain Behav Evol 15: 185–234, 1978.
549 550 551
Lancaster JL, Woldorff MG, Parsons LM, Liotti M, Freitas CS, Rainey L, Kochunov P V, Nickerson D, Mikiten SA, Fox PT. Automated Talairach atlas labels for functional brain mapping. Hum Brain Mapp 10: 120–31, 2000.
552 553 554
Leichnetz GR. Afferent and efferent connections of the dorsolateral precentral gyrus (area 4, hand/arm region) in the macaque monkey, with comparisons to area 8. J Comp Neurol 254: 460– 92, 1986.
555 556
Lu MT, Preston JB, Strick PL. Interconnections between the prefrontal cortex and the premotor areas in the frontal lobe. J Comp Neurol 341: 375–92, 1994.
557 558
Makeig S, Bell AJ, Jung T-P, Sejnowski TJ. Independent component analysis of electroencephalographic data. Adv Neural Inf Process Syst 8: 145–151, 1996.
559 560
Makeig S, Gramann K, Jung T-P, Sejnowski TJ, Poizner H. Linking brain, mind and behavior. Int J Psychophysiol 73: 95–100, 2009.
561 562 563
Münchau A, Bloem BR, Irlbacher K, Trimble MR, Rothwell JC. Functional connectivity of human premotor and motor cortex explored with repetitive transcranial magnetic stimulation. J Neurosci 22: 554–61, 2002.
564 565
Neuper C, Pfurtscheller G. Event-related dynamics of cortical rhythms: frequency-specific features and functional correlates. Int J Psychophysiol 43: 41–58, 2001.
566 567 568
Oostenveld R, Oostendorp TF. Validating the boundary element method for forward and inverse EEG computations in the presence of a hole in the skull. Hum Brain Mapp 17: 179–92, 2002.
569 570
Petersen TH, Willerslev-Olsen M, Conway B, Bo Nielsen J. The motor cortex drives the muscles during walking in human subjects. J Physiol 10: 2443–2452, 2012.
571 572
Pfurtscheller G, Berghold A. Patterns of cortical activation during planning of voluntary movement. Electroencephalogr Clin Neurophysiol 72: 250–8, 1989.
573 574
Pfurtscheller G, Lopes da Silva FH. Event-related EEG/MEG synchronization and desynchronization: basic principles. Clin Neurophysiol 110: 1842–57, 1999.
575 576
Pfurtscheller G, Neuper C, Krausz G. Functional dissociation of lower and upper frequency mu rhythms in relation to voluntary limb movement. Clin Neurophysiol 111: 1873–1879, 2000.
577 578
Presacco A, Goodman R, Forrester L, Contreras-Vidal J. Neural decoding of treadmill walking from noninvasive electroencephalographic signals. J Neurophysiol 106: 1875–87, 2011.
579 580
Quallo MM, Kraskov A, Lemon RN. The activity of primary motor cortex corticospinal neurons during tool use by macaque monkeys. J Neurosci 32: 17351–64, 2012.
581
Rizzolatti G, Luppino G. The cortical motor system. Neuron 31: 889–901, 2001.
582 583
Ryan MM, Gregor RJ. EMG profiles of lower extremity muscles during cycling at constant workload and cadence. J Electromyogr Kinesiol 2: 69–80, 1992.
584 585 586
De Sanctis P, Butler JS, Green JM, Snyder AC, Foxe JJ. Mobile brain/body imaging (MoBI): High-density electrical mapping of inhibitory processes during walking. Proc. Annu. Int. Conf. IEEE Eng. Med. Biol. Soc. EMBS (2012). doi: 10.1109/EMBC.2012.6346236.
587 588 589
Schneider S, Brümmer V, Abel T, Askew C, Strüder H. Changes in brain cortical activity measured by EEG are related to individual exercise preferences. Physiol Behav 98: 447–52, 2009.
590 591 592
Seeber M, Scherer R, Wagner J, Solis-Escalante T, Muller-Putz GR. EEG beta suppression and low gamma modulation are different elements of human upright walking. Front Hum Neurosci 8: 1–9, 2014.
593 594 595
Seeber M, Scherer R, Wagner J, Solis-Escalante T, Müller-Putz GR. High and low gamma EEG oscillations in central sensorimotor areas are conversely modulated during the human gait cycle. Neuroimage 112: 318–326, 2015.
596 597 598
Stepniewska I, Preuss TM, Kaas JH. Architectonics, somatotopic organization, and ipsilateral cortical connections of the primary motor area (M1) of owl monkeys. J Comp Neurol 330: 238– 71, 1993.
599 600
Taira M, Mine S, Georgopoulos AP, Murata A, Sakata H. Parietal cortex neurons of the monkey related to the visual guidance of hand movement. Exp Brain Res 83, 1990.
601 602 603 604
Taniguchi M, Kato A, Fujita N, Hirata M, Tanaka H, Kihara T, Ninomiya H, Hirabuki N, Nakamura H, Robinson SE, Cheyne D, Yoshimine T. Movement-related desynchronization of the cerebral cortex studied with spatially filtered magnetoencephalography. Neuroimage 12: 298–306, 2000.
605 606
Tanji J, Hoshi E. Role of the lateral prefrontal cortex in executive behavioral control. Physiol Rev 88: 37–57, 2008.
607 608
Tanji J. Sequential organization of multiple movements: involvement of cortical motor areas. Annu Rev Neurosci 24: 631–51, 2001.
609 610
Von Tscharner V. Intensity analysis in time-frequency space of surface myoelectric signals by wavelets of specified resolution. J Electromyogr Kinesiol 10: 433–45, 2000.
611 612 613
Wagner J, Solis-Escalante T, Grieshofer P, Neuper C, Müller-Putz G, Scherer R. Level of participation in robotic-assisted treadmill walking modulates midline sensorimotor EEG rhythms in able-bodied subjects. Neuroimage 63: 1203–11, 2012.
614 615 616
Wagner J, Solis-Escalante T, Scherer R, Neuper C, Müller-Putz G. It’s how you get there: walking down a virtual alley activates premotor and parietal areas. Front Hum Neurosci 8: 93, 2014.
617 618
Waldert S, Preissl H, Demandt E, Braun C, Birbaumer N, Aertsen A, Mehring C. Hand movement direction decoded from MEG and EEG. J Neurosci 28: 1000–8, 2008.
619 620
Wise SP, Boussaoud D, Johnson PB, Caminiti R. Premotor and parietal cortex: corticocortical connectivity and combinatorial computations. Annu Rev Neurosci 20: 25–42, 1997.
621
Figure captions
622
Figure 1: The top panel (A) shows a schematic drawing of the experimental setup. Participants
623
were equipped with a tight fitting 64 channel EEG system and several electrodes to record lower
624
limb electromyography (not used in this study). Cables were securely attached to the subjects
625
back to avoid artifacts due to cable movement. The middle panel (B) shows the cadence of a
626
representative subject for a complete time-to-exhaustion trial (gray line). The black line shows
627
the cadence when applying a moving average with a window size of 50 pedal revolution. The
628
dashed black line shows the cadence that the participant voluntarily chose to start the trial. The
629
dotted line shows the corresponding cut off cadence to define task failure. The areas with the
630
horizontal lines depict the time frame corresponding to the fresh and fatigue condition used for
631
the analysis. The shaded gray area corresponds to the last 5% of the trial ultimately prior to task
632
failure. The bottom panel (C) shows the average cadence across all subjects for the fresh and
633
fatigue phase as well as immediately before task failure. * denotes a significant difference
634
compared to the fresh condition (p < 0.01). † denotes a significant difference compared to the
635
fresh and fatigue condition (p < 0.01).
636
Figure 2: Scalp map projections from the dipole clusters and power spectra for electrocortical
637
sources corresponding to frontal cortex (BA 8, top left), SMA (BA 6, bottom left), right parietal
638
cortex (BA 7, top right), and left parietal cortex (BA 39, bottom right). The color scale indicates
639
the strength of the cluster average scalp projection to each electrode. The traces indicate the fresh
640
(black) and fatigue (gray) conditions with shaded gray areas indicating the 95% confidence
641
intervals. The MRI images show the dipole clusters in the sagittal and coronal planes.
642
Figure 3: ERSP plot (A) of BA 6 showing modulations in EEG spectral power throughout the
643
pedalling cycle (x-axis) as a function of frequency (y-axis, log scaled) with warm and cool colors
644
indicating increased and decreased activity with respect to the average spectrogram, respectively.
645
The left column shows data for the fresh condition (first 20%) and the middle column for the
646
fatigue condition (last 20%). The right column shows the difference (fresh minus fatigue)
647
between the two conditions with any non-significant differences set to zero. The plots below (B)
648
show the time-dependent brain activation across the pedaling cycle in the α, β and low γ
649
frequency bands. The shaded pink and white areas indicate the down stroke (power production
650
phase) of the right and left leg, respectively. The onset of right down stroke (RDS) and left down
651
stroke (LDS) corresponds to the right and left leg passing the top dead centre and is indicated by
652
the arrows below the x-axis. (For interpretation of the references to color in this figure legend,
653
the reader is referred to the web version of this article.)
654
Figure 4: Grand average EMG waveforms in the fresh (black) and fatigue (gray) condition for
655
six lower limb muscles. The x-axis shows the crank cycle with zero referring to the top dead
656
centre (TDC). The y-axis shows the EMG activity normalized to the mean activity of each
657
muscle across the pedalling cycle in the fresh condition.
658
Tables Table 1: Subjects’ characteristics including lactate thresholds (LT 1 and LT2) and 85% of maximum aerobic power (MAP) that was used for the time-to-exhaustion trial. Subject ID Age Weight MAP MAP LT 1 LT 2 85% MAP # [years] [kg] [W] [W/kg] [W] [W] [W] 35 77.4 425 5.5 250 300 360 1 23 82.5 386 4.7 200 250 330 2 22 68.2 437 6.4 250 300 370 3 29 71.4 438 6.1 250 300 370 4 30 62.2 360 5.8 175 225 310 5 31 72.4 351 4.8 150 225 300 6 35 67.0 450 6.7 275 325 380 7 18 75.1 430 5.7 225 300 365 8 25 65.8 413 6.3 225 275 350 9 27 75.5 413 5.5 225 275 350 10 27.5 71.8 410.3 5.8 222.5 277.5 348.5 Mean 5.6 6.1 33.8 0.7 38.1 34.3 26.9 SD
659
Table 2: Talairach coordinates for identified clusters (geometric mean) Cluster Nearest Gray Talairach Matter coordinates 1 Brodmann Area Frontal cortex, BA 8 (-21, 48, 52) superior frontal gyrus Frontal cortex, BA 6 (-23,-19,70) precentral gyrus Parietal cortex, BA 7 (38, -60, 60) superior parietal lobule Parietal cortex, BA 39 (-48, -65, 46) inferior parietal lobule 1 Determined using the Talairach client (www.talairach.org) (Lancaster et al., 2000) 660
A
EEG cap Pedal switch
EMG electrodes
Controlled power output Data acquisition boxes (EMG & EEG)
B
Cadence [RPM]
100 80 60 40 20 0 0
Fresh 1
Fatigue 2
3 4 Time [min]
5
6
7
C Cadence [RPM]
100
*
†
80 60 40 20 0
Fresh
Fatigue
Task Failure
(last 5% of pedal revolutions)
Frontal (BA 8)
9 Subjects (10 ICs)
7 Subjects (9 ICs)
Parietal Right (BA 7)
Fresh Fatigue
15
Power(10*log10(μV2))
Power(10*log10(μV2))
20
15 10
10 5 0 -5 5
0 -5
10 15 20 25 30 35 40 45 Frequency (Hz)
5
6 Subjects (8 ICs)
10 15 20 25 30 35 40 45 Frequency (Hz)
Parietal Left (BA 39) Power(10*log10(μV2))
Supplementary Motor Area (BA 6) Power(10*log10(μV2))
8 Subjects (9 ICs)
5
15
15
10
10
5
0 5
10 15 20 25 30 35 40 45 Frequency (Hz)
-3
-2
-1
0
1
2
3
5 0 -5 5
10 15 20 25 30 35 40 45 Frequency (Hz)
Frontal (BA 8)
9 Subjects (10 ICs)
7 Subjects (9 ICs)
Parietal Right (BA 7)
Fresh Fatigue
15
15
Power(10*log10(μV2))
Power(10*log10(μV2))
20
10
10 5 0 -5 5
-5 5
10 15 20 25 30 35 40 45 Frequency (Hz)
Parietal Left (BA 39)
6 Subjects (8 ICs) Power(10*log10(μV2))
Power(10*log10(μV2))
0
10 15 20 25 30 35 40 45 Frequency (Hz)
Supplementary Motor Area (BA 6)
8 Subjects (9 ICs)
5
15
15
10
10 5 0 5
10 15 20 25 30 35 40 45 Frequency (Hz)
-3
-2
-1
0
1
2
3
5 0 -5 5
10 15 20 25 30 35 40 45 Frequency (Hz)
Fresh Frequency (Hz)
12 8
3
25 50 75 100 Pedal Cycle (%) -0.5
α
0
B
Alpha 8−12Hz
0.5
0.5
0
0
−0.5
−0.5
−1
−1
ERSP (dB)
1
Beta 12−30Hz
Beta 12−30Hz
1
1
0.5
0.5
0
0
−0.5
−0.5
−1
−1 Low Gamma 30−50Hz
Low Gamma 30−50Hz 1
1
0.5
0.5
0
0
−0.5
−0.5
−1
−1
RDS
50
75 100 LDS
Pedal Cycle (%)
25 50 75 100 Pedal Cycle (%)
-0.5
1
25
0
0.5
Alpha 8−12Hz
0
12 8 3
25 50 75 100 Pedal Cycle (%)
0 ERSP (dB)
ERSP (dB)
0
30
β 12 8
Difference
50
γ
30
ERSP (dB)
Frequency (Hz)
30
3
Fatigue
50
Frequency (Hz)
A 50
0
25 RDS
50
75 100 LDS
Pedal Cycle (%)
0 ERSP (dB)
0.5
Fresh Frequency (Hz)
12 8
3
25 50 75 100 Pedal Cycle (%) -0.5
α
0
B
Alpha 8−12Hz
0.5
0.5
0
0
−0.5
−0.5
−1
−1
ERSP (dB)
1
Beta 12−30Hz
Beta 12−30Hz
1
1
0.5
0.5
0
0
−0.5
−0.5
−1
−1 Low Gamma 30−50Hz
Low Gamma 30−50Hz 1
1
0.5
0.5
0
0
−0.5
−0.5
−1
−1
RDS
50
75 100 LDS
Pedal Cycle (%)
25 50 75 100 Pedal Cycle (%)
-0.5
1
25
0
0.5
Alpha 8−12Hz
0
12 8 3
25 50 75 100 Pedal Cycle (%)
0 ERSP (dB)
ERSP (dB)
0
30
β 12 8
Difference
50
γ
30
ERSP (dB)
Frequency (Hz)
30
3
Fatigue
50
Frequency (Hz)
A 50
0
25 RDS
50
75 100 LDS
Pedal Cycle (%)
0 ERSP (dB)
0.5
EMG Activity (rel)
3 2 1 0
0
90
0 (TDC) 90 Crank Angle (degrees)
180
EMG Activity (rel)
Vastus Medialis EMG Activity (rel)
3 2 1 0 Biceps Femoris
2
3 2 1 0 180
Soleus
Vastus Lateralis
4
EMG Activity (rel)
EMG Activity (rel)
EMG Activity (rel)
Med. Gastrocnemius
2 1 0 Rectus Femoris
3
Fresh Fatigue
2 1 0 180
90
0 (TDC) 90 Crank Angle (degrees)
180