Hemodynamic Response Alteration As a Function of Task Complexity and Expertise—An fNIRS Study in Jugglers

Detailed knowledge about online brain processing during the execution of complex motor tasks with a high motion range still remains elusive. The aim of the present study was to investigate the hemodynamic responses within sensorimotor networks as well as in visual motion area during the execution of a complex visuomotor task such as juggling. More specifically, we were interested in how far the hemodynamic response as measured with functional near infrared spectroscopy (fNIRS) adapts as a function of task complexity and the level of the juggling expertise. We asked expert jugglers to perform different juggling tasks with different levels of complexity such as a 2-ball juggling, 3- and 5-ball juggling cascades. We here demonstrate that expert jugglers show an altered neurovascular response with increasing task complexity, since a 5-ball juggling cascade showed enhanced hemodynamic responses for oxygenated hemoglobin as compared to less complex tasks such as a 3- or 2-ball juggling pattern. Moreover, correlations between the hemodynamic response and the level of the juggling expertise during the 5-ball juggling cascade, acquired by cinematographic video analysis, revealed only a non-significant trend in primary motor cortex, indicating that a higher level of expertise might be associated with lower hemodynamic responses.

The aim of the present study was to investigate the hemodynamic response within sensorimotor networks [bilateral primary motor cortex (M1), somatosensory cortex (S1) and premotor cortex (PMC)] as well as in visual motion area (MT, V5) during the execution of a complex visuomotor task such as juggling. More specifically, we were interested in how far the hemodynamic response is altered by task complexity. Based on previous findings, higher task complexity in at least simple motor tasks seems to be related with altered hemodynamic response function (Holper et al., 2009). In the present study, we aimed at differentiating between task complexity and altered brain response during the execution of a unilateral 2-ball juggling pattern (right and left hand separately), bilateral 3-and 5-ball cascades.
Furthermore, we were interested in the amount of change of the hemodynamic response during juggling as a function of juggling expertise, an index for the individual learning stage. Short-term learning (over weeks) of simple juggling patterns is associated with structural cerebral changes (Draganski et al., 2004;Boyke et al., 2008;Driemeyer et al., 2008;Scholz et al., 2009). Further studies considering different complex motor tasks report similar findings for structural plasticity associated with short-term motor learning (Taubert et al., 2010;Gryga et al., 2012). There is also strong evidence from studies using fNIRS (Hatakenaka et al., 2007;Ikegami and Taga, 2008;Leff et al., 2008;Gentili et al., 2013) that shortterm learning is associated with functional cerebral changes. Since structural MRI studies provided first evidence that the individual juggling expertise (long-term motor learning over years) is associated with larger gray matter density in MT/V5 (Gerber et al., 2014), we also investigated if a similar relationship could be observed in experienced jugglers on a functional level.
In summary, we hypothesized, according to previous reports on simple motor tasks (Holper et al., 2009), that the hemodynamic response in sensorimotor networks as well as visual motion areas increases as a function of task complexity. Since this is the first study investigating neurovascular effects during the execution of a complex visuomotor task (juggling), we obviously cannot make direct inferences about the regional specific modulatory effects in sensorimotor and/ or visual motion areas. Apart from this limitation, we further hypothesized that task-complexity dependent hemodynamic response alterations are reliant on the individual juggling expertise. More specifically, highly experienced jugglers should show less neurovascular alterations in highly complex juggling tasks (e.g., 5-ball cascade) as compared to less experienced peers.

Subjects
Twenty-three healthy male jugglers (2 left-handed) were enrolled in the present study. Jugglers had to fulfill the following criteria in order to be eligible for the study: Participants should be able to maintain a bimanual 5-ball cascade for at least 20 s in eight consecutive trials. Those who fulfilled this criterion were regarded as expert jugglers. Apart from this, all participants had to be free of any neurological and/or psychological disease. Three participants were not included in the final analyses since they were not able to maintain the requested 5-ball cascade. Furthermore, five participants had to be removed due to poor fNIRS signal quality. Hence, a total number of 15 expert jugglers were finally analyzed in the present study. The average age of the 15 expert jugglers was 26.3 ± 5.2 years (ranging from 17 to 38 years). We used a questionnaire regarding their level of expertise (e.g., hours of training/week and years of juggling experience), which is similar to Gerber et al. (2014). Handedness was acquired by the Edinburgh handedness questionnaire (Oldfield, 1971).

