Young female participants show blunted placebo effects associated with blunted responses to a cue predicting a safe stimulus in the right dorsolateral prefrontal cortex.

Discrimination of cues predicting non-nociceptive/nociceptive stimuli is essential for predicting whether a non-painful or painful stimulus will be administered and for eliciting placebo/nocebo (pain reduction/pain enhancement) effects. Dysfunction of the neural system involved in placebo effects has been implicated in the pathology of chronic pain, while female sex is one of the important risk factors for development of chronic pain in young adults. The dorsolateral prefrontal cortex (dl-PFC) is suggested to be involved in placebo effects and is sensitive to sex and age. In this study, to examine the neural mechanisms by which sex and age alter placebo and nocebo effects, we analyzed cerebral hemodynamic activities in the dl-PFC in different sex and age groups during a differential conditioning task. During the training session, two different sounds were followed by low- and high-intensity electrical shocks. In the following recording session, electrical shocks, the intensity of which was mismatched to the sounds, were occasionally administered to elicit placebo and nocebo effects. In young female participants, both placebo effects and hemodynamic responses to the conditioned sounds in the right dl-PFC were significantly lower than those in elderly female participants, while there were no age differences in male participants. The hemodynamic responses to the sound paired with the safe stimulus in the right dl-PFC were significantly correlated with placebo effects, except in the young female group. These results suggest that blunted placebo effects in the young female participants are ascribed to blunted responses to the sound associated with the safe stimulus in the right dl-PFC, and that sex- and age-related factors may alter the responsiveness of the right dl-PFC to associative cues predicting a safe stimulus.

Time-course of physical properties of the psoas major muscle after exercise as assessed by MR elastography.

This study aimed to analyze the time-course of the physical properties of the psoas major muscle (PM) before and after exercise using magnetic resonance elastography (MRE). Muscle stiffness is one of the important properties associated with muscle function. However, there was no research on the stiffness of the PM after exercise. In this study, we investigated time-course changes of the shear modulus of the PM after exercise. Furthermore, T2 values and apparent diffusion coefficient (ADC), as the additional information associated with muscular physical properties, were also measured simultaneously. Healthy young male volunteers were recruited in this study (n = 9) and they were required to perform a hand-to-knee isometric and unilateral exercise (left side). At each time-point before and after exercise, a set of 3 types of MR scans to measure multiple physical properties of the PM [shear modulus (MRE), T2 values, and ADC] were repeatedly taken. On day 1, a single set MR scan was taken before exercise (pre-exercise MR scan), and 6 sets MR scans were taken (5.5 to 38.0 min after exercise). After about 10-min rest (46.0 to 56.0 min after exercise), 4 sets MR scans were taken (57.5 to 77.0 min after exercise). About 10-min rest was taken again (85.0-95.0 min after exercise), 4 sets MR scans were taken (96.5 to 116.0 min after exercise). On days 2 and 7, a single set MR scan (MRE, T2 value, and ADC) was taken on each experimental day. The data were analyzed as relative changes (%) of the given parameters to the pre-exercise values. The results indicated significant decreases in PM shear modulus up to about 30 min after exercise. Then, it gradually increased and showed significant increases at about 100 min after exercise compared to that before exercise. T2 values and ADC showed significant increases up to about 65 min after exercise compared to those before exercise, and then returned to the pre-exercise values. On days 2 and 7, all values showed no significant changes compared to the pre-exercise values. This study is the first to report the time-course of the physical properties of the PM after exercise.

Hamstring strains in professional rugby players result in increased fascial stiffness without muscle quality changes as assessed using shear wave elastography.

