Encoding cortical dynamics in sparse features.

Distributed cortical solutions of magnetoencephalography (MEG) and electroencephalography (EEG) exhibit complex spatial and temporal dynamics. The extraction of patterns of interest and dynamic features from these cortical signals has so far relied on the expertise of investigators. There is a definite need in both clinical and neuroscience research for a method that will extract critical features from high-dimensional neuroimaging data in an automatic fashion. We have previously demonstrated the use of optical flow techniques for evaluating the kinematic properties of motion field projected on non-flat manifolds like in a cortical surface. We have further extended this framework to automatically detect features in the optical flow vector field by using the modified and extended 2-Riemannian Helmholtz-Hodge decomposition (HHD). Here, we applied these mathematical models on simulation and MEG data recorded from a healthy individual during a somatosensory experiment and an epilepsy pediatric patient during sleep. We tested whether our technique can automatically extract salient dynamical features of cortical activity. Simulation results indicated that we can precisely reproduce the simulated cortical dynamics with HHD; encode them in sparse features and represent the propagation of brain activity between distinct cortical areas. Using HHD, we decoded the somatosensory N20 component into two HHD features and represented the dynamics of brain activity as a traveling source between two primary somatosensory regions. In the epilepsy patient, we displayed the propagation of the epileptic activity around the margins of a brain lesion. Our findings indicate that HHD measures computed from cortical dynamics can: (i) quantitatively access the cortical dynamics in both healthy and disease brain in terms of sparse features and dynamic brain activity propagation between distinct cortical areas, and (ii) facilitate a reproducible, automated analysis of experimental and clinical MEG/EEG source imaging data

Rapidly re-computable EEG (electroencephalography) forward models for realistic head shapes

Solution of the EEG source localization (inverse) problem utilizing model-based methods typically requires a significant number of forward model evaluations. For subspace based inverse methods like MUSIC [6], the total number of forward model evaluations can often approach an order of 10{sup 3} or 10{sup 4}. Techniques based on least-squares minimization may require significantly more evaluations. The observed set of measurements over an M-sensor array is often expressed as a linear forward spatio-temporal model of the form: F = GQ + N (1) where the observed forward field F (M-sensors x N-time samples) can be expressed in terms of the forward model G, a set of dipole moment(s) Q (3xP-dipoles x N-time samples) and additive noise N. Because of their simplicity, ease of computation, and relatively good accuracy, multi-layer spherical models [7] (or fast approximations described in [1], [7]) have traditionally been the 'forward model of choice' for approximating the human head. However, approximation of the human head via a spherical model does have several key drawbacks. By its very shape, the use of a spherical model distorts the true distribution of passive currents in the skull cavity. Spherical models also require that the sensor positions be projected onto the fitted sphere (Fig. 1), resulting in a distortion of the true sensor-dipole spatial geometry (and ultimately the computed surface potential). The use of a single 'best-fitted' sphere has the added drawback of incomplete coverage of the inner skull region, often ignoring areas such as the frontal cortex. In practice, this problem is typically countered by fitting additional sphere(s) to those region(s) not covered by the primary sphere. The use of these additional spheres results in added complication to the forward model. Using high-resolution spatial information obtained via X-ray CT or MR imaging, a realistic head model can be formed by tessellating the head into a set of contiguous regions (typically the scalp, outer skull, and inner skull surfaces). Since accurate in vivo determination of internal conductivities is currently not currently possible, the head is typically assumed to consist of a set of contiguous isotropic regions, each with constant conductivity

The past, present, and future of the brain imaging data structure (BIDS)

The Brain Imaging Data Structure (BIDS) is a community-driven standard for the organization of data and metadata from a growing range of neuroscience modalities. This paper is meant as a history of how the standard has developed and grown over time. We outline the principles behind the project, the mechanisms by which it has been extended, and some of the challenges being addressed as it evolves. We also discuss the lessons learned through the project, with the aim of enabling researchers in other domains to learn from the success of BIDS

Mapping neurotransmitter systems to the structural and functional organization of the human neocortex

Neurotransmitter receptors support the propagation of signals in the human brain. How receptor systems are situated within macro-scale neuroanatomy and how they shape emergent function remain poorly understood, and there exists no comprehensive atlas of receptors. Here we collate positron emission tomography data from more than 1,200 healthy individuals to construct a whole-brain three-dimensional normative atlas of 19 receptors and transporters across nine different neurotransmitter systems. We found that receptor profiles align with structural connectivity and mediate function, including neurophysiological oscillatory dynamics and resting-state hemodynamic functional connectivity. Using the Neurosynth cognitive atlas, we uncovered a topographic gradient of overlapping receptor distributions that separates extrinsic and intrinsic psychological processes. Finally, we found both expected and novel associations between receptor distributions and cortical abnormality patterns across 13 disorders. We replicated all findings in an independently collected autoradiography dataset. This work demonstrates how chemoarchitecture shapes brain structure and function, providing a new direction for studying multi-scale brain organization.</p

The past, present, and future of the Brain Imaging Data Structure (BIDS)

