Mapping Interictal activity in epilepsy using a hidden Markov model: A magnetoencephalography study.

Epilepsy is a highly heterogeneous neurological disorder with variable etiology, manifestation, and response to treatment. It is imperative that new models of epileptiform brain activity account for this variability, to identify individual needs and allow clinicians to curate personalized care. Here, we use a hidden Markov model (HMM) to create a unique statistical model of interictal brain activity for 10 pediatric patients. We use magnetoencephalography (MEG) data acquired as part of standard clinical care for patients at the Children's Hospital of Philadelphia. These data are routinely analyzed using excess kurtosis mapping (EKM); however, as cases become more complex (extreme multifocal and/or polymorphic activity), they become harder to interpret with EKM. We assessed the performance of the HMM against EKM for three patient groups, with increasingly complicated presentation. The difference in localization of epileptogenic foci for the two methods was 7 ± 2 mm (mean ± SD over all 10 patients); and 94% ± 13% of EKM temporal markers were matched by an HMM state visit. The HMM localizes epileptogenic areas (in agreement with EKM) and provides additional information about the relationship between those areas. A key advantage over current methods is that the HMM is a data-driven model, so the output is tuned to each individual. Finally, the model output is intuitive, allowing a user (clinician) to review the result and manually select the HMM epileptiform state, offering multiple advantages over previous methods and allowing for broader implementation of MEG epileptiform analysis in surgical decision-making for patients with intractable epilepsy.

Upper and lower limb motor axons demonstrate differential excitability and accommodation to strong hyperpolarizing currents during induced hyperthermia.

Length-dependent peripheral neuropathy typically involves the insidious onset of sensory loss in the lower limbs before later progressing proximally. Recent evidence proposes hyperpolarization-activated cyclic nucleotide-gated (HCN) channels as dysfunctional in rodent models of peripheral neuropathy, and therefore differential expression of HCN channels in the lower limbs was hypothesized as a pathophysiological mechanism accounting for the pattern of symptomatology within this study. We studied six healthy participants, using motor axon excitability including strong and long [-70% and -100% hyperpolarizing threshold electrotonus (TEh)] hyperpolarizing currents to preferably study HCN channel function from the median and tibial nerves from high (40%) and low (20%) threshold. This was recorded at normothermia (~32°C) and then repeated during hyperthermia (~40°C) as an artificial hyperpolarizing axon stress. Significant differences between recovery cycle, superexcitability, accommodation to small depolarizing currents, and alterations in late stages of the inward-rectifying currents of strongest (-70% and -100% TEh) currents were observed in the lower limbs during hyperthermia. We demonstrate differences in late IH current flow, which implies higher expression of HCN channel isoforms. The findings also indicate their potential inference in the symptomatology of length-dependent peripheral neuropathies and may be a unique target for minimizing symptomatology and pathogenesis in acquired disease. NEW & NOTEWORTHY This study demonstrates nerve excitability differences between the upper and lower limbs during hyperthermia, an experimentally induced axonal stress. The findings indicate that there is differential expression of slow hyperpolarization-activated cyclic nucleotide-gated (HCN) channel isoforms between the upper and lower limbs, which was demonstrated through strong, long hyperpolarizing currents during hyperthermia. Such mechanisms may underlie postural control but render the lower limbs susceptible to dysfunction in disease states.

Oscillatory beta activity mediates neuroplastic effects of motor cortex stimulation in humans.

Continuous theta burst stimulation (cTBS) is a repetitive transcranial magnetic stimulation protocol that can inhibit human motor cortex (M1) excitability and impair movement for ≤ 1 h. While offering valuable insights into brain function and potential therapeutic benefits, these neuroplastic effects are highly variable between individuals. The source of this variability, and the electrophysiological mechanisms underlying the inhibitory after-effects, are largely unknown. In this regard, oscillatory activity at beta frequency (15-35 Hz) is of particular interest as it is elevated in motor disorders such as Parkinson's disease and modulated during the generation of movements. Here, we used a source-level magnetoencephalography approach to investigate the hypothesis that the presence of neuroplastic effects following cTBS is associated with concurrent changes in oscillatory M1 beta activity. M1 cortices were localized with a synthetic aperture magnetometry beamforming analysis of visually cued index finger movements. Virtual electrode analysis was used to reconstruct the spontaneous and movement-related oscillatory activity in bilateral M1 cortices, before and from 10 to 45 min after cTBS. We demonstrate that 40 s of cTBS applied over left M1 reduced corticospinal excitability in the right index finger of 8/16 participants. In these responder participants only, cTBS increased the power of the spontaneous beta oscillations in stimulated M1 and delayed reaction times in the contralateral index finger. No further changes were observed in the latency or power of movement-related beta oscillations. These data provide insights into the electrophysiological mechanisms underlying cTBS-mediated impairment of motor function and demonstrate the association between spontaneous oscillatory beta activity in M1 and the inhibition of motor function.

