Comparing localization of conventional functional magnetic resonance imaging and magnetoencephalography

Enhanced magnetic transduction of neuronal activity by nanofabricated inductors quantified via finite element analysis

Objective.Methods for the detection of neural signals involve a compromise between invasiveness, spatiotemporal resolution, and the number of neurons or brain regions recorded. Electrode-based probes provide excellent response but usually require transcranial wiring and capture activity from limited neuronal populations. Noninvasive methods such as electroencephalography and magnetoencephalography offer fast readouts of field potentials or biomagnetic signals, respectively, but have spatial constraints that prohibit recording from single neurons. A cell-sized device that enhances neurogenic magnetic fields can be used as anin situsensor for magnetic-based modalities and increase the ability to detect diverse signals across multiple brain regions.Approach.We designed and modeled a device capable of forming a tight electromagnetic junction with single neurons, thereby transducing changes in cellular potential to magnetic field perturbations by driving current through a nanofabricated inductor element.Main results.We present detailed quantification of the device performance using realistic finite element simulations with signals and geometries acquired from patch-clamped neuronsin vitroand demonstrate the capability of the device to produce magnetic signals readable via existing modalities. We compare the magnetic output of the device to intrinsic neuronal magnetic fields (NMFs) and show that the transduced magnetic field intensity from a single neuron is more than three-fold higher at its peak (1.62 nT vs 0.51 nT). Importantly, we report on a large spatial enhancement of the transduced magnetic field output within a typical voxel (40 × 40 × 10µm) over 250 times higher than the intrinsic NMF strength (0.64 nT vs 2.5 pT). We use this framework to perform optimizations of device performance based on nanofabrication constraints and material choices.Significance.Our quantifications institute a foundation for synthesizing and applying electromagnetic sensors for detecting brain activity and can serve as a general method for quantifying recording devices at the single cell level.

Sensorimotor Mapping With MEG: An Update on the Current State of Clinical Research and Practice With Considerations for Clinical Practice Guidelines

In this article, we present the clinical indications and advances in the use of magnetoencephalography to map the primary sensorimotor (SM1) cortex in neurosurgical patients noninvasively. We emphasize the advantages of magnetoencephalography over sensorimotor mapping using functional magnetic resonance imaging. Recommendations to the referring physicians and the clinical magnetoencephalographers to achieve appropriate sensorimotor cortex mapping using magnetoencephalography are proposed. We finally provide some practical advice for the use of corticomuscular coherence, cortico-kinematic coherence, and mu rhythm suppression in this indication. Magnetoencephalography should now be considered as a method of reference for presurgical functional mapping of the sensorimotor cortex.

Presurgical Functional Cortical Mapping Using Electromagnetic Source Imaging

Preoperative localization of functionally eloquent cortex (functional cortical mapping) is common clinical practice in order to avoid or reduce postoperative morbidity. This review aims at providing a general overview of magnetoencephalography (MEG) and high-density electroencephalography (hdEEG) based methods and their clinical role as compared to common alternatives for functional cortical mapping of (1) verbal language function, (2) sensorimotor cortex, (3) memory, (4) visual, and (5) auditory cortex. We highlight strengths, weaknesses and limitations of these functional cortical mapping modalities based on findings in the recent literature. We also compare their performance relative to other non-invasive functional cortical mapping methods, such as functional Magnetic Resonance Imaging (fMRI), Transcranial Magnetic Stimulation (TMS), and to invasive methods like the intracarotid Amobarbital Test (WADA-Test) or intracranial investigations.

Comparing the Microvascular Specificity of the 3- and 7-T BOLD Response Using ICA and Susceptibility-Weighted Imaging

In functional MRI it is desirable for the blood-oxygenation level dependent (BOLD) signal to be localized to the tissue containing activated neurons rather than the veins draining that tissue. This study addresses the dependence of the specificity of the BOLD signal – the relative contribution of the BOLD signal arising from tissue compared to venous vessels – on magnetic field strength. To date, studies of specificity have been based on models or indirect measures of BOLD sensitivity such as signal to noise ratio and relaxation rates, and assessment has been made in isolated vein and tissue voxels. The consensus has been that ultra-high field systems not only significantly increase BOLD sensitivity but also specificity, that is, there is a proportionately reduced signal contribution from draining veins. Specificity was not quantified in prior studies, however, due to the difficulty of establishing a reliable network of veins in the activated volume. In this study we use a map of venous vessel networks extracted from 7 T high resolution Susceptibility-Weighted Images to quantify the relative contributions of micro- and macro-vasculature to functional MRI results obtained at 3 and 7 T. High resolution measurements made here minimize the contribution of physiological noise and Independent Component Analysis (ICA) is used to separate activation from technical, physiological, and motion artifacts. ICA also avoids the possibility of timing-dependent bias from different micro- and macro-vasculature responses. We find a significant increase in the number of activated voxels at 7 T in both the veins and the microvasculature – a BOLD sensitivity increase – with the increase in the microvasculature being higher. However, the small increase in sensitivity at 7 T was not significant. For the experimental conditions of this study, our findings do not support the hypothesis of an increased specificity of the BOLD response at ultra-high field.

