Electro/magnetoencephalography beamforming with spatial restrictions based on sparsity

We present a source localization method for electroencephalographic (EEG) and magnetoencephalographic (MEG) data which is based on an estimate of the sparsity obtained through the eigencanceler (EIG), which is a spatial filter whose weights are constrained to lie in the noise subspace. The EIG provides rejection of directional interferences while minimizing noise contributions and maintaining specified beam pattern constraints. In our case, the EIG is used to estimate the sparsity of the signal as a function of the position, then we use this information to spatially restrict the neural sources to locations out of the sparsity maxima. As proof of the concept, we incorporate this restriction in the “classical” linearly constrained minimum variance (LCMV) source localization approach in order to enhance its performance. We present numerical examples to evaluate the proposed method using realistically simulated EEG/MEG data for different signal-to-noise (SNR) conditions and various levels of correlation between sources, as well as real EEG/MEG measurements of median nerve stimulation. Our results show that the proposed method has the potential of reducing the bias on the search of neural sources in the classical approach, as well as making it more effective in localizing correlated sources.     

EEG-MEG Integration Enhances the Characterization of Functional and Effective Connectivity in the Resting State Network

At the sensor level many aspects, such as spectral power, functional and effective connectivity as well as relative-power-ratio ratio (RPR) and spatial resolution have been comprehensively investigated through both electroencephalography (EEG) and magnetoencephalography (MEG). Despite this, differences between both modalities have not yet been systematically studied by direct comparison. It remains an open question as to whether the integration of EEG and MEG data would improve the information obtained from the above mentioned parameters. Here, EEG (64-channel system) and MEG (275 sensor system) were recorded simultaneously in conditions with eyes open (EO) and eyes closed (EC) in 29 healthy adults. Spectral power, functional and effective connectivity, RPR, and spatial resolution were analyzed at five different frequency bands (delta, theta, alpha, beta and gamma). Networks of functional and effective connectivity were described using a spatial filter approach called the dynamic imaging of coherent sources (DICS) followed by the renormalized partial directed coherence (RPDC). Absolute mean power at the sensor level was significantly higher in EEG than in MEG data in both EO and EC conditions. At the source level, there was a trend towards a better performance of the combined EEG+MEG analysis compared with separate EEG or MEG analyses for the source mean power, functional correlation, effective connectivity for both EO and EC. The network of coherent sources and the spatial resolution were similar for both the EEG and MEG data if they were analyzed separately. Results indicate that the combined approach has several advantages over the separate analyses of both EEG and MEG. Moreover, by a direct comparison of EEG and MEG, EEG was characterized by significantly higher values in all measured parameters in both sensor and source level. All the above conclusions are specific to the resting state task and the specific analysis used in this study to have general conclusion multi-center studies would be helpful.     

基于粒子群优化算法的脑磁图源定位

脑磁图作为一种新型的脑探测技术,具有较高定位精度和毫秒级时间分辨率的特点。快速准确地利用脑磁图技术对三维空间中的脑神经活动源进行定位,对于脑功能研究和医学临床应用都具有重要的应用价值。可是,目前的脑磁图源定位广泛采用了多信号分类方法,它要求对三维大脑空间进行全局扫描,需要大量的计算,存在速度慢的缺点。针对这一问题,提出了一种基于粒子群优化算法的脑磁图源定位方法。先利用粒子群优化算法全局搜索能力强的特点寻找出目标函数的全局最优值,进行初步的脑磁图源定位;然后,再在小范围内进行小网格的搜索,进一步实现精确的定位。实验结果表明,基于粒子群优化算法的脑磁图源定位能够很好地解决上述问题,具有计算速度快、定位精度高的特点。     

Reconstructing Coherent Networks from Electroencephalography and Magnetoencephalography with Reduced Contamination from Volume Conduction or Magnetic Field Spread

