Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Effects of Increasing Neuromuscular Electrical Stimulation Current Intensity on Cortical Sensorimotor Network Activation: A Time Domain fNIRS Study

  • Makii Muthalib ,

    Contributed equally to this work with: Makii Muthalib, Rebecca Re

    makii.muthalib@univ-montp1.fr

    Affiliations Movement To Health (M2H), EuroMov, University of Montpellier, Montpellier, France, Movement Neuroscience, IHBI, Queensland University of Technology, Brisbane, Australia

  • Rebecca Re ,

    Contributed equally to this work with: Makii Muthalib, Rebecca Re

    Affiliation Dipartimento di Fisica, Politecnico di Milano, Milan, Italy

  • Lucia Zucchelli,

    Affiliation Dipartimento di Fisica, Politecnico di Milano, Milan, Italy

  • Stephane Perrey,

    Affiliation Movement To Health (M2H), EuroMov, University of Montpellier, Montpellier, France

  • Davide Contini,

    Affiliation Dipartimento di Fisica, Politecnico di Milano, Milan, Italy

  • Matteo Caffini,

    Affiliation Dipartimento di Elettronica, Informazione e Bioingegneria, Politecnico di Milano, Milan, Italy

  • Lorenzo Spinelli,

    Affiliation Istituto di Fotonica e Nanotecnologie, CNR, Politecnico di Milano, Milan, Italy

  • Graham Kerr,

    Affiliation Movement Neuroscience, IHBI, Queensland University of Technology, Brisbane, Australia

  • Valentina Quaresima,

    Affiliation Dipartimento di Medicina Clinica, Sanità Pubblica, Scienze della Vita e dell'Ambiente, Università degli Studi dell'Aquila, L’Aquila, Italy

  • Marco Ferrari,

    Affiliation Dipartimento di Medicina Clinica, Sanità Pubblica, Scienze della Vita e dell'Ambiente, Università degli Studi dell'Aquila, L’Aquila, Italy

  • Alessandro Torricelli

    Affiliation Dipartimento di Fisica, Politecnico di Milano, Milan, Italy

Abstract

Neuroimaging studies have shown neuromuscular electrical stimulation (NMES)-evoked movements activate regions of the cortical sensorimotor network, including the primary sensorimotor cortex (SMC), premotor cortex (PMC), supplementary motor area (SMA), and secondary somatosensory area (S2), as well as regions of the prefrontal cortex (PFC) known to be involved in pain processing. The aim of this study, on nine healthy subjects, was to compare the cortical network activation profile and pain ratings during NMES of the right forearm wrist extensor muscles at increasing current intensities up to and slightly over the individual maximal tolerated intensity (MTI), and with reference to voluntary (VOL) wrist extension movements. By exploiting the capability of the multi-channel time domain functional near-infrared spectroscopy technique to relate depth information to the photon time-of-flight, the cortical and superficial oxygenated (O2Hb) and deoxygenated (HHb) hemoglobin concentrations were estimated. The O2Hb and HHb maps obtained using the General Linear Model (NIRS-SPM) analysis method, showed that the VOL and NMES-evoked movements significantly increased activation (i.e., increase in O2Hb and corresponding decrease in HHb) in the cortical layer of the contralateral sensorimotor network (SMC, PMC/SMA, and S2). However, the level and area of contralateral sensorimotor network (including PFC) activation was significantly greater for NMES than VOL. Furthermore, there was greater bilateral sensorimotor network activation with the high NMES current intensities which corresponded with increased pain ratings. In conclusion, our findings suggest that greater bilateral sensorimotor network activation profile with high NMES current intensities could be in part attributable to increased attentional/pain processing and to increased bilateral sensorimotor integration in these cortical regions.

Introduction

Neuromuscular electrical stimulation (NMES) has been clinically demonstrated to improve movements of neurological populations with motor disabilities, such as post-stroke and incomplete spinal cord injury [1]. The principle of NMES is to apply repeated electrical currents to the peripheral motoneuronal axonal branches overlying the muscle of interest at a stimulation intensity that elicits muscle contractions. This in turn leads to the well characterised peripheral neuromuscular adaptations such as increased muscle strength and oxidative capacity [25]. More recently, it has been reported that rapid supraspinal central nervous system (CNS) reorganization/neuroplastic mechanisms are implicated during NMES [68]. This is because NMES activates peripheral sensory neuronal axons that send proprioceptive (and pain) afferent signals from the stimulated muscle to the CNS leading to cortical neuroplasticity and improvements in voluntary activation [6, 7] through a process of sensorimotor integration [7]. In addition, NMES current intensity above motor threshold is an important factor in determining the rapid cortical neuroplastic changes. In particular, NMES current intensities above motor threshold increases corticospinal excitability while NMES current intensities at sensory threshold decreases corticospinal excitability [7, 9].

Cortical correlates of NMES-evoked movements have been demonstrated in a few functional magnetic resonance imaging (fMRI) studies [1014]. For instance, Blickenstorfer et al. (2009) observed a cortical activation pattern including the contralateral primary motor (M1) and sensory (S1) cortex (i.e., sensorimotor cortex-SMC) and bilateral secondary somatosensory area (S2), supplementary motor area (SMA), premotor cortex (PMC), and regions of the prefrontal cortex (PFC) during NMES-evoked wrist extension/flexion movements at current intensities just above the individual motor threshold (9–23mA) but minimising pain discomfort. Smith et al. [12] observed a positive relationship between increases in sensorimotor network activation and increases in NMES current intensity up to the motor threshold. Since NMES current intensities well above the motor threshold and up to the maximum tolerated current intensity (MTI) have been commonly utilized in strength and clinical neurorehabilitation programmes [5, 6], neuroimaging studies are necessary to investigate the cortical correlates of increasing NMES current intensity at the MTI on healthy subjects and patients.

A neuroimaging modality applicable in these studies is represented by functional near-infrared spectroscopy (fNIRS). This technique, exploiting the hemoglobin absorption spectra in the NIR range, gives a non-invasive measure of the hemodynamic changes in tissues such as the cortical microcirculation blood vessels (for review see [15, 16]). If compared with the other neuroimaging methods, fNIRS represents an optimal brain imaging monitoring tool given that it does not require stringent physical and motor constraints compared to fMRI. The principle of fNIRS lies in the mechanism of neurovascular coupling [17]. Cerebral blood flow adequate for brain activity and metabolic demand is maintained through feedforward and feedback processes of autoregulation and neurovascular coupling [17]. Thus when a specific brain region is activated, cerebral blood flow increases in a temporally and spatially coordinated manner tightly linked to changes in neural activity through a complex sequence of coordinated events involving neurons, glia, arteries/arterioles, and signalling molecules [17]. Therefore, fNIRS infers changes in neural activity that is mirrored by changes in blood oxygenation in the region of the activated cortical area (i.e., increase in O2Hb and decrease in HHb) [18].