Clinical magnetoencephalography for neurosurgery

Noninvasive neuroimaging aids in surgical planning and in counseling patients about possible risks of surgery. Magnetoencephalography (MEG) performs the most common types of surgical planning that the neurosurgeon faces, including localization of epileptic discharges, determination of the hemispheric dominance of verbal processing, and the ability to locate eloquent cortex. MEG is most useful when it is combined with structural imaging, most commonly with structural magnetic resonance (MR) imaging and MR diffusion imaging. This article reviews the history of clinical MEG, introduces the basic concepts about the biophysics of MEG, and outlines the basic neurosurgical applications of MEG.

Spontaneous locally restricted EEG alpha activity determines cortical excitability in the motor cortex

There is growing interest in the functional meaning of rhythmical brain activity. For oscillatory brain activity around 10 Hz in the human electroencephalogram (EEG) it is discussed whether it is associated with the level of cortical excitation. However, it is not clear whether the relation between 10 Hz EEG oscillatory activity and cortical excitability is a global, locally very unspecific phenomenon or whether focal 10 Hz oscillations in the human brain are a highly specific correlate of the cortical excitation level. To determine this open question, multichannel EEG was combined with transcranial magnetic stimulation (TMS) applied to the primary motor cortex in this study. The present data showed that a motor evoked potential was elicited more easily when alpha power immediately preceding the magnetic pulse was low, and vice versa. Interestingly, this effect was only found for very local EEG alpha activity at sites overlying the cortical motor areas to which the TMS pulses were applied. This was verified using source localization in 3D space. These data provide evidence that the magnitude of motor cortical excitability is determined by the amount of topographically specific alpha oscillations in the sensorimotor cortex.

Cardiac responses induced during thought-based control of a virtual environment

Cardiac responses induced by motor imagery were investigated in 3 subjects in a series of experiments with a synchronous (cue-based) Brain-Computer Interface (BCI). The cue specified right hand vs. leg/foot motor imagery. After a number of BCI training sessions reaching a classification accuracy of at least 80%, the BCI experiments were carried out in an immersive virtual environment (VE), commonly referred as a "CAVE". In this VE, the subjects were able to move along a virtual street by motor imagery alone. The thought-based control of VE resulted in an acceleration of the heart rate in 2 subjects and a heart rate deceleration in the other subject. In control experiments in front of a PC, all 3 subjects displayed a significant heart rate deceleration of the order of about 3-5%. This heart rate decrease during motor imagery in a normal environment is similar to that observed during preparation for a voluntary movement. The heart rate acceleration in the VE is interpreted as effect of an increased mental effort to walk as far as possible in VE.

Spatial resolution of fMRI in the human parasylvian cortex: Comparison of somatosensory and auditory activation

In spite of its outstanding spatial resolution, the biological resolution of functional MRI may be worse because it depends on the vascular architecture of the brain. Here, we compared the activation patterns of the secondary somatosensory and parietal ventral cortex (SII/PV) with that of the primary auditory cortex and adjacent areas (AI/AII). These two brain regions are located immediately adjacent to each other on opposite banks of the Sylvian fissure, and are anatomically and functionally distinct. In 12 healthy subjects, SII/PV was activated by pneumatic tactile stimuli applied to the index finger (0.5 cm2 contact area, 4 bar pressure), and AI/AII by amplitude-modulated tones (800 Hz carrier frequency, modulated at 24-36 Hz). Functional images were obtained with a 1.5-T scanner and were evaluated using SPM99. Sensitivity of fMRI activation in this unselected sample was 71% for tactile and 83% for auditory stimulation. Group analysis showed activation of SII/PV by tactile and activation of three locations in AI/AII by auditory stimuli. Distributions extended to the opposite side of the fissure (19-58% after tactile and 13-14% after auditory stimulation, depending on the side of stimulation/hemisphere). Morphometry of individual sulcal anatomy revealed that the course of the Sylvian fissure varied by 5.3 mm (SD) in vertical direction. Taking this into account, SII/PV was located 5.8 +/- 2.7 mm above the Sylvian fissure, whereas AI/AII was located 6.3 +/- 1.7 mm below the Sylvian fissure. Even in individual analysis, the most significant voxel after tactile stimuli in one subject was found on the "wrong" side of the fissure; this error could be ascribed to the spatial normalization procedure. These data show that fMRI signals may overlap substantially, even if the activated regions are separated by 12 mm across a major sulcus. Spatial normalization to an atlas template can introduce additional variance. Individual sulcal anatomy should be preferred over mean atlas locations.

Contralateral and ipsilateral responses in primary somatosensory cortex following electrical median nerve stimulation—an fMRI study

OBJECTIVE

Ten healthy adult subjects were examined using functional magnetic resonance imaging (fMRI) to investigate responses in the contralateral and ipsilateral primary somatosensory cortex (SI) following electrical stimulation of the median nerve.

METHODS

