A new approach to neuroimaging with magnetoencephalography

An extended motor network generates beta and gamma oscillatory perturbations during development

This study examines the time course and neural generators of oscillatory beta and gamma motor responses in typically-developing children. Participants completed a unilateral flexion-extension task using each index finger as whole-head magnetoencephalography (MEG) data were acquired. These MEG data were imaged in the frequency-domain using spatial filtering and the resulting event-related synchronizations and desynchronizations (ERS/ERD) were subjected to voxel-wise statistical analyses to illuminate time-frequency specific activation patterns. Consistent with adult data, these children exhibited a pre-movement ERD that was strongest over the contralateral post-central gyrus, and a post-movement ERS response with the most prominent peak being in the contralateral precentral gyrus near premotor cortices. We also observed a high-frequency (approximately 80 Hz) ERS response that coincided with movement onset and was centered on the contralateral precentral gyrus, slightly superior and posterior to the beta ERS. In addition to pre- and post-central gyri activations, these children exhibited beta and gamma activity in supplementary motor areas (SMA) before and during movement, and beta activation in cerebellar cortices before and after movement. We believe the gamma synchronization may be an excellent candidate signal of basic cortical motor control, as the spatiotemporal dynamics indicate the primary motor cortex generates this response (and not the beta oscillations) which is closely yoked to the initial muscle activation. Lastly, these data suggest several additional neural regions including the SMA and cerebellum are involved in basic movements during development.

GABA(A) alpha-1 subunit mediated desynchronization of elevated low frequency oscillations alleviates specific dysfunction in stroke--a case report

OBJECTIVE

The paradoxical effects of the hypnotic imidazopyridine zolpidem, widely reported in persistent vegetative state, have been replicated recently in brain-injured and cognitively impaired patients. However, the neuronal mechanisms underlying these benefits are yet to be demonstrated. We implemented contemporary neuroimaging methods to investigate sensorimotor and cognitive improvements, observed in stroke patient JP following zolpidem administration.

METHODS

We used Magnetic-Resonance-Imaging (MRI) and Magnetic-Resonance-Spectroscopy (MRS) to anatomically and chemically characterize stroke damage. Single-photon-emission-computed-tomography (SPECT) and magnetoencephalography (MEG) were used to identify changes in cerebrovascular perfusion and neuronal network activity in response to sub-sedative doses of zolpidem, zopiclone and placebo. Cognitive improvements were measured using the WAIS-III and auditory-verbal tasks.

RESULTS

MRI and MRS revealed a lesion with complete loss of neuronal viability in the left temporal-parietal region; whilst SPECT indicated improved perfusion in the affected hemisphere following zolpidem. MEG demonstrated high-amplitude theta (4-10 Hz) and beta (15-30 Hz) oscillations within the peri-infarct region, which reduced in power coincident with zolpidem uptake and improvements in cognitive and motor function.

CONCLUSIONS

In JP, functional deficits and pathological oscillations appear coincidentally reduced following administration of low-dose zolpidem.

SIGNIFICANCE

GABA(A) alpha-1 sensitive desynchronisation of pathological oscillations may represent a biomarker and potential therapeutic target in brain injury.

Visual gamma oscillations and evoked responses: variability, repeatability and structural MRI correlates

There is increasing interest in the role gamma oscillations ( approximately 40 Hz) play in visual information processing. Despite this interest, and in contrast to the classically studied visual evoked potential, surprisingly little is known about the intra-individual repeatability of induced gamma oscillations. Similarly, little is known about inter-individual variability in terms of gamma oscillation frequency, bandwidth and amplitude with no extant normative data for these parameters. The purpose of the current study was therefore to examine the repeatability of visual gamma oscillations and to provide the first normative data on them. Our results demonstrate that evoked responses were highly repeatable across recording sessions whereas for induced visual gamma oscillations a large amount of inter-individual variability existed in terms of frequency, bandwidth and amplitude. However, these parameters and the general morphology of the gamma band response were stable within the same individuals for at least 4 weeks. The high degree of individual variability in gamma oscillations for gamma amplitude, bandwidth and frequency suggests that between-group studies on gamma oscillations will be difficult, requiring relatively large amounts of data to detect differences. However, the high degree of individual repeatability for gamma oscillation frequency, bandwidth and amplitude suggests that these dependent variables will be well suited for repeated-measure designs such as pharmacological studies. A number of individuals are described which show clear evoked responses yet a near absence of gamma oscillations and vice versa suggesting dissociations between the generative mechanisms of these responses. Our results also demonstrate that gamma frequency tends to decline with age and is positively correlated with the thickness of the pericalcarine cortex.

Aberrant high-frequency desynchronization of cerebellar cortices in early-onset psychosis

Sensorimotor integration deficits are routinely observed in both schizophreniform and mood-disordered psychoses. Neurobiological theories of schizophrenia and related psychoses have proposed that aberrations in large-scale cortico-thalamic-cerebellar-thalamic-cortical loops may underlie integration abnormalities, and that such dysfunctional connectivity may be central to the pathophysiology. In this study, we utilized a basic mechanoreception task to probe cortical-cerebellar circuitry in early-onset psychosis. Ten adolescents with psychosis and 10 controls completed unilateral tactile stimulation of the right and left index finger, as whole-head magnetoencephalography (MEG) data were acquired. MEG data were imaged in the frequency domain, using spatial filtering, and the resulting event-related synchronizations and desynchronizations (ERS/ERD) were subjected to voxel-wise analyses of group and task effects using statistical parametric mapping. Our results indicated bilateral ERD activation of cerebellar regions and postcentral gyri in both groups during stimulation of either hand. Interestingly, during left finger stimulations, adolescents with psychosis exhibited greater alpha and gamma ERD activity in right cerebellar cortices relative to controls. Subjects with psychosis also showed greater ERD in bilateral cerebellum and the right postcentral gyrus during right finger stimulation, and these differences were statistically stronger for higher frequency bins. Lastly, controls exhibited greater alpha ERS of the right postcentral gyrus during right finger stimulation. These findings provide new data on the neurodevelopmental trajectory of basic mechanoreception in adolescents, and also indicate aberrant cerebellar functioning in early-onset psychoses, especially in the right cerebellum, which may be the crucial dysfunctional node in cortico-thalamic-cerebellar-thalamic-cortical circuits.