Experimental Procedure
Participants filled out the Edinburgh questionnaire and a questionnaire on their level of expertise during the preparation of the fNIRS setup. Subsequently, participants completed a 5 min standardized warm-up session to familiarize themselves with specific juggling balls (SIL-X balls, diameter 78 mm, mass 150 g). Cinematographic and hemodynamic measurements (fNIRS) were performed during the execution of a unimanual 2-ball juggling pattern with either the left or right hand, a bimanual 3-ball and 5-ball cascade. Since task complexity during juggling is associated with higher movement speed, which may influence hemodynamic responses (Obrig et al., 1996), two control tasks were additionally tested. Here, participants had to move their left and right arm in alternating order with a frequency of 1-2 Hz while holding juggling balls in their hands (which was controlled by a metronome).
During each of the afore mentioned tasks, participants performed the specific juggling task in a block design for 8 × 20 s with 20-25 s rest in between. Between each juggling task, there was a rest period of ∼60 s. The order of juggling tasks was randomized for each individual (see Figure 1A). At the end of the experiment, participants rated their perceived level of complexity of each task on a scale from 1 to 10. Additionally, as a qualitative measure, an expert rater ranked the participants (1 being the best, 15 the weakest), by judging their performance (i.e., number of balls and clubs and variability with the execution of the 5-ball cascade; similar to Gerber et al., 2014). This ranking is presented in Table 1.

Cinematographic Analysis
For a quantitative rating and ranking of the individual juggling expertise, participants were video-taped using a HD-Camcorder placed on a tripod (Panasonic HDC-SD909, CMOS camera systems, Kadoma, Japan) during the entire experimental session.   (Oldfield, 1971); R: right-handed; L: left-handed; RL: no preference. *The jugglers 6, 7, and 8, as well as 9, 10, and 11, as well as 13 and 14 were ranked as equally good.
The video was synchronized with the fNIRS system. To enhance the contrast for the cinematographic analysis, the jugglers used one orange ball, which allows the number of loops to be counted, and the remaining balls were white in front of a black carpet (see Figure 1B, Supplementary Video). Videos were taken with a sampling rate of 50 frames per second and a short shutter (1/500 s). For a detailed description of the analyses steps see Section Data Analyses.

Functional Near-Infrared Spectroscopy (fNIRS)
We recorded hemodynamic responses in sensorimotor (M1, S1, PMC) and visual motion area (V5/MT) on both hemispheres (Hashimoto et al., 2006;Mehnert et al., 2013) using a portable fNIRS system (NIRSport 88 by NIRx Medizintechnik GmbH, Berlin, Germany). fNIRS optode placement was performed using an fNIRS cap (with different sizes) which ensures standardized sensor placement according to the well establishes 10-20 system, which is also used for EEG measurements. The fNIRS setup was similar to Piper et al. (2014) (see Figure 1C) using eight light sources and eight detectors with an inter-optode distance of 2.5-3 cm, which form 16 actual measurement channels. Data was acquired with a sampling frequency of 7.81 Hz.