The right and left median nerves were stimulated alternately at the wrist in the different sessions. First, the location of the response in contralateral SI was identified following median nerve stimulation, and then, a spherical search volume with a 10mm radius centered on the region of the contralateral response was determined. Whether or not fMRI activation occurred within this sphere following ipsilateral stimulation was examined using a 3T MR imager.

RESULTS

A response in contralateral SI was observed in 8 of the 10 subjects in right and left hemisphere. Responses in ipsilateral SI were observed in 6 of 8 subjects in right hemisphere, and the region of the response tended to be posterior to the contralateral region. On the other hand, in left hemisphere, the ipsilateral responses were found in three.

CONCLUSIONS

In the present study, not only contralateral SI but also ipsilateral SI was activated following median nerve. The location of the ipsilateral activation was significantly more posterior than the contralateral one in right hemisphere.

SIGNIFICANCE

The region of activation in ipsilateral SI was located in the posterior portion of post central gyrus, corresponding to around BA2 and 5 in human.

Magnetic source imaging supports clinical decision making in glioma patients

OBJECTIVE

This study addresses the potential utility of preoperative functional imaging with magnetoencephalography (MEG) for the selection of glioma patients who are likely to benefit from resective surgical treatment regarding postoperative morbidity.

METHODS

One hundred and nineteen patients with gliomas adjacent to sensorimotor, visual and speech related brain areas were investigated preoperatively with a MAGNES II biomagnetometer. In each patient the pre-surgical evaluation was focussed on the visual, sensorimotor cortex and/or of the speech related brain areas. A grading system was then used according to the distance of the MEG activation sources to the nearest tumour border to determine the further treatment. The therapeutic options consisted in conservative treatment, stereotactic biopsy and/or a radiation and chemotherapy, substantial cytoreduction and the gross total removal of the lesion.

RESULTS

From 119 investigated patients, 55 patients (46.2%) were not considered for surgery due to tumour invasion to functional cortex. Sixty four patients (53.8%) were chosen for resective surgery. In the surgical group only four patients (6.2%) suffered from neurological deterioration.

CONCLUSIONS

Magnetic source imaging (MSI) proved to be a valuable help in the clinical decision making process of lesions adjacent to functional important brain areas. The relative high number of patients in whom MSI warns of the postoperative crippling sequelae may lead to a better selection of patients who benefit from resective surgery. This method may help to find the patients for whom conservative treatment seems to be more favourable concerning quality of life in the surviving time.

Functional neuronavigation and intraoperative MRI

Our concept of computer assisted surgery is based on the combination of intraoperative magnetic resonance (MR) imaging with microscope-based neuronavigation, providing anatomical and functional guidance simultaneously. Intraoperative imaging evaluates the extent of a resection, while the additional use of functional neuronavigation, which displays the position of eloquent brain areas in the operative field, prevents increasing neurological deficits, which would otherwise result from extended resections. Up to mid 2001 we performed intraoperative MR imaging using a low-field 0.2 Tesla scanner in 330 patients. The main indications were the evaluation of the extent of resection in gliomas, pituitary tumours, and in epilepsy surgery. Intraoperative MR imaging proved to serve as intraoperative quality control with the possibility of an immediate modification of the surgical strategy, i.e. extension of the resection. Integrated use of functional neuronavigation prevented increased neurological deficits. Compared to routine pre- or postoperative imaging being performed with high-Tesla machines, intraoperative image quality and sequence spectrum could not compete. This led to the development of the concept to adapt a high-field MR scanner to the operating environment, preserving the benefits of using standard microsurgical equipment and microscope-based neuronavigational guidance with integrated functional data, which was successfully implemented by April 2002. Up to the end of 2002, 95 patients were investigated with the new setup. Improved image quality, intraoperative workflow, as well as enhanced sophisticated intraoperative imaging possibilities are the major benefits of the high-field setup.

Comparison of memory fMRI response among normal, MCI, and Alzheimer's patients

OBJECTIVE

To determine whether an fMRI memory encoding task distinguishes among cognitively normal elderly individuals, patients with mild cognitive impairment (MCI), and patients with early Alzheimer's disease (AD).

METHODS

Twenty-nine subjects (11 normal, 9 MCI, 9 AD) were studied with an fMRI memory encoding task. A passive sensory task was also performed to assess potential intergroup differences in fMRI responsiveness. Activation in the medial temporal lobe for the memory task and in the anatomic rolandic area for the sensory task was studied. Intergroup comparisons were performed using receiver operating characteristic (ROC) analyses. The ROC method provides rigorous control of artifactual false-positive "activation." Subjects were tested for recall and recognition of the encoding task stimuli following the fMRI study.

RESULTS

Medial temporal lobe activation was greater in normal subjects than MCI and AD patients (p = 0.03 and p = 0.04). There was no difference between AD and MCI patients in fMRI memory performance [corrected]. There was an association between fMRI memory activation (area under the ROC curve) and post-fMRI performance on recognition and free recall. There was no difference among the three groups on the sensory task.

CONCLUSIONS

