The development of constrained source localisation algorithms for human brain imaging

This work sets out to evaluate the potential benefits and pit-falls in using a priori information to help solve the Magnetoencephalographic (MEG) inverse problem. In chapter one the forward problem in MEG is introduced, together with a scheme that demonstrates how a priori information can be incorporated into the inverse problem. Chapter two contains a literature review of techniques currently used to solve the inverse problem. Emphasis is put on the kind of a priori information that is used by each of these techniques and the ease with which additional constraints can be applied. The formalism of the FOCUSS algorithm is shown to allow for the incorporation of a priori information in an insightful and straightforward manner. In chapter three it is described how anatomical constraints, in the form of a realistically shaped source space, can be extracted from a subject’s Magnetic Resonance Image (MRI). The use of such constraints relies on accurate co-registration of the MEG and MRI co-ordinate systems. Variations of the two main co-registration approaches, based on fiducial markers or on surface matching, are described and the accuracy and robustness of a surface matching algorithm is evaluated. Figures of merit introduced in chapter four are shown to given insight into the limitations of a typical measurement set-up and potential value of a priori information. It is shown in chapter five that constrained dipole fitting and FOCUSS outperform unconstrained dipole fitting when data with low SNR is used. However, the effect of errors in the constraints can reduce this advantage. Finally, it is demonstrated in chapter six that the results of different localisation techniques give corroborative evidence about the location and activation sequence of the human visual cortical areas underlying the first 125ms of the visual magnetic evoked response recorded with a whole head neuromagnetometer

Combining temporal and spectral information with spatial mapping to identify differences between phonological and semantic networks: a Magnetoencephalographic (MEG) approach

Early, lesion-based models of language processing suggested that semantic and phonological processes are associated with distinct temporal and parietal regions respectively, with frontal areas more indirectly involved. Contemporary spatial brain mapping techniques have not supported such clear-cut segregation, with strong evidence of activation in left temporal areas by both processes and disputed evidence of involvement of frontal areas in both processes. We suggest that combining spatial information with temporal and spectral data may allow a closer scrutiny of the differential involvement of closely overlapping cortical areas in language processing. Using beamforming techniques to analyze magnetoencephalography data, we localized the neuronal substrates underlying primed responses to nouns requiring either phonological or semantic processing, and examined the associated measures of time and frequency in those areas where activation was common to both tasks. Power changes in the beta (14–30 Hz) and gamma (30–50 Hz) frequency bands were analyzed in pre-selected time windows of 350–550 and 500–700 ms In left temporal regions, both tasks elicited power changes in the same time window (350–550 ms), but with different spectral characteristics, low beta (14–20 Hz) for the phonological task and high beta (20–30 Hz) for the semantic task. In frontal areas (BA10), both tasks elicited power changes in the gamma band (30–50 Hz), but in different time windows, 500–700 ms for the phonological task and 350–550 ms for the semantic task. In the left inferior parietal area (BA40), both tasks elicited changes in the 20–30 Hz beta frequency band but in different time windows, 350–550 ms for the phonological task and 500–700 ms for the semantic task. Our findings suggest that, where spatial measures may indicate overlapping areas of involvement, additional beamforming techniques can demonstrate differential activation in time and frequency domains

Emergence of Modular Structure in a Large-Scale Brain Network with Interactions between Dynamics and Connectivity

A network of 32 or 64 connected neural masses, each representing a large population of interacting excitatory and inhibitory neurons and generating an electroencephalography/magnetoencephalography like output signal, was used to demonstrate how an interaction between dynamics and connectivity might explain the emergence of complex network features, in particular modularity. Network evolution was modeled by two processes: (i) synchronization dependent plasticity (SDP) and (ii) growth dependent plasticity (GDP). In the case of SDP, connections between neural masses were strengthened when they were strongly synchronized, and were weakened when they were not. GDP was modeled as a homeostatic process with random, distance dependent outgrowth of new connections between neural masses. GDP alone resulted in stable networks with distance dependent connection strengths, typical small-world features, but no degree correlations and only weak modularity. SDP applied to random networks induced clustering, but no clear modules. Stronger modularity evolved only through an interaction of SDP and GDP, with the number and size of the modules depending on the relative strength of both processes, as well as on the size of the network. Lesioning part of the network, after a stable state was achieved, resulted in a temporary disruption of the network structure. The model gives a possible scenario to explain how modularity can arise in developing brain networks, and makes predictions about the time course of network changes during development and following acute lesions

Induced gamma activity in primary visual cortex is related to luminance and not color contrast: an MEG study

