Event-related functional MRI: Past, present, and future

The past two decades have seen an enormous growth in the field of human brain mapping. Investigators have extensively exploited techniques such as positron emission tomography and MRI to map patterns of brain activity based on changes in cerebral hemodynamics. However, until recently, most studies have investigated equilibrium changes in blood flow measured over time periods upward of 1 min. The advent of high-speed MRI methods, capable of imaging the entire brain with a temporal resolution of a few seconds, allows for brain mapping based on more transient aspects of the hemodynamic response. Today it is now possible to map changes in cerebrovascular parameters essentially in real time, conferring the ability to observe changes in brain state that occur over time periods of seconds. Furthermore, because robust hemodynamic alterations are detectable after neuronal stimuli lasting only a few tens of milliseconds, a new class of task paradigms designed to measure regional responses to single sensory or cognitive events can now be studied. Such “event related” functional MRI should provide for fundamentally new ways to interrogate brain function, and allow for the direct comparison and ultimately integration of data acquired by using more traditional behavioral and electrophysiological methods.

Effects of anesthesia on functional activation of cerebral blood flow and metabolism

Functional brain mapping based on changes in local cerebral blood flow (lCBF) or glucose utilization (lCMRglc) induced by functional activation is generally carried out in animals under anesthesia, usually α-chloralose because of its lesser effects on cardiovascular, respiratory, and reflex functions. Results of studies on the role of nitric oxide (NO) in the mechanism of functional activation of lCBF have differed in unanesthetized and anesthetized animals. NO synthase inhibition markedly attenuates or eliminates the lCBF responses in anesthetized animals but not in unanesthetized animals. The present study examines in conscious rats and rats anesthetized with α-chloralose the effects of vibrissal stimulation on lCMRglc and lCBF in the whisker-to-barrel cortex pathway and on the effects of NO synthase inhibition with NG-nitro-l-arginine methyl ester (l-NAME) on the magnitude of the responses. Anesthesia markedly reduced the lCBF and lCMRglc responses in the ventral posteromedial thalamic nucleus and barrel cortex but not in the spinal and principal trigeminal nuclei. l-NAME did not alter the lCBF responses in any of the structures of the pathway in the unanesthetized rats and also not in the trigeminal nuclei of the anesthetized rats. In the thalamus and sensory cortex of the anesthetized rats, where the lCBF responses to stimulation had already been drastically diminished by the anesthesia, l-NAME treatment resulted in loss of statistically significant activation of lCBF by vibrissal stimulation. These results indicate that NO does not mediate functional activation of lCBF under physiological conditions.

Functional source separation applied to induced visual gamma activity

Objective of this work was to explore the performance of a recently introduced source extraction method, FSS (Functional Source Separation), in recovering induced oscillatory change responses from extra-cephalic magnetoencephalographic (MEG) signals. Unlike algorithms used to solve the inverse problem, FSS does not make any assumption about the underlying biophysical source model; instead, it makes use of task-related features (functional constraints) to estimate source/s of interest. FSS was compared with blind source separation (BSS) approaches such as Principal and Independent Component Analysis, PCA and ICA, which are not subject to any explicit forward solution or functional constraint, but require source uncorrelatedness (PCA), or independence (ICA). A visual MEG experiment with signals recorded from six subjects viewing a set of static horizontal black/white square-wave grating patterns at different spatial frequencies was analyzed. The beamforming technique Synthetic Aperture Magnetometry (SAM) was applied to localize task-related sources; obtained spatial filters were used to automatically select BSS and FSS components in the spatial area of interest. Source spectral properties were investigated by using Morlet-wavelet time-frequency representations and significant task-induced changes were evaluated by means of a resampling technique; the resulting spectral behaviours in the gamma frequency band of interest (20-70 Hz), as well as the spatial frequency-dependent gamma reactivity, were quantified and compared among methods. Among the tested approaches, only FSS was able to estimate the expected sustained gamma activity enhancement in primary visual cortex, throughout the whole duration of the stimulus presentation for all subjects, and to obtain sources comparable to invasively recorded data.

A new approach to neuroimaging with magnetoencephalography