MCI and AD patients had less medial temporal lobe activation on the memory task than the normal subjects but similar activation as normal subjects on the sensory task. These findings suggest decreased medial temporal activation may be a specific marker of limbic dysfunction due to the neurodegenerative changes of AD. In addition, fMRI is sufficiently sensitive to detect changes in the prodromal, MCI, phase of the disease.

Combining MEG and MRI with neuronavigation for treatment of an epileptiform spike focus in the precentral region: a technical case report

BACKGROUND

Epileptic foci are often located in the vicinity but not necessarily within the boundaries of intra-axial brain tumors. Resection of these tumors is based on two major goals: first, maximizing tumor removal without provoking new neurologic deficits, and second, minimizing epileptic seizure activity. Magnetic source imaging (MSI) depicts the generators of magnetic fields overlaid on individual magnetic resonance (MR) images. Established application areas are lesions located adjacent to or partly within the sensory and motor cortex, or in the depth of the brain, necessitating a surgical approach through functionally highly relevant cortical regions. Magnetoencephalography (MEG) is also applicable for epileptiform spike foci recording during interictal activity.

CASE DESCRIPTION

A patient with a recurrent glioma close to the Rolandic cortex scheduled for epilepsy and tumor surgery was investigated with MSI. The MSI data showed an epileptiform spike focus outside the tumor boundaries. The resulting MSI images were integrated into our neuronavigation system. This procedure allowed for the preoperative identification of the sensory and motor cortex, the precise localization of the epileptiform spike focus, and careful planning of the surgical procedure. In this case, we were able to safely resect the recurrent tumor and the epileptiform spike focus under general anesthesia using MSI-based neuronavigational guidance but no conventional intraoperative mapping techniques.

CONCLUSION

Magnetic source imaging can be a valuable, noninvasive method for planning and performing tumor resections in high-risk brain regions, especially if an epileptiform spike focus has to be localized and included into the resection strategy.

Functional magnetic resonance imaging evidence for a representation of the ear in human primary somatosensory cortex: comparison with magnetoencephalography study

Our previous study (T. Nihashi et al., 2001, Neuro- Image 13: 295-304), using magnetoencephalography (MEG), revealed somatotopy of the ear in the human primary somatosensory cortex (SI); that is, the signals following stimulation of the ear reach both the neck and face areas of the SI. However, since this was the first report on somatotopy of the ear in humans, we performed an fMRI activation study to confirm the somatotopic representation, and compared the electrical activity by MEG and the cerebral blood flow change by fMRI. We studied eight healthy subjects using 3-T MRI. We stimulated three parts of the left ear: the helix, the lobulus, and the tragus. First, we identified the location of the ear area in the SI based on our previous MEG study, in which equivalent current dipoles (ECDs) were located in the neck and/or face areas of the SI. Then, we determined the search volume as a sphere with a 15-mm radius, which was placed in the neck and/or face area. We analyzed whether or not fMRI activation occurred inside such spheres. Stimulation of the helix activated the neck area of the SI in four of eight subjects, and both the neck and face areas in two. No activation was observed in two subjects. Stimulation of the lobulus activated the neck area in one subject, the face area in two, both in four, and neither in one. Stimulation of the tragus activated the face in four, both in three, and neither in one. These fMRI findings confirm the result of MEG that the representation of the ear in the SI is separated into neck and face areas.

Relationship of spikes, synaptic activity, and local changes of cerebral blood flow

The coupling of electrical activity in the brain to changes in cerebral blood flow (CBF) is of interest because hemodynamic changes are used to track brain function. Recent studies, especially those investigating the cerebellar cortex, have shown that the spike rate in the principal target cell of a brain region (i.e. the efferent cell) does not affect vascular response amplitude. Subthreshold integrative synaptic processes trigger changes in the local microcirculation and local glucose consumption. The spatial specificity of the vascular response on the brain surface is limited because of the functional anatomy of the pial vessels. Within the cortex there is a characteristic laminar flow distribution, the largest changes of which are observed at the depth of maximal synaptic activity (i.e. layer IV) for an afferent input system. Under most conditions, increases in CBF are explained by activity in postsynaptic neurons, but presynaptic elements can contribute. Neurotransmitters do not mediate increases in CBF that are triggered by the concerted action of several second messenger molecules. It is important to distinguish between effective synaptic inhibition and deactivation that increase and decrease CBF and glucose consumption, respectively. In summary, hemodynamic changes evoked by neuronal activity depend on the afferent input function (i.e. all aspects of presynaptic and postsynaptic processing), but are totally independent of the efferent function (i.e., the spike rate of the same region). Thus, it is not possible to conclude whether the output level of activity of a region is increased based on brain maps that use blood-flow changes as markers.

Correlation of sensorimotor activation with functional magnetic resonance imaging and magnetoencephalography in presurgical functional imaging: a spatial analysis