The value of Magnetoencephalography to guide electrode implantation in epilepsy.

To investigate if Magnetoencephalography (MEG) can add non-redundant information to guide implantation sites for intracranial recordings (IR). The contribution of MEG to intracranial recording planning was evaluated in 12 consecutive patients assessed pre-surgically with MEG followed by IR. Primary outcome measures were the identification of focal seizure onset in IR and favorable surgical outcome. Outcome measures were compared to those of 12 patients matched for implantation type in whom non-invasive pre-surgical assessment suggested clear hypotheses for implantation (non-MEG group). In the MEG group, non-invasive assessment without MEG was inconclusive, and MEG was then used to further help identify implantation sites. In all MEG patients, at least one virtual MEG electrode generated suitable hypotheses for the location of implantations. No differences in outcome measures were found between non-MEG and MEG groups. Although the MEG group included more complex patients, it showed similar percentage of successful implantations as the non-MEG group. This suggests that MEG can contribute to identify implantation sites where standard methods failed.

A novel binaural pitch elicited by phase-modulated noise: MEG and psychophysical observations.

Binaural pitches are auditory percepts that emerge from combined inputs to the ears but that cannot be heard if the stimulus is presented to either ear alone. Here, we describe a binaural pitch that is not easily accommodated within current models of binaural processing. Convergent magnetoencephalography (MEG) and psychophysical measurements were used to characterize the pitch, heard when band-limited noise had a rapidly changing interaural phase difference. Several interesting features emerged: First, the pitch was perceptually lateralized, in agreement with the lateralization of the evoked changes in MEG spectral power, and its salience depended on dichotic binaural presentation. Second, the frequency of the pure tone that matched the binaural pitch lay within a lower spectral sideband of the phase-modulated noise and followed the frequency of that sideband when the modulation frequency or center frequency and bandwidth of the noise changed. Thus, the binaural pitch depended on the processing of binaural information in that lower sideband.

Primary and secondary somatosensory cortex responses to anticipation and pain: a magnetoencephalography study.

Several brain regions, including the primary and secondary somatosensory cortices (SI and SII, respectively), are functionally active during the pain experience. Both of these regions are thought to be involved in the sensory-discriminative processing of pain and recent evidence suggests that SI in particular may also be involved in more affective processing. In this study we used MEG to investigate the hypothesis that frequency-specific oscillatory activity may be differentially associated with the sensory and affective components of pain. In eight healthy participants (four male), MEG was recorded during a visceral pain experiment comprising baseline, anticipation, pain and post-pain phases. Pain was delivered via intraluminal oesophageal balloon distension (four stimuli at 1 Hz). Significant bilateral but asymmetrical changes in neural activity occurred in the ß-band within SI and SII. In SI, a continuous increase in neural activity occurred during the anticipation phase (20-30 Hz), which continued during the pain phase but at a lower frequency (10-15 Hz). In SII, oscillatory changes only occurred during the pain phase, predominantly in the 20-30 Hz ß band, and were coincident with the stimulus. These data provide novel evidence of functional diversity within SI, indicating a role in attentional and sensory aspects of pain processing. In SII, oscillatory changes were predominantly stimulus-related, indicating a role in encoding the characteristics of the stimulus. We therefore provide objective evidence of functional heterogeneity within SI and functional segregation between SI and SII, and suggest that the temporal and frequency dynamics within cortical regions may offer valuable insights into pain processing.

Magnetoencephalography abnormalities in adult mild traumatic brain injury: A systematic review.

