Signal Space Separation Research Articles

Optimizing magnetometers arrays and analysis pipelines for multivariate pattern analysis

BackgroundMultivariate pattern analysis (MVPA) has proven an excellent tool in cognitive neuroscience. It also holds a strong promise when applied to optically-pumped magnetometer-based magnetoencephalography. New methodTo optimize OPM-MEG systems for MVPA experiments this study examines data from a conventional MEG magnetometer array, focusing on appropriate noise reduction techniques for magnetometers. We determined the least required number of sensors needed for robust MVPA for image categorization experiments. ResultsWe found that the use of signal space separation (SSS) without a proper regularization significantly lowered the classification accuracy considering a sub-array of 102 magnetometers or a sub-array of 204 gradiometers. We also found that classification accuracy did not improve when going beyond 30 sensors irrespective of whether SSS has been applied. Comparison with existing methodsThe power spectra of data filtered with SSS has a substantially higher noise floor that data cleaned with SSP or HFC. Consequently, MVPA decoding results obtained from the SSS-filtered data are significantly lower compared to all other methods employed. ConclusionsWhen designing MEG system based on SQUID magnetometers optimized for multivariate analysis for image categorization experiments, about 30 magnetometers are sufficient. We advise against applying SSS filters without a proper regularization to data from MEG and OPM systems prior to performing MVPA as this method, albeit reducing low-frequency external noise contributions, also introduces an increase in broadband noise. We recommend employing noise reduction techniques that either decrease or maintain the noise floor of the data like signal-space projection, homogeneous field correction and gradient noise reduction.

Origin coordinate influence on performance of temporally extended signal space separation in magnetoencephalography

ObjectiveTemporally extended signal space separation (tSSS) is a powerful method for artifact suppression in magnetoencephalography (MEG). Because tSSS first separates MEG signals coming from inside and outside a certain sphere, definition of the sphere origin is important. For this study, we explored the influence of origin choice on tSSS performance in spontaneous and evoked activity from epilepsy patients. MethodsInterictal epileptiform discharges (IEDs) and somatosensory evoked fields (SEFs) were processed with two tSSSs: one with the default origin of (0, 0, 40 mm) in the head coordinate, and the other with an individual origin estimated using each patient’s anatomical magnetic resonance imaging (MRI). Equivalent current dipoles (ECDs) were calculated for the data. The ECD location and quality of estimation were compared across conditions. ResultsMEG data from 21 patients revealed marginal differences in ECD location, but the estimation quality inferred from goodness of fit (GOF) and confidence volume (CV) was better for the tSSS with individual origins. This choice affected IEDs more than it affected SEFs. ConclusionsIndividual sphere model resulted in better GOF and CV. SignificanceApplication of tSSS using an individual origin would be more desirable when available. This parameter might influence spontaneous activity more strongly.

Adaptive multipole models of optically pumped magnetometer data.

Multipole expansions have been used extensively in the Magnetoencephalography (MEG) literature for mitigating environmental interference and modelling brain signal. However, their application to Optically Pumped Magnetometer (OPM) data is challenging due to the wide variety of existing OPM sensor and array designs. We therefore explore how such multipole models can be adapted to provide stable models of brain signal and interference across OPM systems. Firstly, we demonstrate how prolate spheroidal (rather than spherical) harmonics can provide a compact representation of brain signal when sampling on the scalp surface with as few as 100 channels. We then introduce a type of orthogonal projection incorporating this basis set. The Adaptive Multipole Models (AMM), which provides robust interference rejection across systems, even in the presence of spatially structured nonlinearity errors (shielding factor is the reciprocal of the maximum fractional nonlinearity error). Furthermore, this projection is always stable, as it is an orthogonal projection, and will only ever decrease the white noise in the data. However, for array designs that are suboptimal for spatially separating brain signal and interference, this method can remove brain signal components. We contrast these properties with the more typically used multipole expansion, Signal Space Separation (SSS), which never reduces brain signal amplitude but is less robust to the effect of sensor nonlinearity errors on interference rejection and can increase noise in the data if the system is sub-optimally designed (as it is an oblique projection). We conclude with an empirical example utilizing AMM to maximize signal to noise ratio (SNR) for the stimulus locked neuronal response to a flickering visual checkerboard in a 128-channel OPM system and demonstrate up to 40 dB software shielding in real data.