In this study we investigated the spatial heterotopy of MEG and fMRI localizations after sensory and motor stimulation tasks. Both methods are frequently used to study the topology of the primary and secondary motor cortex, as well as a tool for presurgical brain mapping. fMRI was performed with a 1.5T MR system, using echo-planar imaging with a motor and a sensory task. Somatosensory and motor evoked fields were recorded with a biomagnetometer. fMRI activation was determined with a cross-correlation analysis. MEG source localization was performed with a single equivalent current dipole model and a current density localization approach. Distances between MEG and fMRI activation sites were measured within the same anatomical 3-D-MR image set. The central region could be identified by MEG and fMRI in 33 of 34 cases. However, MEG and fMRI localization results showed significantly different activation sites for the motor and sensory task with a distance of 10 and 15 mm, respectively. This reflects the different neurophysiological mechanisms: direct neuronal current flow (MEG) and secondary changes in cerebral blood flow and oxygenation level of activated versus non activated brain structures (fMRI). The result of our study has clinical implications when MEG and fMRI localizations are used for pre- and intraoperative brain mapping. Although both modalities are useful for the estimation of the motor cortex, a single modality may err in the exact topographical labeling of the motor cortex. In some unclear cases a combination of both methods should be used in order to avoid neurological deficits.

Spatial localization and resolution of BOLD fMRI

It has been demonstrated that the blood-oxygenation-level-dependent (BOLD) fMRI initial dip allows us to resolve (without differential subtraction) structures of the order of 0.5 mm. However, recent results support the proposition that even the later, positive BOLD fMRI signal component can allow us to resolve structures less than 1 mm in size by using differential subtraction when the signal-to-noise ratio is high. So, with a sufficient signal-to-noise ratio, the later, positive component should be useable as a probe for testing cognitive neuroscientific hypotheses that predict neuroanatomical dissociations of less than 1mm.

Quantification of fMRI artifact reduction by a novel plaster cast head holder

In light of artifact-induced high variability of activation in fMRI repeat studies, we developed and tested a clinically useful plaster cast head holder (PCH) with improved immobilization, repositioning, and comfort. With PCH, there were considerably lower levels of translational and rotational head motion components compared to head fixation with conventional restraining straps (CRS). Rotational components cannot be fully compensated by realignment and lead to "false activations." In addition, task-correlated head motion, which highly increases the risk of artifacts, was considerably reduced with PCH, especially in a motion prone subject. Compared with PCH, head motion was 133% larger with CRS in a highly cooperative subject. With a motion prone subject, head motion range was increased by 769% (PCH: 0.9 mm, CRS: 7.8 mm), which may indicate the usefulness of PCH for restless patients. In functional activation maps, PCH alone yielded fewer residual motion artifacts than CRS + image registration. Subject tolerance of the head holder during the long measurement times of up to 2.5 hr was good, and slice orientation on different days confirmed the quality of repositioning.

Spatiotemporal activity of a cortical network for processing visual motion revealed by MEG and fMRI

A sudden change in the direction of motion is a particularly salient and relevant feature of visual information. Extensive research has identified cortical areas responsive to visual motion and characterized their sensitivity to different features of motion, such as directional specificity. However, relatively little is known about responses to sudden changes in direction. Electrophysiological data from animals and functional imaging data from humans suggest a number of brain areas responsive to motion, presumably working as a network. Temporal patterns of activity allow the same network to process information in different ways. The present study in humans sought to determine which motion-sensitive areas are involved in processing changes in the direction of motion and to characterize the temporal patterns of processing within this network of brain regions. To accomplish this, we used both magnetoencephalography (MEG) and functional magnetic resonance imaging (fMRI). The fMRI data were used as supplementary information in the localization of MEG sources. The change in the direction of visual motion was found to activate a number of areas, each displaying a different temporal behavior. The fMRI revealed motion-related activity in areas MT+ (the human homologue of monkey middle temporal area and possibly also other motion sensitive areas next to MT), a region near the posterior end of the superior temporal sulcus (pSTS), V3A, and V1/V2. The MEG data suggested additional frontal sources. An equivalent dipole model for the generators of MEG signals indicated activity in MT+, starting at 130 ms and peaking at 170 ms after the reversal of the direction of motion, and then again at approximately 260 ms. Frontal activity began 0-20 ms later than in MT+, and peaked approximately 180 ms. Both pSTS and FEF+ showed long-duration activity continuing over the latency range of 200-400 ms. MEG responses in the region of V3A and V1/V2 were relatively small, and peaked at longer latencies than the initial peak in MT+. These data revealed characteristic patterns of activity in this cortical network for processing sudden changes in the direction of visual motion.

Activation of multiple cortical areas in response to somatosensory stimulation: combined magnetoencephalographic and functional magnetic resonance imaging