Multichannel fNIRS neuroimaging has shown unilateral hand motor tasks activate regions of the cortical sensorimotor network, including the contralateral SMC, PMC/SMA, and PFC [19], which is consistent with fMRI neuroimaging [2022]. Multichannel fNIRS neuroimaging has also shown an increase in activation of SMC and S2 regions during median nerve electrical stimulation just below the motor threshold current intensity (3–10mA) but minimising eliciting painful sensations [23, 24]. To the best of our knowledge, only our previous fNIRS study [25] has investigated the contralateral PFC activation induced by NMES of the elbow flexors (biceps bachii and brachioradials muscles) at increasing current intensities from well above motor threshold (45–55mA), which evoked strong muscle contractions with some pain perception, to MTI levels (55–70mA). That study found a relationship between the NMES current intensity to the MTI levels and the levels of PFC activation, which was most likely due to increased attentional demands and pain processing. However, the effects of increasing NMES current intensity to MTI levels on activation of the cortical sensorimotor network and perceived pain levels associated to the evoked-sensory/systemic physiological responses are still unknown.

In fNIRS neuroimaging, before reaching the deeper cortical region and exiting from the head, NIR photons must travel through the overlaying superficial layers (i.e. skin, subcutaneous tissue, aponeurosis, connective tissue, periosteum, cranium, meninges, cerebro-spinal fluid) [26]. Thus, the measured O2Hb and HHb concentration changes are a mixture of hemodynamic responses occurring in the cortical region and interfering responses occurring in these overlaying superficial (extra-cerebral) layers [27]. Recent reports have raised a question against the assumption that O2Hb/HHb changes measured by continuous wave (CW) fNIRS actually originated only from the cortical hemodynamic response (for review see [2729]). Although several methods have been proposed to separate cortical and extra-cranial components in CW fNIRS signals (for review see [28, 30]), no consensus has been reached yet on the best approach. Since blood flow increase in the scalp can be expected during NMES we used the time domain (TD) fNIRS technique to exploit its ability in separating the superficial from the cortical layer fNIRS signals (for review see [31]). The ability of TD fNIRS to discriminate hemodynamic signals from intra- and extra-cerebral (i.e. cortical and superficial) layers has been already demonstrated through numerical simulations, experiments on tissue phantoms, and in vivo measurements [28, 3241].

Although there have been several neuroimaging studies investigating neural activity accompanying NMES-evoked movements, these previous studies adopted low electrical stimulation not exceeding motor threshold, and hence the effects of high NMES intensities on neural activity, that is used practically in clinical neurorehabilitation programs, remains to be clarified. Such study will contribute to identify more effective neurorehabilitation programs. Therefore, the aim of the present study was to utilise a multi-channel TD fNIRS system to map the bilateral cortical sensorimotor network activation profile with increasing NMES current intensities from above motor threshold to the MTI, and with reference to the activation profile during VOL wrist extension movements. Considering the SMC, PMC/SMA, S2 involvement in sensorimotor integration, we hypothesized that greater cortical sensorimotor network (SMC, PMC, SMA, S2) activation would be induced with increased NMES current intensities; and considering the PFC involvement in attention and pain processing, we further hypothesized that the PFC would also show greater pain related activation at NMES current intensities at MTI levels.

Materials and Methods

Nine male healthy volunteers (age: 39.2±13.0 y, height: 179.1±7.4 cm, weight: 81.7±15.7 kg) participated in this study. All subjects had no known health problems, no history of neurological disorders, and no upper extremity muscle or joint injuries. All subjects were right handed according to the Edinburgh handedness inventory [42].

Ethics Statement

The study was in accordance with the Declaration of Helsinki and was approved by the Institutional Review Board at the Department of Physics, Politecnico di Milano. Written informed consent was obtained from all subjects.

Experimental procedure

The study was conducted in a quiet and dimly lit room. Subjects were seated in a comfortable chair with both arms on a table. The right wrist extensor muscles (i.e., extensor carpi radialis longus, extensor carpi radialis brevis, extensor digitorum communi, and extensor carpi ulnaris) were stimulated (biphasic symmetrical rectangular pulse shapes at 30 Hz and 200 μs pulse width) with a pair of 5x5 cm self-adhesive electrodes (CEFAR Physio 5, DJO France SAS, Mouguerre, France). One electrode was located on the motor point of the right wrist extensor muscles, the other one on the distal end of the muscle near the wrist. Each subject’s MTI was first determined by using a series of 6–8 brief (3–5 s) electrically stimulated contractions with increasing current intensity. After each increase in current intensity, the participants were asked to report their tolerance to further increases in current intensity. Then, MTI was defined as the intensity of stimulation received when the subjects were unable to tolerate an increase in current intensity. Three attempts were made at an individual level in order to ensure a consistent MTI determination and the group-range of MTIs was between 26–50 mA with a group-average MTI of 35.4±6.8 mA. The fNIRS experimental session started at least 10 minutes after determining the individual MTI.

Each subject was asked to perform 10 blocks of voluntary (VOL) contractions, each consisting in a series of ten 70° wrist extensions (visually monitored by the experimenter) 1-s duration interleaved with 1-s rest (0° wrist extension). Consecutive blocks were separated by a 20-s interval. Five min after the end of the VOL contractions, the subjects underwent NMES at four increasing levels of current intensity based on their individual MTI: 10%MTI, 50%MTI, MTI, and above MTI (MTI+). The NMES conditions consisted of 10 blocks of 10 stimulations (1-s stimulation, 1-s rest) with a 20-s interval between blocks. A 2-min rest period separated each of the 4 NMES conditions. For the MTI+ condition, the first block was performed at the MTI current level determined prior to the start of the experiment. During each of the subsequent 9 blocks, each subject was asked after the 5th stimulation if a further current increase could be tolerated. In the case of an affirmative response, the current was increased by 1 mA increments until the subject indicated their new tolerance level.

Pain/discomfort ratings

At the end of each experimental condition subjective pain was measured using a pain rating scale (PRS) [43] and discomfort ratings using a visual analogue scale (VAS) [44]. Participants were asked to verbally rate the level of pain on the PRS from 0 (no pain) to 12 (extremely painful), and to manually mark the level of discomfort on the VAS, consisting of a 10-cm line with”no discomfort” on one end and”intolerable sensation” on the other end.

Physiological measurements

During the experiment, physiological and autonomic parameters such as heart rate (HR), respiratory rate (RR) and skin conductance (SC) were continuously monitored at a 256 Hz sampling rate with a commercial physiological sensing system (FlexComp Infiniti System T7555M, Thought Technology Ltd., Montreal, Canada). The electrocardiogram electrodes for HR measurements were positioned with a standard configuration with the negative electrode placed on the right shoulder, the positive electrode placed on the lower center or left side of the chest (xiphoid process) and the ground electrode on the left shoulder. A respiration band for RR measurements was placed over the subject’s clothing at the mid chest level. SC was measured between the middle and the ring fingers of the left hand. Subjects were asked to keep their hand as still as possible to minimize artefacts in the SC readings. To observe variations of physiological parameters across conditions, HR, RR and SC measurements were averaged at an individual level within each experimental block.

TD fNIRS instrumentation

The fNIRS measurements were performed by using a multi-channel dual wavelength TD fNIRS system developed at the Department of Physics, Politecnico di Milano, Italy [45, 46]. Two pulsed diode lasers (operating at 690 and 829 nm, with 80 MHz repetition rate) were used as the NIR light source. The NIR light pulses were time multiplexed over 10 locations, using optical fiber switches, each delivering an average power of ~0.2 mW. Diffused light was collected at 16 locations using single photon counting detectors based on multianode photomultipliers and processed by time-correlated single photon counting boards [47]. The fNIRS optodes were arranged symmetrically (5 sources and 8 detectors on each hemisphere) over the entire brain obtaining 32 measurements points (i.e. mid-point between source-detector pairs, see Fig 1A).

