A Reproducible MEG/EEG Group Study With the MNE Software: Recommendations, Quality Assessments, and Good Practices

Limitations of traditional magnetoencephalography (MEG) exclude some important patient groups from MEG examinations, such as epilepsy patients with a vagus nerve stimulator, patients with magnetic particles on the head or having magnetic dental materials that cause severe movement-related artefact signals. Conventional interference rejection methods are not able to remove the artefacts originating this close to the MEG sensor array. For example, the reference array method is unable to suppress interference generated by sources closer to the sensors than the reference array, about 20-40 cm. The spatiotemporal signal space separation method proposed in this paper recognizes and removes both external interference and the artefacts produced by these nearby sources, even on the scalp. First, the basic separation into brain-related and external interference signals is accomplished with signal space separation based on sensor geometry and Maxwell's equations only. After this, the artefacts from nearby sources are extracted by a simple statistical analysis in the time domain, and projected out. Practical examples with artificial current dipoles and interference sources as well as data from real patients demonstrate that the method removes the artefacts without altering the field patterns of the brain signals.

Functional connectivity between cortical areas may appear as correlated time behavior of neural activity. It has been suggested that merging of separate features into a single percept ("binding") is associated with coherent gamma band activity across the cortical areas involved. Therefore, it would be of utmost interest to image cortico-cortical coherence in the working human brain. The frequency specificity and transient nature of these interactions requires time-sensitive tools such as magneto- or electroencephalography (MEG/EEG). Coherence between signals of sensors covering different scalp areas is commonly taken as a measure of functional coupling. However, this approach provides vague information on the actual cortical areas involved, owing to the complex relation between the active brain areas and the sensor recordings. We propose a solution to the crucial issue of proceeding beyond the MEG sensor level to estimate coherences between cortical areas. Dynamic imaging of coherent sources (DICS) uses a spatial filter to localize coherent brain regions and provides the time courses of their activity. Reference points for the computation of neural coupling may be based on brain areas of maximum power or other physiologically meaningful information, or they may be estimated starting from sensor coherences. The performance of DICS is evaluated with simulated data and illustrated with recordings of spontaneous activity in a healthy subject and a parkinsonian patient. Methods for estimating functional connectivities between brain areas will facilitate characterization of cortical networks involved in sensory, motor, or cognitive tasks and will allow investigation of pathological connectivities in neurological disorders.

The increase in spatiotemporal resolution of neuroimaging devices is accompanied by a trend towards more powerful multivariate analysis methods. Often it is desired to interpret the outcome of these methods with respect to the cognitive processes under study. Here we discuss which methods allow for such interpretations, and provide guidelines for choosing an appropriate analysis for a given experimental goal: For a surgeon who needs to decide where to remove brain tissue it is most important to determine the origin of cognitive functions and associated neural processes. In contrast, when communicating with paralyzed or comatose patients via brain-computer interfaces, it is most important to accurately extract the neural processes specific to a certain mental state. These equally important but complementary objectives require different analysis methods. Determining the origin of neural processes in time or space from the parameters of a data-driven model requires what we call a forward model of the data; such a model explains how the measured data was generated from the neural sources. Examples are general linear models (GLMs). Methods for the extraction of neural information from data can be considered as backward models, as they attempt to reverse the data generating process. Examples are multivariate classifiers. Here we demonstrate that the parameters of forward models are neurophysiologically interpretable in the sense that significant nonzero weights are only observed at channels the activity of which is related to the brain process under study. In contrast, the interpretation of backward model parameters can lead to wrong conclusions regarding the spatial or temporal origin of the neural signals of interest, since significant nonzero weights may also be observed at channels the activity of which is statistically independent of the brain process under study. As a remedy for the linear case, we propose a procedure for transforming backward models into forward models. This procedure enables the neurophysiological interpretation of the parameters of linear backward models. We hope that this work raises awareness for an often encountered problem and provides a theoretical basis for conducting better interpretable multivariate neuroimaging analyses. Copyright © 2013 The Authors. Published by Elsevier Inc. All rights reserved.