We combined information from functional magnetic resonance imaging (fMRI) and magnetoencephalography (MEG) to assess which cortical areas and in which temporal order show macroscopic activation after right median nerve stimulation. Five healthy subjects were studied with the two imaging modalities, which both revealed significant activation in the contra- and ipsilateral primary somatosensory cortex (SI), the contra- and ipsilateral opercular areas, the walls of the contralateral postcentral sulcus (PoCS), and the contralateral supplementary motor area (SMA). In fMRI, two separate foci of activation in the opercular cortex were discerned, one posteriorly in the parietal operculum (PO), and one anteriorly near the insula or frontal operculum (anterior operculum, AO). The activation sites from fMRI were used to constrain the solution of the inverse problem of MEG, which allowed us to construct a model of the temporal sequence of activation of the different sites. According to this model, the mean onset latency for significant activation at the contralateral SI was 20 msec (range, 17-22 msec), followed by activation of PoCS at 23 msec (range, 21-25 msec). The contralateral PO was activated at 26 msec (range, 19-32 msec) and AO at 33 msec (range, 22-51 msec). The contralateral SMA became active at 36 msec (range, 24-48 msec). The ipsilateral SI, PO, and AO became activated at 54-67 msec. We conclude that fMRI provides a useful means to constrain the inverse problem of MEG, allowing the construction of spatiotemporal models of cortical activation, which may have significant implications for the understanding of cortical network functioning.

Functional neuronavigation with magnetoencephalography: outcome in 50 patients with lesions around the motor cortex

OBJECT

The authors conducted a study to evaluate the clinical outcome in 50 patients with lesions around the motor cortex who underwent surgery in which functional neuronavigation was performed.

METHODS

The sensorimotor cortex was identified in all patients with the use of magnetoencephalography (MEG). The MEG-source localizations were superimposed onto a three-dimensional magnetic resonance image and the image data set was implemented into a neuronavigation system. Based on this setup, the surgeon chose the best surgical strategy. During surgery, the pre- and postcentral gyri were identified by neuronavigation and, in addition, the central sulcus was localized using intraoperative recording of somatosensory evoked potentials. In all cases MEG localizations of the sensory or motor cortex were correct. In 30% of the patients preoperative paresis improved, in 66% no additional deficits occurred, and in only 4% (two patients) deterioration of neurological function occurred. In one of these patients the deterioration was not related to the procedure.

CONCLUSIONS

The method of incorporating functional data into neuronavigation systems is a promising tool that can be used in more radical surgery to lessen morbidity around eloquent brain areas.

Integration of functional magnetic resonance imaging supported by magnetoencephalography in functional neuronavigation

OBJECTIVE: In this study, the intraoperative visualization of functional data provided by functional magnetic resonance imaging (fMRI) and magnetoencephalography (MEG) leading to functional neuronavigation is demonstrated in surgery around the motor strip. METHODS: In seven patients with lesions adjacent to the central region, fMRI was performed with a 1.5-Tesla magnetic resonance system, using axial echo-planar imaging with a motor and a sensory task. Somatosensory and motor evoked fields were recorded with a biomagnetometer. fMRI and MEG were matched to an anatomic three-dimensional magnetic resonance image set by a contour fit. Then this three-dimensional image data set was transferred to the navigation microscope and displayed in the eyepieces of the microscope during surgery. Additionally, intraoperative recording of somatosensory evoked potentials was performed for verification of the central sulcus. RESULTS: In all cases, the projection of fMRI and MEG data into the operating viewing field allowed easy identification of the central region, which was confirmed by phase reversal of somatosensory evoked potentials in each case. fMRI and MEG measurements yielded corresponding results in each patient. CONCLUSION: Functional neuronavigation with integration of fMRI and MEG allows the fast identification of eloquent brain areas. The widespread availability of fMRI will result in a broad availability of functional neuronavigation, which will, in turn, contribute to the successful surgery of lesions in eloquent brain areas with lower morbidity.

Functional neuronavigation with magnetoencephalography: outcome in 50 patients with lesions around the motor cortex

The authors conducted a study to evaluate the clinical outcome in 50 patients with lesions around the motor cortex who underwent surgery in which functional neuronavigation was performed.

The sensorimotor cortex was identified in all patients with the use of magnetoencephalography (MEG). The MEG-source localizations were superimposed onto a three-dimensional magnetic resonance image, and the image data set was then implemented into a neuronavigation system. Based on this setup, the surgeon chose the best surgical strategy. During surgery, the pre- and postcentral gyrus were identified by neuronavigation, and in addition, the central sulcus was localized using intraoperative recording of somatosensory evoked potentials. In all cases MEG localizations of the sensory or motor cortex were correct. In 30% of the patients preoperative paresis improved, in 66% no additional deficits occurred, and in only 4% (two patients) deterioration of neurological function occurred. In one of these patients the deterioration was not related to the method.

The method of incorporating functional data into neuronavigation systems is a promising tool that can be used in more radical surgery to cause less morbidity around eloquent brain areas.

The influence of stimulus properties, complexity, and contingency on the stability and variability of ongoing and evoked activity in human auditory cortex

The real-time, single-trial activity in the human auditory cortex was extracted from magnetoencephalographic signals. A predictor of single-trial activity was defined as the sum of the average response and a mean-free base level computed over a range of base times. For simple stimuli the residual (predicted-actual) activity had a stimulus-independent oscillatory (10 Hz) component. This component was larger and more durable in trained subjects, reaching saturation only in the most trained of the five subjects studied (S1). Changes in variability and associated reduction of the absolute value and duration of the oscillations were evident in experiments with stimuli loaded with information, saliency, or task contingency. Repetition reintroduces stimulus-independent oscillations very slowly. For S1, after training, the stimulus-independent oscillations were reestablished in the auditory cortices to the level seen for simple stimuli, except for the time periods and in the hemisphere associated with the combination of task demands and stimulus processing.