Volume conduction (VC) and magnetic field spread (MFS) induce spurious correlations between EEG/MEG sensors, such that the estimation of functional networks from scalp recordings is inaccurate. Imaginary coherency [1] reduces VC/MFS artefacts between sensors by assuming that instantaneous interactions are caused predominantly by VC/MFS and do not contribute to the imaginary part of the cross-spectral densities (CSDs). We propose an adaptation of the dynamic imaging of coherent sources (DICS) [2] - a method for reconstructing the CSDs between sources, and subsequently inferring functional connectivity based on coherences between those sources. Firstly, we reformulate the principle of imaginary coherency by performing an eigenvector decomposition of the imaginary part of the CSD to estimate the power that only contributes to the non-zero phase-lagged (NZPL) interactions. Secondly, we construct an NZPL-optimised spatial filter with two a priori assumptions: (1) that only NZPL interactions exist at the source level and (2) the NZPL CSD at the sensor level is a good approximation of the projected source NZPL CSDs. We compare the performance of the NZPL method to the standard method by reconstructing a coherent network from simulated EEG/MEG recordings. We demonstrate that, as long as there are phase differences between the sources, the NZPL method reliably detects the underlying networks from EEG and MEG. We show that the method is also robust to very small phase lags, noise from phase jitter, and is less sensitive to regularisation parameters. The method is applied to a human dataset to infer parts of a coherent network underpinning face recognition.     

Investigation of a two-point maximum entropy regularization method for signal enhancement applied to magnetoencephalography data

A two-point maximum entropy method (TPMEM) was investigated for post-acquisition signal recovery in magnetoencephalography (MEG) data, as a potential replacement of a low-pass (LP) filtering technique currently in use. We first applied TPMEM and the LP filter for signal recovery of synthetically noise corrupted MEG “phantom” data sets in which the true underlying signal was known. Results were quantified with the use of visual plots, percent error histograms, and the statistical parameters root mean squared error and Pearson’s correlation coefficient. Synthetically noise corrupted data from a simulated magnetic dipole was used to quantify the improvements gained in using TPMEM over LP filters in reconstructing known dipole parameters such as position, orientation, and magnitude. Finally, we applied TPMEM and LP filters to a sample MEG patient data set. Our results show that TPMEM has improved noise-reduction and signal recovery capabilities than those of the LP filter, and furthermore data processed with TPMEM shows less error in the reconstructed dipole parameters. We propose that TPMEM can be used for MEG signal processing, resulting in improved MEG source characterization.     

Combined mapping of human auditory EEG and MEG responses

Auditory electric and magnetic P50(m), N1(m) and MMN(m) responses to standard, deviant and novel sounds were studied by recording brain electrical activity with 25 EEG electrodes simultaneously with the corresponding magnetic signals measured with 122 MEG gradiometer coils. The sources of these responses were located on the basis of the MEG responses; all were found to be in the supratemporal plane. The goal of the present paper was to investigate to what degree the source locations and orientations determined from the magnetic data account for the measured EEG signals. It was found that the electric P50, N1 and MMN responses can to a considerable degree be explained by the sources of the corresponding magnetic responses. In addition, source-current components not detectable by MEG were shown to contribute to the measured EEG signals.     

Accounting for Linear Transformations of EEG and MEG Data in Source Analysis

Analyses of electro- and magnetoencephalography (EEG, MEG) data often involve a linear modification of signals at the sensor level. Examples include re-referencing of the EEG, computation of synthetic gradiometer in MEG, or the removal of artifactual independent components to clean EEG and MEG data. A question of practical relevance is, if such modifications must be accounted for by adapting the physical forward model (leadfield) before subsequent source analysis. Here, we show that two scenarios need to be differentiated. In the first scenario, which corresponds to re-referencing the EEG and synthetic gradiometer computation in MEG, the leadfield must be adapted before source analysis. In the second scenario, which corresponds to removing artifactual components to ‘clean’ the data, the leadfield must not be changed. We demonstrate and discuss the consequences of wrongly modifying the leadfield in the latter case for an example. Future EEG and MEG studies employing source analyses should carefully consider whether and, if so, how the leadfield must be modified as explicated here.     

High-frequency (600 Hz) SEP activities originating in the subcortical and cortical human somatosensory system