INTRODUCTION: Hamstring strain injury is common among sports injuries. A previous history of this injury is considered a strong predictor of recurrent hamstring strain injury. Fascial tissue reportedly becomes stiffer after hamstring strain injury. However, the association between fascial stiffness and previous hamstring strain injury has not been investigated in clinical studies. We aimed to determine whether a previous history of hamstring strain injury affects fascial tissue and muscle tissues using shear wave elastography. METHOD: In eleven male professional rugby players, the stiffness as a shear modulus (kPa) of fascial tissue and muscle was measured on the specific injured area measured by magnetic resonance imaging (MRI) at resting position by using shear wave elastography. The side-to-side differences between the injured and the uninjured side were analyzed. The length and area of the muscle scar tissue were evaluated by MRI in relation to fascial stiffness. RESULTS: The shear elastic modulus of fascia was stiffer in the injured vs. the uninjured side; however, no difference was observed in the muscle. No significant relationship was detected between the length and area of the muscle scar tissue (all P > 0.05). DISCUSSION: Rugby players with a previous history of hamstring strain injury exhibited passive stiffness of fascial tissues in the injured leg, regardless of the length or area of the muscle scar tissue. However, the passive stiffness of muscles was same between the injured and the uninjured leg. CONCLUSION: The results can be beneficial to consider future risk for hamstring strain injuries.

Previous hamstring muscle strain injury alters passive tissue stiffness and vibration sense.

BACKGROUND: Hamstring strain is one of the most common among sports injuries. A previous history of this injury is considered a strong predictor of recurrent hamstring strain injury. It has been suggested that fascial tissue alters after muscle strain injury. However, the association between previous hamstring strain injury and tissue stiffness and vibration sense detection has not been investigated. OBJECTIVES: We aimed to determine whether a previous history of hamstring strain injury affects tissue stiffness and vibration sense in professional soccer players. METHOD: The stiffness (MyotonPRO®) and vibration disappearance threshold (tuning fork) were measured in eight professional soccer players with previous history of hamstring strain and eight uninjured players. The differences between two groups' means were analyzed. Side-to-side differences between injured and uninjured legs were also analyzed. RESULTS: The tissue stiffness was higher, and the vibration disappearance threshold was lower, in previously injured players when compared to uninjured players. Similar differences were found between injured and uninjured legs. No significant relationship was detected between the age or body mass index (BMI) for both tissue stiffness and vibration disappearance threshold (all P < 0.05). CONCLUSIONS: Soccer players with a previous history of hamstring strain injury exhibited higher tissue stiffness and lower vibration sensitivity in the injured leg, regardless of the age and BMI. The results that players who have a previous hamstring strain injury with altered tissue stiffness and vibration sense will be useful and feasible evaluation for chronic muscle strain condition.

Magnetic resonance elastography of the supraspinatus muscle: A preliminary study on test-retest repeatability and wave quality with different frequencies and image filtering.

