Whole-head MEG analysis of cortical spatial organization from unilateral stimulation of median nerve in both hands: No complete hemispheric homology

The effect of stimulation type, head modeling, and combined EEG and MEG on the source reconstruction of the somatosensory P20/N20 component

Modeling and experimental parameters influence the Electro- (EEG) and Magnetoencephalography (MEG) source analysis of the somatosensory P20/N20 component. In a sensitivity group study, we compare P20/N20 source analysis due to different stimulation type (Electric-Wrist [EW], Braille-Tactile [BT], or Pneumato-Tactile [PT]), measurement modality (combined EEG/MEG - EMEG, EEG, or MEG) and head model (standard or individually skull-conductivity calibrated including brain anisotropic conductivity). Considerable differences between pairs of stimulation types occurred (EW-BT: 8.7 ± 3.3 mm/27.1° ± 16.4°, BT-PT: 9 ± 5 mm/29.9° ± 17.3°, and EW-PT: 9.8 ± 7.4 mm/15.9° ± 16.5° and 75% strength reduction of BT or PT when compared to EW) regardless of the head model used. EMEG has nearly no localization differences to MEG, but large ones to EEG (16.1 ± 4.9 mm), while source orientation differences are non-negligible to both EEG (14° ± 3.7°) and MEG (12.5° ± 10.9°). Our calibration results show a considerable inter-subject variability (3.1-14 mS/m) for skull conductivity. The comparison due to different head model show localization differences smaller for EMEG (EW: 3.4 ± 2.4 mm, BT: 3.7 ± 3.4 mm, and PT: 5.9 ± 6.8 mm) than for EEG (EW: 8.6 ± 8.3 mm, BT: 11.8 ± 6.2 mm, and PT: 10.5 ± 5.3 mm), while source orientation differences for EMEG (EW: 15.4° ± 6.3°, BT: 25.7° ± 15.2° and PT: 14° ± 11.5°) and EEG (EW: 14.6° ± 9.5°, BT: 16.3° ± 11.1° and PT: 12.9° ± 8.9°) are in the same range. Our results show that stimulation type, modality and head modeling all have a non-negligible influence on the source reconstruction of the P20/N20 component. The complementary information of both modalities in EMEG can be exploited on the basis of detailed and individualized head models.

Localized N20 Component of Somatosensory Evoked Magnetic Fields in Frontoparietal Brain Tumor Patients Using Noise-Normalized Approaches

PURPOSE

To localize sensorimotor cortical activation in 10 patients with frontoparietal tumors using quantitative magnetoencephalography (MEG) with noise-normalized approaches.

MATERIAL AND METHODS

Somatosensory evoked magnetic fields (SEFs) were elicited in 10 patients with somatosensory tumors and in 10 control participants using electrical stimulation of the median nerve via the right and left wrists. We localized the N20m component of the SEFs using dynamic statistical parametric mapping (dSPM) and standardized low-resolution brain electromagnetic tomography (sLORETA) combined with 3D magnetic resonance imaging (MRI). The obtained coordinates were compared between groups. Finally, we statistically evaluated the N20m parameters across hemispheres using non-parametric statistical tests.

RESULTS

The N20m sources were accurately localized to Brodmann area 3b in all members of the control group and in seven of the patients; however, the sources were shifted in three patients relative to locations outside the primary somatosensory cortex (SI). Compared with the affected (tumor) hemispheres in the patient group, N20m amplitudes and the strengths of the current sources were significantly lower in the unaffected hemispheres and in both hemispheres of the control group. These results were consistent for both dSPM and sLORETA approaches.

CONCLUSION

Tumors in the sensorimotor cortex lead to cortical functional reorganization and an increase in N20m amplitude and current-source strengths. Noise-normalized approaches for MEG analysis that are integrated with MRI show accurate and reliable localization of sensorimotor function.

Invariance in current dipole moment density across brain structures and species: physiological constraint for neuroimaging