Magnetoencephalography (MEG) Data Processing in Epilepsy Patients with Implanted Responsive Neurostimulation (RNS) Devices.

Drug-resistant epilepsy (DRE) is often treated with surgery or neuromodulation. Specifically, responsive neurostimulation (RNS) is a widely used therapy that is programmed to detect abnormal brain activity and intervene with tailored stimulation. Despite the success of RNS, some patients require further interventions. However, having an RNS device in situ is a hindrance to the performance of neuroimaging techniques. Magnetoencephalography (MEG), a non-invasive neurophysiologic and functional imaging technique, aids epilepsy assessment and surgery planning. MEG performed post-RNS is complicated by signal distortions. This study proposes an independent component analysis (ICA)-based approach to enhance MEG signal quality, facilitating improved assessment for epilepsy patients with implanted RNS devices. Three epilepsy patients, two with RNS implants and one without, underwent MEG scans. Preprocessing included temporal signal space separation (tSSS) and an automated ICA-based approach with MNE-Python. Power spectral density (PSD) and signal-to-noise ratio (SNR) were analyzed, and MEG dipole analysis was conducted using single equivalent current dipole (SECD) modeling. The ICA-based noise removal preprocessing method substantially improved the signal-to-noise ratio (SNR) for MEG data from epilepsy patients with implanted RNS devices. Qualitative assessment confirmed enhanced signal readability and improved MEG dipole analysis. ICA-based processing markedly enhanced MEG data quality in RNS patients, emphasizing its clinical relevance.

Spatiotemporal signal space separation for regions of interest: Application for extracting neuromagnetic responses evoked by deep brain stimulation.

Magnetoencephalography (MEG) recordings are often contaminated by interference that can exceed the amplitude of physiological brain activity by several orders of magnitude. Furthermore, the activity of interference sources may spatially extend (known as source leakage) into the activity of brain signals of interest, resulting in source estimation inaccuracies. This problem is particularly apparent when using MEG to interrogate the effects of brain stimulation on large-scale cortical networks. In this technical report, we develop a novel denoising approach for suppressing the leakage of interference source activity into the activity representing a brain region of interest. This approach leverages spatial and temporal domain projectors for signal arising from prespecified anatomical regions of interest. We apply this denoising approach to reconstruct simulated evoked response topographies to deep brain stimulation (DBS) in a phantom recording. We highlight the advantages of our approach compared to the benchmark-spatiotemporal signal space separation-and show that it can more accurately reveal brain stimulation-evoked response topographies. Finally, we apply our method to MEG recordings from a single patient with Parkinson's disease, to reveal early cortical-evoked responses to DBS of the subthalamic nucleus.

OPM-MEG bad channel identification method based on the improved box-isolation forest algorithm

Bad channels result from the inherent instability of wearable optically pumped magnetometers (OPMs). These unreliable channels significantly impede subsequent signal analysis. The signal space separation method for OPMs, which relies on time-domain features such as the z-score, cannot accurately identify these unreliable channels. We introduce a novel approach for automatic bad channel identification in OPMs, utilizing frequency domain features. Box-plot technique enables predicting the number of unreliable channels. The predicted ratio of unreliable channels is incorporated as a crucial parameter in the Isolation Forest algorithm to achieve the automated identification of bad channels. To evaluate the effectiveness of our method, we employed three accuracy metrics and conducted comparative assessments via the following methods: box-plot, density-based spatial clustering of applications with noise, and local outlier factor. Our approach is validated via simulation experiments and real measurements using OPM-magnetoencephalography data. The efficiency of subsequent OPM-MEG preprocessing and source location analysis is improved.

Investigating the technical feasibility of magnetoencephalography during transcranial direct current stimulation.