Digitally high-pass filtered median nerve SEP show an oscillatory burst of low-amplitude high-frequency (600 Hz) wavelets superimposed on the N20 component which itself is generated by excitatory postsynaptic potentials of area 3b pyramidal cells. Prior studies using magnetoencephalography (MEG) localized one wavelet generator close to the primary somatosensory hand cortex. Since MEG recordings are biased towards tangentially oriented and superficial generators, a dipole source analysis of 32-channel electric SEP recordings was employed here to test for the possibility of deep and/or radially oriented burst generators: in 10 normal subjects low noise (16 000 averages) median nerve SEP were evaluated using dipole source analysis before and after applying a digital 475 Hz high-pass filter. Two main oscillatory 600 Hz burst sources were modeled; (i) a deep burst source close to the thalamus, most active in a time window between the brain-stem P14 and the cortical N20 sources, detectable in 7 of 10 subjects; most probably, this activity originates from deep axon segments of thalamocortical fibers; and (ii) a subsequent burst source timed around the N20 and located in the vicinity of the primary somatosensory hand cortex in all subjects, which was already known from MEG data. This superficial oscillatory source may be dominated by repetitive activity conducted in the terminal segments of the thalamocortical projection fibers initiated by the thalamic burst generator.     

Divergent Cortical Generators of MEG and EEG during Human Sleep Spindles Suggested by Distributed Source Modeling

Background

Sleep spindles are ∼1-second bursts of 10–15 Hz activity, occurring during normal stage 2 sleep. In animals, sleep spindles can be synchronous across multiple cortical and thalamic locations, suggesting a distributed stable phase-locked generating system. The high synchrony of spindles across scalp EEG sites suggests that this may also be true in humans. However, prior MEG studies suggest multiple and varying generators.

Methodology/Principal Findings

We recorded 306 channels of MEG simultaneously with 60 channels of EEG during naturally occurring spindles of stage 2 sleep in 7 healthy subjects. High-resolution structural MRI was obtained in each subject, to define the shells for a boundary element forward solution and to reconstruct the cortex providing the solution space for a noise-normalized minimum norm source estimation procedure. Integrated across the entire duration of all spindles, sources estimated from EEG and MEG are similar, diffuse and widespread, including all lobes from both hemispheres. However, the locations, phase and amplitude of sources simultaneously estimated from MEG versus EEG are highly distinct during the same spindles. Specifically, the sources estimated from EEG are highly synchronous across the cortex, whereas those from MEG rapidly shift in phase, hemisphere, and the location within the hemisphere.

Conclusions/Significance

The heterogeneity of MEG sources implies that multiple generators are active during human sleep spindles. If the source modeling is correct, then EEG spindles are generated by a different, diffusely synchronous system. Animal studies have identified two thalamo-cortical systems, core and matrix, that produce focal or diffuse activation and thus could underlie MEG and EEG spindles, respectively. Alternatively, EEG spindles could reflect overlap at the sensors of the same sources as are seen from the MEG. Although our results generally match human intracranial recordings, additional improvements are possible and simultaneous intra- and extra-cranial measures are needed to test their accuracy.     

Replicability of MEG and EEG measures of the auditory N1/N1m-response

We investigated the replicability of the source location, amplitude and latency measures of the auditory evoked N1 (EEG) and N1m (MEG) responses. Each of the 5 subjects was measured 6 times in two recording sessions. Responses to monaural stimuli were recorded from 122 MEG and 64 EEG channels simultaneously. The EEG data were modeled with a symmetrically-located dipole pair. For the MEG data, one dipole in each hemisphere was located independently using a subset of channels. Standard deviation (SD) was used as a measure for replicability. The average SD of the x, y and z coordinates of the contralateral N1m dipole was about 2 mm, whereas the corresponding figures for the ipsilateral N1m and the contra- and ipsilateral N1 were about twice as large. The SDs of the dipole amplitudes and latencies were almost equal with MEG and EEG. The amplitude and latency measures of the MEG field gradient waveforms were almost as replicable as those of the dipole models. The results suggest that both MEG and EEG can be used for investigating the simultaneous activity of the left and right auditory cortices independently, MEG being superior in certain experimental setups.     

Fermentation data analysis and state estimation in the presence of incomplete mass balance

A method is developed for identifying measurement errors and estimating fermentation states in the presence of unidentified reactant or product. Unlike conventional approaches using elemental balances, this method employs an empirically determined basis, which can tolerate unidentified reaction species. The essence of this approach is derived from the concept of reaction subspace and the technique of singular value decomposition. It is shown that the subspace determined via singular value decomposition of multiple experimental data provides an empirical basis for identifying measurement errors. The same approach is applied to fermentation state estimation. Via the formulation of the reaction subspace, the sensitivity of state estimates to measurement errors is quantified in terms of a dimensionless quantity, maximum error gain (MEG). It is shown that using the empirically determined subspace, one can circumvent the problem of unidentified reaction species, meanwhile reducing the sensitivity of the estimates.     