Although anatomical constraints have been shown to be effective for MEG and EEG inverse solutions, there are still no effective physiological constraints. Strength of the current generator is normally described by the moment of an equivalent current dipole Q. This value is quite variable since it depends on size of active tissue. In contrast, the current dipole moment density q, defined as Q per surface area of active cortex, is independent of size of active tissue. Here we studied whether the value of q has a maximum in physiological conditions across brain structures and species. We determined the value due to the primary neuronal current (q primary) alone, correcting for distortions due to measurement conditions and secondary current sources at boundaries separating regions of differing electrical conductivities. The values were in the same range for turtle cerebellum (0.56-1.48 nAm/mm(2)), guinea pig hippocampus (0.30-1.34 nAm/mm(2)), and swine neocortex (0.18-1.63 nAm/mm(2)), rat neocortex (~2.2 nAm/mm(2)), monkey neocortex (~0.40 nAm/mm(2)) and human neocortex (0.16-0.77 nAm/mm(2)). Thus, there appears to be a maximum value across the brain structures and species (1-2 nAm/mm(2)). The empirical values closely matched the theoretical values obtained with our independently validated neural network model (1.6-2.8 nAm/mm(2) for initial spike and 0.7-3.1 nAm/mm(2) for burst), indicating that the apparent invariance is not coincidental. Our model study shows that a single maximum value may exist across a wide range of brain structures and species, varying in neuron density, due to fundamental electrical properties of neurons. The maximum value of q primary may serve as an effective physiological constraint for MEG/EEG inverse solutions.

The effect of anodal transcranial direct current stimulation over the primary motor or somatosensory cortices on somatosensory evoked magnetic fields

OBJECTIVES

The purpose of this study was to investigate the effect of anodal transcranial direct-current stimulation (tDCS) applied over the primary motor (M1) or the primary somatosensory (S1) cortices on somatosensory evoked magnetic fields (SEFs) following median nerve stimulation.

METHODS

Anodal tDCS was applied for 15min on the left motor or somatosensory cortices at 1mA. SEFs were recorded following right median nerve stimulation using a magnetoencephalography (MEG) system before and after the application of tDCS. SEFs was measured and compared before and after tDCS was applied over M1 or S1.

RESULTS

The source strengths for the P35m and P60m increased after tDCS was applied over M1 and that for the P60m increased after tDCS was applied over S1. The mean equivalent current dipole (ECD) location for the P35m was located significantly anterior to that of the N20m, but only during post 1 (10-20min after tDCS was applied over M1).

CONCLUSION

Our results indicated that the anodal tDCS applied over M1 affected the P35m and P60m sources on SEF components, while that applied over S1 influenced the P60m source.

SIGNIFICANCE

We demonstrated anodal tDCS applied over M1 or S1 can modulate somatosensory processing and components of SEFs, confirming the hypothesis for locally distinct generators of the P35m and P60m sources.

Sensory handedness is not reflected in cortical responses after basic nerve stimulation: a MEG study

Motor dominance is well established, but sensory dominance is much less clear. We therefore studied the cortical evoked magnetic fields using magnetoencephalography (MEG) in a group of 20 healthy right handed subjects in order to examine whether standard electrical stimulation of the median and ulnar nerve demonstrated sensory lateralization. The global field power (GFP) curves, as an indication of cortical activation, did not depict sensory lateralization to the dominant left hemisphere. Comparison of the M20, M30, and M70 peak latencies and GFP values exhibited no statistical differences between the hemispheres, indicating no sensory hemispherical dominance at these latencies for each nerve. Field maps at these latencies presented a first and second polarity reversal for both median and ulnar stimulation. Spatial dipole position parameters did not reveal statistical left-right differences at the M20, M30 and M70 peaks for both nerves. Neither did the dipolar strengths at M20, M30 and M70 show a statistical left-right difference for both nerves. Finally, the Laterality Indices of the M20, M30 and M70 strengths did not indicate complete lateralization to one of the hemispheres. After electrical median and ulnar nerve stimulation no evidence was found for sensory hand dominance in brain responses of either hand, as measured by MEG. The results can provide a new assessment of patients with sensory dysfunctions or perceptual distortion when sensory dominance occurs way beyond the estimated norm.

Primary sensory and motor cortex activities during voluntary and passive ankle mobilization by the SHADE orthosis

This study investigates cortical involvement during ankle passive mobilization in healthy subjects, and is part of a pilot study on stroke patient rehabilitation. Magnetoencephalographic signals from the primary sensorimotor areas devoted to the lower limb were collected together with simultaneous electromyographic activities from tibialis anterior (TA). This was done bilaterally, on seven healthy subjects (aged 29 ± 7), during rest, left and right passive ankle dorsiflexion (imparted through the SHADE orthosis, O-PM, or neuromuscular electrical stimulation, NMES-PM), and during active isometric contraction (IC-AM). The effects of focussing attention on ankle passive movements were considered. Primary sensory (FS(S1)) and motor (FS(M1)) area activities were discriminated by the Functional Source Separation algorithm. Only contralateral FS(S1) was recruited by common peroneal nerve stimulation and only contralateral FS(M1) displayed coherence with TA muscular activity. FS(M1) showed higher power of gamma rhythms (33-90 Hz) than FS(S1). Both sources displayed higher beta (14-32 Hz) and gamma powers in the left than in the right hemisphere. Both sources displayed a bilateral reduction of beta power during IC-AM with respect to rest. Only FS(S1) beta band power reduced during O-PM. No beta band modulation was observed of either source during NMES-PM. Mutual FS(S1)-FS(M1) coherence in gamma2 band (61-90 Hz) showed a slight trend towards an increase when focussing attention during O-PM. Somatosensory and motor counterparts of lower limb cortical representations were discriminated in both hemispheres. SHADE was effective in generating repeatable dorsiflexion and inducing primary sensory involvement similarly to voluntary movement.