Magnetoencephalography (MEG) can measure weak magnetic fields produced by electrical brain activity. Transcranial direct current stimulation (tDCS) can affect such brain activities. The concurrent application of both, however, is challenging because tDCS presents artifacts on the MEG signal. If brain activity during tDCS can be elucidated by MEG, mechanisms of plasticity-inducing and other effects of tDCS would be more comprehensively understood. We tested the technical feasibility of MEG during tDCS using a phantom that produces an artificial current dipole simulating focal brain activity. An earlier study investigated estimation of a single oscillating phantom dipole during tDCS, and we systematically tested multiple dipole locations with a different MEG device. A phantom provided by the manufacturer was used to produce current dipoles from 32 locations. For the 32 dipoles, MEG was recorded with and without tDCS. Temporally extended signal space separation (tSSS) was applied for artifact rejection. Current dipole sources were estimated as equivalent current dipoles (ECDs). The ECD modeling quality was assessed using localization error, amplitude error, and goodness of fit (GOF). The ECD modeling performance with and without tDCS, and with and without tSSS was assessed. Mean localization errors of the 32 dipoles were 1.70 ± 0.72 mm (tDCS off, tSSS off, mean ± standard deviation), 6.13 ± 3.32 mm (tDCS on, tSSS off), 1.78 ± 0.83 mm (tDCS off, tSSS on), and 5.73 ± 1.60 mm (tDCS on, tSSS on). Mean GOF findings were, respectively, 92.3, 87.4, 97.5, and 96.7%. Modeling was affected by tDCS and restored by tSSS, but improvement of the localization error was marginal, even with tSSS. Also, the quality was dependent on the dipole location. Concurrent tDCS-MEG recording is feasible, especially when tSSS is applied for artifact rejection and when the assumed location of the source of activity is favorable for modeling. More technical studies must be conducted to confirm its feasibility with different source modeling methods and stimulation protocols. Recovery of single-trial activity under tDCS warrants further research.

An Iterative Implementation of the Signal Space Separation Method for Magnetoencephalography Systems with Low Channel Counts.

The signal space separation (SSS) method is routinely employed in the analysis of multichannel magnetic field recordings (such as magnetoencephalography (MEG) data). In the SSS method, signal vectors are posed as a multipole expansion of the magnetic field, allowing contributions from sources internal and external to a sensor array to be separated via computation of the pseudo-inverse of a matrix of the basis vectors. Although powerful, the standard implementation of the SSS method on MEG systems based on optically pumped magnetometers (OPMs) is unstable due to the approximate parity of the required number of dimensions of the SSS basis and the number of channels in the data. Here we exploit the hierarchical nature of the multipole expansion to perform a stable, iterative implementation of the SSS method. We describe the method and investigate its performance via a simulation study on a 192-channel OPM-MEG helmet. We assess performance for different levels of truncation of the SSS basis and a varying number of iterations. Results show that the iterative method provides stable performance, with a clear separation of internal and external sources.

Optimization of Signal Space Separation for Optically Pumped Magnetometer in Magnetoencephalography.

Magnetoencephalography (MEG) is a noninvasive functional neuroimaging modality but highly susceptible to environmental interference. Signal space separation (SSS) is a method for improving the SNR to separate the MEG signals from external interference. The origin and truncation values of SSS significantly affect the SSS performance. The origin value fluctuates with respect to the helmet array, and determining the truncation values using the traversal method is time-consuming; thus, this method is inappropriate for optically pumped magnetometer (OPM) systems with flexible array designs. Herein, an automatic optimization method for the SSS parameters is proposed. Virtual sources are set inside and outside the brain to simulate the signals of interest and interference, respectively, via forward model, with the sensor array as prior information. The objective function is determined as the error between the signals from simulated sources inside the brain and the SSS reconstructed signals; thus, the optimized parameters are solved inversely by minimizing the objective function. To validate the proposed method, a simulation analysis and MEG auditory-evoked experiments were conducted. For an OPM sensor array, this method can precisely determine the optimized origin and truncation values of the SSS simultaneously, and the auditory-evoked component, for example, N100, can be accurately located in the temporal cortex. The proposed optimization procedure outperforms the traditional method with regard to the computation time and accuracy, simplifying the SSS process in signal preprocessing and enhancing the performance of SSS denoising.

Development of Magnetocardiograph without Magnetically Shielded Room Using High-Detectivity TMR Sensors.

