Investigating the functional properties of the somatosensory cortex during experimental visceral pain using magnetoencephalography

Background The somatosensory cortex has been inconsistently activated in pain studies and the functional properties of subregions within this cortical area are poorly understood. To address this we used magnetoencephalography (MEG), a brain imaging technique capable of recording changes in cortical neural activity in real-time, to investigate the functional properties of the somatosensory cortex during different phases of the visceral pain experience. Methods In eight participants (4 male), 151-channel whole cortex MEG was used to detect cortical neural activity during 25 trials lasting 20 seconds each. Each trial comprised four separate periods of 5 seconds in duration. During each of the periods, different visual cues were presented, indicating that period 1=rest, period 2=anticipation, period 3=pain and period 4=post pain. During period 3, participants received painful oesophageal balloon distensions (four at 1 Hz). Regions of cortical activity were identified using Synthetic Aperture Magnetometry (SAM) and by the placement of virtual electrodes in regions of interest within the somatosensory cortex, time-frequency wavelet plots were generated. Results SAM analysis revealed significant activation with the primary (S1) and secondary (S2) somatosensory cortices. The time-frequency wavelet spectrograms showed that activation in S1 increased during the anticipation phase and continued during the presentation of the stimulus. In S2, activation was tightly time and phase-locked to the stimulus within the pain period. Activations in both regions predominantly occurred within the 10–15 Hz and 20–30 Hz frequency bandwidths. Discussion These data are consistent with the role of S1 and S2 in the sensory discriminatory aspects of pain processing. Activation of S1 during anticipation and then pain may be linked to its proposed role in attentional as well as sensory processing. The stimulus-related phasic activity seen in S2 demonstrates that this region predominantly encodes information pertaining to the nature and intensity of the stimulus.

Dissociating the spatio-temporal characteristics of cortical neuronal activity associated with human volitional swallowing in the healthy adult brain

Human swallowing represents a complex highly coordinated sensorimotor function whose functional neuroanatomy remains incompletely understood. Specifically, previous studies have failed to delineate the temporo-spatial sequence of those cerebral loci active during the differing phases of swallowing. We therefore sought to define the temporal characteristics of cortical activity associated with human swallowing behaviour using a novel application of magnetoencephalography (MEG). In healthy volunteers (n = 8, aged 28-45), 151-channel whole cortex MEG was recorded during the conditions of oral water infusion, volitional wet swallowing (5 ml bolus), tongue thrust or rest. Each condition lasted for 5 s and was repeated 20 times. Synthetic aperture magnetometry (SAM) analysis was performed on each active epoch and compared to rest. Temporal sequencing of brain activations utilised time-frequency wavelet plots of regions selected using virtual electrodes. Following SAM analysis, water infusion preferentially activated the caudolateral sensorimotor cortex, whereas during volitional swallowing and tongue movement, the superior sensorimotor cortex was more strongly active. Time-frequency wavelet analysis indicated that sensory input from the tongue simultaneously activated caudolateral sensorimotor and primary gustatory cortex, which appeared to prime the superior sensory and motor cortical areas, involved in the volitional phase of swallowing. Our data support the existence of a temporal synchrony across the whole cortical swallowing network, with sensory input from the tongue being critical. Thus, the ability to non-invasively image this network, with intra-individual and high temporal resolution, provides new insights into the brain processing of human swallowing. © 2004 Elsevier Inc. All rights reserved.

Spatio-temporal imaging of corticaldesynchronization in migraine visual aura:a magnetoencephalography case study

Objective.-To determine cortical oscillatory changes involved in migraine visual aura using magnetoencephalography (MEG). Background.-Visual aura in the form of scintillating scotoma precedes migraine in many cases. The involvement of cortical spreading depression within striate and extra-striate cortical areas is implicated in the generation of the disturbance, but the details of its progression, the effects on cortical oscillations, and the mechanisms of aura generation are unclear. Methods.-We used MEG to directly image changes in cortical oscillatory power during an episode of scintillating scotoma in a patient who experiences aura without subsequent migraine headache. Using the synthetic aperture magnetometry method of MEG source imaging, focal changes in cortical oscillatory power were observed over a 20-minute period and visualized in coregistration with the patient's magnetic resonance image. Results.-Alpha band desynchronization in both the left extra-striate and temporal cortex persisted for the duration of reported visual disturbance, terminating abruptly upon disappearance of scintillations. Gamma frequency desynchronization in the left temporal lobe continued for 8 to 10 minutes following the reported end of aura. Conclusions.-Observations implicate the extra-striate and temporal cortex in migraine visual aura and suggest involvement of alpha desynchronization in generation of phosphenes and gamma desynchronization in sustained inhibition of visual function.

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. Copyright © 2006 Taylor & Francis Group, LLC.

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.