Dipole source analyses of early median nerve SEP components obtained from subdural grid recordings

The median nerve N20 and P22 SEP components constitute the initial response of the primary somatosensory cortex to somatosensory stimulation of the upper extremity. Knowledge of the underlying generators is important both for basic understanding of the initial sequence of cortical activation and to identify landmarks for eloquent areas to spare in resection planning of cortex in epilepsy surgery. We now set out to localize the N20 and P22 using subdural grid recording with special emphasis on the question of the origin of P22: Brodmann area 4 versus area 1. Electroencephalographic dipole source analysis of the N20 and P22 responses obtained from subdural grids over the primary somatosensory cortex after median nerve stimulation was performed in four patients undergoing epilepsy surgery. Based on anatomical landmarks, equivalent current dipoles of N20 and P22 were localized posterior to (n = 2) or on the central sulcus (n = 2). In three patients, the P22 dipole was located posterior to the N20 dipole, whereas in one patient, the P22 dipole was located on the same coordinate in anterior-posterior direction. On average, P22 sources were found to be 6.6 mm posterior [and 1 mm more superficial] compared with the N20 sources. These data strongly suggest a postcentral origin of the P22 SEP component in Brodmann area 1 and render a major precentral contribution to the earliest stages of processing from the primary motor cortex less likely.

Structural and functional asymmetry in the human parietal opercular cortex

In this combined electroencephalographic and magnetic resonance imaging (MRI) study, the asymmetry of functional and structural measures in the human parietal operculum (PO) were investigated. Median nerve somatosensory evoked potential recordings showed maximum scalp potentials over contralateral (N80, N110) and ipsilateral (N100, N130) temporal electrode positions. In accordance, MRI-coregistered source analysis revealed two electrical sources in the contralateral (N80, N110) and two in the ipsilateral (N100, N130) PO. The dipole orientations of the contra- and ipsilateral sources with earlier peak activation, N80 and N100, were more tangential than those of the later peaking N110 and N130 sources. The most prominent contralateral N110 source exhibited pronounced left lateralized dipole strengths in the 80- to 120-ms latency range, in contrast to symmetrical N80 and ipsilateral source responses. The asymmetry of the N110 source activity explained both the asymmetry of N110 and N100 scalp potentials. Morphometric analysis demonstrated no interhemispheric differences in the sizes of the anterior PO (aPO), containing the cytoarchitectonic areas OP3 and OP4, but left lateralized sizes of the posterior PO (pPO), which encompasses the anatomically defined areas OP1 and OP2. The N110 source was located in the pPO and its asymmetry was significantly correlated with the structural pPO asymmetry but not with handedness and auditory lateralization. Thus both structural and functional asymmetries exist in the human PO and they are closely related to each other but not to measures of brain asymmetry in other functional systems, i.e., auditory lateralization and handedness.

Time-dependent hemispheric shift of the cortical control of volitional swallowing

An important part of the cortical processing of swallowing takes place in the sensorimotor cortex, predominantly in the left hemisphere. However, until now, only deglutition related brain activation with low time resolution exceeding a time interval of 1 s has been reported. In this study, we have examined the chronological sequence of cortical swallowing processing in humans by means of high temporal resolution magnetoencephalography (MEG). The cortical MEG activity was recorded during self-paced volitional swallowing in 10 healthy subjects. Data were analyzed using synthetic aperture magnetometry and the group analysis was performed using a permutation test. Swallowing-related muscle activity was recorded by electromyography. Within the time interval of 1 s of the most pronounced muscular swallowing execution, the MEG analysis revealed neural activation in the primary sensorimotor cortex. During the first 600 ms, only left hemispheric activation was found, bihemispheric activation during the next 200 ms and a right hemispheric activation during the last 200 ms. Thus, our results demonstrate a time-dependent shift of neural activation from left to right sensorimotor cortex during deglutition with left hemispheric dominance in the early stage of volitional swallowing and right hemispheric dominance during its later part.