A magnetocardiograph that enables the clear observation of heart magnetic field mappings without magnetically shielded rooms at room temperatures has been successfully manufactured. Compared to widespread electrocardiographs, magnetocardiographs commonly have a higher spatial resolution, which is expected to lead to early diagnoses of ischemic heart disease and high diagnostic accuracy of ventricular arrhythmia, which involves the risk of sudden death. However, as the conventional superconducting quantum interference device (SQUID) magnetocardiographs require large magnetically shielded rooms and huge running costs to cool the SQUID sensors, magnetocardiography is still unfamiliar technology. Here, in order to achieve the heart field detectivity of 1.0 pT without magnetically shielded rooms and enough magnetocardiography accuracy, we aimed to improve the detectivity of tunneling magnetoresistance (TMR) sensors and to decrease the environmental and sensor noises with a mathematical algorithm. The magnetic detectivity of the TMR sensors was confirmed to be 14.1 pTrms on average in the frequency band between 0.2 and 100 Hz in uncooled states, thanks to the original multilayer structure and the innovative pattern of free layers. By constructing a sensor array using 288 TMR sensors and applying the mathematical magnetic shield technology of signal space separation (SSS), we confirmed that SSS reduces the environmental magnetic noise by -73 dB, which overtakes the general triple magnetically shielded rooms. Moreover, applying digital processing that combined the signal average of heart magnetic fields for one minute and the projection operation, we succeeded in reducing the sensor noise by about -23 dB. The heart magnetic field resolution measured on a subject in a laboratory in an office building was 0.99 pTrms and obtained magnetocardiograms and current arrow maps as clear as the SQUID magnetocardiograph does in the QRS and ST segments. Upon utilizing its superior spatial resolution, this magnetocardiograph has the potential to be an important tool for the early diagnosis of ischemic heart disease and the risk management of sudden death triggered by ventricular arrhythmia.

Environmental Magnetic Noise Reduction of Flat Multichannel Optically Pumped Magnetometer with Signal Space Separation Method together with Auxiliary Sensors

In this study, we carried out computer simulations to apply the SSS method to flat multi-channel OPMs with auxiliary sensors and examined the suppression effect of environmental magnetic noise. We assumed that auxiliary sensors are added so that the environmental magnetic noise can be effectively suppressed, and assessed the appropriate position of the sensors. When the auxiliary sensors were circularly aligned with a radius of 50 cm in the plane 40 cm or 30 cm above from the flat multi-channel OPMs with the origin of the SSS method at 10 cm above from the OPMs, the target signal was extracted as an internal component. In addition, in the case of these sensor arrangements, the S/N of 4.3 could be improved to 259.

Extended Signal-Space Separation Method for Improved Interference Suppression in MEG.

Objective:Magnetoencephalography (MEG) signals typically reflect a mixture of neuromagnetic fields, subject-related artifacts, external interference and sensor noise. Even inside a magnetically shielded room, external interference can be significantly stronger than brain signals. Methods such as signal-space projection (SSP) and signal-space separation (SSS) have been developed to suppress this residual interference, but their performance might not be sufficient in cases of strong interference or when the sources of interference change over time.Methods:Here we suggest a new method, extended signal-space separation (eSSS), which combines a physical model of the magnetic fields (as in SSS) with a statistical description of the interference (as in SSP). We demonstrate the performance of this method via simulations and experimental MEG data.Results:The eSSS method clearly outperforms SSS and SSP in interference suppression regardless of the extent of a priori information available on the interference sources. We also show that the method does not cause location or amplitude bias in dipole modeling.Conclusion:Our eSSS method provides better data quality than SSP or SSS and can be readily combined with other SSS-based methods, such as spatiotemporal SSS or head movement compensation. Thus, eSSS extends and complements the interference suppression techniques currently available for MEG.Significance:Due to its ability to suppress external interference to the level of sensor noise, eSSS can facilitate single-trial data analysis, exemplified in automated analysis of epileptic data. Such an enhanced suppression is especially important in environments with large interference fields.

The comparative performance of DBS artefact rejection methods for MEG recordings

Deep brain stimulation (DBS) can be a very efficient treatment option for movement disorders and psychiatric diseases. To better understand DBS mechanisms, brain activity can be recorded using magnetoencephalography (MEG) with the stimulator turned on. However, DBS produces large artefacts compromising MEG data quality due to both the applied current and the movement of wires connecting the stimulator with the electrode. To filter out these artefacts, several methods to suppress the DBS artefact have been proposed in the literature. A comparative study evaluating each method’s effectiveness, however, is missing so far.In this study, we evaluate the performance of four artefact rejection methods on MEG data from phantom recordings with DBS acquired with an Elekta Neuromag and a CTF system: (i) Hampel-filter, (ii) spectral signal space projection (S3P), (iii) independent component analysis with mutual information (ICA-MI), and (iv) temporal signal space separation (tSSS). In the sensor space, the largest increase in signal-to-noise (SNR) ratio was achieved by ICA-MI, while the best correspondence in terms of source activations was obtained by tSSS. LCMV beamforming alone was not sufficient to suppress the DBS-induced artefacts.