Neural responses to rigidly moving faces displaying shifts in social attention investigated with fMRI and MEG

A widely adopted neural model of face perception (Haxby, Hoffman, & Gobbini, 2000) proposes that the posterior superior temporal sulcus (STS) represents the changeable features of a face, while the face-responsive fusiform gyrus (FFA) encodes invariant aspects of facial structure. 'Changeable features' of a face can include rigid and non-rigid movements. The current study investigated neural responses to rigid, moving faces displaying shifts in social attention. Both functional magnetic resonance imaging (fMRI) and magnetoencephalography (MEG) were used to investigate neural responses elicited when participants viewed video clips in which actors made a rigid shift of attention, signalled congruently from both the eyes and head. These responses were compared to those elicited by viewing static faces displaying stationary social attention information or a scrambled video displaying directional motion. Both the fMRI and MEG analyses demonstrated heightened responses along the STS to turning heads compared to static faces or scrambled movement conditions. The FFA responded to both turning heads and static faces, showing only a slight increase in response to the dynamic stimuli. These results establish the applicability of the Haxby model to the perception of rigid face motions expressing changes in social attention direction. Furthermore, the MEG beamforming analyses found an STS response in an upper frequency band (30-80 Hz) which peaked in the right anterior region. These findings, derived from two complementary neuroimaging techniques, clarify the contribution of the STS during the encoding of rigid facial action patterns of social attention, emphasising the role of anterior sulcal regions alongside previously observed posterior areas.

Spatial detection of multiple movement intentions from SAM-filtered single-trial MEG signals

OBJECTIVE

To test whether human intentions to sustain or cease movements in right and left hands can be decoded reliably from spatially filtered single-trial magnetoencephalographic (MEG) signals for motor execution and motor imagery.

METHODS

Seven healthy volunteers, naïve to BCI technology, participated in this study. Signals were recorded from 275-channel MEG, and synthetic aperture magnetometry (SAM) was employed as the spatial filter. The four-class classification was performed offline. Genetic algorithm based Mahalanobis linear distance (GA-MLD) and direct-decision tree classifier (DTC) techniques were adopted for the classification through 10-fold cross-validation.

RESULTS

Through SAM imaging, strong and distinct event-related desynchronization (ERD) associated with sustaining, and event-related synchronization (ERS) patterns associated with ceasing of right and left hand movements were observed in the beta band (15-30Hz) on the contralateral hemispheres for motor execution and motor imagery sessions. Virtual channels were selected from these areas of high activity for the corresponding events as per the paradigm of the study. Through a statistical comparison between SAM-filtered virtual channels from single-trial MEG signals and basic MEG sensors, it was found that SAM-filtered virtual channels significantly increased the classification accuracy for motor execution (GA-MLD: 96.51+/-2.43%) as well as motor imagery sessions (GA-MLD: 89.69+/-3.34%).

CONCLUSION

Multiple movement intentions can be reliably detected from SAM-based spatially filtered single-trial MEG signals.

SIGNIFICANCE

MEG signals associated with natural motor behavior may be utilized for a reliable high-performance brain-computer interface (BCI) and may reduce long-term training compared with conventional BCI methods using rhythm control.

Functional decoupling of BOLD and gamma-band amplitudes in human primary visual cortex

Although functional magnetic resonance imaging is an important tool for measuring brain activity, the hemodynamic blood oxygenation level dependent (BOLD) response is only an indirect measure of neuronal activity. Converging evidence obtained from simultaneous recording of hemodynamic and electrical measures suggest that the best correlate of the BOLD response in primary visual cortex is gamma-band oscillations ( approximately 40 Hz). Here, we examined the coupling between BOLD and gamma-band amplitudes measured with magntoencephalography (MEG) in human primary visual cortex in 10 participants. In Experiment A, participants were exposed to grating stimuli at two contrast levels and two spatial frequencies and in Experiment B square and sine wave stimuli at two spatial frequencies. The amplitudes of both gamma-band oscillations and BOLD showed tuning with stimulus contrast and stimulus type; however, gamma-band oscillations showed a 300% increase across two spatial frequencies, whereas BOLD exhibited no change. This functional decoupling demonstrates that increased amplitude of gamma-band oscillations as measured with MEG is not sufficient to drive the subsequent BOLD response.

Neuroimaging in Human Amblyopia

Functional magnetic resonance imaging (fMRI), positron emission tomography (PET) and magnetoencephalography (MEG) have been the principal neuroimaging tools used to assess the site and nature of cortical deficits in human amblyopia. A review of this growing body of work is presented here with particular reference to various controversial issues, including whether or not the primary visual cortex is dysfunctional, the involvement of higher-order visual areas, neural differences between strabismic and anisometropic amblyopes, and the effects of modern-day drug treatments. We also present our own recent MEG work in which we used the analysis technique of synthetic aperture magnetometry (SAM) to examine the effects of strabismic amblyopia on cortical function. Our results provide evidence that the neuronal assembly associated with form perception in the extrastriate cortex may be dysfunctional in amblyopia, and that the nature of this dysfunction may relate to a change in the normal temporal pattern of neuronal discharges. Based on these results and existing literature, we conclude that a number of cortical areas show reduced levels of activation in amblyopia, including primary and secondary visual areas and regions within the parieto-occipital cortex and ventral temporal cortex.

Spectral properties of induced and evoked gamma oscillations in human early visual cortex to moving and stationary stimuli

In two experiments, magnetoencephalography (MEG) was used to investigate the effects of motion on gamma oscillations in human early visual cortex. When presented centrally, but not peripherally, stationary and moving gratings elicited several evoked and induced response components in early visual cortex. Time-frequency analysis revealed two nonphase locked gamma power increases-an initial, rapidly adapting response and one sustained throughout stimulus presentation and varying in frequency across observers from 28 to 64 Hz. Stimulus motion raised the sustained gamma oscillation frequency by a mean of approximately 10 Hz. The largest motion-induced frequency increases were in those observers with the lowest gamma response frequencies for stationary stimuli, suggesting a possible saturation mechanism. Moderate gamma amplitude increases to moving versus stationary stimuli were also observed but were not correlated with the magnitude of the frequency increase. At the same site in visual cortex, sustained alpha/beta power reductions and an onset evoked response were observed, but these effects did not change significantly with the presence of motion and did not correlate with the magnitude of gamma power changes. These findings suggest that early visual areas encode moving and stationary percepts via activity at higher and lower gamma frequencies, respectively.