The Brain Imaging Data Structure (BIDS) is a community-driven standard for the organization of data and metadata from a growing range of neuroscience modalities. This paper is meant as a history of how the standard has developed and grown over time. We outline the principles behind the project, the mechanisms by which it has been extended, and some of the challenges being addressed as it evolves. We also discuss the lessons learned through the project, with the aim of enabling researchers in other domains to learn from the success of BIDS

Age-Related Reduced Somatosensory Gating Is Associated with Altered Alpha Frequency Desynchronization

Sensory gating (SG), referring to an attenuated neural response to the second identical stimulus, is considered as preattentive processing in the central nervous system to filter redundant sensory inputs. Insufficient somatosensory SG has been found in the aged adults, particularly in the secondary somatosensory cortex (SII). However, it remains unclear which variables leading to the age-related somatosensory SG decline. There has been evidence showing a relationship between brain oscillations and cortical evoked excitability. Thus, this study used whole-head magnetoencephalography to record responses to paired-pulse electrical stimulation to the left median nerve in healthy young and elderly participants to test whether insufficient stimulus 1- (S1-) induced event-related desynchronization (ERD) contributes to a less-suppressed stimulus 2- (S2-) evoked response. Our analysis revealed that the minimum norm estimates showed age-related reduction of SG in the bilateral SII regions. Spectral power analysis showed that the elderly demonstrated significantly reduced alpha ERD in the contralateral SII (SIIc). Moreover, it was striking to note that lower S1-induced alpha ERD was associated with higher S2-evoked amplitudes in the SIIc among the aged adults. Conclusively, our findings suggest that age-related decline of somatosensory SG is partially attributed to the altered S1-induced oscillatory activity

Hybrid MEG (Magnetoencephalography) source characterization by cortical remapping and imaging of parametric source models

Reliable estimation of the local spatial extent of neural activity is a key to the quantitative analysis of MEG sources across subjects and conditions. In association with an understanding of the temporal dynamics among multiple areas, this would represent a major advance in electrophysiological source imaging. Parametric current dipole approaches to MEG (and EEG) source localization can rapidly generate a physical model of neural current generators using a limited number of parameters. However, physiological interpretation of these models is often difficult, especially in terms of the spatial extent of the true cortical activity. In new approaches using multipolar source models [3, 5], similar problems remain in the analysis of the higher-order source moments as parameters of cortical extent. Image-based approaches to the inverse problem provide a direct estimate of cortical current generators, but computationally expensive nonlinear methods are required to produce focal sources [1,4]. Recent efforts describe how a cortical patch can be grown until a best fit to the data is reached in the least-squares sense [6], but computational considerations necessitate that the growth be seeded in predefined regions of interest. In a previous study [2], a source obtained using a parametric model was remapped onto the cortex by growing a patch of cortical dipoles in the vicinity of the parametric source until the forward MEG or EEG fields of the parametric and cortical sources matched. The source models were dipoles and first-order multipoles. We propose to combine the parametric and imaging methods for MEG source characterization to take advantage of (i) the parsimonious and computationally efficient nature of parametric source localization methods and (ii) the anatomical and physiological consistency of imaging techniques that use relevant a priori information. By performing the cortical remapping imaging step by matching the multipole expansions of the original parametric source and the equivalent cortical patch, rather than their forward fields, we achieve significant reductions in computational complexity

Identification of growth seeds in the neonate brain through surfacic Helmholtz decomposition

We report on a new framework to investigate the rapid brain development of newborns. It is based on the analysis of depth maps of the cortical surface through the study of a displacement field estimated by surfacic optical flow methods. This displacement field shows local evolution of sulci directly on the cortical surface. Detection of its critical points is performed with the Helmholtz decomposition which allows us to identify sources of the developmental process. They can be viewed as growth seeds or in other terms points around which the sulcal growth organizes itself. We show the reproducibility of such growth seeds across 4 neonates and make a link of this new concept to the "sulcal roots" one proposed to explain the variability of human brain anatomy

Simultaneous MEG and intracranial EEG recordings during attentive reading.

International audienceThe relationship between neural oscillations recorded at various spatial scales remains poorly understood partly due to an overall dearth of studies utilizing simultaneous measurements. In an effort to study quantitative markers of attention during reading, we performed simultaneous magnetoencephalography (MEG) and intracranial electroencephalography (iEEG) recordings in four epileptic patients. Patients were asked to attend to a specific color when presented with an intermixed series of red words and green words, with words of a given color forming a cohesive story. We analyzed alpha, beta, and gamma band oscillatory responses to the word presentation and compared the strength and spatial organization of those responses in both electrophysiological recordings. Time-frequency analysis of iEEG revealed a network of clear attention-modulated high gamma band (50-150 Hz) power increases and alpha/beta (9-25 Hz) suppressions in response to the words. In addition to analyses at the sensor level, MEG time-frequency analysis was performed at the source level using a sliding window beamformer technique. Strong alpha/beta suppressions were observed in MEG reconstructions, in tandem with iEEG effects. While the MEG counterpart of high gamma band enhancement was difficult to interpret at the sensor level in two patients, MEG time-frequency source reconstruction revealed additional activation patterns in accordance with iEEG results. Importantly, iEEG allowed us to confirm that several sources of gamma band modulation observed with MEG were indeed of cortical origin rather than EMG muscular or ocular artifact