thumbnail
Fig 1. (a) Positions of the 32 fNIRS channels projected over a three-dimensional reconstruction of the brain cortex (L: Left hemisphere, R: Right hemisphere).

The main cortical regions of interest are shown in the green circles: SMC- sensorimotor cortex (S1/M1); PMC/SMA- premotor cortex and supplementary motor area; PFC- prefrontal cortex; S2- secondary somatosensory area (See Table 1 for specific Brodmann areas representing each cortical region of interest). (b) Scheme illustrating the concept of early and late photons in TD fNIRS measurements (N is the number of photons, t the photon time-of-flight, and ρ the source detector distance).

https://doi.org/10.1371/journal.pone.0131951.g001

fNIRS optodes were positioned over the subjects scalp by means of a flexible electroencephalography (EEG) cap (g.EEGcap, g.tec medical engineering GmbH, Austria), presenting holes in correspondence to the positions of the International 10–20 system for EEG electrode placement (~3 cm source-detector distance). Home-made spring-loaded holders assured that optodes fitted into the holes and were firmly coupled to the skin. Moreover, the large holes on the cap allowed for good hair removal from under the optodes.

Since sufficient separation existed between optical channels in opposite hemispheres, a pair of light sources, one in each hemisphere, can be simultaneously activated without the risk of signal contamination. The sequential activation of the five pairs of sources every 0.2 s determined an overall acquisition time of 1 s for a full topographic map from 32 optical channels. Before the measurements, the count rate was set to ~50,000 ph/s for each wavelength which allowed for an acceptable value for the signal-to-noise-ratio during the measurements.

Cortical (intra-cerebral) and superficial (extra-cerebral) layer NIR light absorption changes, at experiment time t for wavelength λ, were calculated by means of the time-resolved Beer-Lambert formula in turbid media [48], eventually corrected to enhance the contribution from cortical (intra) layer and to remove possible disturbances caused by superficial (extra) layers [49]. The following formulas were used and a detailed description of the analysis can be found in [31]: (1) (2) where N(d,w; t;λ) is the number of photons collected in a time window with delay d and width w, and N0(d,w;λ) is the number of photons, collected in a window with the same width, at the same delay and for the same wavelength, averaged over the resting period of the protocol. Parameters dE, wE (dL, wL) represent delay and width of the time window corresponding to early (late) arriving photons, whereas LE (LL) represent the pathlength of the early (late) photons (see Fig 1B). As in our previous studies [41, 50, 51], we have used dE = 0 ps, wE = 800 ps, dL = 2,400 ps and wL = 1,200 ps. A rough assumption for the path-length is L = , where v is the speed of light in the medium and τ is the average photon time-of-flight. Therefore, in Eq.(1) and (2) we have used LE = E, and LL = L, with τE = 400 ps and τL = 3,000 ps calculated as average time-of-flight in the early and late time window, respectively. The number of received photons in the late time window were typically 6,000, with a coefficient of variation of 1.3% [45]. Making the assumption that hemoglobin is the only chromophore contributing to absorption, the time courses of O2Hb and HHb concentration changes were then derived by the Beer-Lambert law using the hemoglobin absorption spectra [27]. It should be noted that the thickness of the superficial layer was assumed to be constant across channels and subjects. However, choosing a time window of analysis with appropriate delay and width (e.g. dL = 2,400 ps, wL = 1200 ps) allows optimal representation of the photon migration through deeper (cortical) regions that accounts for inter-subject anatomical variability [52, 53].

Statistical parametric mapping of TD-fNIRS signals

General Linear Model (GLM) analysis was applied to the time courses of O2Hb and HHb concentration changes (in superficial and cortical layers) by using NIRS-SPM software [54, 55] both at the single subject level and for the whole group. The positions of the 32 channels were registered over a reference MRI atlas in the Montreal Neurological Institute (MNI) coordinates system [56], and the points on the scalp were projected over a three-dimensional reconstruction of the brain cortex (see Fig 1A). Cortical regions corresponding to each of the 32 channels were extrapolated using the Anatomy 1.8 toolbox for SPM [57] (see Table 1).

thumbnail
Table 1. Cortical regions (Brodmann area, BA) corresponding to the 32 fNIRS channels on the left (LH) and right (RH) hemisphere (see Fig 1 for the location of individual channels).

PFC: prefrontal cortex (BA 9, 10, 11, 46); SMC: primary sensorimotor cortex (S1: BA 1, 2, 3; M1: BA 4); PMC/SMA: premotor cortex and supplementary motor area (BA 6); S2: secondary somatosensory area (BA 43, 40); AG: angular gyrus (BA 39).

https://doi.org/10.1371/journal.pone.0131951.t001

Pre-processing algorithms were applied using internal routines of the NIRS-SPM software: a detrending filter (Wavelet-MDL) was used to remove physiological components from the NIRS signal due to breathing, cardiac, and vasomotion [54]. Time series (O2Hb or HHb time course) were modeled as a linear combination of L regressors (known functions) plus an error term ε for each of the J channels: (3) where X is the design matrix (TxL) containing the L regressors (columns) defined over T time points (rows), and Beta is the matrix (LxJ) representing the weight parameter of each regressor fitting the original data by minimization of the least-squares error For each experiment, a design matrix for O2Hb was built by two regressors (rest and task) obtained through a convolution between the canonical hemodynamic response function (HRF, modeled as a linear combination of gamma functions) [58] and a step-function equal to 1 during the task periods, 0 elsewhere; the rest regressor was modelled with a step-function being 1 during the rest and 0 during the task, convolved with the HRF. Conversely, given that HHb is expected to decrease during the task, the two regressors for HHb were built as a convolution between the HRF and a step-function equal to -1 during the task periods, 0 elsewhere; the opposite step function was used for the rest regressor.

GLM analyses were performed on the single experimental condition (VOL and different NMES conditions) in order to investigate if task-related hemodynamic responses (e.g. changes of O2Hb and HHb signals) were significantly different from changes measured during resting conditions. For each condition, a T-test was conducted for each channel k on the expected value of . Furthermore the SPM method was used to compare the activations registered during different conditions. Thus a new design matrix was built composed only by the five task regressors for all the conditions. T-tests were performed on a new linear combination of coefficients, contrasting the regressors related to the task performed in different experiments. The different NMES conditions were compared firstly with the task period during the VOL condition (chosen as reference), and then between each of the NMES conditions, by following the hypothesis that O2Hb and HHb variations should increase when increasing the current intensity.

The T statistic was calculated for every channel, and sparse T-values associated to each optical channel were interpolated over the whole probe extension with a cubic interpolation. The p-values were finally inferred from the calculated T-values and graphically represented on a map in a logarithmic scale. The same visualization method of the statistical analyses has been previously used by Fazli et al. [59].

Statistical analysis of current intensity, pain/discomfort and physiological measurements.