基于独立元分析的脑磁图源定位

在脑磁图源定位问题中,通常感兴趣的是脑内众多神经活动源中的一个或几个,而传统源定位方法,如多信号分类方法,需要将所有源的位置都确定后,通过重组信号波形才能获得所感兴趣源的位置信息。为了提高定位效率,文章作者提出了一种结合独立元分析的脑磁图源定位方法。实验结果表明,该方法能加快定位的速度,同时能够在一定程度上克服噪声的影响,具有更强的抗噪声能力。     

Population Level Inference for Multivariate MEG Analysis

Multivariate analysis is a very general and powerful technique for analysing Magnetoencephalography (MEG) data. An outstanding problem however is how to make inferences that are consistent over a group of subjects as to whether there are condition-specific differences in data features, and what are those features that maximise these differences. Here we propose a solution based on Canonical Variates Analysis (CVA) model scoring at the subject level and random effects Bayesian model selection at the group level. We apply this approach to beamformer reconstructed MEG data in source space. CVA estimates those multivariate patterns of activation that correlate most highly with the experimental design; the order of a CVA model is then determined by the number of significant canonical vectors. Random effects Bayesian model comparison then provides machinery for inferring the optimal order over the group of subjects. Absence of a multivariate dependence is indicated by the null model being the most likely. This approach can also be applied to CVA models with a fixed number of canonical vectors but supplied with different feature sets. We illustrate the method by identifying feature sets based on variable-dimension MEG power spectra in the primary visual cortex and fusiform gyrus that are maximally discriminative of data epochs before versus after visual stimulation.     

Assessment of Subcortical Source Localization Using Deep Brain Activity Imaging Model with Minimum Norm Operators: A MEG Study

Subcortical structures are involved in many healthy and pathological brain processes. It is crucial for many studies to use magnetoencephalography (MEG) to assess the ability to detect subcortical generators. This study aims to assess the source localization accuracy and to compare the characteristics of three inverse operators in the specific case of subcortical generators. MEG has a low sensitivity to subcortical sources mainly because of their distance from sensors and their complex cyto-architecture. However, we show that using a realistic anatomical and electrophysiological model of deep brain activity (DBA), the sources make measurable contributions to MEG sensors signals. Furthermore, we study the point-spread and cross-talk functions of the wMNE, sLORETA and dSPM inverse operators to characterize distortions in cortical and subcortical regions and to study how noise-normalization methods can improve or bias accuracy. We then run Monte Carlo simulations with neocortical and subcortical activations. In the case of single hippocampus patch activations, the results indicate that MEG can indeed localize the generators in the head and the body of the hippocampus with good accuracy. We then tackle the question of simultaneous cortical and subcortical activations. wMNE can detect hippocampal activations that are embedded in cortical activations that have less than double their amplitude, but it does not completely correct the bias to more superficial sources. dSPM and sLORETA can still detect hippocampal activity above this threshold, but such detection might include the creation of ghost deeper sources. Finally, using the DBA model, we showed that the detection of weak thalamic modulations of ongoing brain activity is possible.     

Modelling and detecting deep brain activity with MEG and EEG

We introduce an anatomical and electrophysiological model of deep brain structures dedicated to magnetoencephalography (MEG) and electroencephalography (EEG) source imaging. So far, most imaging inverse models considered that MEG/EEG surface signals were predominantly produced by cortical, hence superficial, neural currents. Here we question whether crucial deep brain structures such as the basal ganglia and the hippocampus may also contribute to distant, scalp MEG and EEG measurements. We first design a realistic anatomical and electrophysiological model of these structures and subsequently run Monte-Carlo experiments to evaluate the respective sensitivity of the MEG and EEG to signals from deeper origins. Results indicate that MEG/EEG may indeed localize these deeper generators, which is confirmed here from experimental MEG data reporting on the modulation of alpha (10–12 Hz) brain waves.     

脑磁源的定位研究

在脑电脑磁研究中,常将大脑视为一电磁系统,利用准静态近似下的麦克斯韦方程,可以发现代表脑神经元活动的原在电流密度与脑外磁场之间呈线性关系。在高斯哭喊杨存在的情况下,采用极大似然估计理论,讨论了一种基于脑磁场时域－罕注解磁源定位问题的一般方法。球对称导体模型下的模拟计算表明,这一方法是有效的。对于考虑真实头模型下的磁源定位问题求解磁源定位,提出了一种联合使用脑磁脑电数据的近似方法。     