The purpose of this study was to determine an optimal condition (vibration frequency and image filtering) for stiffness estimation with high accuracy and stiffness measurement with high repeatability in magnetic resonance elastography (MRE) of the supraspinatus muscle. Nine healthy volunteers underwent two MRE exams separated by at least a 30 min break, on the same day. MRE acquisitions were performed with a gradient-echo type multi-echo MR sequence at 75, 100, and 125 Hz pneumatic vibration. Wave images were processed by a bandpass filter or filter combining bandpass and directional filters (bandpass-directional filter). An observer specified the region of interest (ROI) on clear wave propagation in the supraspinatus muscle, within which the observer measured the stiffness. This study assessed wave image quality according to two indices, as a substitute for the assessment of the accuracy of the stiffness estimation. One is the size of the clear wave propagation area (ROI size used to measure the stiffness) and the other is the qualitative stiffness resolution score in that area. These measurements made by the observer were repeated twice at least one month apart after each MRE exam. This study assessed the intra-examiner and observer repeatability of the stiffness value, ROI size and resolution score in each combination of vibration frequency and image filter. Repeatability of the data was analyzed using the intraclass correlation coefficient (ICC) and 95% limits-of-agreement (LOA) in Bland-Altman analysis. The analyses on intra-examiner and observer repeatability of stiffness indicated that the ICC and 95% LOA were not varied greatly depending on vibration frequency and image filter (intra-examiner repeatability, ICC range, 0.79 to 0.88; 95% LOA range, ±23.95 to ±32.42%, intra-observer repeatability, ICC range, 0.98 to 1.00; 95% LOA range, ±5.10 to ±10.99%). In the analyses on intra-examiner repeatability of ROI size, ICCs were rather low (ranging from: 0.03 to 0.69) while 95% LOA was large in all the combinations of vibration frequency and image filter (ranging from: ±62.66 to ±83.33%). In the analyses on intra-observer repeatability of ROI size, ICCs were sufficiently high in the total combination of vibration frequency and image filter (ranging from 0.80 to 0.87) while the 95% LOAs were better (lower) in the bandpass-directional filter than the bandpass filter (bandpass directional filter vs. bandpass filter, ±28.81 vs. ±54.83% at 75 Hz; ±25.63 vs. ±37.83% at 100 Hz; ±34.51 vs. ±43.36% at 125 Hz). In the analyses on intra-examiner and observer repeatability of resolution score, the mean difference (bias) between the two exams (or observations) was significantly low and there was almost no difference across all the combinations of vibration frequency and image filter (range of bias: -0.11-0.11 and -0.17-0.00, respectively). Additionally, effects of vibration frequency and image filter on wave image quality (ROI size and resolution score) were assessed separately in each exam. Both mean ROI size and resolution score in the bandpass-directional filter were larger than those in the bandpass filter. Among the data in the bandpass-directional filter, mean ROI size was larger at 75 and 100 Hz, and mean resolution score was larger at 100 and 125 Hz. Taking into consideration with the results of repeatability and wave image quality, the present results suggest that optimal vibration frequency and image filter for MRE of the supraspinatus muscles is 100 Hz and bandpass-directional filter, respectively.

Motor Imagery Training With Neurofeedback From the Frontal Pole Facilitated Sensorimotor Cortical Activity and Improved Hand Dexterity.

To develop a real-time neurofeedback system from the anterior prefrontal cortex (aPFC) using functional near-infrared spectroscopy (fNIRS) for motor rehabilitation, we investigated the effects of motor imagery training with neurofeedback from the aPFC on hand dexterity and cerebral hemodynamic activity during a motor rehabilitation task. Thirty-one right-handed healthy subjects participated in this study. They received motor imagery training six times for 2 weeks under fNIRS neurofeedback from the aPFC, in which they were instructed to increase aPFC activity. The real group subjects (n = 16) were shown real fNIRS neurofeedback signals from the aPFC, whereas the sham group subjects (n = 15) were shown irrelevant randomized signals during neurofeedback training. Before and after the training, hand dexterity was assessed by a motor rehabilitation task, during which cerebral hemodynamic activity was also measured. The results indicated that aPFC activity was increased during the training, and performance improvement rates in the rehabilitation task after the training was increased in the real group when compared with the sham group. Improvement rates of mean aPFC activity across the training were positively correlated with performance improvement rates in the motor rehabilitation task. During the motor rehabilitation task after the training, the hemodynamic activity in the left somatosensory motor-related areas [premotor area (PM), primary motor area (M1), and primary somatosensory area (S1)] was increased in the real group, whereas the hemodynamic activity was increased in the supplementary motor area in the sham group. This hemodynamic activity increases in the somatosensory motor-related areas after the training correlated with aPFC activity during the last 2 days of motor imagery training. Furthermore, improvement rates of M1 hemodynamic activity after the training was positively correlated with performance improvement rates in the motor rehabilitation task. The results suggest that the aPFC might shape activity in the somatosensory motor-related areas to improve hand dexterity. These findings further suggest that the motor imagery training using neurofeedback signals from the aPFC might be useful to patients with motor disability.