Activation of the left inferior frontal gyrus in the first 200 ms of reading: evidence from magnetoencephalography (MEG)

BACKGROUND

It is well established that the left inferior frontal gyrus plays a key role in the cerebral cortical network that supports reading and visual word recognition. Less clear is when in time this contribution begins. We used magnetoencephalography (MEG), which has both good spatial and excellent temporal resolution, to address this question.

METHODOLOGY/PRINCIPAL FINDINGS

MEG data were recorded during a passive viewing paradigm, chosen to emphasize the stimulus-driven component of the cortical response, in which right-handed participants were presented words, consonant strings, and unfamiliar faces to central vision. Time-frequency analyses showed a left-lateralized inferior frontal gyrus (pars opercularis) response to words between 100-250 ms in the beta frequency band that was significantly stronger than the response to consonant strings or faces. The left inferior frontal gyrus response to words peaked at approximately 130 ms. This response was significantly later in time than the left middle occipital gyrus, which peaked at approximately 115 ms, but not significantly different from the peak response in the left mid fusiform gyrus, which peaked at approximately 140 ms, at a location coincident with the fMRI-defined visual word form area (VWFA). Significant responses were also detected to words in other parts of the reading network, including the anterior middle temporal gyrus, the left posterior middle temporal gyrus, the angular and supramarginal gyri, and the left superior temporal gyrus.

CONCLUSIONS/SIGNIFICANCE

These findings suggest very early interactions between the vision and language domains during visual word recognition, with speech motor areas being activated at the same time as the orthographic word-form is being resolved within the fusiform gyrus. This challenges the conventional view of a temporally serial processing sequence for visual word recognition in which letter forms are initially decoded, interact with their phonological and semantic representations, and only then gain access to a speech code.

Cortical oscillatory activity associated with the perception of illusory and real visual contours

We used magnetoencephalography (MEG) to examine the nature of oscillatory brain rhythms when passively viewing both illusory and real visual contours. Three stimuli were employed: a Kanizsa triangle; a Kanizsa triangle with a real triangular contour superimposed; and a control figure in which the corner elements used to form the Kanizsa triangle were rotated to negate the formation of illusory contours. The MEG data were analysed using synthetic aperture magnetometry (SAM) to enable the spatial localisation of task-related oscillatory power changes within specific frequency bands, and the time-course of activity within given locations-of-interest was determined by calculating time-frequency plots using a Morlet wavelet transform. In contrast to earlier studies, we did not find increases in gamma activity (>30 Hz) to illusory shapes, but instead a decrease in 10-30 Hz activity approximately 200 ms after stimulus presentation. The reduction in oscillatory activity was primarily evident within extrastriate areas, including the lateral occipital complex (LOC). Importantly, this same pattern of results was evident for each stimulus type. Our results further highlight the importance of the LOC and a network of posterior brain regions in processing visual contours, be they illusory or real in nature. The similarity of the results for both real and illusory contours, however, leads us to conclude that the broadband (<30 Hz) decrease in power we observed is more likely to reflect general changes in visual attention than neural computations specific to processing visual contours.

Functional neuroimaging: a brief overview and feasibility for use in chiropractic research

There is a need to further our understanding of the neurophysiological effects of chiropractic spinal manipulation on brain activity as it pertains to both musculoskeletal and non-musculoskeletal complaints. This paper aims to provide a basic overview of the most commonly utilised techniques in the neurosciences for functional imaging the brain (positron emission tomography, single-photon emission computerised tomography, functional magnetic resonance imaging, electroencephalography, and magnetoencephalography), and discuss their applicability in future chiropractic research. Functional neuroimaging modalities are used in a wide range of different research and clinical settings, and are powerful tools in the investigation of neuronal activity in the human brain. There are many potential applications for functional neuroimaging in future chiropractic research, but there are some feasibility issues, mainly pertaining to access and funding. We strongly encourage the use of functional neuroimaging in future investigations of the effects of chiropractic spinal manipulation on brain function.

Increased anterior cingulate cortical activity in response to fearful faces: a neurophysiological biomarker that predicts rapid antidepressant response to ketamine

BACKGROUND

Most patients with major depressive disorder (MDD) experience a period of lengthy trial and error when trying to find optimal antidepressant treatment; identifying biomarkers that could predict response to antidepressant treatment would be of enormous benefit. We tested the hypothesis that pretreatment anterior cingulate cortex (ACC) activity could be a putative biomarker of rapid antidepressant response to ketamine, in line with previous findings that investigated the effects of conventional antidepressants. We also investigated patterns of ACC activity to rapid presentation of fearful faces compared with the normal habituation observed in healthy subjects.

METHODS

We elicited ACC activity in drug-free patients with MDD (n = 11) and healthy control subjects (n = 11) by rapidly presenting fearful faces, a paradigm known to activate rostral regions of the ACC. Spatial-filtering analyses were performed on magnetoencephalographic (MEG) recordings, which offer the temporal precision necessary to estimate ACC activity elicited by the rapid presentation of stimuli. Magnetoencephalographic recordings were obtained only once for both patients and control subjects. Patients were subsequently administered a single ketamine infusion followed by assessment of depressive symptoms 4 hours later.

RESULTS

Although healthy subjects had decreased neuromagnetic activity in the rostral ACC across repeated exposures, patients with MDD showed robust increases in pretreatment ACC activity. Notably, this increase was positively correlated with subsequent rapid antidepressant response to ketamine. Exploratory analyses showed that pretreatment amygdala activity was negatively correlated with change in depressive symptoms.

CONCLUSIONS

Pretreatment rostral ACC activation may be a useful biomarker that identifies a subgroup of patients who will respond favorably to ketamine's antidepressant effects.