BACKGROUND: The global incidence of traumatic brain injuries is rising, with at least 80% being classified as mild. These mild injuries are not visible on routine clinical imaging. The potential clinical role of a specific imaging biomarker be it diagnostic, prognostic or directing and monitoring progress of personalised treatment and rehabilitation has driven the exploration of several new neuroimaging modalities. This systematic review examined the evidence for magnetoencephalography (MEG) to provide an imaging biomarker in mild traumatic brain injury (mTBI). METHODS: Our review was prospectively registered on PROSPERO: CRD42019151387. We searched EMBASE, MEDLINE, trial registers, PsycINFO, Cochrane Library and conference abstracts and identified 37 papers describing MEG changes in mTBI eligible for inclusion. Since meta-analysis was not possible, based on the heterogeneity of reported outcomes, we provide a narrative synthesis of results. RESULTS: The two most promising MEG biomarkers are excess resting state low frequency power, and widespread connectivity changes in all frequency bands. These may represent biomarkers with potential for diagnostic application, which reflect time sensitive changes, or may be capable of offering clinically relevant prognostic information. In addition, the rich data that MEG produces are well-suited to new methods of machine learning analysis, which is now being actively explored. INTERPRETATION: MEG reveals several promising biomarkers, in the absence of structural abnormalities demonstrable with either computerised tomography or magnetic resonance imaging. This review has not identified sufficient evidence to support routine clinical use of MEG in mTBI currently. However, verifying MEG's potential would help meet an urgent clinical need within civilian, sports and military medicine.

Mild traumatic brain injury impairs the coordination of intrinsic and motor-related neural dynamics.

Mild traumatic brain injury (mTBI) poses a considerable burden on healthcare systems. Whilst most patients recover quickly, a significant number suffer from sequelae that are not accompanied by measurable structural damage. Understanding the neural underpinnings of these debilitating effects and developing a means to detect injury, would address an important unmet clinical need. It could inform interventions and help predict prognosis. Magnetoencephalography (MEG) affords excellent sensitivity in probing neural function and presents significant promise for assessing mTBI, with abnormal neural oscillations being a potential specific biomarker. However, growing evidence suggests that neural dynamics are (at least in part) driven by transient, pan-spectral bursting and in this paper, we employ this model to investigate mTBI. We applied a Hidden Markov Model to MEG data recorded during resting state and a motor task and show that previous findings of diminished intrinsic beta amplitude in individuals with mTBI are largely due to the reduced beta band spectral content of bursts, and that diminished beta connectivity results from a loss in the temporal coincidence of burst states. In a motor task, mTBI results in diminished burst amplitude, altered modulation of burst probability during movement, and a loss in connectivity in motor networks. These results suggest that, mechanistically, mTBI disrupts the structural framework underlying neural synchrony, which impairs network function. Whilst the damage may be too subtle for structural imaging to see, the functional consequences are detectable and persist after injury. Our work shows that mTBI impairs the dynamic coordination of neural network activity and proposes a potent new method for understanding mTBI.

Neuronal network pharmacodynamics of GABAergic modulation in the human cortex determined using pharmaco-magnetoencephalography.

Neuronal network oscillations are a unifying phenomenon in neuroscience research, with comparable measurements across scales and species. Cortical oscillations are of central importance in the characterization of neuronal network function in health and disease and are influential in effective drug development. Whilst animal in vitro and in vivo electrophysiology is able to characterize pharmacologically induced modulations in neuronal activity, present human counterparts have spatial and temporal limitations. Consequently, the potential applications for a human equivalent are extensive. Here, we demonstrate a novel implementation of contemporary neuroimaging methods called pharmaco-magnetoencephalography. This approach determines the spatial profile of neuronal network oscillatory power change across the cortex following drug administration and reconstructs the time course of these modulations at focal regions of interest. As a proof of concept, we characterize the nonspecific GABAergic modulator diazepam, which has a broad range of therapeutic applications. We demonstrate that diazepam variously modulates theta (4-7 Hz), alpha (7-14 Hz), beta (15-25 Hz), and gamma (30-80 Hz) frequency oscillations in specific regions of the cortex, with a pharmacodynamic profile consistent with that of drug uptake. We examine the relevance of these results with regard to the spatial and temporal observations from other modalities and the various therapeutic consequences of diazepam and discuss the potential applications of such an approach in terms of drug development and translational neuroscience.