Statistical analyses were conducted with SigmaStat software (Systat Software Inc., Erkrath, Germany). All data were first examined for normality and homogeneity using Skewness-Kurtosis and Levene tests, respectively. A Friedman repeated measures ANOVA was used to determine whether current intensity amplitude level (mA) was different between the 4 NMES conditions (10%MTI, 50%MTI, MTI and MTI+). A one-way repeated measures ANOVA was used to determine whether each physiological measurement (HR, RR, and SC) and pain (PRS) rating was different across the experimental conditions (VOL, 10%MTI, 50%MTI, MTI and MTI+). A Friedman repeated measures ANOVA was used to determine whether discomfort (VAS) ratings were different between the experimental conditions (VOL, 10%MTI, 50%MTI, MTI and MTI+). When appropriate, the Tukey's HSD post-hoc test was used to find significant differences between pairs of conditions. Data are presented as mean±SD. The significance level was set at p<0.05.

Results

Current intensity, pain/discomfort ratings and physiological parameters

Table 2 shows the current intensities and pain/discomfort ratings for the experimental conditions. No subjective indications of pain or discomfort during the VOL condition were reported. For the 10%MTI condition, the current intensity was significantly lower compared to the other NMES conditions and below the motor threshold, such that wrist movements were not observed, and no or minimal pain/discomfort were reported by the subjects. As expected, at current intensities greater than the motor threshold (i.e. the 50%MTI, MTI and MTI+ NMES conditions) wrist extension movements were produced. The PRS (and VAS) scores were found to progressively increase from the 10%MTI (not painful/discomforting) to the 50%MTI (moderately painful/discomfort), and up to the MTI and MTI+ conditions (extremely painful, discomforting), which corresponded with the increase in current intensities (see Table 2). Although the VAS score and current intensity for MTI+ were significantly greater than the corresponding values at 50%MTI, they were both not significantly different from MTI.

thumbnail
Table 2. Grand average (n = 9, mean±SD) of the pain/discomfort and physiological parameters during the voluntary (VOL) and neuromuscular electrical stimulation (NMES) conditions at percentages of the individual maximal tolerated current intensity (MTI).

https://doi.org/10.1371/journal.pone.0131951.t002

The physiological parameters (SC, HR and RR) during the VOL and NMES conditions are shown in Table 2. Although SC during the NMES conditions with current intensity at 50%MTI and more (50%MTI, MTI and MTI+) was significantly greater than compared to VOL and 10%MTI, there were no significant differences between 50%MTI, MTI and MTI+. HR and RR remained relatively unchanged over the experimental conditions.

Cortical and superficial O2Hb and HHb maps during VOL and NMES conditions

Fig 2 shows the O2Hb and HHb maps for the cortical and superficial layers during the VOL and NMES conditions. For the VOL condition, the cortical layer O2Hb and HHb maps showed significant task related activations (i.e., O2Hb increase and HHb decrease) in the contralateral (left) hemisphere, primarily in the SMC and PMC/SMA regions; while the ipsilateral hemisphere showed no significant O2Hb increases or HHb decreases (see Fig 2A). The superficial layer O2Hb maps also showed a task related activation in the sensorimotor network area of the contralateral hemisphere (see Fig 2B), but these areas were narrower and characterized by a lower statistical significance than the corresponding area in the cortical layer. This was even more pronounced in the superficial layer HHb maps, where very narrow and not co-located areas with lower statistical significance were found compared to the cortical layer.

thumbnail
Fig 2. Group mean (n = 9) cortical (a) and superficial (b) oxygenated (O2Hb) and deoxygenated (HHb) hemoglobin maps during the voluntary (VOL) and NMES conditions at percentages of the individual maximal tolerated current intensity (MTI).

The vertical colored scale indicates the significance level based on the uncorrected p-values (p<0.05-p<0.0001).

https://doi.org/10.1371/journal.pone.0131951.g002

The GLM analysis indicated no significant superficial or cortical layer O2Hb increases or HHb decreases for the 10%MTI condition (data not shown in Fig 2). The O2Hb and HHb maps both showed that the NMES conditions that evoked wrist extension movements (50%MTI, MTI and MTI+) induced a significant area of activation in the cortical layer of the contralateral (left) sensorimotor network (primarily SMC, PMC/SMA, and S2) and PFC regions, and a smaller area of significant O2Hb increase in the ipsilateral hemisphere (see Fig 2A). Although, the NMES conditions also induced significant O2Hb increases in the superficial layer of the contralateral sensorimotor network (see Fig 2B), these areas were narrower (in the case of 50%MTI and MTI) or more spread out (in the case of MTI+) than in the cortical layer. Similarly, the superficial layer HHb maps showed a minimal area of HHb decrease with low significance in the contralateral sensorimotor network regions for the NMES conditions.

Comparison of cortical O2Hb and HHb maps between NMES and VOL conditions

The cortical layer O2Hb and HHb contrast maps (i.e. the representations of the statistical significance resulting from the comparison of different NMES conditions) between the three NMES conditions (50%MTI, MTI and MTI+) are shown in Fig 3. When contrasting MTI/MTI+ with 50%MTI, although the O2Hb maps showed no (MTI+) or minimal (MTI) greater activation of the contralateral sensorimotor network compared to 50%MTI, the HHb maps revealed greater activation in the contralateral sensorimotor network corresponding to the SMC and PMC/SMA regions. Furthermore, the O2Hb maps revealed a significantly greater area of activation in the ipsilateral sensorimotor network for both the MTI and MTI+ conditions compared to 50%MTI. For the contrast between MTI and MTI+, only the O2Hb maps showed a significantly greater area of activation in the ipsilateral sensorimotor network and PFC regions for the MTI+ than MTI condition, and no significant difference in the area of activation in the contralateral hemisphere.

thumbnail
Fig 3. Group mean (n = 9) cortical oxygenated (O2Hb) and deoxygenated (HHb) hemoglobin contrast maps between the NMES conditions at percentages of the individual maximal tolerated current intensity (MTI).

The vertical colored scale indicates the significance level based on the uncorrected p-values (p<0.05-p<0.0001).

https://doi.org/10.1371/journal.pone.0131951.g003

The cortical layer O2Hb and HHb contrast maps between the VOL and NMES conditions are shown in Fig 4. The O2Hb and HHb contrast maps indicated that 50%MTI, MTI and MTI+ activated a significantly greater area of the contralateral sensorimotor network regions compared with the VOL condition (see Fig 4). Furthermore, the O2Hb maps indicated that compared with the VOL condition, the 50%MTI, MTI and MTI+ conditions significantly activated also the ipsilateral sensorimotor network regions, including the PFC regions for the MTI+ condition.

thumbnail
Fig 4. Group mean (n = 9) cortical oxygenated (O2Hb) and deoxygenated (HHb) hemoglobin contrast maps between the voluntary (VOL) and NMES conditions at percentages of the individual maximal tolerated current intensity (MTI).

The vertical colored scale indicates the significance level based on the uncorrected p-values (p<0.05-p<0.0001).

https://doi.org/10.1371/journal.pone.0131951.g004

Discussion