Effectively combining temporal projection noise suppression methods in magnetoencephalography

BackgroundMagnetoencephalography (MEG) is an excellent non-invasive tool to study the brain. However, measurements often suffer from the contribution of external interference, including noise from the sensors. Suppression of noise from the data is critical for an accurate representation of brain signals. Due to MEG’s limited spatial resolution and superior temporal resolution, noise suppression methods that operate in the temporal domain can be favorable. New MethodWe examined the independent and joint effects of two temporal projection noise suppression algorithms for MEG measurements: One commonly used algorithm which suppresses correlated noise; temporal signal space separation (tSSS) and one new method which suppresses uncorrelated sensor noise; oversampled temporal projection (OTP). ResultsWe found that both OTP and tSSS effectively suppress noise in raw MEG data and have the greatest effect of joint operation in cases where SNR is low, or when detecting higher SNR single-trial responses from raw data. We additionally demonstrate how the combination of OTP and tSSS is useful for the detectability of high-frequency brain oscillations (HFO). Comparison with existing MethodsAlthough the mathematical description of OTP has been described before (Larson and Taulu, 2017), OTP’s effect on HFOs in MEG data is novel. Additionally, the combination of OTP and commonly used temporal noise suppression algorithms (i.e., tSSS) has not been shown. ConclusionsThis finding is applicable to clinical populations such as epilepsy, where HFO signals are thought to be important markers for areas of seizure onset and are typically difficult to detect with non-invasive neuroimaging methods.

Comparison of beamformer implementations for MEG source localization

Beamformers are applied for estimating spatiotemporal characteristics of neuronal sources underlying measured MEG/EEG signals. Several MEG analysis toolboxes include an implementation of a linearly constrained minimum-variance (LCMV) beamformer. However, differences in implementations and in their results complicate the selection and application of beamformers and may hinder their wider adoption in research and clinical use. Additionally, combinations of different MEG sensor types (such as magnetometers and planar gradiometers) and application of preprocessing methods for interference suppression, such as signal space separation (SSS), can affect the results in different ways for different implementations. So far, a systematic evaluation of the different implementations has not been performed. Here, we compared the localization performance of the LCMV beamformer pipelines in four widely used open-source toolboxes (MNE-Python, FieldTrip, DAiSS (SPM12), and Brainstorm) using datasets both with and without SSS interference suppression.We analyzed MEG data that were i) simulated, ii) recorded from a static and moving phantom, and iii) recorded from a healthy volunteer receiving auditory, visual, and somatosensory stimulation. We also investigated the effects of SSS and the combination of the magnetometer and gradiometer signals. We quantified how localization error and point-spread volume vary with the signal-to-noise ratio (SNR) in all four toolboxes.When applied carefully to MEG data with a typical SNR (3–15 ​dB), all four toolboxes localized the sources reliably; however, they differed in their sensitivity to preprocessing parameters. As expected, localizations were highly unreliable at very low SNR, but we found high localization error also at very high SNRs for the first three toolboxes while Brainstorm showed greater robustness but with lower spatial resolution. We also found that the SNR improvement offered by SSS led to more accurate localization.

Comparison of DSSP and tSSS algorithms for removing artifacts from vagus nerve stimulators in magnetoencephalography data

Objective. Large-amplitude artifacts from vagus nerve stimulator (VNS) implants for refractory epilepsy affect magnetoencephalography (MEG) recordings and are difficult to reject, resulting in unusable data from this important population of patients who are frequently evaluated for surgical treatment of epilepsy. Here we compare the performance of two artifact removal algorithms for MEG data: dual signal subspace projection (DSSP) and temporally extended signal space separation (tSSS). Approach. Each algorithm’s performance was first evaluated in simulations. We then tested the performance of each algorithm on resting-state MEG data from patients with VNS implants. We also examined how each algorithm improved source localization of somatosensory evoked fields in patients with VNS implants. Main results. DSSP and tSSS algorithms have a similar ability to reject interference in both simulated and real MEG data if the origin location for tSSS is appropriately set. If the origin set for tSSS is inappropriate, the signal after tSSS can be distorted due to a mismatch between the internal region and the actual source space. Both DSSP and tSSS are able to remove large-amplitude artifacts from outside the brain. DSSP might be a better choice than tSSS when the choice of origin location for tSSS is difficult. Significance. Both DSSP and tSSS algorithms can recover distorted MEG recordings from people with intractable epilepsy and VNS implants, improving epileptic spike identification and source localization of both functional activity and epileptiform activity.