Anatomical congruence of metabolic and electromagnetic activation signals during a self-paced motor task: a combined PET-MEG study

We have investigated the degree of spatial correlation between the cerebral blood flow variations measured by positron emission tomography (PET) and the electromagnetic sources as measured by magnetoencephalography (MEG) in five subjects while performing a self-paced right index finger tapping task. Data were processed independently for each technique using both single-case and intersubject analysis. PET and MEG were coregistered with anatomical magnetic resonance images for each subject. Both extension and flexion motor-related fields were extracted from the MEG signal. Using the single dipole model we identified the motor evoked field 1 (MEF1) in all subjects and the motor field (MF) in three subjects. Individual and intersubject averaged PET data showed consistent contralateral primary sensorimotor (PSM) hand area and bilateral supplementary motor area activation. MEG individual and intersubject averaged results demonstrated that both MEF1 and MF dipoles were localized within the PSM PET activated area. Individual PSM mass center to dipole distance was 12 and 15.3 mm on average for the MEF1 and the MF component, respectively. For the same components, the intersubject averaged analysis shows distances between the PET Z-score maximum and the dipole locations of 6.3 and 15.0 mm, respectively. These results show that PET and MEG MEF1 activation signals spatially coincide within instrumental, registration, and modeling errors.

A comparison between electric source localisation and fMRI during somatosensory stimulation

The present study investigates the electric source localisation of somatosensory evoked potentials (SEPs) using high-resolution EEG (61 scalp electrodes) considering the individual brain morphology as obtained from magnetic resonance images (MRI). A comparison with the activation maps in fMRI under the same somatosensory stimulation paradigm was done. The somatosensory evoked potentials (SEPs) to electrical stimulation of the right median nerve were collected from the scalp of 8 healthy right-handed subjects. The source reconstruction for the 20 ms SEP component was performed by using a single moving dipole model as a source model and a spherical three-shell model as a head model. In 6 of the subjects fMRI was performed using the same electric stimulation of the right median nerve. The source location of the 20 ms SEP component was found to be within the postcentral gyrus. The fMRI activation maps were also located in the postcentral gyrus when using the same somatosensory stimulation paradigm. The appropriateness of using high-resolution EEG and fMRI in the functional localisation of the primary somatosensory cortex is discussed.

Localizing the central sulcus by functional magnetic resonance imaging and magnetoencephalography

To further validate the potential of functional magnetic resonance imaging (fMRI) for localization of the sensorimotor cortex, fMRI was compared with somatosensory evoked fields (SEFs) in eight normal volunteers. A conventional 1.5 T MRI scanner and an MRI-linked 66-channel whole head magnetoencephalography system were used. fMRI activated by unilateral hand squeeze movement indicated the highest activation on the central sulci that were localized by SEFs in all 16 contralateral hemispheres. This indicates that although the fMRI signal activation may originate from a vein running along the central sulcus, fMRI is reliable to detect the central sulcus. The pre-central gyrus also indicated some signal activation on fMRI implying better visualization of spatial distribution of activation. fMRI and SEFs are complementary methods for localizing the central sulcus.

Fuzzy clustering of gradient-echo functional MRI in the human visual cortex. Part II: quantification

Fuzzy cluster analysis (FCA) is a new exploratory method for analyzing fMRI data. Using simulated functional MRI (fMRI) data, the performance of FCA, as implemented in the software package Evident, was tested and a quantitative comparison with correlation analysis is presented. Furthermore, the fMRI model fit allows separation and quantification of flow and blood oxygen level dependent (BOLD) contributions in the human visual cortex. In gradient-recalled echo fMRI at 1.5 T (TR = 60 ms, TE = 42 ms, radiofrequency excitation flip angle [theta] = 10 degrees-60 degrees) total signal enhancement in the human visual cortex, ie, flow-enhanced BOLD plus inflow contributions, on average varies from 5% to 10% in or close to the visual cortex (average cerebral blood volume [CBV] = 4%) and from 100% to 20% in areas containing medium-sized vessels (ie, average CBV = 12% per voxel), respectively. Inflow enhancement, however, is restricted to intravascular space (= CBV) and increases with increasing radiofrequency (RF) flip angle, whereas BOLD contributions may be obtained from a region up to three times larger and, applying an unspoiled gradient-echo (GRE) sequence, also show a flip angle dependency with a minimum at approximately 30 degrees. This result suggests that a localized hemodynamic response from the microvasculature at 1.5 T may be extracted via fuzzy clustering. In summary, fuzzy clustering of fMRI data, as realized in the Evident software, is a robust and efficient method to (a) separate functional brain activation from noise or other sources resulting in time-dependent signal changes as proven by simulated fMRI data analysis and in vivo data from the visual cortex, and (b) allows separation of different levels of activation even if the temporal pattern is indistinguishable. Combining fuzzy cluster separation of brain activation with appropriate model calculations allows quantification of flow and (flow-enhanced) BOLD contributions in areas with different vascularization.

