Serial activation of distinct cytoarchitectonic areas of the human S1 cortex after posterior tibial nerve stimulation

Suppression of somatosensory stimuli during motor planning may explain levels of balance and mobility after stroke

The ability to actively suppress, or gate, irrelevant sensory information is required for safe and efficient walking in sensory-rich environments. Both motor attention and motor planning alter somatosensory evoked potentials (SEPs) in healthy adults. This study's aim was to examine the effect of motor attention on processing of irrelevant somatosensory information during plantar flexion motor planning after stroke. Thirteen healthy older adults and 11 individuals with stroke participated. Irrelevant tibial nerve stimulation was delivered while SEPs were recorded over Cz, overlaying the leg portion of the sensorimotor cortex at the vertex of the head. Three conditions were tested in both legs: (1) Rest, (2) Attend To, and (3) Attend Away from the stimulated limb. In conditions 2 and 3, relevant vibration cued voluntary plantar flexion movements of the stimulated (Attend To) or non-stimulated (Attend Away) leg. SEP amplitudes were averaged during motor planning per condition. Individuals with stroke did not show attention-mediated gating of the N40 component associated with irrelevant somatosensory information during motor planning. It may be that dysfunction in pathways connecting to area 3b explains the lack of attention-mediated gating of the N40. Also, attention-mediated gating during motor planning explained significant and unique variance in a measure of community balance and mobility combined with response time. Thus, the ability to gate irrelevant somatosensory information appears important for stepping in both older adults and after stroke. Our data suggest that therapies that direct motor attention could positively impact walking after stroke.

Dipole orientation of receptive fields in the somatosensory cortex after stimulation of the posterior tibial nerve in humans

The origins of the earliest evoked potentials and magnetic fields after tibial nerve electrical stimulation are still controversial. We focused on the initial activity from the gyrus area and analyzed the components for the coronal and sagittal planes. In 12 healthy right-handed subjects, electrical stimuli were delivered to the left posterior tibial nerve at the ankle. The cortical somatosensory evoked fields were recorded, and the equivalent current dipoles were calculated and separated into the sagittal plane (y-components) and coronal plane (x-components) components. In nine subjects, we observed two deflections (y1 and y2) in the y-component waveform and two deflections (x1 and x2) in the x-component waveform over 60 milliseconds; y1 was directed anteriorly, y2 posteriorly, x1 to the left, and x2 to the right. The y1 was originated in the anterior wall of the central sulcus. By focusing on the y-component, we elucidated the existence of the posteroanterior component, which may originate from the area 3b (gyrus) in tibial nerve somatosensory evoked fields and has the same quality as N20m for median nerve somatosensory evoked fields. This is the first report to suggest that the posteroanterior component in the tibial nerve is analogous to N20m in the median nerve using magnetoencephalography.

Human Tonotopic Maps and their Rapid Task-Related Changes Studied by Magnetic Source Imaging

A brief review of previous studies is presented on tonotopic organization of primary auditory cortex (AI) in humans. Based on the place theory for pitch perception, in which place information from the cochlea is used to derive pitch, a well-organized layout of tonotopic map is likely in human AI. The conventional view of tonotopy in human AI is a layout inwhich the medial-to-lateral portion of Heschl's gyrus represents high-to-low frequency tones. However, we have shown that the equivalent current dipole (BCD) in auditory evoked magnetic fields in the rising phase of N100m response dynamically moves along the long axis of Heschl's gyrus. Based on analyses of the current sources for high-pitched and low-pitched tones in the right and left hemispheres, we propose an alternative tonotopic map in human AI. In the right AI, isofrequency bands for each tone frequency are parallell to the first transverse sulcus; on the other hand, the layout for tonotopy in the left AI seems poorly organized. The validity of single dipole modelling in the calculation of a moving source and the discrepancy as to tonotopic maps in the results between auditory evoked fields or intracerebral recordings and neuroimaging studies also are discussed. The difference in the layout of isofrequency bands between the right and left auditory cortices may reflect distinct functional roles in auditory information processing such as pitch versus phonetic analysis.

A new approach to neuroimaging with magnetoencephalography

We discuss the application of beamforming techniques to the field of magnetoencephalography (MEG). We argue that beamformers have given us an insight into the dynamics of oscillatory changes across the cortex not explored previously with traditional analysis techniques that rely on averaged evoked responses. We review several experiments that have used beamformers, with special emphasis on those in which the results have been compared to those observed in functional magnetic resonance imaging (fMRI) and on those studying induced phenomena. We suggest that the success of the beamformer technique, despite the assumption that there are no linear interactions between the mesoscopic local field potentials across distinct cortical areas, may tell us something of the balance between functional integration and segregation in the human brain. What is more, MEG beamformer analysis facilitates the study of these complex interactions within cortical networks that are involved in both sensory-motor and cognitive processes.