Effective electromagnetic noise cancellation with beamformers and synthetic gradiometry in shielded and partly shielded environments

The major challenge of MEG, the inverse problem, is to estimate the very weak primary neuronal currents from the measurements of extracranial magnetic fields. The non-uniqueness of this inverse solution is compounded by the fact that MEG signals contain large environmental and physiological noise that further complicates the problem. In this paper, we evaluate the effectiveness of magnetic noise cancellation by synthetic gradiometers and the beamformer analysis method of synthetic aperture magnetometry (SAM) for source localisation in the presence of large stimulus-generated noise. We demonstrate that activation of primary somatosensory cortex can be accurately identified using SAM despite the presence of significant stimulus-related magnetic interference. This interference was generated by a contact heat evoked potential stimulator (CHEPS), recently developed for thermal pain research, but which to date has not been used in a MEG environment. We also show that in a reduced shielding environment the use of higher order synthetic gradiometry is sufficient to obtain signal-to-noise ratios (SNRs) that allow for accurate localisation of cortical sensory function.

Visual awareness, emotion, and gamma band synchronization

What makes us become aware? A popular hypothesis is that if cortical neurons fire in synchrony at a certain frequency band (gamma), we become aware of what they are representing. We tested this hypothesis adopting brain-imaging techniques with good spatiotemporal resolution and frequency-specific information. Specifically, we examined the degree to which increases in event-related synchronization (ERS) in the gamma band were associated with awareness of a stimulus (its detectability) and/or the emotional content of the stimulus. We observed increases in gamma band ERS within prefrontal-anterior cingulate, visual, parietal, posterior cingulate, and superior temporal cortices to stimuli available to conscious awareness. However, we also observed increases in gamma band ERS within the amygdala, visual, prefrontal, parietal, and posterior cingulate cortices to emotional relative to neutral stimuli, irrespective of their availability to conscious access. This suggests that increased gamma band ERS is related to, but not sufficient for, consciousness.

Magnetoencephalography (MEG)

Magnetoencephalography (MEG) encompasses a family of non-contact, non-invasive techniques for detecting the magnetic field generated by the electrical activity of the brain, for analyzing this MEG signal and for using the results to study brain function. The overall purpose of MEG is to extract estimates of the spatiotemporal patterns of electrical activity in the brain from the measured magnetic field outside the head. The electrical activity in the brain is a manifestation of collective neuronal activity and, to a large extent, the currency of brain function. The estimates of brain activity derived from MEG can therefore be used to study mechanisms and processes that support normal brain function in humans and help us understand why, when and how they fail.

Evoked amygdala responses to negative faces revealed by adaptive MEG beamformers

Adaptive beamformer analyses of magnetoencephalograms (MEG) have shown promise as a method for functional imaging of cortical processes. Although recent evidence is encouraging, it is unclear whether these methods can both localize and reconstruct the time course of activity in subcortical structures such as the amygdala. Fourteen healthy participants (7 women) performed a perceptual matching task of negative emotional faces (angry and fearful) and geometric shapes that was designed for functional magnetic resonance imaging (fMRI) studies to maximize amygdala activation. Neuromagnetic data were collected with a 275-channel whole-head magnetometer, and event-related adaptive beamformer analyses were conducted to estimate broadband evoked responses to faces and shapes across the whole brain in 7 mm steps. Group analyses revealed greater left amygdala activity to faces over shapes, both when face-matching and shape-matching trials were presented in separate blocks and when they were randomly intermixed. This finding was replicated in a second experiment with 7 new participants (3 women). Virtual sensor time series showed clear evoked responses in the left amygdala and left fusiform gyrus in both runs and experiments. We conclude that amygdala activity can be resolved from MEGs with adaptive beamformers with temporal resolution superior to other neuroimaging modalities. This demonstration should encourage the use of MEG for elucidating functional networks mediating fear-related neural phenomena that likely unfold rapidly in time across cortical and subcortical structures.

Human hippocampal and parahippocampal theta during goal-directed spatial navigation predicts performance on a virtual Morris water maze

The hippocampus and parahippocampal cortices exhibit theta oscillations during spatial navigation in animals and humans, and in the former are thought to mediate spatial memory formation. Functional specificity of human hippocampal theta, however, is unclear. Neuromagnetic activity was recorded with a whole-head 275-channel magnetoencephalographic (MEG) system as healthy participants navigated to a hidden platform in a virtual reality Morris water maze. MEG data were analyzed for underlying oscillatory sources in the 4-8 Hz band using a spatial filtering technique (i.e., synthetic aperture magnetometry). Source analyses revealed greater theta activity in the left anterior hippocampus and parahippocampal cortices during goal-directed navigation relative to aimless movements in a sensorimotor control condition. Additional analyses showed that left anterior hippocampal activity was predominantly observed during the first one-half of training, pointing to a role for this region in early learning. Moreover, posterior hippocampal theta was highly correlated with navigation performance, with the former accounting for 76% of the variance of the latter. Our findings suggest human spatial learning is dependent on hippocampal and parahippocampal theta oscillations, extending to humans a significant body of research demonstrating such a pivotal role for hippocampal theta in animal navigation.

Interhemispheric differences of spectral power in expressive language: a MEG study with clinical applications

In the last decade we have seen an exponential growth of functional imaging studies investigating multiple aspects of language processing. These studies have sparked an interest in applying some of the paradigms to various clinically relevant questions, such as the identification of the cortical regions mediating language function in surgical candidates for refractory epilepsy. Here we present data from a group of adult control participants in order to investigate the potential of using frequency specific spectral power changes in MEG activation patterns to establish lateralisation of language function using expressive language tasks. In addition, we report on a paediatric patient whose language function was assessed before and after a left hemisphere amygdalo-hippocampectomy. Our verb generation task produced left hemisphere decreases in beta-band power accompanied by right hemisphere increases in low beta-band power in the majority of the control group, a previously unreported phenomenon. This pattern of spectral power was also found in the patient's post-surgery data, though not her pre-surgery data. Comparison of pre and post-operative results also provided some evidence of reorganisation in language related cortex both inter- and intra-hemispherically following surgery. The differences were not limited to changes in localisation of language specific cortex but also changes in the spectral and temporal profile of frontal brain regions during verb generation. While further investigation is required to establish concordance with invasive measures, our data suggest that the methods described may serve as a reliable lateralisation marker for clinical assessment. Furthermore, our findings highlight the potential utility of MEG for the investigation of cortical language functioning in both healthy development and pathology.

