The prefrontal cortex (PFC), often described as the center of cognitive function, is involved in Executive Functions (EFs) that include cognitive processes such as problem solving, decision making, working memory, planning, inhibition of responses, and cognitive flexibility. In addition, the PFC also contributes to motor functions. The activation in the PFC has been linked to motor timing, movement selection, and control of gait    . Motor timing involves determining the appropriate time for movement initiation  . Movement selection helps to suppress automatically triggered responses and prepare for a proper reaction to a task  . Examining the activity in the PFC can help us to better understand its role in the neural basis of motor control and this knowledge can be applied to populations with motor dysfunctions. In the current study, task-induced activity in the PFC was measured during a simple motor task and an EF task to examine the role of the PFC in each task.
Many tasks have been used to assess EFs of the PFC. The Tower of Hanoi is a commonly used tool to assess EFs  . The game consists of three vertical pegs and a number of disks graduated in size with the largest disk at the bottom. The object of the game is to move all the disks from a start state to a goal state while following two rules: 1) only one disk may be moved at a time, and 2) a disk may never be placed on top of a disk smaller than itself. The completion of the Tower of Hanoi requires complex EFs including inhibition, working memory, and planning   .
Due to the EF elements required in completion of the Tower of Hanoi, it has been used in patients with PFC lesions  . In the study, participants were asked to move disks to a goal state, in which all disks were stacked in descending order on a specific peg, from different start states. A performance score was calculated by a computer based on movement accuracy, movement speed, and problem difficulty. Patients with PFC lesions had lower score than the control group but the strategy they used to solve the puzzle was similar to the one used by the controls. The finding suggested that these patients had deficits in inhibition control rather than planning. Nonetheless, the contribution of the PFC during the Tower of Hanoi has been verified by studies using functional Magnetic Resonance Imaging (fMRI)  and electroencephalogram (EEG)   . Increased prefrontal-parietal intrahemispheric correlation was observed during the Tower of Hanoi task in both left and right hemispheres  . Changes in absolute power of EEG activity in the PFC were also observed during the Tower of Hanoi  . However, limited resources can be found on using functional near-infrared (fNIR) spectroscopy to examine brain activity during the Tower of Hanoi.
fNIR measures the hemodynamics changes in cortical function based on the optical properties of hemoglobin. It provides a balance between spatial and temporal resolution for monitoring local and instance changes of oxygenation in the cerebral cortex     . In addition, fNIR is more tolerant to motion artifacts than other non-invasive neuroimaging technologies making its feasibility enhanced for clinical studies, especially in motor tasks   .
The overall goal of our research is to understand brain/behavior relationships in motor tasks requiring executive function. In this study, we wanted to examine PFC activation in a computerized Tower of Hanoi task and in a simple tapping task with similar motor and spatial requirements but no problem-solving component. By controlling for movement type and amplitude, we could compare activation patterns which may inform us on the relative contributions of the PFC to each task. We had two purposes: first, we wanted to compare levels of activation among TOH, tapping and a resting condition. We hypothesized that TOH would have significantly higher activation levels than the other two conditions. Second, we wanted to compare the sensitivity of four biomarkers in detecting differences among conditions. We made no a priori prediction about sensitivity.
Prior to initiation of this study, the protocol was approved by the Institutional Review Board of the institution. A total of 18 participants participated in this study. All participants were students at the University of Delaware recruited through flyers posted around the campus and through word of mouth. Prior to participation, they were screened for suitability for this research and all were healthy and without family history of cognitive or mental disorders. Data of two participants was removed because one fell asleep during the study and high motion artifacts in another participant’s data, leaving a total of 16 adult participants (7 females and 9 males, 22.3 ± 3.2 years old) in the study. There was no a priori power analysis conducted for this work as the dependent measures would provide sufficient information on variability to calculate effect sizes for follow-up studies. Data were collected between April and September, 2014.
Upon entering the laboratory, participants received a description of the research and were asked to provide written informed consent. The principal investigator instructed participants how to perform the tasks and participants practiced the Tower of Hanoi on an iPad until they understood the rules and were able to complete the 3-disk task once. Next, the fNIR sensor pad was applied on the forehead of participants. During the resting period, the participants were instructed to “Try to relax and do nothing” and “Do not fall asleep” and participants stared at a “+” sign (1 cm * 1 cm) in the middle of a monitor on a “15.6” laptop (Satellite L755-S5351, Toshiba). This was followed by performance of the tasks in a randomly assigned order (Figure 1). An auditory sound was provided to notify participants of the start and end of each condition. Participants were instructed to keep quiet and minimize their movements during data collection. Data collection was completed within 40 minutes for all participants.
2.3.1. Tower of Hanoi
The Tower of Hanoi is a disk-transfer task. The goal of the task is to move disks from a
Figure 1. Data collection procedure.
start state to a goal state using the fewest number of moves possible while following two rules: 1) only one disk may be moved at a time, and 2) a disk may never be placed on top of a disk smaller than itself. Tower Hanoi task in our study required participants to move disks from the leftmost peg to the rightmost peg. Participants performed the task on an iPad using Tap Towers (Madcap Studios, Inc.), starting from 3 disks (Figure 2) then progressed to 4 disks, 5 disks, and 6 disks if time allowed. The level participants reached and the number of moves to achieve goal states of each level were recorded.
The tapping task was designed to be a simple motor task with minimal EF demands. Participants tapped two sides of iPad (about 5 cm in width and 13 cm apart, index of difficulty = 2.4  ) back and forth on an iPad at a self-selected pace for two minutes. The number of taps was recorded.
2.4. Functional Near Infrared (fNIR) Spectroscopy
The activity in the prefrontal cortex was measured by an fNIR system (fNIR100A, Biopac System Inc.). fNIR detects hemodynamic changes resulted from brain activation based on the general transparency of biological tissues to infrared light as well as optical properties of oxygenated-hemoglobin (oxy-Hb) and deoxygenated-hemoglobin (deoxy-Hb)  . The fNIR system has a sensor pad, a control box, and a computer installed with Cognitive Optical Brain Imaging (COBI) Studio software and fnirSoft  for data collection and analysis. The sensor pad contains 4 light sources (LED) and 10 light detectors creating an array of 16 channels (Figure 3). The source-detector distance is 2.5 cm and the depth of penetration is approximately 1.5 cm beneath the scalp   . The light sources and detectors are embedded in a flexible pad made with medical grade silicone. The pad can be worn as a headband and match the contour of the foreheads. The light sources emit infrared light at 730 nm and 850 nm which penetrates the scalp. Part of light is absorbed by hemoglobin and the rest of light reaches detectors in a banana-shaped path. Oxygenated-hemoglobin and deoxygenated-he- moglobin has different absorption coefficients. Concentrations of hemoglobin can then be calculated by the ratio of light absorbed at different wavelengths. LED current and detector gain were adjusted before prior to the study to prevent signal saturation Base-
Figure 2. The Tower of Hanoi task using three disks: (a) starting position (b) ending position.
Figure 3. The sensor pad with 8 light detectors and 4 light sources (left, top) creates 16 channels (left, bottom). The sensor pad positioned on a participant (right).
line measurements were obtained with the fNIR signal recorded at 2 Hz.
2.5. Data Processing
Raw data was visually inspected and channels with saturated or very low signals were removed. Optodes on the sensor pad that may not contact completely with the skin resulted in saturated channels. Low signal values were due to interference with light transmission from the hair. Motion artifacts were excluded manually. Raw light intensity was filtered by a low-pass finite impulse response (FIR) filter with an order of 20 and cutoff frequency at 0.1 Hz to attenuate the high frequency noise, respiration, and cardiac cycle effects  . Relative changes of concentrations in oxygenated hemoglobin (Δoxy-Hb), deoxygenated hemoglobin (Δdeoxy-Hb), blood oxygenation (Δoxy; calculated from Δoxy-Hb-Δdeoxy-Hb) and blood volume changes (ΔHbT; calculated as ∑(Δoxy-Hb + Δdeoxy-Hb) were calculated using modified Beer-Lambert law, with a 10- second baseline recorded at the beginning of data collection  . Data was processed using fnir Soft  . The first and the last 15 seconds of each condition were removed to prevent an unstable state at the beginning of the test and the effect from previous activity  . The remaining data points were averaged in each channel for each task. Z scores were calculated using the mean and standard deviation of all 16 channels in each task of each participant using z = (x-mean)/SD. Channels with extreme values (z scores above 2 or less than-2) were removed as outliers  . Remaining channels were averaged for each condition. Four regions of interest (ROIs) were created for the prefrontal cortex which included averages of 4 channels for each biomarker: left dorsolateral prefrontal cortex (LDLPFC) which involved channels 1 - 4; left medial prefrontal cortex (LMPFC) channel 5 - 8 were used; right medial prefrontal cortex (RMPFC) used channels 9 - 12 and right dorsolateral prefrontal cortex (RDLPFC) involved channels 13 - 16.
A one-way repeated measures ANOVA was used to examine the differences in the four biomarkers (i.e., Δoxy-Hb, Δdeoxy-Hb, Δoxy, and ΔHbT) between the eyes open resting, tapping, and the Tower of Hanoi tasks. Tukey’s post-hoc comparisons were calculated to determine specific differences. Effect sizes (Cohen’s d) were calculated to aid in the interpretation of the results. We interpreted the values of d = 0.20 as considered small, d = 0.50 as medium, d = 0.80 as large and d = 1.3 as very large according to Cohen  . The level of significance was set at 0.05. Number Cruncher Statistical Systems (NCSS 9-www.ncss.com) was used for the analyses.
3.1. Tower of Hanoi Behavioral Data
All participants completed 3-disk level in 2 minutes (9.9 ± 4.3 moves). Eleven participants finished at the 4-disk level (25.1 ± 9.3 moves) and 3 participants completed the 5-disk level (39.3 ± 4.9 moves) within the 2-minute time limit.
3.2. PFC Activation
The results of repeated measures ANOVA showed significant differences in all 4 biomarkers in RDLPFC among the three conditions. Tukey’s post-hoc comparisons revealed significant increases in both Δoxy and Δoxy-Hb from eyes open resting condition and tapping condition to the Tower of Hanoi condition. A significant increase of ΔHbT was found from tapping to the Tower of Hanoi condition. These changes along with a significant decrease in Δdeoxy-Hb from eyes open resting to the Tower of Hanoi condition indicated that there was a higher activation in RDLPFC when participants were performing the Tower of Hanoi task. In addition, an increased activity in LDLPFC was observed with significant higher Δoxy and Δoxy-Hb during the Tower of Hanoi task. A significant change in Δoxy was also found in RMPFC among 3 conditions. Figure 4 shows the average areas of activation for the three conditions. Note that in the TOH task, there are high levels of activity in the DLPFC on both the left and right sides.
Detailed results were presented in Table 1.
The aim of current study was to examine the differences in prefrontal cortex activity
Figure 4. Comparison of areas of activation associated with a) resting, b) tapping and c) tower of hanoi condition. Areas of significant activation (t (15) = 1.756, one-tailed) are denoted in red and yellow, and areas with limited or no activation are denoted in cyan and blue.
Table 1. Four biomarkers descriptive and inferential statistics (N = 16).
between an executive function task (Tower of Hanoi) and a tapping task with identical movement and spatial requirements using a fNIR system. Four fNIR biomarkers were measured during a simple motor task and Tower of Hanoi in adults. Our results found increased activities in bilateral DLPFC during Tower of Hanoi task. In contrast, there were no significant changes in brain activity during the tapping task as compared to the resting condition. The finding corresponded with previous studies exploring the involvement of the PFC in problem solving tasks using fMRI, in which activations of the DLPFC have been observed during Tower of London task     . The involvement of the DLPFC in executive functions including working memory, inhibition, and planning has been long established   . It is thought that the left DLPFC supports cognitive processes of recognizing specific features, reasoning, and problem solving, while the right DLPFC supports cognitive processes of planning and decision making    . Our results showing bilateral activation of DLPFC matched the signatures of Tower of Hanoi that require functions of DLPFC from both hemispheres.
In the current study, bilateral activation of DLPFC during Tower of Hanoi task was reflected by two biomarkers, Δoxy and Δoxy-Hb. This was not the case with Δdeoxy- Hb or ΔHbT which only had statistically significant changes in the right hemisphere. Furthermore, the associated effect sizes for Δoxy and Δoxy-Hb ranged from d = 0.44 - 0.76 which indicate a medium effect. In contrast, effect sizes for the measures of Ddeoxy were small. Our results suggested that Δoxy and Δoxy-Hb were more sensitive in event-related response than Δdeoxy-Hb or ΔHbT. To date, most commonly reported fNIR biomarkers are oxy-Hb, deoxy-Hb, and total-Hb. Limited studies exist comparing the sensitivity of Δoxy to other biomarkers. To our knowledge, the changes of Δoxy has not been compared to other neuroimaging methods (e.g., EEG, fMRI). Prior studies examined the correlation between fNIR biomarkers (oxy-Hb, deoxy-Hb, and HbT) and blood oxygen level dependent (BOLD) response measured by fMRI had provided conflicting results. Although significant correlations exist between the BOLD signal and all fNIR biomarkers, different biomarkers were reported to have the strongest correlation to BOLD signals   . Strangman and colleagues correlated the event-related changes of BOLD signals to quantitative changes of oxy-Hb, deoxy-Hb, and total-Hb measured by fNIR during a motor task  . They found that oxy-Hb had the strongest correlation to BOLD signal and deoxy-Hb had the weakest correlation to BOLD signal. The result supported part of our finding that oxy-Hb was more sensitive than deoxy- Hb in response to an executive function task. Conversely, Huppert and colleagues examined the temporal correlation between the BOLDresponse measured by fMRI and fNIR biomarkers during a motor activity in healthy adults  . The result showed an improved correlation between deoxy-Hb and a BOLD signal than the correlation between oxy-Hb and a BOLD signal. Combining the results from both studies, it suggested that sensitivity of each fNIR biomarker varies in different dimensions of the hemodynamic response. Deoxy-Hb is more sensitive to temporal changes whereas oxy- Hb is more sensitive in quantitative measurements. Our results were closer to Strangman et al.  since quantitative changes were examined in our study.
A limitation of the current study was that we did not analyze temporal changes of fNIR data. We were not able to explore the time course of hemodynamic changes in response to the Tower of Hanoi task. In addition, cortical activity was not recorded during the orientation session nor was a control baseline session recorded. Assessing during an orientation and during a control session might have provided valuable information about cognitive processing and neural activation.
5. Conclusions and Future Directions
In the current study, we found a significant increase in bilateral DLPFC activation during the Tower of Hanoi, an EF task, over a simple tapping task and a rest condition in our sample of adults. The results confirmed the involvement of the PFC in the EF task. In addition, we found that Δoxy and Δoxy-Hb were more sensitive than Δdeoxy-Hb and ΔHbT. Further, Δoxy may be a more valuable measure when an event-related response is examined.
The fact that we could detect activation differences between TOH and tapping suggests that this methodology can be used in the future with clinical populations who have movement difficulties. In particular, we plan to replicate this study by using a population of children with Developmental Coordination Disorder (DCD), which is a neuromotor disability that affects motor coordination and the capacity to learn and acquire motor skills. While the etiology of DCD has yet to be fully understood, previous research has suggested associations between motor deficits and cognitive processing (e.g.  ). The movement requirements in our existing protocol are relatively simple and have few motor planning requirements, thus enabling participants with poor motor coordination to successfully perform the task. Our proposed future work will allow us to look for EF differences in populations with and without DCD.
 Bortoletto, M. and Cunnington, R. (2010) Motor Timing and Motor Sequencing Contribute Differently to the Preparation for Voluntary Movement. Neuroimage, 49, 3338-3348.
 Koenraadt, K.L., Roelofsen, E.G., Duysens, J. and Keijsers, N.L. (2014) Cortical Control of Normal Gait and precision Stepping: An fNIRS Study. Neuroimage, 85, 415-422.
 Krams, M., Rushworth, M.F., Deiber, M.P., Frackowiak, R.S. and Passingham, R.E. (1998) The Preparation, Execution and Suppression of Copied Movements in the Human Brain. Experimental Brain Research, 120, 386-398.
 Goel, V. and Grafman, J. (1995) Are the Frontal Lobes Implicated in “Planning” Functions? Interpreting Data from the Tower of Hanoi. Neuropsychologia, 33, 623-642.
 Welsh, M.C. and Huizinga, M. (2005) Tower of Hanoi Disk-Transfer Task: Influences of Strategy Knowledge and Learning on Performance. Learning and Individual Differences, 15, 283-298.
 Welsh, M.C., Satterlee-Cartmell, T. and Stine, M. (1999) Towers of Hanoi and London: Contribution of Working Memory and Inhibition to Performance. Brain and Cognition, 41, 231-242.
 Fincham, J.M., Carter, C.S., van Veen, V., Stenger, V.A. and Anderson, J.R. (2002) Neural Mechanisms of Planning: A Computational Analysis Using Event-Related fMRI. Proceedings of the National Academy of Sciences of the United States of America, 99, 3346-3351.
 Guevara, M.A., Hernandez Gonzalez, M., Rizo Martinez, L.E. and Robles Aguirre, F.A. (2013) Performance of the Towers of Hanoi Task and Cortical Electroencephalographic Power Changes Associated with Infancy, Adolescence, and Early Adulthood. Experimental Brain Research, 231, 315-324.
 Guevara, M.A., Rizo Martinez, L.E., Robles Aguirre, F.A. and Hernandez Gonzalez, M. (2012) Prefrontal-Parietal Correlation during Performance of the Towers of Hanoi Task in Male Children, Adolescents and Young Adults. Developmental Cognitive Neuroscience, 2, 129-138.
 Babiloni, C., Pizzella, V., Gratta, C.D., Ferretti, A. and Romani, G.L. (2009) Fundamentals of Electroencefalography, Magnetoencefalography, and Functional Magnetic Resonance Imaging. International Review of Neurobiology, 86, 67-80.
 Chance, B., Anday, E., Nioka, S., Zhou, S., Hong, L., Worden, K., Li, C., Murray, T., Ovetsky, Y., Pidikiti, D. and Thomas, R. (1998) A Novel Method for Fast Imaging of Brain Function, Non-Invasively, with Light. Optics Express, 2, 411-423.
 Dieler, A.C., Tupak, S.V. and Fallgatter, A.J. (2012) Functional Near-Infrared Spectroscopy for the Assessment of Speech Related Tasks. Brain and Language, 121, 90-109.
 Lloyd-Fox, S., Blasi, A. and Elwell, C.E. (2010) Illuminating the Developing Brain: The Past, Present and Future of Functional Near Infrared Spectroscopy. Neuroscience & Biobehavioral Reviews, 34, 269-284.
 Ayaz, H., Izzetoglu, M., Shewokis, P.A. and Onaral, B. (2010) Sliding-Window Motion Artifact Rejection for Functional Near-Infrared Spectroscopy. 2010 Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Buenos Aires, 31 August-4 September 2010, 6567-6570.
 Izzetoglu, M., Chitrapu, P., Bunce, S. and Onaral, B. (2010) Motion Artifact Cancellation in NIR Spectroscopy Using Discrete Kalman Filtering. Biomedical Engineering Online, 9, 16.
 Fitts, P.M. (1954) The Information Capacity of the Human Motor System in Controlling the Amplitude of Movement. Journal of Experimental Psychology, 47, 381-391.
 Okada, E. and Delpy, D.T. (2003) Near-Infrared Light Propagation in an Adult Head Model. II. Effect of Superficial Tissue Thickness on the Sensitivity of the Near-Infrared Spectroscopy Signal. Applied Optics, 42, 2915-2922.
 Okada, E. and Delpy, D.T. (2003) Near-Infrared Light Propagation in an Adult Head Model. I. Modeling of Low-Level Scattering in the Cerebrospinal Fluid Layer. Applied Optics, 42, 2906-2914.
 Izzetoglu, M., Izzetoglu, K., Bunce, S., Ayaz, H., Devaraj, A., Onaral, B. and Pourrezaei, K. (2005) Functional Near-Infrared Neuroimaging. IEEE Transactions on Neural Systems and Rehabilitation Engineering, 13, 153-159.
 Chance, B., Zhuang, Z., UnAh, C., Alter, C. and Lipton, L. (1993) Cognition-Activated Low-Frequency Modulation of Light Absorption in Human Brain. Proceedings of the National Academy of Sciences of the United States of America, 90, 3770-3774.
 Izzetoglu, M., Bunce, S.C., Izzetoglu, K., Onaral, B. and Pourrezaei, K. (2007) Functional Brain Imaging Using Near-Infrared Technology. IEEE Engineering in Medicine and Biology Magazine, 26, 38-46.
 Van Selst, M. and Jolicoeur, P. (1994) A Solution to the Effect of Sample-Size on Outlier Elimination. The Quarterly Journal of Experimental Psychology Section A-Human Experimental Psychology, 47, 631-650.
 Boghi, A., Rampado, O., Bergui, M., Avidano, F., Manzone, C., Coriasco, M., Mortara, P., Orsi, L., Ropolo, R. and Bradac, G.B. (2006) Functional MR Study of a Motor Task and the Tower of London Task at 1.0 T. Neuroradiology, 48, 763-771.
 Newman, S.D., Carpenter, P.A., Varma, S. and Just, M.A. (2003) Frontal and Parietal Participation in Problem Solving in the Tower of London: fMRI and Computational Modeling of Planning and High-Level Perception. Neuropsychologia, 41, 1668-1682.
 Ruh, N., Rahm, B., Unterrainer, J.M., Weiller, C. and Kaller, C.P. (2012) Dissociable Stages of Problem Solving (II): First Evidence for Process-Contingent Temporal Order of Activation in Dorsolateral Prefrontal Cortex. Brain and Cognition, 80, 170-176.
 Wagner, G., Koch, K., Reichenbach, J.R., Sauer, H. and Schlosser, R.G. (2006) The Special Involvement of the Rostrolateral Prefrontal Cortex in Planning Abilities: An Event-Related fMRI Study with the Tower of London Paradigm. Neuropsychologia, 44, 2337-2347.
 Moriguchi, Y. and Hiraki, K. (2013) Prefrontal Cortex and Executive Function in Young Children: A Review of NIRS Studies. Frontiers in Human Neuroscience, 7, 867.
 Barbey, A.K., Krueger, F. and Grafman, J. (2009) An Evolutionarily Adaptive Neural Architecture for Social Reasoning. Trends in Neurosciences, 32, 603-610.
 Huey, E.D., Krueger, F. and Grafman, J. (2006) Representations in the Human Prefrontal Cortex. Current Directions in Psychological Science, 15, 167-171.
 Huppert, T.J., Hoge, R.D., Diamond, S.G., Franceschini, M.A. and Boas, D.A. (2006) A Temporal Comparison of BOLD, ASL, and NIRS Hemodynamic Responses to Motor Stimuli in Adult Humans. Neuroimage, 29, 368-382.
 Strangman, G., Culver, J.P., Thompson, J.H. and Boas, D.A. (2002) A Quantitative Comparison of Simultaneous BOLD fMRI and NIRS Recordings during Functional Brain Activation. Neuroimage, 17, 719-731.
 Asonitou, K., Koutsouki, D., Kourtessis, T. and Charitou, S. (2012) Motor and Cognitive Performance Differences between Children with and without Developmental Coordination Disorder (DCD). Research in Developmental Disabilities, 33, 996-1005.