To the best of our knowledge, the present TD fNIRS study was the first to determine the cortical sensorimotor network activation profile and pain responses induced by NMES-evoked wrist extension movements at increasing current intensities relative to the individual maximal tolerated current intensity (MTI), and with reference to those obtained in response to voluntary (VOL) wrist extension movements. The cortical O2Hb and HHb maps indicated that NMES-evoked wrist extension movements (50%MTI, MTI, and MTI+) activated a contralateral sensorimotor network (SMC, PMC/SMA, S2) as that during VOL wrist extension movements, which is consistent with previous fMRI studies that compared between NMES-evoked and VOL movements [10, 13]. The novel findings were that i) when the NMES current intensity was increased from 50%MTI (mild pain/discomfort) to MTI/MTI+ (extremely painful/discomfortable), distinct patterns of bilateral sensorimotor network (including PFC) activation were evident (Figs 2A and 3); ii) compared to the VOL condition, NMES at high current intensities activated the contralateral sensorimotor network in a larger area and with higher significance, and activations were also seen in the ipsilateral sensorimotor and bilateral PFC regions (Figs 2A and 4).

Cortical activation during NMES and VOL conditions

NMES-evoked movements are induced by applying to the skin overlying the muscles of interest repeated electrical currents at sufficient intensity to activate the peripheral motoneuronal axons that elicit the excitation-contraction coupling process for muscle force production. Although less well characterised, the NMES electrical field also activate peripheral sensory neuronal axons that send proprioceptive (and nociceptive) signals from the stimulated muscle to the CNS leading to supraspinal neural adaptations [57]. We have shown that unilateral NMES-evoked wrist extension movements (50%MTI, MTI and MTI+) activated primarily a contralateral sensorimotor network comprising the SMC, PMC/SMA and S2 regions, which are in agreement with previous fMRI findings of NMES-evoked movements [10, 11, 13, 60]

As expected, the unilateral VOL wrist extension movements were characterized by activation of contralateral sensorimotor network regions (SMC, PMC/SMA, S2) that are well known to be activated to perform hand motor tasks [19, 20, 22]. Furthermore, the regions of the contralateral sensorimotor network activated during VOL wrist extension movements were also activated during NMES-evoked wrist extension movements (Fig 2A). These findings are also in agreement with fMRI studies that compared the cortical sensorimotor network activation profile between VOL and NMES-evoked movements at current intensities at the motor threshold [10, 13]. A new finding was that NMES at painful/discomfortable current intensities (50%MTI, MTI, and MTI+) activated a greater and wider area of the contralateral sensorimotor network compared to VOL. It also activated regions of the ipsilateral sensorimotor network including bilateral PFC, which was particularly evident for the extremely painful/discomfortable NMES conditions (MTI and MTI+ conditions, Fig 3). It is possible to hypothesize that the greater bilateral sensorimotor network and PFC activation profile for the NMES conditions at high current intensities are related to greater pain processing in these cortical regions. These observations support our previous study [25] that showed an increased PFC activation with increasing NMES current intensities at MTI levels, and add to the previous studies that have alluded to the activation of the attention and pain processing network with NMES current intensities above the motor threshold [11, 14].

The pain processing network comprising the S1/S2 (pain intensity) and PFC (attention and evaluation of an emotional response to pain) regions receives parallel inputs from multiple nociceptive pathways, suggesting that pain is processed in a distributed fashion [61]. If a stimulus is intense enough to activate peripheral nociceptors, multiple brain areas of the pain matrix respond to this input in a correlated manner with perceived pain intensity [62, 63]. This schema is supported by the increased pain/discomfort scales, SC responses (a physiological measure of subjective pain/arousal [64]) and bilateral sensorimotor network activation profile during high NMES current intensities (Table 2 and Fig 2A). The S1, S2 and PFC are amongst those brain regions that consistently respond to painful/uncomfortable stimulations [61, 63, 6570], and the PFC is densely interconnected with sensory areas [71], which integrates motor with sensory/proprioceptive information [72]. The M1 is also known to be implicated in pain processing [73]. The exercise physiology and clinical benefits of NMES have been proposed to be via a sensorimotor integration mechanism; increased proprioceptive signals from NMES-evoked movements activate the sensorimotor network, particularly the SMC, thereby increasing corticospinal excitability, and facilitating greater voluntary activation of the relevant neuronal network [6, 74]. Accordingly, it can be suggested that NMES-evoked wrist extension movements at high stimulation intensities activate brain regions related to sensorimotor integration, which provides evidence for the concept that NMES-evoked somatosensory and nociceptive inputs lead to changes in M1 excitability [6, 7], which in turn can cause functional improvements, including muscle activation and strength [6].

The detection of bilateral sensorimotor network activation elicited by unilateral NMES of the wrist extensors at high current intensities is consistent with the presence of transcallosal interactions between the two hemispheres [75]. These observations signify the interhemispheric integration of the somatosensory and nociceptive information primarily at the level of the SMC (but also PMC/SMA, S2 and PFC), and provide support for the concept of distributed (bilateral) processing in cortical sensorimotor networks during unilateral NMES of the wrist extensors as previously shown for NMES of the wrist [11] and finger [13] extensors/flexors, unilateral median [24, 76] and peroneal [14] nerve stimulation, and unilateral voluntary movements with increased task demands [7779]. A scheme for the bilateral activation of cortical networks with unilateral VOL and NMES is proposed: During VOL wrist extension movements, the ipsilateral sensorimotor network processing of sensory and motor information was inhibited prior to and during movement execution through interhemispheric inhibitory processes [80] to allow the contralateral sensorimotor network to perform the right arm wrist extension movements, which was evidenced with robust activation of the contralateral sensorimotor network with no obvious ipsilateral activation. The ipsilateral inhibition may be aimed at reducing interference, aiding focus and thereby achieving high manual dexterity in unilateral voluntary motor tasks [81]. During NMES-evoked wrist extension movements, since it has been previously shown that unilateral median nerve stimulation at low current intensities activates contralateral sensorimotor regions while actively inhibiting ipsilateral sensorimotor regions [82], NMES at high stimulation currents may disinhibit ipsilateral sensorimotor network to process somatosensory and nociceptive information, which was evidenced in the present and a previous [14] study with bilateral activation of the sensorimotor network. Therefore we hypothesize that NMES at high current intensities is able to modify the excitability of interhemispheric connections and perhaps the balance between interhemispheric excitation and inhibition [6]. Further neurophysiological research is necessary to confirm this hypothesis.

Cortical hemodynamic measurements using TD fNIRS

The O2Hb and HHb maps of the deeper layers (including mainly the cerebral cortex) revealed a greater area of activation (i.e. significant increase of O2Hb and a corresponding significant decrease of HHb) during the experimental conditions (Fig 2A) which was not found in the superficial layers (including mainly the scalp skin) where mainly O2Hb maps showed significant statistical effects (Fig 2B). This behavior of O2Hb is typical of a greater skin blood efflux probably due to increasing pain/discomfort sensation by the subjects as the NMES current intensity was increased. This emotional state of the subjects was confirmed by the subjective pain/discomfort scales and SC responses, which increased over the NMES current levels, while HR and RR remained at a similar level of that one observed in the VOL condition (Table 2). These findings suggest that attentional/pain/arousal processing in the sensorimotor network and PFC rather than systemic physiological mechanisms were primarily responsible for the changes in cortical O2Hb and HHb maps.