Human brain mapping: hemodynamic response and electrophysiology

In view of the recent advance in functional neuroimaging, the current status of non-invasive techniques applied for human brain mapping was reviewed by integrating two principles: hemodynamic and electrophysiological, from the viewpoint of clinical neurophysiology. The currently available functional neuroimaging techniques based on hemodynamic principles are functional magnetic resonance imaging (fMRI), positron emission tomography (PET) or single-photon emission computed tomography (SPECT), and near-infrared spectroscopy (NIRS). Electrophysiological techniques include electroencephalography (EEG), magnetoencephalography (MEG), and transcranial magnetic stimulation (TMS). As for the coupling between hemodynamic response and neuronal activity (neurovascular coupling), experimental studies suggest that the hemodynamic response is significantly correlated to neuronal activity, especially local field potential (synaptic activity) rather than spiking activity, within a certain range. The hemodynamic response tends to be more widespread in space and lasts longer in time as compared with the neuronal activity. Since each technique has its own characteristic features especially in terms of spatial and temporal resolution, it is important to adopt the most appropriate technique for solving each specific question, and it is useful to combine two techniques either simultaneously or in separate sessions. As for the multi-modal approach, the combined use of EEG and MEG, EEG and PET, or EEG and fMRI is applied for the simultaneous studies, and for the separate use of two different techniques, the information obtained from fMRI is used for estimating the generator source from EEG or MEG data (fMRI-constrained source estimation). Functional connectivity among different brain areas can be studied by using a single technique such as the EEG coherence or the correlation analysis of fMRI or PET data, or by combining the stimulation technique such as TMS with neuroimaging. Further advance of each technology and improvement in the analysis method will promote the understanding of precise functional specialization and inter-areal coupling, and will contribute to the increased efficacy of rapidly developing physiological treatments of neurological and psychiatric disorders.

Population-level inferences for distributed MEG source localization under multiple constraints: Application to face-evoked fields

We address some key issues entailed by population inference about responses evoked in distributed brain systems using magnetoencephalography (MEG). In particular, we look at model selection issues at the within-subject level and feature selection issues at the between-subject level, using responses evoked by intact and scrambled faces around 170 ms (M170). We compared the face validity of subject-specific forward models and their summary statistics in terms of how estimated responses reproduced over subjects. At the within-subject level, we focused on the use of multiple constraints, or priors, for inverting distributed source models. We used restricted maximum likelihood (ReML) estimates of prior covariance components (in both sensor and source space) and show that their relative importance is conserved over subjects. At the between-subject level, we used standard anatomical normalization methods to create posterior probability maps that furnish inference about regionally specific population responses. We used these to compare different summary statistics, namely; (i) whether to test for differences between condition-specific source estimates, or whether to test the source estimate of differences between conditions, and (ii) whether to accommodate differences in source orientation by using signed or unsigned (absolute) estimates of source activity.

Optimising experimental design for MEG beamformer imaging

In recent years, the use of beamformers for source localisation has significantly improved the spatial accuracy of magnetoencephalography. In this paper, we examine techniques by which to optimise experimental design, and ensure that the application of beamformers yields accurate results. We show that variation in the experimental duration, or variation in the bandwidth of a signal of interest, can significantly affect the accuracy of a beamformer reconstruction of source power. Specifically, power will usually be underestimated if covariance windows are made too short, or bandwidths too narrow. The accuracy of spatial localisation may also be reduced. We conclude that for optimum accuracy, experimenters should aim to collect as much data as possible, and use a bandwidth spanning the entire frequency distribution of the signal of interest. This minimises distortion to reconstructed source images, time courses and power estimation. In the case where experimental duration is short, and small covariance windows are therefore used, we show that accurate power estimation can be achieved by matrix regularisation. However, large amounts of regularisation cause a loss in the spatial resolution of the MEG beamformer, hence regularisation should be used carefully, particularly if multiple sources in close proximity are expected.

The temporal sequence of evoked and induced cortical responses to implied-motion processing in human motion area V5/MT+

Previous functional magnetic resonance imaging (fMRI) studies have demonstrated that the human visual motion area V5/MT+ is differentially activated by stimuli in which the presence of motion is implied by the content of static photographs, compared with similar static scenes in which no motion is implied. Here, using a group magnetoencephalography study, we confirm the role of V5/MT+ in the perception of implied motion (IM) by the measurement and localization of task-related evoked and induced oscillatory responses, and demonstrate the temporal sequence of these responses. Within the lateral occipital complex, including V5/MT+, statistically significant differential oscillatory responses to IM and implied-static (IS) stimuli were only found in the beta band (15-20 Hz). An early, evoked, beta power increase (IM>IS) occurred at about 150 ms, whilst a power decrease (IM<IS) occurred at approximately 700 ms. We additionally show that the difference between the early evoked beta responses for the IM and IS conditions only occurs when participants are able to ascertain whether an IM or IS photograph is expected. We therefore hypothesize that the enhanced early response represents covert priming of V5/MT+ for the arrival of a motion stimuli, whilst the late-induced beta power increase for IS stimuli indicates suppression of activity within V5/MT+.

Anticipation of affective images and event-related desynchronization (ERD) of alpha activity: An MEG study

We investigated the event-related power decrease (event-related desynchronization: ERD) of the alpha bands associated with the anticipation of affective images. Participants (n=19) were presented with emotionally positive or negative images under different anticipatory conditions, and their brain responses were recorded using magnetoencephalography (MEG). In the Affective Cue conditions, the cue stimulus indicated the emotional valence (positive or negative) of the image. In the Null Cue condition, the cue stimulus did not include any information about the valence of the image, and in the No Cue condition, the affective image was presented without a preceding cue. The cues in the affective and null conditions were followed by emotional images. During the anticipation period for the affective image, the alpha ERD preceding an anticipated negative image was larger than that preceding an anticipated positive image; this effect had an occipital dominance. Furthermore, during the anticipation period, the lower-2-alpha ERD of the right frontal area showed the same result. These results demonstrate that anticipation of negative stimuli induced alpha ERD in both the visual and the right frontal cortex, indicating that top-down modulation may be provided by the right frontal cortex to the visual cortex.