Gamma activity in the visual cortex has been reported in numerous EEG studies of coherent and illusory figures. A dominant theme of many such findings has been that temporal synchronization in the gamma band in response to these identifiable percepts is related to perceptual binding of the common features of the stimulus. In two recent studies using magnetoencephalography (MEG) and the beamformer analysis technique, we have shown that the magnitude of induced gamma activity in visual cortex is dependent upon independent stimulus features such as spatial frequency and contrast. In particular, we showed that induced gamma activity is maximal in response to gratings of 3 cycles per degree (3 cpd) of high luminance contrast. In this work, we set out to examine stimulus contrast further by using isoluminant red/green gratings that possess color but not luminance contrast using the same cohort of subjects. We found no induced gamma activity in V1 or visual cortex in response to the isoluminant gratings in these subjects who had previously shown strong induced gamma activity in V1 for luminance contrast gratings

Retinotopic mapping of the primary visual cortex – a challenge for MEG imaging of the human cortex

Magnetoencephalography (MEG) can be used to reconstruct neuronal activity with high spatial and temporal resolution. However, this reconstruction problem is ill-posed, and requires the use of prior constraints in order to produce a unique solution. At present there are a multitude of inversion algorithms, each employing different assumptions, but one major problem when comparing the accuracy of these different approaches is that often the true underlying electrical state of the brain is unknown. In this study, we explore one paradigm, retinotopic mapping in the primary visual cortex (V1), for which the ground truth is known to a reasonable degree of accuracy, enabling the comparison of MEG source reconstructions with the true electrical state of the brain. Specifically, we attempted to localize, using a beanforming method, the induced responses in the visual cortex generated by a high contrast, retinotopically varying stimulus. Although well described in primate studies, it has been an open question whether the induced gamma power in humans due to high contrast gratings derives from V1 rather than the prestriate cortex (V2). We show that the beanformer source estimate in the gamma and theta bands does vary in a manner consistent with the known retinotopy of V1. However, these peak locations, although retinotopically organized, did not accurately localize to the cortical surface. We considered possible causes for this discrepancy and suggest that improved MEG/magnetic resonance imaging co-registration and the use of more accurate source models that take into account the spatial extent and shape of the active cortex may, in future, improve the accuracy of the source reconstructions

Can you tell your clunis from your cubitus? A benchmark for functional imaging

Advances in functional brain imaging have allowed the development of new investigative techniques with clinical application—ranging from presurgical mapping of eloquent cortex to identifying cortical regions involved in religious experiences. Similarly a variety of methods are available to referring physicians, ranging from metabolic measures such as functional magnetic resonance imaging and positron emission tomography to measurements based on electrical activity such as electroencephalography and magnetoencephalography. However, there are no universal benchmarks by which to judge between these methods. In this study we attempt to develop a standard for functional localisation, based on the known functional organisation of somatosensory cortex. Studies have shown spatially distinct sites of brain activity in response to stimulation of various body parts. Generally these studies have focused on areas with large cortical representations, such as the index finger and face. We tested the limits of magnetoencephalography source localisation by stimulation of body parts, namely the clunis and the cubitus, that map to proximal and relatively poorly represented regions of somatosensory cortex

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

Does function fit structure? A ground truth for non-invasive neuroimaging.

There are now a number of non-invasive methods to image human brain function in-vivo. However, the accuracy of these images remains unknown and can currently only be estimated through the use of invasive recordings to generate a functional ground truth. Neuronal activity follows grey matter structure and accurate estimates of neuronal activity will have stronger support from accurate generative models of anatomy. Here we introduce a general framework that, for the first time, enables the spatial distortion of a functional brain image to be estimated empirically. We use a spherical harmonic decomposition to modulate each cortical hemisphere from its original form towards progressively simpler structures, ending in an ellipsoid. Functional estimates that are not supported by the simpler cortical structures have less inherent spatial distortion. This method allows us to compare directly between magnetoencephalography (MEG) source reconstructions based upon different assumption sets without recourse to functional ground truth

Localized energy for wave equations with degenerate trapping

Localized energy estimates have become a fundamental tool when studying wave equations in the presence of asymptotically at background geometry. Trapped rays necessitate a loss when compared to the estimate on Minkowski space. A loss of regularity is a common way to incorporate such. When trapping is sufficiently weak, a logarithmic loss of regularity suffices. Here, by studying a warped product manifold introduced by Christianson and Wunsch, we encounter the first explicit example of a situation where an estimate with an algebraic loss of regularity exists and this loss is sharp. Due to the global-in-time nature of the estimate for the wave equation, the situation is more complicated than for the Schr\"{o}dinger equation. An initial estimate with sub-optimal loss is first obtained, where extra care is required due to the low frequency contributions. An improved estimate is then established using energy functionals that are inspired by WKB analysis. Finally, it is shown that the loss cannot be improved by any power by saturating the estimate with a quasimode.Comment: 18 page

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