Localization bias and spatial resolution of adaptive and non-adaptive spatial filters for MEG source reconstruction

This paper discusses the location bias and the spatial resolution in the reconstruction of a single dipole source by various spatial filtering techniques used for neuromagnetic imaging. We first analyze the location bias for several representative adaptive and non-adaptive spatial filters using their resolution kernels. This analysis theoretically validates previously reported empirical findings that standardized low-resolution electromagnetic tomography (sLORETA) has no location bias. We also find that the minimum-variance spatial filter does exhibit bias in the reconstructed location of a single source, but that this bias is eliminated by using the normalized lead field. We then focus on the comparison of sLORETA and the lead-field normalized minimum-variance spatial filter, and analyze the effect of noise on source location bias. We find that the signal-to-noise ratio (SNR) in the measurements determines whether the sLORETA reconstruction has source location bias, while the lead-field normalized minimum-variance spatial filter has no location bias even in the presence of noise. Finally, we compare the spatial resolution for sLORETA and the minimum-variance filter, and show that the minimum-variance filter attains much higher resolution than sLORETA does. The results of these analyses are validated by numerical experiments as well as by reconstructions based on two sets of evoked magnetic responses.

Asymptotic SNR of scalar and vector minimum-variance beamformers for neuromagnetic source reconstruction

To reconstruct neuromagnetic sources, the minimum-variance beamformer has been extended to incorporate the three-dimensional vector nature of the sources, and two types of extensions-the scalar- and vector-type extensions-have been proposed. This paper discusses the asymptotic signal-to-noise ratio (SNR) of the outputs of these two types of beamformers. We first show that these two types of beamformers give exactly the same output power and output SNR if the beamformer pointing direction is optimized. We then compare the output SNR of the beamformer with optimum direction to that of the conventional vector beamformer formulation where the beamformer pointing direction is not optimized. The comparison shows that the beamformer with optimum direction gives an output SNR superior to that of the conventional vector beamformer. Numerical examples validating the results of the analysis are presented.

A general linear model for MEG beamformer imaging

A new general linear model (GLM) beamformer method is described for processing magnetoencephalography (MEG) data. A standard nonlinear beamformer is used to determine the time course of neuronal activation for each point in a predefined source space. A Hilbert transform gives the envelope of oscillatory activity at each location in any chosen frequency band (not necessary in the case of sustained (DC) fields), enabling the general linear model to be applied and a volumetric T statistic image to be determined. The new method is illustrated by a two-source simulation (sustained field and 20 Hz) and is shown to provide accurate localization. The method is also shown to locate accurately the increasing and decreasing gamma activities to the temporal and frontal lobes, respectively, in the case of a scintillating scotoma. The new method brings the advantages of the general linear model to the analysis of MEG data and should prove useful for the localization of changing patterns of activity across all frequency ranges including DC (sustained fields).

Performance of an MEG Adaptive-Beamformer Source Reconstruction Technique in the Presence of Additive Low-Rank Interference

The influence of external interference on neuromagnetic source reconstruction by adaptive beamformer techniques was investigated. In our analysis, we assume that the interference has the following two properties: First, it is additive and uncorrelated with brain activity. Second, its temporal behavior can be characterized by a few distinct activities, and as a result, the spatio-temporal matrix of the interference has a few distinctly large singular values. Namely, the interference can be modeled as a low-rank signal. Under these assumptions, our analysis shows that the adaptive beamformer techniques are insensitive to interference when its spatial singular vectors are so different from a lead field vector of a brain source that the generalized cosine between these two vectors is much smaller than unity. Four types of numerical examples verifying this conclusion are presented.

Muscle afferent inputs from the hand activate human cerebellum sequentially through parallel and climbing fiber systems

OBJECTIVE

Spatio-temporal response characteristics of the human cerebellum to median nerve stimulation (MNS) were studied with the use of a whole-head magnetoencephalographic (MEG) system covering the cerebellum and upper cervical spine.

METHODS

Neuromagnetic responses from the cerebellum were recorded following electric stimulation of the right median nerve in 12 subjects. In 6 out of 12 subjects, the responses to the left median nerve and to the right index or middle finger stimulation were also recorded.

RESULTS

The medial part of the cerebellum (spinocerebellum) was activated by MNS. In contrast, there were no responses from the cerebellum to the finger stimulation, suggesting that muscle afferent inputs are the source of cerebellar activation for MNS. The cerebellar responses consisted of 3 or 4 components of alternating polarity within 90 ms post-stimulus: the current direction for the first component was from the depth to the surface of the anterior lobe.

CONCLUSIONS