Robust source analysis of oscillatory motor cortex activity with inherently variable phase delay

This study evaluated quantitatively the synchronization between the magnetoencephalography (MEG) and electromyography (EMG) signals and developed a novel method for the determination of the synchronization in order to increase the reliability of the source analysis of the oscillatory motor cortex activity. The new method is based on our observation that there are large variances in the time lag due to relatively low muscle-cortex synchronization which reduces the signal-to-noise ratio of the MEG signal when averaged in direct synchrony with the rectified EMG peaks. To improve the localization of the motor cortex activity, time-frequency analysis was performed for each epoch coinciding with an EMG peak to reject the weak oscillatory activity and artifacts. In addition, the MEG signals were shifted to maximize the coherence between MEG and rectified EMG by determining for each accepted epoch the time lag resulting in a maximum cross-correlation. An experiment was carried out using 30 subjects in order to determine the applicability of this method to a real situation. The synchronization and the results of the corresponding source analysis based on the novel method were compared with the data obtained using the non-phase-shift method and Hilbert approach detecting EMG phase. The results showed that the synchronization was significantly enhanced and the signal-to-noise ratio of the MEG signals improved, and that the localized dipoles of all subjects were well clustered at the motor cortex. This method, based on shifting the MEG epochs according to the simultaneously measured time lag, considerably improves the performance of the averaging and localization of the rhythmic activity of the motor cortex.

Neural responses to auditory stimulus deviance under threat of electric shock revealed by spatially-filtered magnetoencephalography

Stimulus novelty or deviance may be especially salient in anxiety-related states due to sensitization to environmental change, a key symptom of anxiety disorders such as posttraumatic stress disorder (PTSD). We aimed to identify human brain regions that show potentiated responses to stimulus deviance during anticipatory anxiety. Twenty participants (14 men) were presented a passive oddball auditory task in which they were exposed to uniform auditory stimulation of tones with occasional deviations in tone frequency, a procedure that elicits the mismatch negativity (MMN) and its magnetic counterpart (MMNm). These stimuli were presented during threat periods when participants anticipated unpleasant electric shocks, and safe periods when no shocks were anticipated. Neuromagnetic data were collected with a 275-channel whole-head MEG system and event-related beamformer analyses were conducted to estimate source power across the brain in response to stimulus deviance. Source analyses revealed greater right auditory and inferior parietal activity to stimulus deviance under threat relative to safe conditions, consistent with locations of MMN and MMNm sources identified in other studies. Structures related to evaluation of threat, left amygdala and right insula, also showed increased activity to stimulus deviance under threat. As anxiety level increased across participants, right and left auditory cortical as well as right amygdala activity increased to stimulus deviance. These findings fit with evidence of a potentiated MMN in PTSD relative to healthy controls, and warrant closer evaluation of how these structures might form a functional network mediating sensitization to stimulus deviance during anticipatory anxiety.

Beamformer reconstruction of correlated sources using a modified source model

This paper introduces a lead field formulation for use in beamformer analysis of MEG data. This 'dual source beamformer' is a technique to image two temporally correlated sources using beamformer methodology. We show that while the standard, single source beamformer suppresses the reconstructed power of two spatially separate but temporally correlated sources, the dual source beamformer allows for their accurate reconstruction. The technique is proven to be accurate using simulations. We also show that it can be used to image accurately the auditory steady state response, which is correlated between the left and right auditory cortices. We suggest that this technique represents a useful way of locating correlated sources, particularly if a seed location can be defined a priori for one of the two sources. Such a priori information could be based on previous studies using similar paradigms, or from other functional neuroimaging techniques.

Stimuli of varying spatial scale induce gamma activity with distinct temporal characteristics in human visual cortex

Gamma activity to stationary grating stimuli was studied non-invasively using MEG recordings in humans. Using a spatial filtering technique, we localized gamma activity to primary visual cortex. We tested the hypothesis that spatial frequency properties of visual stimuli may be related to the temporal frequency characteristics of the associated cortical responses. We devised a method to assess temporal frequency differences between stimulus-related responses that typically exhibit complex spectral shapes. We applied this methodology to either single-trial (induced) or time-averaged (evoked) responses in four frequency ranges (0-40, 20-60, 40-80 and 60-100 Hz) and two time windows (either the entire duration of stimulus presentation or the first second following stimulus onset). Our results suggest that stimuli of varying spatial frequency induce responses that exhibit significantly different temporal frequency characteristics. These effects were particularly accentuated for induced responses in the classical gamma frequency band (20-60 Hz) analyzed over the entire duration of stimulus presentation. Strikingly, examining the first second of the responses following stimulus onset resulted in significant loss in stimulus specificity, suggesting that late signal components contain functionally relevant information. These findings advocate a functional role of gamma activity in sensory representation. We suggest that stimulus specific frequency characteristics of MEG signals can be mapped to processes of neuronal synchronization within the framework of coupled dynamical systems.

Semantic and phonological task-set priming and stimulus processing investigated using magnetoencephalography (MEG)

In this study the neural substrates of semantic and phonological task priming and task performance were investigated using single word task-primes. Magnetoencephalography (MEG) data were analysed using Synthetic Aperture Magnetometry (SAM) to determine the spatiotemporal and spectral characteristics of cortical responses. Comparisons were made between the task-prime conditions for evidence of differential effects as a function of the nature of the task being primed, and between the task-prime and the task performance responses for evidence of parallels in activation associated with preparation for and completion of a specific task. Differential priming effects were found. Left middle temporal and inferior frontal voxels showed a statistically significant power decrease associated with the semantic task-prime, and a power increase associated with the phonological task-prime, within beta and gamma frequency bands respectively. Similarities between the task-related differential effects associated with task-prime presentation and those associated with target stimulus presentation were also found. For example, within the semantic task condition, left superior frontal and middle temporal regions showed a significant power decrease within both task-prime and target epochs; within the phonological task condition there were significant parietal and cerebellar power decreases within both types of epoch. In addition there was evidence within the priming epochs of dissociable patterns of activity which could be interpreted as indices of de-activation of task-irrelevant networks. Following a phonological task-prime, significant power increases were observed in those inferior frontal and middle temporal regions in which significant power decreases were associated with semantic task priming and performance.