Mapping language networks and their association with verbal abilities in paediatric epilepsy using MEG and graph analysis.

Recent theoretical models of language have emphasised the importance of integration within distributed networks during language processing. This is particularly relevant to young patients with epilepsy, as the topology of the functional network and its dynamics may be altered by the disease, resulting in reorganisation of functional language networks. Thus, understanding connectivity within the language network in patients with epilepsy could provide valuable insights into healthy and pathological brain function, particularly when combined with clinical correlates. The objective of this study was to investigate interactions within the language network in a paediatric population of epilepsy patients using measures of MEG phase synchronisation and graph-theoretical analysis, and to examine their association with language abilities. Task dependent increases in connectivity were observed in fronto-temporal networks during verb generation across a group of 22 paediatric patients (9 males and 13 females; mean age 14 years). Differences in network connectivity were observed between patients with typical and atypical language representation and between patients with good and poor language abilities. In addition, node centrality in left frontal and temporal regions was significantly associated with language abilities, where patients with good language abilities had significantly higher node centrality within inferior frontal and superior temporal regions of the left hemisphere, compared to patients with poor language abilities. Our study is one of the first to apply task-based measures of MEG network synchronisation in paediatric epilepsy, and we propose that these measures of functional connectivity and node centrality could be used as tools to identify critical regions of the language network prior to epilepsy surgery.

Rogue bioelectrical waves in the brain: the Hurst exponent as a potential measure for presurgical mapping in epilepsy.

OBJECTIVE: Brain electromagnetic activity in patients with epilepsy is characterized by abnormal high-amplitude transient events (spikes) and abnormal patterns of synchronization of brain rhythms that accompany epileptic seizures. With the aim of improving methods for identifying epileptogenic sources in magnetoencephalographic (MEG) recordings of brain data, we applied methods previously used in the study of oceanic 'rogue waves' and other freak events in complex systems. APPROACH: For data from three patients who were awaiting surgical treatment for epilepsy, we used a beamformer source model to produce volumetric maps showing areas with a high proportion of spikes that could be classified as 'rogue waves', and areas with high Hurst exponent (HE). The HE describes the extent to which a system is exhibiting persistent behavior, may predict the likelihood of freak events. These measures were compared with the more standard measure of kurtosis, which has been shown to be a reliable method for localizing interictal spikes. MAIN RESULTS: There was partial concordance between the three different volumetric maps indicating that each measure provides different information about the underlying brain data. The HE, when combined with a simple connectivity analysis based on phase slope index (PSI), was able to identify the probable epileptogenic zone in all three patients, despite very different patterns of abnormal activity. The differences between distributions of high HE and high kurtosis values indicates that while spikes are propagated through cortex from the epileptogenic zone, the persistent dynamical conditions under which the spikes are generated may not be propagated in a similar way. Finally, the patterns of persistent activity, indicating a departure from 'healthy criticality' in brain networks may explain the wide range of social and cognitive impairments that are seen in epilepsy patients. SIGNIFICANCE: The HE is a potentially useful addition to the clinician's battery of measures which may be used convergently to guide surgical intervention.

An evaluation of kurtosis beamforming in magnetoencephalography to localize the epileptogenic zone in drug resistant epilepsy patients.

OBJECTIVE: Kurtosis beamforming is a useful technique for analysing magnetoencephalograpy (MEG) data containing epileptic spikes. However, the implementation varies and few studies measure concordance with subsequently resected areas. We evaluated kurtosis beamforming as a means of localizing spikes in drug-resistant epilepsy patients. METHODS: We retrospectively applied kurtosis beamforming to MEG recordings of 22 epilepsy patients that had previously been analysed using equivalent current dipole (ECD) fitting. Virtual electrodes were placed in the kurtosis volumetric peaks and visually inspected to select a candidate source. The candidate sources were compared to the ECD localizations and resection areas. RESULTS: The kurtosis beamformer produced interpretable localizations in 18/22 patients, of which the candidate source coincided with the resection lobe in 9/13 seizure-free patients and in 3/5 patients with persistent seizures. The sublobar accuracy of the kurtosis beamformer with respect to the resection zone was higher than ECD (56% and 50%, respectively), however, ECD resulted in a higher lobar accuracy (75%, 67%). CONCLUSIONS: Kurtosis beamforming may provide additional value when spikes are not clearly discernible on the sensors and support ECD localizations when dipoles are scattered. SIGNIFICANCE: Kurtosis beamforming should be integrated with existing clinical protocols to assist in localizing the epileptogenic zone.