From the timing and current direction, we speculate that the 4 components reflect, respectively, (1) excitatory postsynaptic potentials (EPSPs) of granule cells, (2) Purkinje cell EPSPs at the distal dendrites driven by parallel fibers, (3) Purkinje cell EPSPs at the soma and the proximal dendrites mediated by climbing fibers and (4) second Purkinje cell EPSPs at the distal dendrites driven by parallel fibers.

SIGNIFICANCE

We first visualized serial activation of the human spinocerebellum following MNS noninvasively with MEG.

Dynamic movement of N100m dipoles in evoked magnetic field reflects sequential activation of isofrequency bands in human auditory cortex

OBJECTIVE

To investigate spatiotemporal features of the isofrequency bands for 400 and 4000 Hz tones in human auditory cortex and on the hemispheric differences in the arrangement of isofrequency bands.

METHODS

We recorded auditory evoked magnetic fields (AEFs) to 400 or 4000 Hz tone pips presented at right or left ear from 31 normal subjects. The dipole location for the N100m sources was successively calculated from the AEFs obtained from the hemisphere contralateral to the stimulated ear.

RESULTS

In the right hemisphere, the current sources for 400 and 4000 Hz moved toward the anterolateral direction before the N100m peak, showing parallel arrangement of the isofrequency bands (4000 Hz in medial location). In the left hemisphere, the movement direction of 400 Hz dipoles was anterolateral, while that of 4000 Hz dipoles was lateral.

CONCLUSIONS

This difference in the organization of isofrequency bands between right and left auditory cortices reflects distinct functional roles in auditory information processing such as pitch vs. language discrimination.

SIGNIFICANCE

This work is the first to disclose isofrequency bands in human auditory cortex based on the analysis of magnetoencephalography.

Performance of an MEG adaptive-beamformer technique in the presence of correlated neural activities: effects on signal intensity and time-course estimates

The influence of temporarily correlated source activities on neuromagnetic reconstruction by adaptive beamformer techniques was investigated. It is known that the spatial filter weight of an adaptive beamformer cannot perfectly block correlated signals. This causes two major influences on the reconstruction results: time course distortions and reductions in reconstructed signal intensities. Our theoretical analysis and numerical experiments both showed that the reduction in signal intensity for sources with a medium degree of correlation is small. The time-course distortion for such sources, however, may be discernible. Our analysis also showed that the magnitude correlation coefficient between two correlated sources can be accurately estimated by using the beamformer outputs. A method of retrieving the original time courses using estimated correlation coefficients was developed. Our numerical experiments demonstrated that reasonably accurate time courses can be retrieved from considerably distorted time courses even when the signal-to-noise ratio is low.

Dipole orientation differs between high frequency oscillations and N20m current sources in human somatosensory evoked magnetic fields to median nerve stimulation

The dipole orientation of equivalent current source for the high frequency oscillations (HFOs) above 300 Hz and that for the underlying N20m were compared. Somatic magnetic fields were recorded over the left hand somatosensory area to right median nerve stimulation at the wrist with a wide-bandpass (0.1-2000 Hz). The HFOs and underlying N20m were extracted by digital filtering of 300-900 Hz and 1-300 Hz, respectively. We found that the orientation of the HFOs and underlying N20m current sources differs and that the HFO source orientation shows a more divergent pattern than the N20m. These results suggest that the somatosensory HFOs are not generated from the pyramidal cell population in area 3b which produces the underlying N20m and that they may reflect activities of the non-pyramidal neuron population.

Dynamic activation of distinct cytoarchitectonic areas of the human SI cortex after median nerve stimulation

MEG recordings visualized non-invasively a serial mediolateral activation of the human somatosensory 3b area followed by a stationary activation of area 1 after median nerve stimulation. Somatosensory evoked fields (SEFs) were recorded over the hand area contralateral to the right median nerve stimulation at the wrist in six normal subjects. A newly developed MEG vector beamformer technique applied to the SEFs revealed two distinct sources (areas 3b and 1) in the primary somatosensory cortex (SI) during the primary N20m-P22m response in all subjects. The first source was located in area 3b, which started to move sequentially toward mediolateral direction 0.7 ms prior to the peak of N20m and ended its movement 1.4 ms after the peak with a total distance of 11.2 mm. We speculate that the movement reflects a sequential mediolateral activation of the pyramidal cells in area 3b, which is mediated by horizontal connections running parallel to the cortical surface. The second source in area 1, located 5.6 mm medial and 4.2 mm posterior to the first source, was active 1.0 ms after the N20m peak. Then, the first source became inactive and the second source was dominant. In sharp contrast with the first source, the second source was stationary. The different behavior of these two components (moving vs stationary) indicates independent parallel inputs to area 3b and area 1 from the thalamus.