Generators of the gamma-band activities in response to rare and novel stimuli during the auditory oddball paradigm

In this study, we report the cortical sources of the gamma-band activity emitted during the auditory oddball paradigm using the adaptive beamformer algorithm and non-parametric permutation test and then compare them with those of the P3a and P3b components. The results of the gamma-band current sources revealed that the same gamma-band activities are in process during the rare target and novelty task. In the low (30-55 Hz) gamma-band activity, the common sources were localized in the (inferior) anterior cingulate and adjacent inferior frontal cortex. In the high (65-85 Hz) gamma-band activity, the generator was represented in the superior frontal cortex. On the other hand, the P3a and P3b generators showed widespread distributions including the well-known fronto-parietal network [J. Polich, Theoretical overview of P3a and P3b, in: J. Polish (Ed.), Detection of Change: Event-Related Potential and fMRI Findings, Kluwer Academic Press, Boston, 2003, pp. 83-98]. In conclusion, the same frontal generators of gamma-band activities in the present study may be associated with the functions of attentional control for the binding of consecutive cognitive stages corresponding to earlier P3a and later P3b components, which have distinct source distributions except for some overlaps.

Neural dynamics for facial threat processing as revealed by gamma band synchronization using MEG

Facial threat conveys important information about imminent environmental danger. The rapid detection of this information is critical for survival and social interaction. However, due to technical and methodological difficulties, the spatiotemporal profile for facial threat processing is unknown. By utilizing magnetoencephalography (MEG), a brain-imaging technique with superb temporal resolution and fairly good spatial resolution, Synthetic Aperture Magnetometry (SAM), a recently developed source analysis technique, and a sliding window analysis, we identified the spatiotemporal development of facial threat processing in the gamma frequency band. We also tested the dual-route hypothesis by LeDoux who proposed, based on animal research, that there are two routes to the amygdala: a quick subcortical route and a slower and cortical route. Direct evidence with humans supporting this model has been lacking. Moreover, it has been unclear whether the subcortical route responds specifically to fearful expressions or to threatening expressions in general. We found early event-related synchronizations (ERS) in response to fearful faces in the hypothalamus/thalamus area (10-20 ms) and then the amygdala (20-30 ms). This was even earlier than the ERS response seen to fearful faces in visual cortex (40-50 ms). These data support LeDoux's suggestion of a quick, subcortical thamalo-amygdala route. Moreover, this route was specific for fear expressions; the ERS response in the amygdala to angry expressions had a late onset (150-160 ms). The ERS onset in prefrontal cortex followed that seen within the amygdala (around 160-210 ms). This is consistent with its role in higher-level emotional/cognitive processing.

The spatial distribution and temporal dynamics of brain regions activated during the perception of object and non-object patterns

Both animal and human studies suggest that the efficiency with which we are able to grasp objects is attributable to a repertoire of motor signals derived directly from vision. This is in general agreement with the long-held belief that the automatic generation of motor signals by the perception of objects is based on the actions they afford. In this study, we used magnetoencephalography (MEG) to determine the spatial distribution and temporal dynamics of brain regions activated during passive viewing of object and non-object targets that varied in the extent to which they afforded a grasping action. Synthetic Aperture Magnetometry (SAM) was used to localize task-related oscillatory power changes within specific frequency bands, and the time course of activity within given regions-of-interest was determined by calculating time-frequency plots using a Morlet wavelet transform. Both single subject and group-averaged data on the spatial distribution of brain activity are presented. We show that: (i) significant reductions in 10-25 Hz activity within extrastriate cortex, occipito-temporal cortex, sensori-motor cortex and cerebellum were evident with passive viewing of both objects and non-objects; and (ii) reductions in oscillatory activity within the posterior part of the superior parietal cortex (area Ba7) were only evident with the perception of objects. Assuming that focal reductions in low-frequency oscillations (<30 Hz) reflect areas of heightened neural activity, we conclude that: (i) activity within a network of brain areas, including the sensori-motor cortex, is not critically dependent on stimulus type and may reflect general changes in visual attention; and (ii) the posterior part of the superior parietal cortex, area Ba7, is activated preferentially by objects and may play a role in computations related to grasping.

Activation of Area MT/V5 and the Right Inferior Parietal Cortex during the Discrimination of Transient Direction Changes in Translational Motion

The perception of changes in the direction of objects that translate in space is an important function of our visual system. Here we investigate the brain electrical phenomena underlying such a function by using a combination of magnetoencephalography (MEG) and magnetic resonance imaging. We recorded MEG-evoked responses in 9 healthy human subjects while they discriminated the direction of a transient change in a translationally moving random dot pattern presented either to the right or to the left of a central fixation point. We found that responses reached their maximum in 2 main regions corresponding to motion processing area middle temporal (MT)/V5 contralateral to the stimulated visual field, and to the right inferior parietal lobe (rIPL). The activation latencies were very similar in both regions ( approximately 135 ms) following the direction change onset. Our findings suggest that area MT/V5 provides the strongest sensory signal in response to changes in the direction of translational motion, whereas area rIPL may be involved either in the sensory processing of transient motion signals or in the processing of signals related to orienting of attention.

Effects of memory load on cortical oscillatory activity during auditory pattern working memory