Impaired hemodynamic activity in the right dorsolateral prefrontal cortex is associated with impairment of placebo analgesia and clinical symptoms in postherpetic neuralgia.

The dorsolateral prefrontal cortex (dlPFC) is functionally linked to the descending pain modulation system and has been implicated in top down pain inhibition, including placebo analgesia. Therefore, functions of the dlPFC may be impaired in patients with chronic pain. Postherpetic neuralgia (PHN) is one of several syndromes with chronic neuropathic pain. In the present study, we investigated possible dysfunction of the dlPFC in chronic pain using patients with PHN. In a conditioning phase, heathy controls (n = 15) and patients with PHN (n = 7) were exposed to low (LF) and high (HF) frequency tones associated with noxious stimuli: weak (WS) and strong (SS) electrical stimulation, respectively. After the conditioning, cerebral hemodynamic activity was recorded from the bilateral dlPFC while the subjects were subjected to the cue tone-noxious electrical stimulation paradigm, in which incorrectly cued noxious stimuli were sometimes delivered to induce placebo and nocebo effects. The results indicated that hemodynamic responses to the LF tone in the right dlPFC was significantly lower in patients with PHN compared to the healthy controls. Furthermore, the same hemodynamic responses in the right dlPFC were correlated with placebo effects. In addition, clinical symptoms of PHN were negatively correlated to cerebral hemodynamic responses in the right dlPFC and magnitudes of the placebo effects. The results suggest that the right dlPFC, which is closely associated with the descending pain modulation system, is disturbed in PHN.

Distinct effects of thermal treatments after lengthening contraction on mechanical hyperalgesia and exercise-induced physiological changes in rat muscle.

Delayed-onset muscle soreness (DOMS) is a common but displeasing event induced by excessive muscle use or unaccustomed exercise and characterized by tenderness and movement-related pain in the exercised muscle. Thermal therapies, either icing or heating applied to muscles immediately after exercise, have been used as therapeutic interventions for DOMS. However, the mechanisms of their analgesic effects, and physiological and metabolic changes in the muscle during thermal therapy, remain unclear. In the present study, we investigated the effects of both thermal treatments on mechanical hyperalgesia of DOMS and physiological and muscle metabolite changes using the rat DOMS model induced by lengthening contraction (LC) to the gastrocnemius muscle. Heating treatment just after LC induced analgesic effects, while rats with icing treatment showed mechanical hyperalgesia similar to that of the LC group. Furthermore, increased physiological responses (e.g., muscle temperature and blood flow) following the LC were significantly kept high only in the rats with heating treatment. In addition, heating treatment increased metabolites involved in the improvement of blood flow and oxidative metabolisms in the exercised muscle. The results indicated that heating treatment just after LC has analgesic effects on DOMS, which might be mediated partly through the improvement of muscle oxidative metabolisms by changes in metabolites and elevated physiological responses.NEW & NOTEWORTHY Physiological effects of thermal therapy in the muscle and its mechanisms of analgesic effects remain unclear. The results indicated that heating, but not icing, treatment just after lengthening contractions induced analgesic effects in the rat muscle. Increases in hemodynamics, muscle temperature, and metabolites such as nicotinamide were more prominent in heating treatment, consistent with improvement of muscle oxidative metabolisms, which might reduce chemical factors to induce mechanical hyperalgesia.

Analgesic Effects of Compression at Trigger Points Are Associated With Reduction of Frontal Polar Cortical Activity as Well as Functional Connectivity Between the Frontal Polar Area and Insula in Patients With Chronic Low Back Pain: A Randomized Trial.