Quantifying interhemispheric symmetry of somatosensory evoked potentials with the intraclass correlation coefficient

Although large intersubject variability is reported for cortical somatosensory evoked potentials (SEPs), variability between hemispheres within one subject is thought to be small. Therefore, interhemispheric comparison of SEP waveforms might be clinically useful to detect unilateral abnormalities in cortical sensory processing. We developed and evaluated a new technique to quantify interhemispheric SEP symmetry that uses a time interval including multiple SEP components, measures similarity of SEP waveforms between both hemispheres and results in high symmetry values even in the presence of small interhemispheric anatomic differences. Median nerve SEPs were recorded in 50 healthy subjects (20-70 years) using 128-channel EEG. Symmetry was quantified by the intraclass correlation coefficient and correlation coefficient between global field power of left and right median nerve SEPs. In 74% of subjects left-right intraclass correlation coefficient was higher than 0.60, implying high SEP hemispheric symmetry in terms of shape and amplitude. Left-right intraclass correlation coefficients lower than 0.60 were due to differences in amplitude, unilateral absence of peaks, or shape differences. We quantified SEP waveform interhemispheric symmetry and found it to be high in most healthy subjects. This technique may therefore be useful for detection of unilateral abnormalities in cortical sensory processing.

MAGNETOENCEPHALOGRAPHIC STUDY OF POSTERIOR TIBIAL NERVE STIMULATION IN PATIENTS WITH INTRACRANIAL LESIONS AROUND THE CENTRAL SULCUS

OBJECTIVE

To study interhemispheric differences of somatosensory evoked field (SEF) characteristics and the spatial distribution of equivalent current dipole sources in patients with unilateral hemispheric lesions around the central sulcus region.

METHODS

In 17 patients with perirolandic lesions, averaged somatosensory responses after posterior tibial nerve stimulation at the ankle were recorded with magnetoencephalography. Dipole source solutions in the affected (AH) and unaffected (UH) hemispheres were analyzed and compared for latency, equivalent current dipole strength, root mean square, and spatial distribution in relation to clinical findings.

RESULTS

Three main SEF components, P45m, N60m, and P75m, were identified in the hemisphere contralateral to the stimulated nerve. Dipole strength for the P45m component was significantly higher in the AH compared with the UH. SEF characteristics in the AH and UH showed no significant differences with respect to component latency or dipole strength of the N60m and P75m components. Interdipole location asymmetries exceeded 1.0 cm in 71% of the patients. Comparison of the posterior tibial nerve evoked responses (P45m and N60m) in patients with motor deficits and patients without deficits showed that these responses are enlarged in the AH when perirolandic lesions are present. Patients with motor deficits also showed an increased response for P45m in the UH.

CONCLUSION

The results of posterior tibial nerve SEFs suggest spatial and functional changes in the somatosensory network as a result of perirolandic lesions with a possible relationship with clinical symptoms. The results can provide further basis for the evaluation of cortical changes in the presence of perirolandic lesions.

Cortical Characterization and Inter-Dipole Distance Between Unilateral Median Versus Ulnar Nerve Stimulation of Both Hands in MEG

Contralateral somatosensory evoked fields (SEF) by whole head MEG after unilateral median and ulnar nerve stimulation of both hands were studied in 10 healthy right-handed subjects. Major parameters describing cortical activity were examined to discriminate median and ulnar nerve evoked responses. Somatic sensitivity showed high similarity in the 4 study conditions for both hand and nerve. The brain SEFs consisted of 7-8 major peak stages with consistent responses in all subjects at M20, M30, M70 and M90. Comparable inter-hemispheric waveform profile but high inter-subject variability was found. Median nerve induced significantly shorter latencies in the early activities than those of the ulnar nerve. The 3D cortical maps in the post stimulus 450 ms timeframe showed for both nerves two polarity reversals, an early and a late one which is a new finding. Dipole characteristics showed differential sites for the M20 and M30 in the respective nerve. Higher dipole moments evoked by the median nerve were noticed when compared to the ulnar. Furthermore, the results of the dipole distances between both nerves for M20 were calculated to be at 11.17 mm +/- 4.93 (LH) and 16.73 mm +/- 5.66 (RH), respectively after right hand versus left hand stimulation. This study showed substantial differences in the cortical responses between median and ulnar nerve. Especially the dipole distance between median and ulnar nerve on the cortex was computed accurately for the first time in MEG. Little is known however of the cortical responses in chronic pain patients and the parameter(s) that may change in an individual patient or a group. These results provide precise basis for further evaluating cortical changes in functional disorders and disease sequelae related to median and ulnar nerves.