Quantification of signal changes in gradient recalled echo FMRI

Understanding and quantifying the various contributions to functional magnetic resonance imaging (FMRI) signal changes in activated cortical areas is paramount for a clinical application of brain mapping by FMRI. Therefore, all significant contributions to FMRI signal changes, both extra- and intravascular, from macrovessels down to the capillary network, should be taken into account. We present a gradient-recalled-echo FMRI model based on in-flow effects described by the Bloch equations, adding susceptibility effects empirically via T2* differences measured in vitro in human blood samples. Results of these calculations (by systematically varying alpha, echo time (TE), repetition time (TR), as well as blood velocity and T2* upon stimulation) may be used to (a) simulate functional MRI experiments with different measurement protocols and (b) estimate realistic values for important anatomical and physiological details that influence local signal changes in FMRI (i.e., size and distribution of vessels, effective relaxation times of blood, etc.). The excellent agreement between our model calculations and experimental results from conventional gradient recalled echo fMRI in vivo suggests a significant contribution from very slow flow and oxygenation changes, predominantly in small vessels (vblood = 1-4 mm/s). The actual contribution of T1- and T2-related effects is strongly dependent on sequence design and actual sequence parameters used. Thus, the model simulations presented may also be used to optimize measurement protocols for investigating various neurophysiological phenomena.

Magnetoencephalography may help to improve functional MRI brain mapping

The validity of functional magnetic resonance imaging (FMRI) brain maps with respect to the sites of neuronal activation is still unknown. One source of localization error may be pixels with large signal amplitudes, since such pixels may be expected to overlie large vessels, running remote from the centre of neuronal activation. In this study, magnetoencephalography was used to determine the centre of neuronal activation in a simple finger tapping task. The localization accuracy of conventional FMRI depending on FMRI signal enhancement was investigated relative to the magnetoencephalography reference. The results show a deterioration of FMRI localization with increasing signal amplitude related to increased contributions from large vessels. We conclude that FMRI data analysis should exclude large signal amplitudes and that magnetoencephalography may help to improve FMRI brain mapping results in a multimethod approach.

Magnetoencephalographic evidence for common sources of long latency fields to rare target and rare novel visual stimuli

This study used magnetoencephalography to examine the possibility that different generators account for the long-latency event-related potential (P300), evoked by rare target and by rare non-target, novel visual stimuli, in a visual oddball counting task performed by seven subjects. As expected, P300 peak latency was longer in response to rare targets compared to novel, non-target stimuli. Two main source regions were found for the Target- as well as for the Novel-P300, one in the temporal and one in the occipital lobe. Centers of neural activity were observed in the vicinity of the superior temporal sulcus, in the hippocampal formation and parahippocampal gyrus and in the occipital extrastriate cortex. It appears that the brain structures which contributed to the generation of the P300 response to both the target and the novel visual stimuli overlapped to a great extent.

Single-epoch neuromagnetic signals during epileptiform activities in guinea pig longitudinal CA3 slices

Neuromagnetic fields with high signal-to-noise ratios can be measured above longitudinal CA3 slices of the guinea pig during single epochs of interictal- and ictal-like synchronized population activities. Technical refinements enabled us to reduce the number of epochs for clear responses from over 5000 needed in an earlier study to single epochs. Simultaneous recording of field potentials revealed that neuromagnetic fields reflect intracellular currents in the pyramidal cells. Intracellular currents could be estimated from the external magnetic fields during paroxysmal depolarization shifts, multiple bursts, and slowly varying potential shifts lasting several seconds.

Reproducibility and postprocessing of gradient-echo functional MRI to improve localization of brain activity in the human visual cortex

High reproducibility of human FMRI studies is imperative for potential clinical applications of this new method for mapping human brain functions. So far, published data are not comparable quantitatively (even at the same field strength) as differences in sequence design and parameters as well as statistical methods applied to enhance function related image contrast, in particular, to extract the size of the "activated areas," are manifold. We present a study on reproducibility of gradient-echo FMRI in the human visual cortex using the different threshold strategies for correlation analysis that shows that, (a) applying adaptive correlation thresholds results in higher reproducibility compared to a fixed (0.5) threshold; (b) highly reproducible data can be obtained on a clinical 1.5 T MRI system, at least for repeated single subject studies (i.e., standard deviation of 2-30% for signal enhancement in 72-94% of the studies and 5-50% for activated area size in 63-88% of the studies, respectively, depending on threshold strategies); however, depending also on subject cooperation; (c) reproducibility across groups (alpha = const.) is worse, i.e., standard deviations are within 33-45% for signal enhancement and 41-74% for activated area size, respectively; (d) SNR is maximum at about 30 degrees flip angle, suggesting significant contributions from T1-effects for larger flip angles. Various technical, methodological, and physiological factors are influencing variability of signal enhancement and apparently activated area size, which should be taken into account if interpreting FMRI data quantitatively.