BACKGROUND: Compression of myofascial trigger points (MTrPs) in muscles is reported to reduce chronic musculoskeletal pain. Although the prefrontal cortex (PFC) is implicated in development of chronic pain, the mechanisms of how MTrP compression at low back regions affects PFC activity remain under debate. In this study, we investigated effects of MTrP compression on brain hemodynamics and EEG oscillation in subjects with chronic low back pain. METHODS: The study was a prospective, randomized, parallel-group trial and an observer and subject-blinded clinical trial. Thirty-two subjects with chronic low back pain were divided into two groups: subjects with compression at MTrPs (n = 16) or those with non-MTrPs (n = 16). Compression at MTrP or non-MTrP for 30 s was applied five times, and hemodynamic activity (near-infrared spectroscopy; NIRS) and EEGs were simultaneously recorded during the experiment. RESULTS: The results indicated that compression at MTrPs significantly (1) reduced subjective pain (P < 0.05) and increased the pressure pain threshold (P < 0.05), (2) decreased the NIRS hemodynamic activity in the frontal polar area (pPFC) (P < 0.05), and (3) increased the current source density (CSD) of EEG theta oscillation in the anterior part of the PFC (P < 0.05). CSD of EEG theta oscillation was negatively correlated with NIRS hemodynamic activity in the pPFC (P < 0.05). Furthermore, functional connectivity in theta bands between the medial pPFC and insula cortex was significantly decreased in the MTrP group (P < 0.05). The functional connectivity between those regions was positively correlated with subjective low back pain (P < 0.05). DISCUSSION: The results suggest that MTrP compression at the lumbar muscle modulates pPFC activity and functional connectivity between the pPFC and insula, which may relieve chronic musculoskeletal pain. TRIAL REGISTRATION: This trial was registered at University Hospital Medical Information Network Clinical Trials Registry (UMIN000033913) on 27 August 2018, at https://upload. umin.ac.jp/cgi-open-bin/ctr/ctr_view.cgi?recptno=R000038660.

Corrigendum: Analgesic Effects of Compression at Trigger Points Are Associated With Reduction of Frontal Polar Cortical Activity as Well as Functional Connectivity Between the Frontal Polar Area and Insula in Patients With Chronic Low Back Pain: A Randomized Trial.

[This corrects the article DOI: 10.3389/fnsys.2019.00068.].

Cerebral Hemodynamics in Speech-Related Cortical Areas: Articulation Learning Involves the Inferior Frontal Gyrus, Ventral Sensory-Motor Cortex, and Parietal-Temporal Sylvian Area.

Although motor training programs have been applied to childhood apraxia of speech (AOS), the neural mechanisms of articulation learning are not well understood. To this aim, we recorded cerebral hemodynamic activity in the left hemisphere of healthy subjects (n = 15) during articulation learning. We used near-infrared spectroscopy (NIRS) while articulated voices were recorded and analyzed using spectrograms. The study consisted of two experimental sessions (modified and control sessions) in which participants were asked to repeat the articulation of the syllables "i-chi-ni" with and without an occlusal splint. This splint was used to increase the vertical dimension of occlusion to mimic conditions of articulation disorder. There were more articulation errors in the modified session, but number of errors were decreased in the final half of the modified session; this suggests that articulation learning took place. The hemodynamic NIRS data revealed significant activation during articulation in the frontal, parietal, and temporal cortices. These areas are involved in phonological processing and articulation planning and execution, and included the following areas: (i) the ventral sensory-motor cortex (vSMC), including the Rolandic operculum, precentral gyrus, and postcentral gyrus, (ii) the dorsal sensory-motor cortex, including the precentral and postcentral gyri, (iii) the opercular part of the inferior frontal gyrus (IFGoperc), (iv) the temporal cortex, including the superior temporal gyrus, and (v) the inferior parietal lobe (IPL), including the supramarginal and angular gyri. The posterior Sylvian fissure at the parietal-temporal boundary (area Spt) was selectively activated in the modified session. Furthermore, hemodynamic activity in the IFGoperc and vSMC was increased in the final half of the modified session compared with its initial half, and negatively correlated with articulation errors during articulation learning in the modified session. The present results suggest an essential role of the frontal regions, including the IFGoperc and vSMC, in articulation learning, with sensory feedback through area Spt and the IPL. The present study provides clues to the underlying pathology and treatment of childhood apraxia of speech.