Assessing interactions of linear and nonlinear neuronal sources using MEG beamformers:a proof of concept

Objective: This study aimed to explore methods of assessing interactions between neuronal sources using MEG beamformers. However, beamformer methodology is based on the assumption of no linear long-term source interdependencies [VanVeen BD, vanDrongelen W, Yuchtman M, Suzuki A. Localization of brain electrical activity via linearly constrained minimum variance spatial filtering. IEEE Trans Biomed Eng 1997;44:867-80; Robinson SE, Vrba J. Functional neuroimaging by synthetic aperture magnetometry (SAM). In: Recent advances in Biomagnetism. Sendai: Tohoku University Press; 1999. p. 302-5]. Although such long-term correlations are not efficient and should not be anticipated in a healthy brain [Friston KJ. The labile brain. I. Neuronal transients and nonlinear coupling. Philos Trans R Soc Lond B Biol Sci 2000;355:215-36], transient correlations seem to underlie functional cortical coordination [Singer W. Neuronal synchrony: a versatile code for the definition of relations? Neuron 1999;49-65; Rodriguez E, George N, Lachaux J, Martinerie J, Renault B, Varela F. Perception's shadow: long-distance synchronization of human brain activity. Nature 1999;397:430-3; Bressler SL, Kelso J. Cortical coordination dynamics and cognition. Trends Cogn Sci 2001;5:26-36]. Methods: Two periodic sources were simulated and the effects of transient source correlation on the spatial and temporal performance of the MEG beamformer were examined. Subsequently, the interdependencies of the reconstructed sources were investigated using coherence and phase synchronization analysis based on Mutual Information. Finally, two interacting nonlinear systems served as neuronal sources and their phase interdependencies were studied under realistic measurement conditions. Results: Both the spatial and the temporal beamformer source reconstructions were accurate as long as the transient source correlation did not exceed 30-40 percent of the duration of beamformer analysis. In addition, the interdependencies of periodic sources were preserved by the beamformer and phase synchronization of interacting nonlinear sources could be detected. Conclusions: MEG beamformer methods in conjunction with analysis of source interdependencies could provide accurate spatial and temporal descriptions of interactions between linear and nonlinear neuronal sources. Significance: The proposed methods can be used for the study of interactions between neuronal sources. © 2005 International Federation of Clinical Neurophysiology. Published by Elsevier Ireland Ltd. All rights reserved.

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. © 2006 Elsevier Inc. All rights reserved.

Development and application of magnetoencephalographic methods in the investigation of the human visual system

This thesis is an exploration of the organisation and functioning of the human visual system using the non-invasive functional imaging modality magnetoencephalography (MEG). Chapters one and two provide an introduction to the ‘human visual system and magnetoencephalographic methodologies. These chapters subsequently describe the methods by which MEG can be used to measure neuronal activity from the visual cortex. Chapter three describes the development and implementation of novel analytical tools; including beamforming based analyses, spectrographic movies and an optimisation of group imaging methods. Chapter four focuses on the use of established and contemporary analytical tools in the investigation of visual function. This is initiated with an investigation of visually evoked and induced responses; covering visual evoked potentials (VEPs) and event related synchronisation/desynchronisation (ERS/ERD). Chapter five describes the employment of novel methods in the investigation of cortical contrast response and demonstrates distinct contrast response functions in striate and extra-striate regions of visual cortex. Chapter six use synthetic aperture magnetometry (SAM) to investigate the phenomena of visual cortical gamma oscillations in response to various visual stimuli; concluding that pattern is central to its generation and that it increases in amplitude linearly as a function of stimulus contrast, consistent with results from invasive electrode studies in the macaque monkey. Chapter seven describes the use of driven visual stimuli and tuned SAM methods in a pilot study of retinotopic mapping using MEG; finding that activity in the primary visual cortex can be distinguished in four quadrants and two eccentricities of the visual field. Chapter eight is a novel implementation of the SAM beamforming method in the investigation of a subject with migraine visual aura; the method reveals desynchronisation of the alpha and gamma frequency bands in occipital and temporal regions contralateral to observed visual abnormalities. The final chapter is a summary of main conclusions and suggested further work.

Functional neuroimaging and behavioural studies on global form processing in the human visual system

Magnetoencephalography (MEG), functional magnetic resonance imaging (fMRI) and behavioural experiments were used to investigate the neural processes underlying global form perception in human vision. Behavioural studies using Glass patterns examined sensitivity for detecting radial, rotational and horizontal structure. Neuroimaging experiments using either Glass patterns or arrays of Gabor patches determined the spatio-temporal neural responseto global form. MEG data were analysed using synthetic aperture magnetometry (SAM) to spatially map event-related cortical oscillatory power changes: the temporal sequencing of activity within a discrete cortical area was determined using a Morlet wavelet transform. A case study was conducted to determine the effects of strbismic amblyopia on global form processing: all other observers were normally-sighted. The main findings from normally-sighted observers were: 1) sensitivity to horizontal structure was less than for radial or rotational structure; 2) the neural response to global structure was a reduction in cortical oscillatory power (10-30 Hz) within a network of extrastriate areas, including V4 and V3a; 3) the extend of reduced cortical power was least for horizontal patters; 4) V1 was not identified as a region of peak activity with either MEG or fMRI. The main findings with the strabismic amblyope were: 1) sensitivity for detection of radial, rotational, and horizontal structure was reduced when viewed with the amblyopic- relative to the fellow- eye; 2) cortical power changes within V4 to the presentation of rotational Glass patterns were less when viewed with the amblyopic- compared with the fellow- eye. The main conclusions are: 1) a network of extrastriate cortical areas are involved in the analysis of global form, with the most prominent change in neural activity being a reduction in oscillatory power within the 10-30 Hz band; 2) in strabismic amblyopia, the neuronal assembly associated with form perception in extrastriate cortex may be dysfunctional, the nature of this dysfunction may be a change in the normal temporal pattern of neuronal discharges; 3) MEG, fMRI and behavioural measures support the notion that different neural processes underlie the perception of horizontal as opposed to radial or rotational structure.