We discuss the application of beamforming techniques to the field of magnetoencephalography (MEG). We argue that beamformers have given us an insight into the dynamics of oscillatory changes across the cortex not explored previously with traditional analysis techniques that rely on averaged evoked responses. We review several experiments that have used beamformers, with special emphasis on those in which the results have been compared to those observed in functional magnetic resonance imaging (fMRI) and on those studying induced phenomena. We suggest that the success of the beamformer technique, despite the assumption that there are no linear interactions between the mesoscopic local field potentials across distinct cortical areas, may tell us something of the balance between functional integration and segregation in the human brain. What is more, MEG beamformer analysis facilitates the study of these complex interactions within cortical networks that are involved in both sensory-motor and cognitive processes. © 2005 Wiley-Liss, Inc.

Accuracy and applications of group MEG studies using cortical source locations estimated from participants' scalp surfaces

We contend that powerful group studies can be conducted using magnetoencephalography (MEG), which can provide useful insights into the approximate distribution of the neural activity detected with MEG without requiring magnetic resonance imaging (MRI) for each participant. Instead, a participant's MRI is approximated with one chosen as a best match on the basis of the scalp surface from a database of available MRIs. Because large inter-individual variability in sulcal and gyral patterns is an inherent source of blurring in studies using grouped functional activity, the additional error introduced by this approximation procedure has little effect on the group results, and offers a sufficiently close approximation to that of the participants to yield a good indication of the true distribution of the grouped neural activity. T1-weighted MRIs of 28 adults were acquired in a variety of MR systems. An artificial functional image was prepared for each person in which eight 5 × 5 × 5 mm regions of brain activation were simulated. Spatial normalisation was applied to each image using transformations calculated using SPM99 with (1) the participant's actual MRI, and (2) the best matched MRI substituted from those of the other 27 participants. The distribution of distances between the locations of points using real and substituted MRIs had a modal value of 6 mm with 90% of cases falling below 12.5 mm. The effects of this -approach on real grouped SAM source imaging of MEG data in a verbal fluency task are also shown. The distribution of MEG activity in the estimated average response is very similar to that produced when using the real MRIs. © 2003 Wiley-Liss, Inc.

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 ? (4–7 Hz), a (7–14 Hz), ß (15–25 Hz), and ? (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.

Topographic studies of scalp potentials evoked by pattern presentation

The waveform and scalp distribution of the visual evoked potentials elicited by stimuli in the foveal and parafoveal regions have been investigated in a group of normal humans using a 16-channel `brain mapping' system. The waveform and topography of the responses to pattern onset and pattern reversal stimulation were investigated, using 4 x 4o full field and 4 x 2o lateral and altitudinal half-field stimuli. The responses were composed of several successive peaks which are in some respects consistent with those demonstrated by other workers using larger field sizes. The differences in the behaviour of these components with respect to the position of the stimulus in the visual field were suggestive of origins in different areas of the visual cortex and/or different visual mechanism. Of particular interest were the major early positive components `P90' and `P95' of the responses to pattern onset and pattern reversal stimulation respectively. More detailed exploration of the behaviour of these major early positive components was carried out using `M-scaled' stimuli selected to activate one square centimetre patches of striate cortex and associated extrastriate re-projections, positioned at different points in the foveal and parafoveal area of the visual field. The inter- and intra-subject variability in amplitude and localisation of the signals elicited by these targets was considered to be a reflection of the individual variations in relationship of visual field projections with the pattern of gyri and fissures on the proximal surface of the occipital lobe. The behaviour of component P90 of the onset response is consistent with a lateral origin in extrastriate visual cortex; that of P95 of the pattern reversal response is consistent in some respects with a striate cortical origin, but in others with a partial origin in extrastriate cortex.

The investigation of the primary response of the flash visual evoked response

The topographical distribution of the early components of the flash visual evoked response (VER) were investigated using a twenty channel brain mapping system. Thirty subjects, ranging in age from 21 to 84 years, had flash VERs recorded using the standard 10-20 electrode system to a balanced non-cephalic reference. The subjects were divided into three age groups: a young group, a middle group and an older group. The P2 component (positive component around 100-120 msec) of the flash VER was recorded consistently over the occipital region throughout the age range, as was a frontal negative component (N120) of about the same latency. Only the young age group had this single negative component on the frontage channels, whilst the middle age group showed an additional negative component at around 75 msec (N75). Neither group had a recordable P1 component (positive component around 60-75 msec) over the occipital region. The older age group showed both P1 and P2 components over the occipital region with the distribution of the P1 component being more widespread anteriorly. The frontal channels showed both the negative N75 and the later N120 components. The frontal negative components were shown not to be related to the electroretinogram or the balanced non-cephalic reference, but were affected by the type of stimulation. Responses recorded to both pattern reversal and onset/offset stimulation did not show the frontal negative potentials seen with flash stimulation. It was shown that the P1 component is more readily recordable in the elderly and is preceded during middle age by the development of a frontal negative component at around the same latency. The changing morphology of the negative activity in the frontal region across the age range suggests that the use of an Fz reference would produce an artificial P1 component in the middle age group and an enhancement of this component in the elderly, as well as enhance the P2 component in all ages.

Hand somatosensory subcortical and cortical sources assessed by functional source separation:an EEG study

We propose a novel electroencephalographic application of a recently developed cerebral source extraction method (Functional Source Separation, FSS), which starts from extracranial signals and adds a functional constraint to the cost function of a basic independent component analysis model without requiring solutions to be independent. Five ad-hoc functional constraints were used to extract the activity reflecting the temporal sequence of sensory information processing along the somatosensory pathway in response to the separate left and right median nerve galvanic stimulation. Constraints required only the maximization of the responsiveness at specific latencies following sensory stimulation, without taking into account that any frequency or spatial information. After source extraction, the reliability of identified FS was assessed based on the position of single dipoles fitted on its retroprojected signals and on a discrepancy measure. The FS positions were consistent with previously reported data (two early subcortical sources localized in the brain stem and thalamus, the three later sources in cortical areas), leaving negligible residual activity at the corresponding latencies. The high-frequency component of the oscillatory activity (HFO) of the extracted component was analyzed. The integrity of the low amplitude HFOs was preserved for each FS. On the basis of our data, we suggest that FSS can be an effective tool to investigate the HFO behavior of the different neuronal pools, recruited at successive times after median nerve galvanic stimulation. As FSs are reconstructed along the entire experimental session, directional and dynamic HFO synchronization phenomena can be studied.

Contradiction in universal and particular reasoning

A wide range of essential reasoning tasks rely on contradiction identification, a cornerstone of human rationality, communication and debate founded on the inversion of the logical operators "Every" and "Some." A high-density electroencephalographic (EEG) study was performed in 11 normal young adults. The cerebral network involved in the identification of contradiction included the orbito-frontal and anterior-cingulate cortices and the temporo-polar cortices. The event-related dynamic of this network showed an early negative deflection lasting 500 ms after sentence presentation. This was followed by a positive deflection lasting 1.5 s, which was different for the two logical operators. A lesser degree of network activation (either in neuron number or their level of phase locking or both) occurred while processing statements with "Some," suggesting that this was a relatively simpler scenario with one example to be figured out, instead of the many examples or the absence of a counterexample searched for while processing statements with "Every." A self-generated reward system seemed to resonate the recruited circuitry when the contradictory task is successfully completed.

A penny for your thoughts! Patterns of fMRI activity reveal the content and the spatial topography of visual mental images

Visual mental imagery is a complex process that may be influenced by the content of mental images. Neuropsychological evidence from patients with hemineglect suggests that in the imagery domain environments and objects may be represented separately and may be selectively affected by brain lesions. In the present study, we used functional magnetic resonance imaging (fMRI) to assess the possibility of neural segregation among mental images depicting parts of an object, of an environment (imagined from a first-person perspective), and of a geographical map, using both a mass univariate and a multivariate approach. Data show that different brain areas are involved in different types of mental images. Imagining an environment relies mainly on regions known to be involved in navigational skills, such as the retrosplenial complex and parahippocampal gyrus, whereas imagining a geographical map mainly requires activation of the left angular gyrus, known to be involved in the representation of categorical relations. Imagining a familiar object mainly requires activation of parietal areas involved in visual space analysis in both the imagery and the perceptual domain. We also found that the pattern of activity in most of these areas specifically codes for the spatial arrangement of the parts of the mental image. Our results clearly demonstrate a functional neural segregation for different contents of mental images and suggest that visuospatial information is coded by different patterns of activity in brain areas involved in visual mental imagery. Hum Brain Mapp 36:945-958, 2015.

Dynamic causal modeling of load-dependent modulation of effective connectivity within the verbal working memory network

Neuroimaging studies have consistently shown that working memory (WM) tasks engage a distributed neural network that primarily includes the dorsolateral prefrontal cortex, the parietal cortex, and the anterior cingulate cortex. The current challenge is to provide a mechanistic account of the changes observed in regional activity. To achieve this, we characterized neuroplastic responses in effective connectivity between these regions at increasing WM loads using dynamic causal modeling of functional magnetic resonance imaging data obtained from healthy individuals during a verbal n-back task. Our data demonstrate that increasing memory load was associated with (a) right-hemisphere dominance, (b) increasing forward (i.e., posterior to anterior) effective connectivity within the WM network, and (c) reduction in individual variability in WM network architecture resulting in the right-hemisphere forward model reaching an exceedance probability of 99% in the most demanding condition. Our results provide direct empirical support that task difficulty, in our case WM load, is a significant moderator of short-term plasticity, complementing existing theories of task-related reduction in variability in neural networks. Hum Brain Mapp, 2013. © 2013 Wiley Periodicals, Inc.

Resting state morphology predicts the effect of theta burst stimulation in false belief reasoning:ventrolateral prefrontal cortex in false belief reasoning

When required to represent a perspective that conflicts with one's own, functional magnetic resonance imaging (fMRI) suggests that the right ventrolateral prefrontal cortex (rvlPFC) supports the inhibition of that conflicting self-perspective. The present task dissociated inhibition of self-perspective from other executive control processes by contrasting belief reasoning-a cognitive state where the presence of conflicting perspectives was manipulated-with a conative desire state wherein no systematic conflict existed. Linear modeling was used to examine the effect of continuous theta burst stimulation (cTBS) to rvlPFC on participants' reaction times in belief and desire reasoning. It was anticipated that cTBS applied to rvlPFC would affect belief but not desire reasoning, by modulating activity in the Ventral Attention System (VAS). We further anticipated that this effect would be mediated by functional connectivity within this network, which was identified using resting state fMRI and an unbiased model-free approach. Simple reaction-time analysis failed to detect an effect of cTBS. However, by additionally modeling individual measures from within the stimulated network, the hypothesized effect of cTBS to belief (but, importantly, not desire) reasoning was demonstrated. Structural morphology within the stimulated region, rvlPFC, and right temporoparietal junction were demonstrated to underlie this effect. These data provide evidence that inconsistencies found with cTBS can be mediated by the composition of the functional network that is being stimulated. We suggest that the common claim that this network constitutes the VAS explains the effect of cTBS to this network on false belief reasoning. Hum Brain Mapp, 2016. © 2016 Wiley Periodicals, Inc.

Identifying spatially overlapping local cortical networks with MEG

Recent modelling studies (Hadjipapas et al. [2009]: Neuroimage 44:1290-1303) have shown that it may be possible to distinguish between different neuronal populations on the basis of their macroscopically measured (EEG/MEG) mean field. We set out to test whether the different orientation columns contributing to a signal at a specific cortical location could be identified based on the measured MEG signal. We used 1.5deg square, static, obliquely oriented grating stimuli to generate sustained gamma oscillations in a focal region of primary visual cortex. We then used multivariate classifier methods to predict the orientation (left or right oblique) of the stimuli based purely on the time-series data from this one location. Both the single trial evoked response (0-300 ms) and induced post-transient power spectra (300-2,300 ms, 20-70 Hz band) due to the different stimuli were classifiable significantly above chance in 11/12 and 10/12 datasets respectively. Interestingly, stimulus-specific information is preserved in the sustained part of the gamma oscillation, long after perception has occurred and all neuronal transients have decayed. Importantly, the classification of this induced oscillation was still possible even when the power spectra were rank-transformed showing that the different underlying networks give rise to different characteristic temporal signatures. © 2009 Wiley-Liss, Inc.

A critical review of both the synthesis approach and the receptor profile of the 8-chloro-1-(2′,4′-dichlorophenyl)-N-piperidin-1-yl-1,4,5,6-tetrahydrobenzo[6,7]cyclohepta[1,2-c]pyrazole-3-carboxamide and analogue derivatives

We thank European Commission (project “PET BRAIN: Mapping the brain with PET radiolabeled cannabinoid CB1 ligands”; FP7-People-2009-IAPP; Grant Agreement N.25142).