The superficial layer O2Hb maps also increased with increase of the stimulation current intensity in good spatial agreement with the cortical O2Hb maps. These findings confirm the results by Kirilina et al. [28] who demonstrated that the spatial agreement between cortical and superficial O2Hb changes are located in the superficial blood vessels supplying/draining the cortical layers. This suggests that the source of task-evoked superficial systemic signals in fNIRS (primarily O2Hb signals) is co-localized with veins draining the scalp. Conversely, the HHb maps of the superficial layer were found unchanged or minimal changes during all experimental conditions. The poor changes in HHb maps of the superficial layer compared with the well localized and significant HHb maps at the cortical layer confirm that the greater activation found in the cortical HHb maps arose primarily from the underlying cortex. Therefore, we are confident that the cortical activation maps found during NMES and VOL conditions arose from the cerebral cortex and were not an artefact of the extracerebral systemic haemodynamic changes.

Conclusion

We have demonstrated that, compared to the VOL wrist extension movements, NMES induces a greater bilateral sensorimotor network activation profile with increasing stimulation intensities and pain ratings. This could partly be explained by the increased attentional/pain processing associated with higher stimulation current intensities and to the increased bilateral sensorimotor integration in these cortical regions. Future studies should apply TD fNIRS neuroimaging in longitudinal patient studies for assessing cortical neuroplastic changes associated with NMES neurorehabilitation treatment.

Author Contributions

Conceived and designed the experiments: MM RR AT MF SP VQ GK. Performed the experiments: MM RR LZ DC MC. Analyzed the data: MM AT RR LZ SP MF LS. Contributed reagents/materials/analysis tools: AT SP MM RR LZ. Wrote the paper: MM AT RR SP MF GK VQ LZ DC LS MC.