The present magnetoencephalography study investigated memory load-dependent changes in cortical oscillatory activity during a modified auditory version of the Sternberg paradigm. Twelve subjects were presented with test stimulus sets of 1-3 syllables spoken in a natural female voice. After an 800-ms delay period, a probe syllable was presented and subjects had to judge whether the probe had been included in the preceding test set. Compared to a non-memory-control task, memory trials elicited an increase of beta activity over right temporal regions and an increase of alpha activity over right middle prefrontal cortex at the end of the delay phase. Monotonic increases in spectral amplitude as a function of memory load were revealed for the beta band over right temporal sensors and the alpha band over right frontal sensors during the delay period. The results demonstrate the relevance of both beta and alpha oscillations for the memorization of multiple stimuli. The former may be associated with the representation of task-relevant stimulus features, while the latter may reflect the top-down control of these representations.

Disordered connectivity in the autistic brain: challenges for the "new psychophysiology"

In 2002, we published a paper [Brock, J., Brown, C., Boucher, J., Rippon, G., 2002. The temporal binding deficit hypothesis of autism. Development and Psychopathology 142, 209-224] highlighting the parallels between the psychological model of 'central coherence' in information processing [Frith, U., 1989. Autism: Explaining the Enigma. Blackwell, Oxford] and the neuroscience model of neural integration or 'temporal binding'. We proposed that autism is associated with abnormalities of information integration that is caused by a reduction in the connectivity between specialised local neural networks in the brain and possible overconnectivity within the isolated individual neural assemblies. The current paper updates this model, providing a summary of theoretical and empirical advances in research implicating disordered connectivity in autism. This is in the context of changes in the approach to the core psychological deficits in autism, of greater emphasis on 'interactive specialisation' and the resultant stress on early and/or low-level deficits and their cascading effects on the developing brain [Johnson, M.H., Halit, H., Grice, S.J., Karmiloff-Smith, A., 2002. Neuroimaging of typical and atypical development: a perspective from multiple levels of analysis. Development and Psychopathology 14, 521-536]. We also highlight recent developments in the measurement and modelling of connectivity, particularly in the emerging ability to track the temporal dynamics of the brain using electroencephalography (EEG) and magnetoencephalography (MEG) and to investigate the signal characteristics of this activity. This advance could be particularly pertinent in testing an emerging model of effective connectivity based on the balance between excitatory and inhibitory cortical activity [Rubenstein, J.L., Merzenich M.M., 2003. Model of autism: increased ratio of excitation/inhibition in key neural systems. Genes, Brain and Behavior 2, 255-267; Brown, C., Gruber, T., Rippon, G., Brock, J., Boucher, J., 2005. Gamma abnormalities during perception of illusory figures in autism. Cortex 41, 364-376]. Finally, we note that the consequence of this convergence of research developments not only enables a greater understanding of autism but also has implications for prevention and remediation.

Preoperative simulation of intracerebral epileptiform discharges: synthetic aperture magnetometry virtual sensor analysis of interictal magnetoencephalography data

OBJECT

Magnetoencephalography (MEG) has been used for the preoperative localization of epileptic equivalent current dipoles (ECDs) in neocortical epilepsy. Spatial filtering can be applied to MEG data by means of synthetic aperture magnetometry (SAM), and SAM virtual sensor analysis can be used to estimate the strength and temporal course of the epileptic source in the region of interest. To evaluate the clinical usefulness of this approach, the authors compare the results of SAM virtual sensor analysis to the results of ECD analysis, subdural electroencephalography (EEG) findings, and surgical outcomes in pediatric patients with neocortical epilepsy.

METHODS

Ten pediatric patients underwent MEG, invasive subdural EEG, and cortical resection for neocortical epilepsy. The authors compared the morphological characteristics, quantity, location, and distribution of the epileptiform discharges assessed using SAM and ECD analysis, and subdural EEG findings (interictal discharges and ictal onset zones). In nine patients, MEG revealed clustered ECDs. The region exhibiting the maximum percentage (> or = 70%) of spikes/sharp waves on SAM was colocalized to clustered ECDs in seven patients. In six patients, SAM demonstrated focal spikes; in two, diffuse spikes; and in two others, focal rhythmic sharp waves. These epileptiform discharges were similar to those recorded on subdural EEG. In nine patients, concordant regions containing the maximum percentage of spikes/sharp waves were revealed by SAM and subdural EEG data. The region of the maximum percentage of spikes/sharp waves as demonstrated by SAM was colocalized to the ictal onset zone identified by subdural EEG findings in seven patients and partially colocalized in two.

CONCLUSIONS

The SAM virtual sensor analysis revealed morphological characteristics, location, and distribution of epileptiform discharges similar to those shown by subdural EEG recordings. By using SAM it is possible to predict intracerebral interictal epileptiform discharges in the region of interest from noninvasively collected preoperative MEG data. The maximum interictal discharge zone identified by SAM virtual sensors correlated to clustered ECDs and the ictal onset zone on subdural EEG findings. Complementary analyses of ECDs and SAM on three-dimensional MR images can improve delineation of epileptogenic zones and lesions in neocortical epilepsy.

A verifiable solution to the MEG inverse problem

Magnetoencephalography (MEG) is a non-invasive brain imaging technique with the potential for very high temporal and spatial resolution of neuronal activity. The main stumbling block for the technique has been that the estimation of a neuronal current distribution, based on sensor data outside the head, is an inverse problem with an infinity of possible solutions. Many inversion techniques exist, all using different a-priori assumptions in order to reduce the number of possible solutions. Although all techniques can be thoroughly tested in simulation, implicit in the simulations are the experimenter's own assumptions about realistic brain function. To date, the only way to test the validity of inversions based on real MEG data has been through direct surgical validation, or through comparison with invasive primate data. In this work, we constructed a null hypothesis that the reconstruction of neuronal activity contains no information on the distribution of the cortical grey matter. To test this, we repeatedly compared rotated sections of grey matter with a beamformer estimate of neuronal activity to generate a distribution of mutual information values. The significance of the comparison between the un-rotated anatomical information and the electrical estimate was subsequently assessed against this distribution. We found that there was significant (P < 0.05) anatomical information contained in the beamformer images across a number of frequency bands. Based on the limited data presented here, we can say that the assumptions behind the beamformer algorithm are not unreasonable for the visual-motor task investigated.

Independent component analysis of the EEG: is this the way forward for understanding abnormalities of brain-gut signalling?

A combination of electroencephalography and independent component analysis has the potential to contribute towards our understanding of brain‐gut signalling