Data Analyses
Quantitative Rating and Ranking For quantitative rating, we analyzed the 5-ball cascade only since this task is obviously the most demanding task within our study design for expert jugglers. Here, 10 s time frames within an activation block (5-ball cascade) were rated for the 3rd-6th block of each juggler. The balls are located in each frame of the video using a Hough Transform (Fernandes and Oliveira, 2012) for circle detection and a limited radius range. Next, associating temporal near ball detections automatically formed the tracks, based on distances in the image plane, color similarity, and temporal distance. During handover, where the hands come to occlude the balls and thereby the circle detection tends to fail, we used a template matching procedure (frame wise comparison) to extend the tracks by the information from previous and later images. For this, an image patch at a track end was used as a template and a similar looking match is searched for within a limited spatial search area around the last detection in the neighboring frame using normed cross correlation. Missing frames were manually tracked. The implementation of the described cinematographic approach uses OpenCV (Bradski, 2000). After this reconstruction we rated the individual trials with three kinematic parameters (Mapelli et al., 2012) to find the level of expertise for each juggler: tosses per second, velocity of the toss, angle of the toss (for left and right hand). For all three parameters, we calculated the average stability, that is, the coefficient of variability independent of the average. The averaged rank of the three behavior measures was calculated. The statistical analysis was performed using SPSS (Version 20.0 for Windows, SPSS Inc., Chicago, IL, USA). Calculating its intra-class correlation (ICC [3,1]) tested the reliability of this procedure with above-mentioned rating of an expert rater. We expect that the quantitative and qualitative rating show comparable ranking positions.
Hemodynamics fNIRS data was analyzed using the software package NIRSlab (v2014.05, NIRx Medical Technologies, New York, United States of America) and Wolfram Mathematica 9.0. Attenuation changes of both wavelengths (850 and 760 nm) were transformed to concentration changes of oxy-and deoxygenated hemoglobin (oxy-Hb and deoxy-Hb, respectively) using the modified Beer-Lambert approach (extinction coefficients: ε oxyHb λ 760 = 1.4865865 [ε oxyHb λ 850 = 2.526391[ l cm mmol ], ε deoxyHb λ 760 = 3.843707[ l cm mmol ], ε deoxyHb λ 850 = 1.798643[ l cm mmol ]; Differential pathlength Factor: 7.25; Cope and Delpy, 1988;Kocsis et al., 2006;Xu et al., 2014). The data was then band-pass filtered to attenuate high-frequency noise and cardiovascular artifacts using 0.01 Hz as high and 0.2 Hz as low pass cutoff frequencies (Huppert et al., 2009). Single trials were baseline corrected (regarding 5 s until stimulus onset). Though the experimental design included long breaks to prevent overlapping of the hemodynamic responses between trials, we directly analyzed the height of amplitude (baseline-corrected average of the temporal window from 5 to 20 s with regard to stimulus onset). Channels were rated as noisy when the first coefficient of the wavelet transform (4th order Daubechies Wavelet) had a higher inter-quartile-range than the 6th one. With this procedure 2.2% of the oxy-Hb and 8.7% of the deoxy-Hb channels were regarded as too noisy and not included in the further analysis (Nakano et al., 2009). Statistical analyses were performed using a one factor (with the six levels: 2-ball, left hand; 2-ball right hand; 3-ball cascade; 5-ball cascade and the two control tasks) repeated measurement ANOVA (within subject design, Greenhouse-Geisser correction; Greenhouse and Geisser, 1959). Partial eta squared of |>0.06| represents a "small" effect size, |>0.01| represents a "medium" effect size and |>0.14| represents a "large" effect size (Cohen, 1992). Post-hoc testing included pair-wise comparisons and Bonferroni correction (Dunn, 1961) for the number of tasks. The level of significance was set to 5% (p ≤ 0.05) for all tests.

Ethics Committee
There were no minors or persons with disabilities involved in the study and no additional considerations had to be made.

Qualitative Rating of the Level of Expertise
The expertise level between jugglers differed individually (see Table 1), especially for the maximum number of reliably handled balls, which fits well with previous studies investigating expert jugglers such as Gerber et al. (2014). Quantitative results provided by an expert rater (also an experienced juggler) are shown in Table 1 as well.

Quantitative Rating of the Level of Expertise
During the execution of the 5-ball cascade, three jugglers had no failure, two jugglers had one failures (dropping balls), five juggler had two failures, two jugglers had three failures and three jugglers had four failures. Finally, we collected eight valid trials for each participant. The cinematographic analysis revealed between 128 and 163 tosses for each participant (see Table 2). The analysis of the tosses per seconds resulted in an inter-subject average frequency of 4.1 ± 0.3 Hz. Average velocities at the toss were 7.6 ± 0.6 ms −1 (left hand) and 7.4 ± 0.6 ms −1 (right hand), and the average angles 83.0 ± 1.1 • (left hand) and 83.2 ± 1.3 • (right hand), respectively.

Subjective Rating of the Complexity of Tasks
Expert jugglers rated the control tasks as the ones with lowest complexity, followed by the 3-ball cascade, 2-ball juggling with the left hand, 2-ball juggling with the right hand, and the 5-ball cascade [Friedman-Test: X 2 (2, 30) = 36.66, p < 0.001, Table 3].

Task Related Differences
A repeated measures ANOVA with factor TASK (3-and 5-ball cascade, one handed 2-ball juggling [left and right hand], 1 and 2 Hz control tasks) revealed a significant increase in oxy-Hb as a function of task complexity for bilateral M1, S1, PMC and MT/V5 (for details see Figures 2, 3 and Table 3). More specifically, those tasks that were ranked as less complex were those that show smallest hemodynamic responses in oxy-Hb and vice versa. On the other hand, no such effects were observed for deoxy-Hb, except for left M1 (see Table 3).

Differences Related to Individual Level of Expertise
For the 5-ball cascade, the level of expertise did not correlate with the height of oxy-nor deoxy-Hb within M1, S1, PMC, and MT/V5 when correcting for multiple comparisons. Furthermore, we did a partial correlation analysis to control for the effect of movement speed. The results of the correlations analysis remained non-significant [ c oxy-Hb: r (9) = −0.37 − 0.43, p = 0.18 − 0.88; c deoxy-Hb: r (9) = −0.01 − 0.74, p = 0.15 − 0.99). However, there was a trend within right M1 toward a negative association between the hemodynamic response and the level of juggling expertise (oxy-Hb: M1 R ρ = 0.694, p = 0.012; deoxy-Hb: M1 R ρ = 0.718, p = 0.009; uncorrected for multiple comparisons, see Figure 4). More specifically, higher level of expertise tends toward lower hemodynamic responses within right M1 only.

DISCUSSION
Based on our findings, we confirmed that fNIRS is a capable tool to investigate online brain processes during the execution of certain motor tasks. Here, we extended previous findings by showing that the execution of a complex, sports-related motor task such as juggling is associated with neurovascular changes in sensorimotor (M1, S1, PMC) as well as visual motion areas (MT/V5). The novelty of the present study is that those brain regions seem to be modulated by task complexity, since a 5-ball cascade showed enhanced hemodynamic responses for oxy-Hb as compared to less complex tasks such as a 3-or 2-ball juggling pattern. However, in contrast to previous findings (Gerber et al., 2014) that showed an association between gray matter density and the level of juggling expertise in MT/V5, we here found only non-significant trend in M1 on a functional level.

Quantitative Rating of the Level of Expertise
The small inter-individual differences in the main tossing parameters release angle and release speed during 5-ball juggling confirms the expert level in our study population. However, the cinematographic analyses revealed individual differences in the variability of the parameters and thereby differences in the level of expertise between the jugglers. As the quantitative analyses reports similar ranking positions as compared to the subjective rating by an expert, our results emphasize the validity of the quantitative approach.

Hemodynamics
For the first time, according to our knowledge, we were able to show the adaptability of the human brain as a consequence of task complexity during the execution of a sports-related motor task. Subjectively rated levels of complexity for the different tasks are associated with the amount of hemodynamic response [oxy-Hb, for all areas p < 0.001; η 2 = 0.60 − 0.84] in sensorimotor and visual motion areas. So far, studies have only reported that higher movement speed when walking and running (Suzuki et al., 2004;Harada et al., 2009) is associated with higher hemodynamic responses (oxy-Hb). We extended these findings by showing that the hemodynamic response seems not only be affected by movement speed since control experiments (2 Hz alternating arm movements) indicated that the hemodynamic response was lower compared to 2-ball juggling. Furthermore, the hemodynamic response of the 3-ball cascade (∼3 Hz movements speed) was lower than 2-ball juggling (∼2 Hz movement speed). In contrast, with regard to deoxy-Hb, the majority of subjects in the present study did not show any significant alteration during the execution of ball juggling with different degrees of complexity. The detailed underlying mechanisms for the lack of effects on deoxy-Hb remain unknown and cannot be answered with our study design. However, it is tempting to speculate FIGURE 4 | Spearmann's rank correlations between the individual level of expertise as rated by the cinematographic video analysis of the 5-ball cascade and the changes in oxy-and deoxygenated hemoglobin (95% confidence intervals in brackets, uncorrected for multiple comparisons), PMC, premotor area; M1, primary motor area; S1, primary somatosensory area; MT/V5, visual area 5. The results in right M1 confirm: the higher the expertise, the lower the hemodynamic response ( 1 too much noisy data).
that the amplitude for deoxy-Hb in general is quite small and or less prone to alterations (induced by task complexity) as compared to oxy-Hb. It is interesting to see, however, that there was an "atypical" increase in deoxy-Hb during 5-ball juggling in the majority of jugglers (see Figure 2 and Table 3). The typical hemodynamic response to motor stimulation would include a reduction in deoxy-Hb concentration. The reason for this response is an increased regional blood flow (rCBF) in order to meet the higher demand in oxygen and thus higher "washout" of deoxy-Hb. In contrast, 5-ball juggling potentially causes a higher extraction of oxygen. This high cerebral metabolic rate of oxygen consumption (CMRO2) compensates the washout effect.
During the execution of 5-ball juggling we observed a trend toward a negative relationship between the 5-ball juggling expertise and alterations in oxy-and deoxy-Hb for right M1. This effect did not reach significance potentially due to the relative small number of participants included in the present study. However, even that trend is in line with previous studies investigating structural correlates of juggling expertise (Gerber et al., 2014). Similar findings of such brain-behavior relationships when learning to perform complex motor tasks have been shown on a structural as well as functional level (Taubert et al., 2010(Taubert et al., , 2011Gryga et al., 2012).

Limitations
Within each of the analyzed juggling and control tasks hemodynamic responses (oxy-Hb) in sensorimotor areas as well as in visual motion area MT/V5 were quite similar (see Table 3 and Figure 2). Hence, it is reasonable to assume that hemodynamic changes in oxy-Hb are not task related but due to a global effect caused by alterations in extra-cerebral blood flow (Kirilina et al., 2012). However, we believe that our findings are potentially not due to an increase in extra-cerebral blood flow, since we were able to show regional specific effects for deoxy-Hb. Deoxy-Hb is consistently less sensitive to alterations in extra-cerebral blood flow (Piper et al., 2014). Furthermore, the execution of unimanual 2-ball juggling in the left hand revealed higher neurovascular changes (see Table 3 and Figure 2, deoxy-Hb) in primary motor cortex on the contralateral hemisphere as compared to the ipsilateral hemisphere [t (14) = 3.19, p < 0. 01]. Testing the same difference during unimanual 2-ball juggling in the right hand fails to attain statistical significance [t (14) = -1.23, p = 0.24].
Since task complexity during juggling is associated with higher movement speed, which may influence hemodynamic responses (Obrig et al., 1996), we performed two control experiments. Here, we were able to show, that hemodynamic responses during the execution of 2-ball juggling compared to a 2 Hz control task, which involves a similar movement frequency, were higher. We interpret this finding as a task complexitydependent effect. Another potential limitation of the present study is that we only controlled for movement frequency for the 2-ball juggling task and not for 3-and 5-ball juggling, which obviously has to be performed with higher movement speed up to 4 Hz for 5-ball juggling. Hence, we cannot totally rule out the possibility that the observed effects (modulation of oxy-Hb with increasing complexity) rely on movement frequency.
Additionally, comparisons in our study design were made between uni-(2-ball juggling) and bimanual juggling tasks (3-and 5-ball cascade), a factor that might influence the interpretation of the observed effects as well. More specifically, one could argue that the observed effects are not related to task complexity but displays a simple difference between uni-and bilateral tasks. However, this is unlikely since 3-ball cascades showed less oxy-Hb alterations as compared to 2 ball (unilateral) juggling.