References

  1. 1. Sheffler LR, Chae J (2007) Neuromuscular electrical stimulation in neurorehabilitation. Muscle Nerve 35: 562–590. pmid:17299744
  2. 2. Gondin J, Brocca L, Bellinzona E, D'Antona G, Maffiuletti NA, Miotti D, et al. (2011) Neuromuscular electrical stimulation training induces atypical adaptations of the human skeletal muscle phenotype: a functional and proteomic analysis. J Appl Physiol 110: 433–450. pmid:21127206
  3. 3. Muthalib M, Jubeau M, Millet GY, Maffiuletti NA, Nosaka K (2009) Comparison between electrically evoked and voluntary isometric contractions for biceps brachii muscle oxidative metabolism using near-infrared spectroscopy. Eur J Appl Physiol 107: 235–241. pmid:19568766
  4. 4. Bax L, Staes F and Verhagen A (2005) Does neuromuscular electrical stimulation strengthen the quadriceps femoris? A systematic review of randomised controlled trials. Sports Med 35: 191–212. pmid:15730336
  5. 5. Vanderthommen M and Duchateau J (2007) Electrical stimulation as a modality to improve performance of the neuromuscular system. Exerc Sport Sci Rev 35: 180–185. pmid:17921786
  6. 6. Hortobagyi T and Maffiuletti NA (2011) Neural adaptations to electrical stimulation strength training. Eur J Appl Physiol 111: 2439–2449. pmid:21643920
  7. 7. Chipchase LS, Schabrun SM and Hodges PW (2011) Peripheral electrical stimulation to induce cortical plasticity: a systematic review of stimulus parameters. Clin Neurophysiol 122: 456–463. pmid:20739217
  8. 8. Papaiordanidou M, Guiraud D and Varray A (2010) Kinetics of neuromuscular changes during low-frequency electrical stimulation. Muscle Nerve 41: 54–62. pmid:19882645
  9. 9. Chipchase LS, Schabrun SM and Hodges PW (2011) Corticospinal excitability is dependent on the parameters of peripheral electric stimulation: a preliminary study. Arch Phys Med Rehabil 92: 1423–1430. pmid:21620374
  10. 10. Joa KL, Han YH, Mun CW, Son BK, Lee CH, Shin YB, et al. (2012) Evaluation of the brain activation induced by functional electrical stimulation and voluntary contraction using functional magnetic resonance imaging. J Neuroeng Rehabil 9: 48. pmid:22828165
  11. 11. Blickenstorfer A, Kleiser R, Keller T, Keisker B, Meyer M, Riener R, et al. (2009) Cortical and subcortical correlates of functional electrical stimulation of wrist extensor and flexor muscles revealed by fMRI. Hum Brain Mapp 30: 963–975. pmid:18344193
  12. 12. Smith GV, Alon G, Roys SR and Gullapalli RP (2003) Functional MRI determination of a dose-response relationship to lower extremity neuromuscular electrical stimulation in healthy subjects. Exp Brain Res 150: 33–39. pmid:12698214
  13. 13. Iftime-Nielsen SD, Christensen MS, Vingborg RJ, Sinkjaer T, Roepstorff A and Grey MJ (2012) Interaction of electrical stimulation and voluntary hand movement in SII and the cerebellum during simulated therapeutic functional electrical stimulation in healthy adults. Hum Brain Mapp 33: 40–49. pmid:21591025
  14. 14. Francis S, Lin X, Aboushoushah S, White TP, Phillips M, Bowtell R, et al. (2009) fMRI analysis of active, passive and electrically stimulated ankle dorsiflexion. Neuroimage 44: 469–479. pmid:18950717
  15. 15. Ferrari M and Quaresima V (2012) A brief review on the history of human functional near-infrared spectroscopy (fNIRS) development and fields of application. Neuroimage 63: 921–935. pmid:22510258
  16. 16. Hoshi Y (2011) Towards the next generation of near-infrared spectroscopy. Philos Trans A Math Phys Eng Sci 369: 4425–4439. pmid:22006899
  17. 17. Girouard H and Iadecola C (2006) Neurovascular coupling in the normal brain and in hypertension, stroke, and Alzheimer disease. J Appl Physiol 100: 328–335. pmid:16357086
  18. 18. Devor A, Sakadzic S, Srinivasan VJ, Yaseen MA, Nizar K, Saisan PA, et al. (2012) Frontiers in optical imaging of cerebral blood flow and metabolism. J Cereb Blood Flow Metab 32: 1259–1276. pmid:22252238
  19. 19. Leff DR, Orihuela-Espina F, Elwell CE, Athanasiou T, Delpy DT, Darzi AW, et al. (2011) Assessment of the cerebral cortex during motor task behaviours in adults: a systematic review of functional near infrared spectroscopy (fNIRS) studies. Neuroimage 54: 2922–2936. pmid:21029781
  20. 20. Ehrsson HH, Fagergren A, Jonsson T, Westling G, Johansson RS and Forssberg H (2000) Cortical activity in precision- versus power-grip tasks: an fMRI study. J Neurophysiol 83: 528–536. pmid:10634893
  21. 21. Muthalib M, Anwar AR, Perrey S, Dat M, Galka A, Wolff S, et al. (2013) Multimodal integration of fNIRS, fMRI and EEG neuroimaging. Clin Neurophysiol 124: 2060–2062. pmid:23648071
  22. 22. Anwar AR, Muthalib M, Perrey S, Galka A, Granert O, Wolff S, et al. (2012) Directionality Analysis on Functional Magnetic Resonance Imaging during Motor Task using Granger Causality. Conf Proc IEEE Eng Med Biol Soc 2012: 2287–2290. pmid:23366380
  23. 23. Tanosaki M, Sato C, Shimada M, Iguchi Y and Hoshi Y (2003) Effect of stimulus frequency on human cerebral hemodynamic responses to electric median nerve stimulation: a near-infrared spectroscopic study. Neurosci Lett 352: 1–4. pmid:14615035
  24. 24. Takeuchi M, Hori E, Takamoto K, Tran AH, Satoru K, Ishikawa A, et al. (2009) Brain cortical mapping by simultaneous recording of functional near infrared spectroscopy and electroencephalograms from the whole brain during right median nerve stimulation. Brain Topogr 22: 197–214. pmid:19705276
  25. 25. Muthalib M, Ferrari M, Quaresima V and Nosaka K (2012) Frontal cortex activation during electrical muscle stimulation as revealed by functional near-infrared spectroscopy (fNIRS). Adv Exp Med Biol 737: 45–49. pmid:22259080
  26. 26. Ellis H (2012) Anatomy of head injury. Surgery 30: 99–101.
  27. 27. Scholkmann F, Kleiser S, Metz AJ, Zimmermann R, Mata Pavia J, Wolf U, et al. (2014) A review on continuous wave functional near-infrared spectroscopy and imaging instrumentation and methodology. Neuroimage 85 6–27 pmid:23684868
  28. 28. Kirilina E, Jelzow A, Heine A, Niessing M, Wabnitz H, Bruhl R, et al. (2012) The physiological origin of task-evoked systemic artefacts in functional near infrared spectroscopy. Neuroimage 61: 70–81. pmid:22426347
  29. 29. Takahashi T, Takikawa Y, Kawagoe R, Shibuya S, Iwano T and Kitazawa S (2011) Influence of skin blood flow on near-infrared spectroscopy signals measured on the forehead during a verbal fluency task. Neuroimage 57: 991–1002. pmid:21600294
  30. 30. Hallacoglu B, Sassaroli A and Fantini S (2013) Optical characterization of two-layered turbid media for non-invasive, absolute oximetry in cerebral and extracerebral tissue. PLoS One 8: e64095. pmid:23724023
  31. 31. Torricelli A, Contini D, Pifferi A, Caffini M, Re R, Zucchelli L, et al. (2014) Time domain functional NIRS imaging for human brain mapping. Neuroimage 85: 28–50. pmid:23747285
  32. 32. Steinbrink J, Wabnitz H, Obrig H, Villringer A and Rinneberg H (2001) Determining changes in NIR absorption using a layered model of the human head. Phys Med Biol 46: 879–896. pmid:11277232
  33. 33. Del Bianco S, Martelli F and Zaccanti G (2002) Penetration depth of light re-emitted by a diffusive medium: theoretical and experimental investigation. Phys Med Biol 47: 4131–4144. pmid:12502038
  34. 34. Selb J, Stott JJ, Franceschini MA, Sorensen AG and Boas DA (2005) Improved sensitivity to cerebral hemodynamics during brain activation with a time-gated optical system: analytical model and experimental validation. J Biomed Opt 10: 11013. pmid:15847579
  35. 35. Selb J, Dale AM and Boas DA (2007) Linear 3D reconstruction of time-domain diffuse optical imaging differential data: improved depth localization and lateral resolution. Opt Express 15: 16400–16412. pmid:19550930
  36. 36. Sato C, Shimada M, Yamada Y and Hoshi Y (2005) Extraction of depth-dependent signals from time-resolved reflectance in layered turbid media. J Biomed Opt 10: 064008. pmid:16409073
  37. 37. Wabnitz H, Jelzow A, Mazurenka M, Steinkellner O, Macdonald R, Milej D, et al. (2014) Performance assessment of time-domain optical brain imagers, part 2: nEUROPt protocol. J Biomed Opt 19: 086012. pmid:25121480
  38. 38. Aletti F, Re R, Pace V, Contini D, Molteni E, Cerutti S, et al. (2012) Deep and surface hemodynamic signal from functional time resolved transcranial near infrared spectroscopy compared to skin flowmotion. Comput Biol Med 42: 282–289. pmid:21742320
  39. 39. Liebert A, Wabnitz H, Steinbrink J, Obrig H, Moller M, Macdonald R, et al. (2004) Time-resolved multidistance near-infrared spectroscopy of the adult head: intracerebral and extracerebral absorption changes from moments of distribution of times of flight of photons. Appl Opt 43: 3037–3047. pmid:15176190
  40. 40. Zhao Q, Spinelli L, Bassi A, Valentini G, Contini D, Torricelli A, et al. (2011) Functional tomography using a time-gated ICCD camera. Biomedical optics express 2: 705–716. pmid:21412474
  41. 41. Molteni E, Contini D, Caffini M, Baselli G, Spinelli L, Cubeddu R, et al. (2012) Load-dependent brain activation assessed by time-domain functional near-infrared spectroscopy during a working memory task with graded levels of difficulty. J Biomed Opt 17: 056005. pmid:22612128
  42. 42. Oldfield RC (1971) The assessment and analysis of handedness: the Edinburgh inventory. Neuropsychologia 9: 97–113. pmid:5146491
  43. 43. Paddon-Jones D, Muthalib M and Jenkins D (2000) The effects of a repeated bout of eccentric exercise on indices of muscle damage and delayed onset muscle soreness. J Sci Med Sport 3: 35–43. pmid:10839227
  44. 44. Peake JM, Nosaka K, Muthalib M and Suzuki K (2006) Systemic inflammatory responses to maximal versus submaximal lengthening contractions of the elbow flexors. Exerc Immunol Rev 12: 72–85. pmid:17201073
  45. 45. Contini D, Torricelli A, Pifferi A, Spinelli L, Paglia F and Cubeddu R (2006) Multi-channel time-resolved system for functional near infrared spectroscopy. Opt Express 14: 5418–5432. pmid:19516708
  46. 46. Contini D, Spinelli L, Caffini M, Cubeddu R and Torricelli A (2009) A multichannel time-domain brain oximeter for clinical studies. Proceedings of SPIE 7369: 73691D.
  47. 47. Becker W (2005) Advanced time-correlated single-photon counting,. Berlin: Springer.
  48. 48. Nomura Y, Hazeki O and Tamura M (1997) Relationship between time-resolved and non-time-resolved Beer-Lambert law in turbid media. Phys Med Biol 42: 1009–1022. pmid:9194125
  49. 49. Contini D, Torricelli A, Pifferi A, Spinelli L and Cubeddu R (2007) Novel method for depth-resolved brain functional imaging by time-domain NIRS. Proceedings of SPIE 6629: 662908.
  50. 50. Quaresima V, Ferrari M, Torricelli A, Spinelli L, Pifferi A and Cubeddu R (2006) Bilateral prefrontal cortex oxygenation responses to a verbal fluency task: a multichannel time-resolved near-infrared topography study. J Biomed Opt 10: 11012.
  51. 51. Butti M, Contini D, Molteni E, Caffini M, Spinelli L, Baselli G, et al. (2009) Effect of prolonged stimulation on cerebral hemodynamic: a time-resolved fNIRS study. Med Phys 36: 4103–4114. pmid:19810483
  52. 52. Haeussinger FB, Heinzel S, Hahn T, Schecklmann M, Ehlis AC and Fallgatter AJ (2012) Simulation of near-infrared light absorption considering individual head and prefrontal cortex anatomy: implications for optical neuroimaging. PLoS One 6: e26377.
  53. 53. Moreira-Gonzalez A, Papay FE and Zins JE (2006) Calvarial thickness and its relation to cranial bone harvest. Plast Reconstr Surg 117: 1964–1971. pmid:16651971
  54. 54. Jang KE, Tak S, Jung J, Jang J, Jeong Y and Ye JC (2009) Wavelet minimum description length detrending for near-infrared spectroscopy. J Biomed Opt 14: 034004. pmid:19566297
  55. 55. Ye JC, Tak S, Jang KE, Jung J and Jang J (2009) NIRS-SPM: statistical parametric mapping for near-infrared spectroscopy. Neuroimage 44: 428–447. pmid:18848897
  56. 56. Singh AK, Okamoto M, Dan H, Jurcak V and Dan I (2005) Spatial registration of multichannel multi-subject fNIRS data to MNI space without MRI. Neuroimage 27: 842–851. pmid:15979346
  57. 57. Eickhoff SB, Stephan KE, Mohlberg H, Grefkes C, Fink GR, Amunts K, et al. (2005) A new SPM toolbox for combining probabilistic cytoarchitectonic maps and functional imaging data. Neuroimage 25: 1325–1335. pmid:15850749
  58. 58. Friston KJ, Ashburner J, Kiebel S, Nichols T and Penny W (2006) Statistical Parametric Mapping: The Analysis of Functional Brain Images. Sandiego: Academic Press
  59. 59. Fazli S, Mehnert J, Steinbrink J, Curio G, Villringer A, Muller KR, et al. (2012) Enhanced performance by a hybrid NIRS-EEG brain computer interface. Neuroimage 59: 519–529. pmid:21840399
  60. 60. Han BS, Jang SH, Chang Y, Byun WM, Lim SK and Kang DS (2003) Functional magnetic resonance image finding of cortical activation by neuromuscular electrical stimulation on wrist extensor muscles. Am J Phys Med Rehabil 82: 17–20. pmid:12510180
  61. 61. Apkarian AV, Bushnell MC, Treede RD and Zubieta JK (2005) Human brain mechanisms of pain perception and regulation in health and disease. Eur J Pain 9: 463–484. pmid:15979027
  62. 62. Lorenz J and Casey KL (2005) Imaging of acute versus pathological pain in humans. Eur J Pain 9: 163–165. pmid:15737808
  63. 63. Lorenz J, Minoshima S and Casey KL (2003) Keeping pain out of mind: the role of the dorsolateral prefrontal cortex in pain modulation. Brain 126: 1079–1091. pmid:12690048
  64. 64. Storm H (2008) Changes in skin conductance as a tool to monitor nociceptive stimulation and pain. Curr Opin Anaesthesiol 21: 796–804. pmid:18997532
  65. 65. Niddam DM, Yeh TC, Wu YT, Lee PL, Ho LT, Arendt-Nielsen L, et al. (2002) Event-related functional MRI study on central representation of acute muscle pain induced by electrical stimulation. Neuroimage 17: 1437–1450. pmid:12414283
  66. 66. Azar AA (2009) Monitoring blood oxygenation changes due to acute pain stimuli using functional near-infrared spectroscopy (fNIRS). Conf Proc IEEE Eng Med Biol Soc 2009: 1375–1379. pmid:19964758
  67. 67. Becerra L, Harris W, Grant M, George E, Boas D and Borsook D (2009) Diffuse optical tomography activation in the somatosensory cortex: specific activation by painful vs. non-painful thermal stimuli. PLoS One 4: e8016. pmid:19956637
  68. 68. Tracey I and Mantyh PW (2007) The cerebral signature for pain perception and its modulation. Neuron 55: 377–391. pmid:17678852
  69. 69. Becerra L, Harris W, Joseph D, Huppert T, Boas DA and Borsook D (2008) Diffuse optical tomography of pain and tactile stimulation: activation in cortical sensory and emotional systems. Neuroimage 41: 252–259. pmid:18394924
  70. 70. Seifert F, Schuberth N, De Col R, Peltz E, Nickel FT and Maihofner C (2012) Brain activity during sympathetic response in anticipation and experience of pain. Hum Brain Mapp 34: 1768–1782. pmid:22438199
  71. 71. Rozzi S, Calzavara R, Belmalih A, Borra E, Gregoriou GG, Matelli M, et al. (2006) Cortical connections of the inferior parietal cortical convexity of the macaque monkey. Cereb Cortex 16: 1389–1417. pmid:16306322
  72. 72. Hinkley LB, Krubitzer LA, Nagarajan SS and Disbrow EA (2007) Sensorimotor integration in S2, PV, and parietal rostroventral areas of the human sylvian fissure. J Neurophysiol 97: 1288–1297. pmid:17122318
  73. 73. Frot M, Magnin M, Mauguiere F and Garcia-Larrea L (2012) Cortical representation of pain in primary sensory-motor areas (S1/M1)—a study using intracortical recordings in humans. Hum Brain Mapp 34: 2655–2668. pmid:22706963
  74. 74. Ridding MC and Rothwell JC (1999) Afferent input and cortical organisation: a study with magnetic stimulation. Exp Brain Res 126: 536–544. pmid:10422717
  75. 75. Iwamura Y (2000) Bilateral receptive field neurons and callosal connections in the somatosensory cortex. Philos Trans R Soc Lond B Biol Sci 355: 267–273. pmid:10724460
  76. 76. Sutherland MT and Tang AC (2006) Reliable detection of bilateral activation in human primary somatosensory cortex by unilateral median nerve stimulation. Neuroimage 33: 1042–1054. pmid:16997579
  77. 77. Stippich C, Blatow M, Durst A, Dreyhaupt J and Sartor K (2007) Global activation of primary motor cortex during voluntary movements in man. Neuroimage 34: 1227–1237. pmid:17137794
  78. 78. Chiou SY, Wang RY, Liao KK, Wu YT, Lu CF and Yang YR (2013) Co-activation of primary motor cortex ipsilateral to muscles contracting in a unilateral motor task. Clin Neurophysiol 124: 1353–1363. pmid:23478202
  79. 79. Derosiere G, Alexandre F, Bourdillon N, Mandrick K, Ward T, Perrey S (2014) Similar scaling of contralateral and ipsilateral cortical responses during graded unimanual force generation. Neuroimage 85: 471–477. pmid:23416251
  80. 80. Ferbert A, Priori A, Rothwell JC, Day BL, Colebatch JG, Marsden CD (1992) Interhemispheric inhibition of the human motor cortex. J Physiol 453: 525–546. pmid:1464843
  81. 81. Stefanovic B, Warnking JM, Pike GB (2004) Hemodynamic and metabolic responses to neuronal inhibition. Neuroimage 22: 771–778. pmid:15193606
  82. 82. Schafer K, Blankenburg F, Kupers R, Gruner JM, Law I, Lauritzen M, et al. (2012) Negative BOLD signal changes in ipsilateral primary somatosensory cortex are associated with perfusion decreases and behavioral evidence for functional inhibition. Neuroimage 59: 3119–3127. pmid:22155327