Multimodal imaging of human brain activity: rational, biophysical aspects and modes of integration

Until relatively recently the vast majority of imaging and electrophysiological studies of human brain activity have relied on single-modality measurements usually correlated with readily observable or experimentally modified behavioural or brain state patterns. Multi-modal imaging is the concept of bringing together observations or measurements from different instruments. We discuss the aims of multi-modal imaging and the ways in which it can be accomplished using representative applications. Given the importance of haemodynamic and electrophysiological signals in current multi-modal imaging applications, we also review some of the basic physiology relevant to understanding their relationship

A method for determining venous contribution to BOLD contrast sensory activation

While BOLD contrast reflects haemodynamic changes within capillaries serving neural tissue, it also has a venous component. Studies that have determined the relation of large blood vessels to the activation map indicate that veins are the source of the largest response, and the most delayed in time. It would be informative if the location of these large veins could be extracted from the properties of the functional responses, since vessels are not visible in BOLD contrast images. The present study describes a method for investigating whether measures taken from the functional response can reliably predict vein location, or at least be useful in down-weighting the venous contribution to the activation response, and illustrates this method using data from one subject. We combined fMRI at 3 Tesla with high-resolution anatomical imaging and MR venography to test whether the intrinsic properties of activation time courses corresponded to tissue type. Measures were taken from a gamma fit to the functional response. Mean magnitude showed a significant effect of tissue type (P veins ≈ grey matter > white matter. Mean delays displayed the same ranking across tissue types (P grey matter. However, measures for all tissue types were distributed across an overlapping range. A logistic regression model correctly discriminated 72% of the veins from grey matter in the absence of independent information of macroscopic vessels (ROC=0.72). Whilst tissue classification was not perfect for this subject, weighting the T contrast by the predicted probabilities materially reduced the venous component to the activation map

The spatio-temporal mapping of epileptic networks: Combination of EEG–fMRI and EEG source imaging

Simultaneous EEG–fMRI acquisitions in patients with epilepsy often reveal distributed patterns of Blood Oxygen Level Dependant (BOLD) change correlated with epileptiform discharges. We investigated if electrical source imaging (ESI) performed on the interictal epileptiform discharges (IED) acquired during fMRI acquisition could be used to study the dynamics of the networks identified by the BOLD effect, thereby avoiding the limitations of combining results from separate recordings. Nine selected patients (13 IED types identified) with focal epilepsy underwent EEG–fMRI. Statistical analysis was performed using SPM5 to create BOLD maps. ESI was performed on the IED recorded during fMRI acquisition using a realistic head model (SMAC) and a distributed linear inverse solution (LAURA). ESI could not be performed in one case. In 10/12 remaining studies, ESI at IED onset (ESIo) was anatomically close to one BOLD cluster. Interestingly, ESIo was closest to the positive BOLD cluster with maximal statistical significance in only 4/12 cases and closest to negative BOLD responses in 4/12 cases. Very small BOLD clusters could also have clinical relevance in some cases. ESI at later time frame (ESIp) showed propagation to remote sources co-localised with other BOLD clusters in half of cases. In concordant cases, the distance between maxima of ESI and the closest EEG–fMRI cluster was less than 33 mm, in agreement with previous studies. We conclude that simultaneous ESI and EEG–fMRI analysis may be able to distinguish areas of BOLD response related to initiation of IED from propagation areas. This combination provides new opportunities for investigating epileptic networks

Structure Learning in Coupled Dynamical Systems and Dynamic Causal Modelling

Identifying a coupled dynamical system out of many plausible candidates, each of which could serve as the underlying generator of some observed measurements, is a profoundly ill posed problem that commonly arises when modelling real world phenomena. In this review, we detail a set of statistical procedures for inferring the structure of nonlinear coupled dynamical systems (structure learning), which has proved useful in neuroscience research. A key focus here is the comparison of competing models of (ie, hypotheses about) network architectures and implicit coupling functions in terms of their Bayesian model evidence. These methods are collectively referred to as dynamical casual modelling (DCM). We focus on a relatively new approach that is proving remarkably useful; namely, Bayesian model reduction (BMR), which enables rapid evaluation and comparison of models that differ in their network architecture. We illustrate the usefulness of these techniques through modelling neurovascular coupling (cellular pathways linking neuronal and vascular systems), whose function is an active focus of research in neurobiology and the imaging of coupled neuronal systems

Imaging haemodynamic changes related to seizures: comparison of EEG-based general linear model, independent component analysis of fMRI and intracranial EEG

Background: Simultaneous EEG-fMRI can reveal haemodynamic changes associated with epileptic activity which may contribute to understanding seizure onset and propagation. Methods: Nine of 83 patients with focal epilepsy undergoing pre-surgical evaluation had seizures during EEG-fMRI and analysed using three approaches, two based on the general linear model (GLM) and one using independent component analysis (ICA): 1. EEGs were divided into up to three phases: early ictal EEG change, clinical seizure onset and late ictal EEG change and convolved with a canonical haemodynamic response function (HRF) (canonical GLM analysis). 2. Seizures lasting three scans or longer were additionally modelled using a Fourier basis set across the entire event (Fourier GLM analysis). 3. Independent component analysis (ICA) was applied to the fMRI data to identify ictal BOLD patterns without EEG. The results were compared with intracranial EEG. Results: The canonical GLM analysis revealed significant BOLD signal changes associated with seizures on EEG in 7/9 patients, concordant with the seizure onset zone in 4/7. The Fourier GLM analysis revealed changes in BOLD signal corresponding with the results of the canonical analysis in two patients. ICA revealed components spatially concordant with the seizure onset zone in all patients (8/9 confirmed by intracranial EEG). Conclusion: Ictal EEG-fMRI visualises plausible seizure related haemodynamic changes. The GLM approach to analysing EEG-fMRI data reveals localised BOLD changes concordant with the ictal onset zone when scalp EEG reflects seizure onset. ICA provides additional information when scalp EEG does not accurately reflect seizures and may give insight into ictal haemodynamics

Simultaneous intracranial EEG and fMRI of interictal epileptic discharges in humans

Simultaneous scalp EEG–fMRI measurements allow the study of epileptic networks and more generally, of the coupling between neuronal activity and haemodynamic changes in the brain. Intracranial EEG (icEEG) has greater sensitivity and spatial specificity than scalp EEG but limited spatial sampling. We performed simultaneous icEEG and functional MRI recordings in epileptic patients to study the haemodynamic correlates of intracranial interictal epileptic discharges (IED). Two patients undergoing icEEG with subdural and depth electrodes as part of the presurgical assessment of their pharmaco-resistant epilepsy participated in the study. They were scanned on a 1.5 T MR scanner following a strict safety protocol. Simultaneous recordings of fMRI and icEEG were obtained at rest. IED were subsequently visually identified on icEEG and their fMRI correlates were mapped using a general linear model (GLM). On scalp EEG–fMRI recordings performed prior to the implantation, no IED were detected. icEEG–fMRI was well tolerated and no adverse health effect was observed. intra-MR icEEG was comparable to that obtained outside the scanner. In both cases, significant haemodynamic changes were revealed in relation to IED, both close to the most active electrode contacts and at distant sites. In one case, results showed an epileptic network including regions that could not be sampled by icEEG, in agreement with findings from magneto-encephalography, offering some explanation for the persistence of seizures after surgery. Hence, icEEG–fMRI allows the study of whole-brain human epileptic networks with unprecedented sensitivity and specificity. This could help improve our understanding of epileptic networks with possible implications for epilepsy surgery