Efficacy and mechanism of sub-sensory sacral (optimised) neuromodulation in adults with faecal incontinence: study protocol for a randomised controlled trial.

BACKGROUND: Faecal incontinence (FI) is a substantial health problem with a prevalence of approximately 8% in community-dwelling populations. Sacral neuromodulation (SNM) is considered the first-line surgical treatment option in adults with FI in whom conservative therapies have failed. The clinical efficacy of SNM has never been rigorously determined in a trial setting and the underlying mechanism of action remains unclear. METHODS/DESIGN: The design encompasses a multicentre, randomised, double-blind crossover trial and cohort follow-up study. Ninety participants will be randomised to one of two groups (SNM/SHAM or SHAM/SNM) in an allocation ratio of 1:1. The main inclusion criteria will be adults aged 18-75 years meeting Rome III and ICI definitions of FI, who have failed non-surgical treatments to the UK standard, who have a minimum of eight FI episodes in a 4-week screening period, and who are clinically suitable for SNM. The primary objective is to estimate the clinical efficacy of sub-sensory SNM vs. SHAM at 32 weeks based on the primary outcome of frequency of FI episodes using a 4-week paper diary, using mixed Poisson regression analysis on the intention-to-treat principle. The study is powered (0.9) to detect a 30% reduction in frequency of FI episodes between sub-sensory SNM and SHAM stimulation over a 32-week crossover period. Secondary objectives include: measurement of established and new clinical outcomes after 1 year of therapy using new (2017 published) optimised therapy (with standardised SNM-lead placement); validation of new electronic outcome measures (events) and a device to record them, and identification of potential biological effects of SNM on underlying anorectal afferent neuronal pathophysiology (hypothesis: SNM leads to increased frequency of perceived transient anal sphincter relaxations; improved conscious sensation of defaecatory urge and cortical/subcortical changes in afferent responses to anorectal electrical stimulation (main techniques: high-resolution anorectal manometry and magnetoencephalography). DISCUSSION: This trial will determine clinical effect size for sub-sensory chronic electrical stimulation of the sacral innervation. It will provide experimental evidence of modifiable afferent neurophysiology that may aid future patient selection as well as a basic understanding of the pathophysiology of FI. TRIAL REGISTRATION: International Standard Randomised Controlled Trial Number: ISRCTN98760715 . Registered on 15 September 2017.

The auditory evoked-gamma response and its relation with the N1m.

This study explored the patterns of oscillatory activity that underpin the N1m auditory evoked response. Evoked gamma activity is a small and relatively rarely-reported component of the auditory evoked response, and the objective of this work was to determine how this component relates to the larger and more prolonged changes in lower frequency bands. An event-related beamformer analysis of MEG data from monaural click stimulation was used to reconstruct volumetric images and virtual electrode time series. Group analysis of localisations showed that activity in the gamma band originated from a source that was more medial than those for activity in the theta-to-beta band, and virtual-electrode analysis showed that the source of the gamma activity could be statistically dissociated from the lower-frequency response. These findings are in accordance with separate functional roles for the activity in each frequency band, and provide evidence that the oscillatory activity that underpins the auditory evoked response may contain important information about the physiological basis of the macroscopic signals recorded by MEG in response to auditory stimulation.

Localising the auditory N1m with event-related beamformers: localisation accuracy following bilateral and unilateral stimulation.

The auditory evoked N1m-P2m response complex presents a challenging case for MEG source-modelling, because symmetrical, phase-locked activity occurs in the hemispheres both contralateral and ipsilateral to stimulation. Beamformer methods, in particular, can be susceptible to localisation bias and spurious sources under these conditions. This study explored the accuracy and efficiency of event-related beamformer source models for auditory MEG data under typical experimental conditions: monaural and diotic stimulation; and whole-head beamformer analysis compared to a half-head analysis using only sensors from the hemisphere contralateral to stimulation. Event-related beamformer localisations were also compared with more traditional single-dipole models. At the group level, the event-related beamformer performed equally well as the single-dipole models in terms of accuracy for both the N1m and the P2m, and in terms of efficiency (number of successful source models) for the N1m. The results yielded by the half-head analysis did not differ significantly from those produced by the traditional whole-head analysis. Any localisation bias caused by the presence of correlated sources is minimal in the context of the inter-individual variability in source localisations. In conclusion, event-related beamformers provide a useful alternative to equivalent-current dipole models in localisation of auditory evoked responses.

Distinct contrast response functions in striate and extra-striate regions of visual cortex revealed with magnetoencephalography (MEG).

OBJECTIVE: To spatially and temporally characterise the cortical contrast response function to pattern onset stimuli in humans. METHODS: Magnetoencephalography (MEG) was used to investigate the human cortical contrast response function to pattern onset stimuli with high temporal and spatial resolution. A beamformer source reconstruction approach was used to spatially localise and identify the time courses of activity at various visual cortical loci. RESULTS: Consistent with the findings of previous studies, MEG beamformer analysis revealed two simultaneous generators of the pattern onset evoked response. These generators arose from anatomically discrete locations in striate and extra-striate visual cortex. Furthermore, these loci demonstrated notably distinct contrast response functions, with striate cortex increasing approximately linearly with contrast, whilst extra-striate visual cortex followed a saturating function. CONCLUSIONS: The generators that underlie the pattern onset visual evoked response arise from two distinct regions in striate and extra-striate visual cortex. SIGNIFICANCE: The spatially, temporally and functionally distinct mechanisms of contrast processing within the visual cortex may account for the disparate results observed across earlier studies and assist in elucidating causal mechanisms of aberrant contrast processing in neurological disorders.

Imaging the dynamics of the auditory steady-state evoked response.

This study used magnetoencephalography (MEG) to examine the dynamic patterns of neural activity underlying the auditory steady-state response. We examined the continuous time-series of responses to a 32-Hz amplitude modulation. Fluctuations in the amplitude of the evoked response were found to be mediated by non-linear interactions with oscillatory processes both at the same source, in the alpha and beta frequency bands, and in the opposite hemisphere.

Psychophysiological responses to pain identify reproducible human clusters.

Pain is a ubiquitous yet highly variable experience. The psychophysiological and genetic factors responsible for this variability remain unresolved. We hypothesised the existence of distinct human pain clusters (PCs) composed of distinct psychophysiological and genetic profiles coupled with differences in the perception and the brain processing of pain. We studied 120 healthy subjects in whom the baseline personality and anxiety traits and the serotonin transporter-linked polymorphic region (5-HTTLPR) genotype were measured. Real-time autonomic nervous system parameters and serum cortisol were measured at baseline and after standardised visceral and somatic pain stimuli. Brain processing reactions to visceral pain were studied in 29 subjects using functional magnetic resonance imaging (fMRI). The reproducibility of the psychophysiological responses to pain was assessed at year. In group analysis, visceral and somatic pain caused an expected increase in sympathetic and cortisol responses and activated the pain matrix according to fMRI studies. However, using cluster analysis, we found 2 reproducible PCs: at baseline, PC1 had higher neuroticism/anxiety scores (P ≤ 0.01); greater sympathetic tone (P<0.05); and higher cortisol levels (P ≤ 0.001). During pain, less stimulus was tolerated (P ≤ 0.01), and there was an increase in parasympathetic tone (P ≤ 0.05). The 5-HTTLPR short allele was over-represented (P ≤ 0.005). PC2 had the converse profile at baseline and during pain. Brain activity differed (P ≤ 0.001); greater activity occurred in the left frontal cortex in PC1, whereas PC2 showed greater activity in the right medial/frontal cortex and right anterior insula. In health, 2 distinct reproducible PCs exist in humans. In the future, PC characterization may help to identify subjects at risk for developing chronic pain and may reduce variability in brain imaging studies.

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.