Quantitatively validating the efficacy of artifact suppression techniques to study the cortical consequences of deep brain stimulation with magnetoencephalography

Deep brain stimulation (DBS) is an established and effective treatment for several movement disorders and is being developed to treat a host of neuropsychiatric disorders including epilepsy, chronic pain, obsessive compulsive disorder, and depression. However, the neural mechanisms through which DBS produces therapeutic benefits, and in some cases unwanted side effects, in these disorders are only partially understood. Non-invasive neuroimaging techniques that can assess the neural effects of active stimulation are important for advancing our understanding of the neural basis of DBS therapy. Magnetoencephalography (MEG) is a safe, passive imaging modality with relatively high spatiotemporal resolution, which makes it a potentially powerful method for examining the cortical network effects of DBS. However, the degree to which magnetic artifacts produced by stimulation and the associated hardware can be suppressed from MEG data, and the comparability between signals measured during DBS-on and DBS-off conditions, have not been fully quantified. The present study used machine learning methods in conjunction with a visual perception task, which should be relatively unaffected by DBS, to quantify how well neural data can be salvaged from artifact contamination introduced by DBS and how comparable DBS-on and DBS-off data are after artifact removal. Machine learning also allowed us to determine whether the spatiotemporal pattern of neural activity recorded during stimulation are comparable to those recorded when stimulation is off. The spatiotemporal patterns of visually evoked neural fields could be accurately classified in all 8 patients with DBS implants during both DBS-on and DBS-off conditions and performed comparably across those two conditions. Further, the classification accuracy for classifiers trained on the spatiotemporal patterns evoked during DBS-on trials and applied to DBS-off trials, and vice versa, were similar to that of the classifiers trained and tested on either trial type, demonstrating the comparability of these patterns across conditions. Together, these results demonstrate the ability of MEG preprocessing techniques, like temporal signal space separation, to salvage neural data from recordings contaminated with DBS artifacts and validate MEG as a powerful tool to study the cortical consequences of DBS.

Tracking numerical and continuous magnitude processing by frequency-tagged neuromagnetic responses