A multimodal investigation of dynamic face perception using functional magnetic resonance imaging and magnetoencephalography

Motion is an important aspect of face perception that has been largely neglected to date. Many of the established findings are based on studies that use static facial images, which do not reflect the unique temporal dynamics available from seeing a moving face. In the present thesis a set of naturalistic dynamic facial emotional expressions was purposely created and used to investigate the neural structures involved in the perception of dynamic facial expressions of emotion, with both functional Magnetic Resonance Imaging (fMRI) and Magnetoencephalography (MEG). Through fMRI and connectivity analysis, a dynamic face perception network was identified, which is demonstrated to extend the distributed neural system for face perception (Haxby et al.,2000). Measures of effective connectivity between these regions revealed that dynamic facial stimuli were associated with specific increases in connectivity between early visual regions, such as inferior occipital gyri and superior temporal sulci, along with coupling between superior temporal sulci and amygdalae, as well as with inferior frontal gyri. MEG and Synthetic Aperture Magnetometry (SAM) were used to examine the spatiotemporal profile of neurophysiological activity within this dynamic face perception network. SAM analysis revealed a number of regions showing differential activation to dynamic versus static faces in the distributed face network, characterised by decreases in cortical oscillatory power in the beta band, which were spatially coincident with those regions that were previously identified with fMRI. These findings support the presence of a distributed network of cortical regions that mediate the perception of dynamic facial expressions, with the fMRI data providing information on the spatial co-ordinates paralleled by the MEG data, which indicate the temporal dynamics within this network. This integrated multimodal approach offers both excellent spatial and temporal resolution, thereby providing an opportunity to explore dynamic brain activity and connectivity during face processing.

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 inhibithumanmotor 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. © 2013 the authors.

Intelligent Computations for Flood Monitoring

Floods represent the most devastating natural hazards in the world, affecting more people and causing more property damage than any other natural phenomena. One of the important problems associated with flood monitoring is flood extent extraction from satellite imagery, since it is impractical to acquire the flood area through field observations. This paper presents a method to flood extent extraction from synthetic-aperture radar (SAR) images that is based on intelligent computations. In particular, we apply artificial neural networks, self-organizing Kohonen’s maps (SOMs), for SAR image segmentation and classification. We tested our approach to process data from three different satellite sensors: ERS-2/SAR (during flooding on Tisza river, Ukraine and Hungary, 2001), ENVISAT/ASAR WSM (Wide Swath Mode) and RADARSAT-1 (during flooding on Huaihe river, China, 2007). Obtained results showed the efficiency of our approach.