The Development of Vibration System for Applying Magnetic Resonance Elastography (MRE) to the Supraspinatus Muscle.

Palpation is a standard clinical tool to diagnose abnormal stiffness changes in soft tissues. However, it is difficult to palpate the supraspinatus muscle because it locates under the trapezius muscle. The magnetic resonance elastography (MRE) uses harmonic mechanical excitation to quantitatively measure the stiffness (shear modulus) of both the superficial and deep tissues. The purpose of this study was to build a vibration system for applying the MRE to the supraspinatus muscle. In this study, a power amplifier and a pneumatic pressure generator were used to supply vibrations to a vibration pad. Six healthy volunteers underwent MRE. We investigated the effects of position (the head of the humerus and the trapezius muscle) of the vibration pad on the patterns of wave propagation (wave image). When the vibration pad was placed in the trapezius muscle, the wave images represented clear wave propagation. On the other hand, when the vibration pad was placed in the head of the humerus, the wave images represented unclear wave propagation. This result might be caused by wave interferences resulting from the vibrations from bones and an intramuscular tendon of the supraspinatus muscle. The mean shear modulus also was 8.12 ± 1.83 (mean ± SD) kPa, when the vibration pad was placed in the trapezius muscle. Our results demonstrated that the vibration pad should be placed in the trapezius muscle in the MRE of the supraspinatus muscle.

Manual therapy ameliorates delayed-onset muscle soreness and alters muscle metabolites in rats.

Delayed-onset muscle soreness (DOMS) can be induced by lengthening contraction (LC); it can be characterized by tenderness and movement-related pain in the exercised muscle. Manual therapy (MT), including compression of exercised muscles, is widely used as physical rehabilitation to reduce pain and promote functional recovery. Although MT is beneficial for reducing musculoskeletal pain (i.e. DOMS), the physiological mechanisms of MT remain unclear. In the present study, we first developed an animal model of MT in DOMS; LC was applied to the rat gastrocnemius muscle under anesthesia, which induced mechanical hyperalgesia 2-4 days after LC. MT (manual compression) ameliorated mechanical hyperalgesia. Then, we used capillary electrophoresis time-of-flight mass spectroscopy (CE-TOFMS) to investigate early effects of MT on the metabolite profiles of the muscle experiencing DOMS. The rats were divided into the following three groups; (1) normal controls, (2) rats with LC application (LC group), and (3) rats undergoing MT after LC (LC + MT group). According to the CE-TOFMS analysis, a total of 171 metabolites were detected among the three groups, and 19 of these metabolites were significant among the groups. Furthermore, the concentrations of eight metabolites, including branched-chain amino acids, carnitine, and malic acid, were significantly different between the LC + MT and LC groups. The results suggest that MT significantly altered metabolite profiles in DOMS. According to our findings and previous data regarding metabolites in mitochondrial metabolism, the ameliorative effects of MT might be mediated partly through alterations in metabolites associated with mitochondrial respiration.

Selective Medial Prefrontal Cortex Responses During Live Mutual Gaze Interactions in Human Infants: An fNIRS Study.