从国际非专利名称纵观全球生物药发展

生物药(bio-therapeutics)是指采用生物技术制备的、临床上用于疾病治疗的大分子生物制品,具有结构复杂、异质性高等特点,科学严谨的生物药通用名命名,是区分生物药物质基础的主要依据,也是药品生命周期管理的重要基础。世界卫生组织(World Health Organization,WHO)协调建立的国际非专利名称(International Nonproprietary Names,INN)是全球药物命名的标准化体系。从INN的起源,以及生物药INN的类别、发生与发展为主线,以较为详实的数据统计和分析,呈现了全球生物药的衍化进程,从不同的角度纵览生物药技术发展历程,对生物药的研发设计、技术标准及监管策略的考量均具有一定的参考意义。     

High-resolution magnetoencephalography--studies with a small-volume phantom]

To investigate the spatiotemporal organisation of neuronal processes in an animal model using magnetoencephalography (MEG), a high temporal resolution (ms) and an appropriate spatial resolution of about 1 mm is necessary. With the aim of determining the localization error and the resolution power of high-resolution MEG systems, we developed a phantom capable of simulating the characteristics of animal models. The phantom enables us to variably position at least two magnetic field sources to within 0.1 mm. For source localization on the basis of the magnetic field data, a spatial filtering algorithm was used. The investigation of a 16-channel micro SQUID-MEG system with a current dipole orientated tangentially to the phantom surface produced the following localization data (min ... max, x, y--horizontal plane, z--depth); systematic localization error e(x) = 1.16 ... 1.67 mm, e(y) = -1.01 ... -1.28 mm, e(z) = -5.22 ... -7.64 mm, standard deviation of the individual measurements perpendicular to the dipole axis s(perp) = 0.05 ... 0.22 mm, along this axis s(long) = 0.20 ... 1.73 mm, in the depths sz = 0.17 ... 3.17 mm. The "goodness of fit" was > 95%. Separation of two dipoles was still possible for parallel dipoles at a distance apart of d(parallel) = 0.03 mm and for those oriented perpendicularly to each other at a distance apart of d(perp) = 0.10 mm. On the basis of these results we conclude that the MEG system can achieve a resolution sufficient to permit the investigation of neuronal microstructures. The spatial errors detected were related to sensor position in the cryostatic vessel as well as to external low-frequency noise.     

A multivariate approach for integrating genome-wide expression data and biological knowledge

MOTIVATION: Several statistical methods that combine analysis of differential gene expression with biological knowledge databases have been proposed for a more rapid interpretation of expression data. However, most such methods are based on a series of univariate statistical tests and do not properly account for the complex structure of gene interactions. RESULTS: We present a simple yet effective multivariate statistical procedure for assessing the correlation between a subspace defined by a group of genes and a binary phenotype. A subspace is deemed significant if the samples corresponding to different phenotypes are well separated in that subspace. The separation is measured using Hotelling's T(2) statistic, which captures the covariance structure of the subspace. When the dimension of the subspace is larger than that of the sample space, we project the original data to a smaller orthonormal subspace. We use this method to search through functional pathway subspaces defined by Reactome, KEGG, BioCarta and Gene Ontology. To demonstrate its performance, we apply this method to the data from two published studies, and visualize the results in the principal component space.     

Harmony: EEG/MEG Linear Inverse Source Reconstruction in the Anatomical Basis of Spherical Harmonics

EEG/MEG source localization based on a “distributed solution” is severely underdetermined, because the number of sources is much larger than the number of measurements. In particular, this makes the solution strongly affected by sensor noise. A new way to constrain the problem is presented. By using the anatomical basis of spherical harmonics (or spherical splines) instead of single dipoles the dimensionality of the inverse solution is greatly reduced without sacrificing the quality of the data fit. The smoothness of the resulting solution reduces the surface bias and scatter of the sources (incoherency) compared to the popular minimum-norm algorithms where single-dipole basis is used (MNE, depth-weighted MNE, dSPM, sLORETA, LORETA, IBF) and allows to efficiently reduce the effect of sensor noise. This approach, termed Harmony, performed well when applied to experimental data (two exemplars of early evoked potentials) and showed better localization precision and solution coherence than the other tested algorithms when applied to realistically simulated data.