Evaluation of TerraSAR-X Observations for Wetland InSAR Application

This paper assesses the potential of using spaceborne X-band synthetic aperture radar (SAR) data for monitoring water-level changes over wetlands. Our analysis is based on three sets of TerraSAR-X (TSX) observations acquired over South Florida's Everglades wetlands during an eight-month period in 2008. The first set was acquired in single HH polarization stripmap mode over our northern study area, consisting of managed wetlands and urban environments. The second set was acquired in dual-polarization stripmap mode over the western half of the same area, consisting mostly of managed wetlands. The third set was also acquired with dual-polarization stripmap mode over our southern study area, consisting of natural flow freshand salt-water wetlands in the southern Everglades. The first data set was used for a proof-of-concept study to verify that X-band data can generate coherent interferograms in wetland areas. Interferometric processing of this data set shows a high level of coherence (> 0.35) over both wetland and urban regions, maintaining interferometric phase in all three interferograms spanning 11 days. Surprisingly, phase is maintained over some of the wetlands even for interferograms spanning 33 days. The other two data sets were used to evaluate interferometric coherence of all four polarization modes and to determine dominant scattering mechanism in each wetland environment. Our results show high coherence values (> 0.4) in all polarization modes, with highest values in HH, then VV, and lowest in HV or VH. Interferograms calculated from multipolarization data show very similar fringe patterns regardless of the polarization type, suggesting that the phase information in all polarization data reflects water-level changes in wetlands and that volume scattering may be less important than commonly believed. We also used the two multipolarization data sets to conduct the Pauli decomposition, finding a strong dependence of scattering mechanism on vegetation t- - ype. The high interferometric coherence level of all polarization data suggests that a significant part of the X-band scattered signal interacts with lower sections of the vegetation (trunks and branches), because scattering from wind-affected canopies cannot support such a high coherence level. The high spatial resolution of TSX, combined with its 11-day repeat orbit, makes this X-band sensor surprisingly suitable for wetland interferometric SAR applications.

Hydrologic Dynamics of the Ground-Water-Dependent Sian Ka’an Wetlands, Mexico, Derived from InSAR and SAR Data

The 5,280 km2 Sian Ka’an Biosphere Reserve includes pristine wetlands fed by ground water from the karst aquifer of the Yucatan Peninsula, Mexico. The inflow through underground karst structures is hard to observe making it difficult to understand, quantify, and predict the wetland dynamics. Remotely sensed Synthetic Aperture Radar (SAR) amplitude and phase observations offer new opportunities to obtain information on hydrologic dynamics useful for wetland management. Backscatter amplitude of SAR data can be used to map flooding extent. Interferometric processing of the backscattered SAR phase data (InSAR) produces temporal phase-changes that can be related to relative water level changes in vegetated wetlands. We used 56 RADARSAT-1 SAR acquisitions to calculate 38 interferograms and 13 flooding maps with 24 day and 48 day time intervals covering July 2006 to March 2008. Flooding extent varied between 1,067 km2 and 2,588 km2 during the study period, and main water input was seen to take place in sloughs during October–December. We propose that main water input areas are associated with water-filled faults that transport ground water from the catchment to the wetlands. InSAR and Landsat data revealed local-scale water divides and surface water flow directions within the wetlands.

Remote-Sensing Monitoring of Tide Propagation Through Coastal Wetlands

Tide propagation through coastal wetlands is a complex phenomenon affected by vegetation, channels, and tidal conditions. Generally, tidal flow is studied using stage (water level) observations, which provide good temporal resolution, but they are acquired in limited locations. Here, a remote-sensing technique, wetland InSAR (interferometric synthetic aperture radar), is used to detect tidal flow in vegetated coastal environments over broad spatial scales. The technique is applied to data sets acquired by three radar satellites over the western Everglades in south Florida. Interferometric analysis of the data shows that the greatest water-level changes occur along tidal channels, reflecting a high velocity gradient between fast horizontal flow in the channel and the slow flow propagation through the vegetation. The high-resolution observations indicate that the tidal flushing zone extends 2–3 km on both sides of tidal channels and can extend 3–4 km inland from the end of the channel. The InSAR observations can also serve as quantitative constraints for detailed coastal wetland flow models.