Event Abstract Back to Event Tracking numerical and continuous magnitude processing by frequency-tagged neuromagnetic responses Amandine Van Rinsveld1*, Vincent Wens2, 3, Mathieu Guillaume1, Anthony Beuel1, Wim Gevers1, Xavier De Tiège2, 3 and Alain Content1 1 ULB Neuroscience Institute, Free University of Brussels, Center for Research in Cognition and Neurosciences , Belgium 2 ULB Neuroscience Institute, Free University of Brussels, Laboratoire de Cartographie fonctionnelle du Cerveau, Belgium 3 Erasmus Hospital, Free University of Brussels, Magnetoencephalography Unit, Department of Functional Neuroimaging, Service of Nuclear Medicine,, Belgium The ability to handle approximate large quantities has been identified as a building block of mathematical skills but the mechanisms allowing to extract numerical magnitudes (e.g., numerosity) from environmental stimuli is still debated. Most authors agree that humans have an approximate number system that specifically processes numerosity. However, a set of objects is not only characterized by its numerosity but also by additional visual information related to its continuous dimensions (e.g., object size), leading to alternative theories to explain numerosity extraction. Until now opposite views could not be tested properly due to the intrinsic correlations between numerosity and continuous dimensions. We isolated the specific responses to numerosity and other continuous dimensions using a 10 Hz-rate presentation of dot patterns that varied randomly in all dimensions but one dimension of interest, which was varied systematically at 1.25 Hz5 (see Fig.1). We recorded the evoked brain activity as frequency-tagged steady-state visual responses measured with magnetoencephalography. Neuromagnetic data were recorded in twenty-one subjects with a 306-channel MEG Triux system. Data were preprocessed by signal space separation (SSS) with movement correction and realignment on default head position. The amplitude spectra for all channels were obtained through Fourier transformation on the average of four 40 s-long epochs (frequency resolution of 0.025 Hz). To assess the neural synchronization at 1.25 Hz, we extracted and summed spectra in 0.6 Hz-wide frequency intervals around the 7 first harmonics of 1.25 Hz. The resulting sum-based amplitudes (SBA) were then converted into statistical maps of individual t scores (df=19) comparing the SBA at 1.25 Hz with the SBA at the surrounding frequencies (i.e., 10 adjacent frequency bins on both sides of the 1.25 Hz bin). At the group level, significance of the t maps was assessed using one-tailed maximum-based permutation tests against the null hypothesis ts = 0. This non-parametric approach automatically corrects for multiple comparisons across all sensors. We found neural synchronization at the presentation rate of 10 Hz in all sequences (average t=10.73, permutation p<.001). Crucially, we observed significant SBA peak at 1.25 Hz, corresponding to the periodic change of one dimension, for three dimensions: total area, number, and convex hull (ps<.01), all emerging at medial occipital sensors (see Fig.2). No significant synchronization occurred for density or dot size. In conclusion, by using an objective measure of implicit discrimination of each magnitude change, we observed neuromagnetic responses synchronized on the periodic change of numerosity, which is coherent with previous EEG results but also of total area and convex hull of the dot collections. In contrast, no synchronization was observed for other continuous dimensions. Numerosity, Area, and Convex hull can thus be rapidly and independently extracted in the visual stream, supporting the idea of an early visual number sense. At this stage of the processing, the discrimination of numerosity changes cannot be reduced as resulting from continuous dimensions’ discrimination. Figure Captions: Fig 1. Sequences of dot patterns were characterized by five dimensions: dot size, total area, convex hull, density, and numerosity. a. Time series of values (200 stimuli) for the five dimensions considered. In this example, the numerosity varies periodically at 1.25 Hz and the other continuous dimensions vary randomly. b. Frequency spectra of the corresponding time series after Fourier transformation. c. Illustration of a subset of stimuli. Fig 2. Topographic maps of group t values for SBA at 1.25 Hz and its harmonics in the conditions that disclosed significant SBA according to the maximum-based permutation test (H0: t=0). The star indicates significant channels, whose SBA spectrum is depicted on the right. Figure 1 Figure 2

Magnetoencephalographic Source Localization of the Eye Area of the Motor Homunculus.

A patient with intractable epilepsy, previous right frontal resection, and active vagus nerve stimulation (VNS) developed new onset quasi-continuous twitching around the left eye. Electroencephalography showed no correlate to the orbicularis oculi twitches apart from myographic potentials at the left supraorbital and anterior frontal electrodes. Magnetoencephalography was performed using spatiotemporal signal space separation to suppress magnetic artifacts associated with the VNS apparatus. Magnetoencephalographic source imaging performed on the data back-averaged from the left supraorbital myographic potentials revealed an intrasulcal cortical generator situated in the posterior wall of the right precentral gyrus representing the eye area of the motor homunculus.

Spontaneous sensorimotor cortical activity is suppressed by deep brain stimulation in patients with advanced Parkinson’s disease

Advanced Parkinson's disease (PD) is characterized by an excessive oscillatory beta band activity in the subthalamic nucleus (STN). Deep brain stimulation (DBS) of STN alleviates motor symptoms in PD and suppresses the STN beta band activity. The effect of DBS on cortical sensorimotor activity is more ambiguous; both increases and decreases of beta band activity have been reported. Non-invasive studies with simultaneous DBS are problematic due to DBS-induced artifacts. We recorded magnetoencephalography (MEG) from 16 advanced PD patients with and without STN DBS during rest and wrist extension. The strong magnetic artifacts related to stimulation were removed by temporal signal space separation. MEG oscillatory activity at 5–25 Hz was suppressed during DBS in a widespread frontoparietal region, including the sensorimotor cortex identified by the cortico-muscular coherence. The strength of suppression did not correlate with clinical improvement. Our results indicate that alpha and beta band oscillations are suppressed at the frontoparietal cortex by STN DBS in PD.