To investigate the role of the prefrontal cortex (PFC) in processing multimodal communicative ostensive signals in infants, we measured cerebral hemodynamic responses by using near-infrared spectroscopy (NIRS) during the social interactive play "peek-a-boo", in which both visual (direct gaze) and auditory (infant-directed speech) stimuli were presented. The infants (mean age, around 7 months) sat on their mother's lap, equipped with an NIRS head cap, and looked at a partner's face during "peek-a-boo". An eye-tracking system simultaneously monitored the infants' visual fixation patterns. The results indicate that, when the partner presented a direct gaze, rather than an averted gaze, toward an infant during social play, the infant fixated on the partner's eye region for a longer duration. Furthermore, hemodynamic activity increased more prominently dorsomedial prefrontal cortex (mPFC) in response to social play with a partner's direct gaze compared to an averted gaze. In contrast, hemodynamic activity increased in the right dorsolateral prefrontal cortex (R-lPFC) regardless of a partner's eye gaze direction. These results indicate that a partner's direct gaze shifts an infant's attention to the partner's eyes for interactive communication, and specifically activates the mPFC. The differences in hemodynamic responses between the mPFC and R-lPFC suggest functional differentiation within the PFC, and a specific role of the mPFC in the perception of face-to-face communication, especially in mutual gaze, which is essential for social interaction.

Cortical hemodynamic responses to intravenous thiamine propyldisulphide administration detected by multichannel near infrared spectroscopy (NIRS) system.

Intravenous injection of thiamine propyldisulphide (TPD), which induces sensation of a garlic-like odor, has been used as a representative subjective olfactory test in Japan. However, cortical loci activated by TPD still remain unclear. We recorded cerebral hemodynamic responses (changes in Oxy-Hb concentrations) induced by TPD administration using whole-head multi-channel near infrared spectroscopy (NIRS) system based on 3D-MRIs. TPD as an odorant and saline as a control were injected from the cephalic vein in the left forearm in ten male normosmic (five young and five elderly) subjects and five dysosmic elderly patients. The all normosmic, but not dysosmic, subjects felt the garlic-like odor in the all TPD trials. There was no significant difference in hemodynamic responses between the young and elderly normosmic subjects. However, TPD injection induced significantly larger hemodynamic responses in the bilateral operculums, bilateral dorsolateral prefrontal cortices (PFC) and anteromedial PFC in the normosmic subjects, compared with saline injection. Onset latencies of these hemodynamic responses were significantly correlated with onset latencies of subjective odor sensation in the normosmic subjects. Comparison of hemodynamic responses between the normosmic and dysosmic subjects indicated a significant difference in the bilateral operculums. The results demonstrated that Oxy-Hb increases in the bilateral operculums reflected olfactory sensation induced by TPD injection. Consideration of a route for intravenous TPD to reach the olfactory mucosa suggests that these hemodynamic responses might be attributed to food-related retronasal olfactory responses to TPD.

Brain cortical mapping by simultaneous recording of functional near infrared spectroscopy and electroencephalograms from the whole brain during right median nerve stimulation.

To investigate relationships between hemodynamic responses and neural activities in the somatosensory cortices, hemodynamic responses by near infrared spectroscopy (NIRS) and electroencephalograms (EEGs) were recorded simultaneously while subjects received electrical stimulation in the right median nerve. The statistical significance of the hemodynamic responses was evaluated by a general linear model (GLM) with the boxcar design matrix convoluted with Gaussian function. The resulting NIRS and EEGs data were stereotaxically superimposed on the reconstructed brain of each subject. The NIRS data indicated that changes in oxy-hemoglobin concentration increased at the contralateral primary somatosensory (SI) area; responses then spread to the more posterior and ipsilateral somatosensory areas. The EEG data indicated that positive somatosensory evoked potentials peaking at 22 ms latency (P22) were recorded from the contralateral SI area. Comparison of these two sets of data indicated that the distance between the dipoles of P22 and NIRS channels with maximum hemodynamic responses was less than 10 mm, and that the two topographical maps of hemodynamic responses and current source density of P22 were significantly correlated. Furthermore, when onset of the boxcar function was delayed 5-15 s (onset delay), hemodynamic responses in the bilateral parietal association cortices posterior to the SI were more strongly correlated to electrical stimulation. This suggests that GLM analysis with onset delay could reveal the temporal ordering of neural activation in the hierarchical somatosensory pathway, consistent with the neurophysiological data. The present results suggest that simultaneous NIRS and EEG recording is useful for correlating hemodynamic responses to neural activity.