Accuracy and precision of navigated transcranial magnetic stimulation.

Objective.Transcranial magnetic stimulation (TMS) induces an electric field (E-field) in the cortex. To facilitate stimulation targeting, image-guided neuronavigation systems have been introduced. Such systems track the placement of the coil with respect to the head and visualize the estimated cortical stimulation location on an anatomical brain image in real time. The accuracy and precision of the neuronavigation is affected by multiple factors. Our aim was to analyze how different factors in TMS neuronavigation affect the accuracy and precision of the coil-head coregistration and the estimated E-field.Approach.By performing simulations, we estimated navigation errors due to distortions in magnetic resonance images (MRIs), head-to-MRI registration (landmark- and surface-based registrations), localization and movement of the head tracker, and localization of the coil tracker. We analyzed the effect of these errors on coil and head coregistration and on the induced E-field as determined with simplistic and realistic head models.Main results.Average total coregistration accuracies were in the range of 2.2-3.6 mm and 1°; precision values were about half of the accuracy values. The coregistration errors were mainly due to head-to-MRI registration with average accuracies 1.5-1.9 mm/0.2-0.4° and precisions 0.5-0.8 mm/0.1-0.2° better with surface-based registration. The other major source of error was the movement of the head tracker with average accuracy of 1.5 mm and precision of 1.1 mm. When assessed within an E-field method, the average accuracies of the peak E-field location, orientation, and magnitude ranged between 1.5 and 5.0 mm, 0.9 and 4.8°, and 4.4 and 8.5% across the E-field models studied. The largest errors were obtained with the landmark-based registration. When computing another accuracy measure with the most realistic E-field model as a reference, the accuracies tended to improve from about 10 mm/15°/25% to about 2 mm/2°/5% when increasing realism of the E-field model.Significance.The results of this comprehensive analysis help TMS operators to recognize the main sources of error in TMS navigation and that the coregistration errors and their effect in the E-field estimation depend on the methods applied. To ensure reliable TMS navigation, we recommend surface-based head-to-MRI registration and realistic models for E-field computations.

µ-rhythm phase from somatosensory but not motor cortex correlates with corticospinal excitability in EEG-triggered TMS.

BACKGROUND: Sensorimotor µ-rhythm phase is correlated with corticospinal excitability. Transcranial magnetic stimulation (TMS) of motor cortex results in larger motor evoked potentials (MEPs) during the negative peak of the EEG oscillation as extracted with a surface Laplacian. However, the anatomical source of the relevant oscillation is not clear and demonstration of the relationship is sensitive to the choice of EEG montage. OBJECTIVE/HYPOTHESIS: Here, we compared two EEG montages preferentially sensitive to oscillations originating from the crown of precentral gyrus (dorsal premotor cortex) vs. postcentral gyrus (secondary somatosensory cortex). We hypothesized that the EEG signal from precentral gyrus would correlate more strongly with MEP amplitude, given that the corticospinal neurons are located in the anterior wall of the sulcus and the corticospinal tract has input from premotor cortex. NEW METHOD: Real-time EEG-triggered TMS of motor cortex was applied in 6 different conditions in randomly interleaved order, 3 phase conditions (positive peak, negative peak, random phase of the ongoing µ-oscillation), and each phase condition for 2 different EEG montages corresponding to oscillations preferentially originating in precentral gyrus (premotor cortex) vs. postcentral gyrus (somatosensory cortex), extracted using FCC3h vs. C3 centered EEG montages. RESULTS: The negative vs. positive peak of sensorimotor µ-rhythm as extracted from the C3 montage (postcentral gyrus, somatosensory cortex) correlated with states of high vs. low corticospinal excitability (p < 0.001), replicating previous findings. However, no significant correlation was found for sensorimotor µ-rhythm as extracted from the neighboring FCC3 montage (precentral gyrus, premotor cortex). This implies that EEG-signals from the somatosensory cortex are better predictors of corticospinal excitability than EEG-signals from the motor areas. CONCLUSIONS: The extraction of a brain oscillation whose phase corresponds to corticospinal excitability is highly sensitive to the selected EEG montage and the location of the EEG sensors on the scalp. Here, the cortical source of EEG oscillations predicting response amplitude does not correspond to the cortical target of the stimulation, indicating that even in this simple case, a specific neuronal pathway from somatosensory cortex to primary motor cortex is involved.

Towards an objective evaluation of EEG/MEG source estimation methods - The linear approach.

The spatial resolution of EEG/MEG source estimates, often described in terms of source leakage in the context of the inverse problem, poses constraints on the inferences that can be drawn from EEG/MEG source estimation results. Software packages for EEG/MEG data analysis offer a large choice of source estimation methods but few tools to experimental researchers for methods evaluation and comparison. Here, we describe a framework and tools for objective and intuitive resolution analysis of EEG/MEG source estimation based on linear systems analysis, and apply those to the most widely used distributed source estimation methods such as L2-minimum-norm estimation (L2-MNE) and linearly constrained minimum variance (LCMV) beamformers. Within this framework it is possible to define resolution metrics that define meaningful aspects of source estimation results (such as localization accuracy in terms of peak localization error, PLE, and spatial extent in terms of spatial deviation, SD) that are relevant to the task at hand and can easily be visualized. At the core of this framework is the resolution matrix, which describes the potential leakage from and into point sources (point-spread and cross-talk functions, or PSFs and CTFs, respectively). Importantly, for linear methods these functions allow generalizations to multiple sources or complex source distributions. This paper provides a tutorial-style introduction into linear EEG/MEG source estimation and resolution analysis aimed at experimental (rather than methods-oriented) researchers. We used this framework to demonstrate how L2-MNE-type as well as LCMV beamforming methods can be evaluated in practice using software tools that have only recently become available for routine use. Our novel methods comparison includes PLE and SD for a larger number of methods than in similar previous studies, such as unweighted, depth-weighted and normalized L2-MNE methods (including dSPM, sLORETA, eLORETA) and two LCMV beamformers. The results demonstrate that some methods can achieve low and even zero PLE for PSFs. However, their SD as well as both PLE and SD for CTFs are far less optimal for all methods, in particular for deep cortical areas. We hope that our paper will encourage EEG/MEG researchers to apply this approach to their own tasks at hand.

Closed-loop optimization of transcranial magnetic stimulation with electroencephalography feedback.

BACKGROUND: Transcranial magnetic stimulation (TMS) is widely used in brain research and treatment of various brain dysfunctions. However, the optimal way to target stimulation and administer TMS therapies, for example, where and in which electric field direction the stimuli should be given, is yet to be determined. OBJECTIVE: To develop an automated closed-loop system for adjusting TMS parameters (in this work, the stimulus orientation) online based on TMS-evoked brain activity measured with electroencephalography (EEG). METHODS: We developed an automated closed-loop TMS-EEG set-up. In this set-up, the stimulus parameters are electronically adjusted with multi-locus TMS. As a proof of concept, we developed an algorithm that automatically optimizes the stimulation orientation based on single-trial EEG responses. We applied the algorithm to determine the electric field orientation that maximizes the amplitude of the TMS-EEG responses. The validation of the algorithm was performed with six healthy volunteers, repeating the search twenty times for each subject. RESULTS: The validation demonstrated that the closed-loop control worked as desired despite the large variation in the single-trial EEG responses. We were often able to get close to the orientation that maximizes the EEG amplitude with only a few tens of pulses. CONCLUSION: Optimizing stimulation with EEG feedback in a closed-loop manner is feasible and enables effective coupling to brain activity.

Prefrontal theta phase-dependent rTMS-induced plasticity of cortical and behavioral responses in human cortex.

BACKGROUND: Prefrontal theta oscillations are involved in neuronal information transfer and retention. Phases along the theta cycle represent varied excitability states, whereby high-excitability states correspond to high-frequency neuronal activity and heightened capacity for plasticity induction, as demonstrated in animal studies. Human studies corroborate this model and suggest a core role of prefrontal theta activity in working memory (WM). OBJECTIVE/HYPOTHESIS: We aimed at modulating prefrontal neuronal excitability and WM performance in healthy humans, using real-time EEG analysis for triggering repetitive transcranial magnetic stimulation (rTMS) theta-phase synchronized to the left dorsomedial prefrontal cortex. METHODS: 16 subjects underwent 3 different rTMS interventions on separate days, with pulses triggered according to the individual's real-time EEG activity: 400 rTMS gamma-frequency (100 Hz) triplet bursts applied during either the negative peak of the prefrontal theta oscillation, the positive peak, or at random phase. Changes in cortical excitability were assessed with EEG responses following single-pulse TMS, and behavioral effects by using a WM task. RESULTS: Negative-peak rTMS increased single-pulse TMS-induced prefrontal theta power and theta-gamma phase-amplitude coupling, and decreased WM response time. In contrast, positive-peak rTMS decreased prefrontal theta power, while no changes were observed after random-phase rTMS. CONCLUSION: Findings point to the feasibility of EEG-TMS technology in a theta-gamma phase-amplitude coupling mode for effectively modifying WM networks in human prefrontal cortex, with potential for therapeutic applications.

Looking through the windows: a study about the dependency of phase-coupling estimates on the data length.

Objective. Being able to characterize functional connectivity (FC) state dynamics in a real-time setting, such as in brain-computer interface, neurofeedback or closed-loop neurostimulation frameworks, requires the rapid detection of the statistical dependencies that quantify FC in short windows of data. The aim of this study is to characterize, through extensive realistic simulations, the reliability of FC estimation as a function of the data length. In particular, we focused on FC as measured by phase-coupling (PC) of neuronal oscillations, one of the most functionally relevant neural coupling modes.Approach. We generated synthetic data corresponding to different scenarios by varying the data length, the signal-to-noise ratio (SNR), the phase difference value, the spectral analysis approach (Hilbert or Fourier) and the fractional bandwidth. We compared seven PC metrics, i.e. imaginary part of phase locking value (iPLV), PLV of orthogonalized signals, phase lag index (PLI), debiased weighted PLI, imaginary part of coherency, coherence of orthogonalized signals and lagged coherence.Main results. Our findings show that, for a SNR of at least 10 dB, a data window that contains 5-8 cycles of the oscillation of interest (e.g. a 500-800 ms window at 10 Hz) is generally required to achieve reliable PC estimates. In general, Hilbert-based approaches were associated with higher performance than Fourier-based approaches. Furthermore, the results suggest that, when the analysis is performed in a narrow frequency range, a larger window is required.Significance. The achieved results pave the way to the introduction of best-practice guidelines to be followed when a real-time frequency-specific PC assessment is at target.

Looking through the windows: a study about the dependency of phase-coupling estimates on the data length.

OBJECTIVE: Being able to characterize functional connectivity (FC) state dynamics in a real-time setting, such as in brain-computer interface, neurofeedback or closed-loop neurostimulation frameworks, requires the rapid detection of the statistical dependencies that quantify FC in short windows of data. The aim of this study is to characterize, through extensive realistic simulations, the reliability of FC estimation as a function of the data length. In particular, we focused on FC as measured by phase-coupling (PC) of neuronal oscillations, one of the most functionally relevant neural coupling modes. APPROACH: We generated synthetic data corresponding to different scenarios by varying the data length, the signal-to-noise ratio, the phase difference value, the spectral analysis approach (Hilbert or Fourier) and the fractional bandwidth. We compared seven PC metrics, i.e. imaginary part of phase locking value (PLV), PLV of orthogonalized signals, phase lag index (PLI), debiased weighted PLI, imaginary part of coherency, coherence of orthogonalized signals and lagged coherence. MAIN RESULTS: Our findings show that, for a signal-to-noise-ratio of at least 10 dB, a data window that contains 5 to 8 cycles of the oscillation of interest (e.g. a 500-800ms window at 10Hz) is generally required to achieve reliable PC estimates. In general, Hilbert-based approaches were associated with higher performance than Fourier-based approaches. Furthermore, the results suggest that, when the analysis is performed in a narrow frequency range, a larger window is required. SIGNIFICANCE: The achieved results pave the way to the introduction of best-practice guidelines to be followed when a real-time frequency-specific PC assessment is at target.

Multi-locus transcranial magnetic stimulation system for electronically targeted brain stimulation.

BACKGROUND: Transcranial magnetic stimulation (TMS) allows non-invasive stimulation of the cortex. In multi-locus TMS (mTMS), the stimulating electric field (E-field) is controlled electronically without coil movement by adjusting currents in the coils of a transducer. OBJECTIVE: To develop an mTMS system that allows adjusting the location and orientation of the E-field maximum within a cortical region. METHODS: We designed and manufactured a planar 5-coil mTMS transducer to allow controlling the maximum of the induced E-field within a cortical region approximately 30 mm in diameter. We developed electronics with a design consisting of independently controlled H-bridge circuits to drive up to six TMS coils. To control the hardware, we programmed software that runs on a field-programmable gate array and a computer. To induce the desired E-field in the cortex, we developed an optimization method to calculate the currents needed in the coils. We characterized the mTMS system and conducted a proof-of-concept motor-mapping experiment on a healthy volunteer. In the motor mapping, we kept the transducer placement fixed while electronically shifting the E-field maximum on the precentral gyrus and measuring electromyography from the contralateral hand. RESULTS: The transducer consists of an oval coil, two figure-of-eight coils, and two four-leaf-clover coils stacked on top of each other. The technical characterization indicated that the mTMS system performs as designed. The measured motor evoked potential amplitudes varied consistently as a function of the location of the E-field maximum. CONCLUSION: The developed mTMS system enables electronically targeted brain stimulation within a cortical region.

Spatial sampling of MEG and EEG based on generalized spatial-frequency analysis and optimal design.

In this paper, we analyze spatial sampling of electro- (EEG) and magnetoencephalography (MEG), where the electric or magnetic field is typically sampled on a curved surface such as the scalp. By simulating fields originating from a representative adult-male head, we study the spatial-frequency content in EEG as well as in on- and off-scalp MEG. This analysis suggests that on-scalp MEG, off-scalp MEG and EEG can benefit from up to 280, 90 and 110 spatial samples, respectively. In addition, we suggest a new approach to obtain sensor locations that are optimal with respect to prior assumptions. The approach also allows to control, e.g., the uniformity of the sensor locations. Based on our simulations, we argue that for a low number of spatial samples, model-informed non-uniform sampling can be beneficial. For a large number of samples, uniform sampling grids yield nearly the same total information as the model-informed grids.

Prefrontal Theta-Phase Synchronized Brain Stimulation With Real-Time EEG-Triggered TMS.

BACKGROUND: Theta-band neuronal oscillations in the prefrontal cortex are associated with several cognitive functions. Oscillatory phase is an important correlate of excitability and phase synchrony mediates information transfer between neuronal populations oscillating at that frequency. The ability to extract and exploit the prefrontal theta rhythm in real time in humans would facilitate insight into neurophysiological mechanisms of cognitive processes involving the prefrontal cortex, and development of brain-state-dependent stimulation for therapeutic applications. OBJECTIVES: We investigate individual source-space beamforming-based estimation of the prefrontal theta oscillation as a method to target specific phases of the ongoing theta oscillations in the human dorsomedial prefrontal cortex (DMPFC) with real-time EEG-triggered transcranial magnetic stimulation (TMS). Different spatial filters for extracting the prefrontal theta oscillation from EEG signals are compared and additional signal quality criteria are assessed to take into account the dynamics of this cortical oscillation. METHODS: Twenty two healthy participants were recruited for anatomical MRI scans and EEG recordings with 18 composing the final analysis. We calculated individual spatial filters based on EEG beamforming in source space. The extracted EEG signal was then used to simulate real-time phase-detection and quantify the accuracy as compared to post-hoc phase estimates. Different spatial filters and triggering parameters were compared. Finally, we validated the feasibility of this approach by actual real-time triggering of TMS pulses at different phases of the prefrontal theta oscillation. RESULTS: Higher phase-detection accuracy was achieved using individualized source-based spatial filters, as compared to an average or standard Laplacian filter, and also by detecting and avoiding periods of low theta amplitude and periods containing a phase reset. Using optimized parameters, prefrontal theta-phase synchronized TMS of DMPFC was achieved with an accuracy of ±55°. CONCLUSION: This study demonstrates the feasibility of triggering TMS pulses during different phases of the ongoing prefrontal theta oscillation in real time. This method is relevant for brain state-dependent stimulation in human studies of cognition. It will also enable new personalized therapeutic repetitive TMS protocols for more effective treatment of neuropsychiatric disorders.

Individual head models for estimating the TMS-induced electric field in rat brain.

In transcranial magnetic stimulation (TMS), the initial cortical activation due to stimulation is determined by the state of the brain and the magnitude, waveform, and direction of the induced electric field (E-field) in the cortex. The E-field distribution depends on the conductivity geometry of the head. The effects of deviations from a spherically symmetric conductivity profile have been studied in detail in humans. In small mammals, such as rats, these effects are more pronounced due to their less spherical head, proportionally much thicker neck region, and overall much smaller size compared to the TMS coils. In this study, we describe a simple method for building individual realistically shaped head models for rats from high-resolution X-ray tomography images. We computed the TMS-induced E-field with the boundary element method and assessed the effect of head-model simplifications on the estimated E-field. The deviations from spherical symmetry have large, non-trivial effects on the E-field distribution: for some coil orientations, the strongest stimulation is in the brainstem even when the coil is over the motor cortex. With modelling prior to an experiment, such problematic coil orientations can be avoided for more accurate targeting.

Concurrent electrophysiological and hemodynamic measurements of evoked neural oscillations in human visual cortex using sparsely interleaved fast fMRI and EEG.

Electroencephalography (EEG) concurrently collected with functional magnetic resonance imaging (fMRI) is heavily distorted by the repetitive gradient coil switching during the fMRI acquisition. The performance of the typical template-based gradient artifact suppression method can be suboptimal because the artifact changes over time. Gradient artifact residuals also impede the subsequent suppression of ballistocardiography artifacts. Here we propose recording continuous EEG with temporally sparse fast fMRI (fast fMRI-EEG) to minimize the EEG artifacts caused by MRI gradient coil switching without significantly compromising the field-of-view and spatiotemporal resolution of fMRI. Using simultaneous multi-slice inverse imaging to achieve whole-brain fMRI with isotropic 5-mm resolution in 0.1 â€‹s, and performing these acquisitions once every 2 â€‹s, we have 95% of the duty cycle available to record EEG with substantially less gradient artifact. We found that the standard deviation of EEG signals over the entire acquisition period in fast fMRI-EEG was reduced to 54% of that in conventional concurrent echo-planar imaging (EPI) and EEG recordings (EPI-EEG) across participants. When measuring 15-Hz steady-state visual evoked potentials (SSVEPs), the baseline-normalized oscillatory neural response in fast fMRI-EEG was 2.5-fold of that in EPI-EEG. The functional MRI responses associated with the SSVEP delineated by EPI and fast fMRI were similar in the spatial distribution, the elicited waveform, and detection power. Sparsely interleaved fast fMRI-EEG provides high-quality EEG without substantially compromising the quality of fMRI in evoked response measurements, and has the potential utility for applications where the onset of the target stimulus cannot be precisely determined, such as epilepsy.

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.

Short-interval intracortical inhibition in human primary motor cortex: A multi-locus transcranial magnetic stimulation study.

Short-interval intracortical inhibition (SICI) has been studied with paired-pulse transcranial magnetic stimulation (TMS) by administering two pulses at a millisecond-scale interstimulus interval (ISI) to a single cortical target. It has, however, been difficult to study the interaction of nearby cortical targets with paired-pulse TMS. To overcome this limitation, we have developed a multi-locus TMS (mTMS) device, which allows controlling the stimulus location electronically. Here, we applied mTMS to study SICI in primary motor cortex with paired pulses targeted to adjacent locations, aiming to quantify the extent of the cortical region producing SICI in the location of a test stimulus. We varied the location and timing of the conditioning stimulus with respect to a test stimulus targeted to the cortical hotspot of the abductor pollicis brevis (APB) in order to study their effects on motor evoked potentials. We further applied a two-coil protocol with the conditioning stimulus given by an oval coil only to the surroundings of the APB hotspot, to which a subsequent test stimulus was administered with a figure-of-eight coil. The strongest SICI occurred at ISIs below 1 ms and at ISIs around 2.5 ms. These ISIs increased when the conditioning stimulus receded from the APB hotspot. Our two-coil paired-pulse TMS study suggests that SICI at ISIs of 0.5 and 2.5 ms originate from different mechanisms or neuronal elements.

Real-time computation of the TMS-induced electric field in a realistic head model.

Transcranial magnetic stimulation (TMS) is often targeted using a model of TMS-induced electric field (E). In such navigated TMS, the E-field models have been based on spherical approximation of the head. Such models omit the effects of cerebrospinal fluid (CSF) and gyral folding, leading to potentially large errors in the computed E-field. So far, realistic models have been too slow for interactive TMS navigation. We present computational methods that enable real-time solving of the E-field in a realistic five-compartment (5-C) head model that contains isotropic white matter, gray matter, CSF, skull and scalp. Using reciprocity and Geselowitz integral equation, we separate the computations to coil-dependent and -independent parts. For the Geselowitz integrals, we present a fast numerical quadrature. Further, we present a moment-matching approach for optimizing dipole-based coil models. We verified and benchmarked the new methods using simulations with over 100 coil locations. The new quadrature introduced a relative error (RE) of 0.3-0.6%. For a coil model with 42 dipoles, the total RE of the quadrature and coil model was 0.44-0.72%. Taking also other model errors into account, the contribution of the new approximations to the RE was 0.1%. For comparison, the RE due to omitting the separation of white and gray matter was >11%, and the RE due to omitting also the CSF was >23%. After the coil-independent part of the model has been built, E-fields can be computed very quickly: Using a standard PC and basic GPU, our solver computed the full E-field in a 5-C model in 9000 points on the cortex in 27 coil positions per second (cps). When the separation of white and gray matter was omitted, the speed was 43-65 cps. Solving only one component of the E-field tripled the speed. The presented methods enable real-time solving of the TMS-induced E-field in a realistic head model that contains the CSF and gyral folding. The new methodology allows more accurate targeting and precise adjustment of stimulation intensity during experimental or clinical TMS mapping.

The impact of improved MEG-MRI co-registration on MEG connectivity analysis.

Co-registration between structural head images and functional MEG data is needed for anatomically-informed MEG data analysis. Despite the efforts to minimize the co-registration error, conventional landmark- and surface-based strategies for co-registering head and MEG device coordinates achieve an accuracy of typically 5-10 mm. Recent advances in instrumentation and technical solutions, such as the development of hybrid ultra-low-field (ULF) MRI-MEG devices or the use of 3D-printed individualized foam head-casts, promise unprecedented co-registration accuracy, i.e., 2 mm or better. In the present study, we assess through simulations the impact of such an improved co-registration on MEG connectivity analysis. We generated synthetic MEG recordings for pairs of connected cortical sources with variable locations. We then assessed the capability to reconstruct source-level connectivity from these recordings for 0-15-mm co-registration error, three levels of head modeling detail (one-, three- and four-compartment models), two source estimation techniques (linearly constrained minimum-variance beamforming and minimum-norm estimation MNE) and five separate connectivity metrics (imaginary coherency, phase-locking value, amplitude-envelope correlation, phase-slope index and frequency-domain Granger causality). We found that beamforming can better take advantage of an accurate co-registration than MNE. Specifically, when the co-registration error was smaller than 3 mm, the relative error in connectivity estimates was down to one-third of that observed with typical co-registration errors. MNE provided stable results for a wide range of co-registration errors, while the performance of beamforming rapidly degraded as the co-registration error increased. Furthermore, we found that even moderate co-registration errors (>6 mm, on average) essentially decrease the difference of four- and three- or one-compartment models. Hence, a precise co-registration is important if one wants to take full advantage of highly accurate head models for connectivity analysis. We conclude that an improved co-registration will be beneficial for reliable connectivity analysis and effective use of highly accurate head models in future MEG connectivity studies.

Individual Activation Patterns After the Stimulation of Different Motor Areas: A Transcranial Magnetic Stimulation-Electroencephalography Study.

The combination of transcranial magnetic stimulation (TMS) and electroencephalography (EEG) enables one to study effective connectivity and activation order in neuronal networks. To characterize effective connectivity originating from the primary motor cortex (M1), dorsal premotor area (PMd), and supplementary motor area (SMA). Three right-handed volunteers (two men, aged 25-30 years) participated in a navigated TMS-EEG experiment. M1, PMd, and SMA over the nondominant hemisphere were stimulated with 150 TMS pulses. Minimum-norm estimates were derived from the EEG data to estimate the spatial spreading of TMS-elicited neuronal activation on an individual level. The activation order of the cortical areas varied depending on the stimulated area. There were similarities and differences in the spatial distribution of the TMS-evoked potentials between subjects. Similarities in cortical activation patterns were seen at short poststimulus latencies and the differences at long latencies. This pilot study suggests that cortical activation patterns and the activation order of motor areas differ interindividually and depend on the stimulated motor area. It further indicates that TMS-activated effective connections or underlying structural connections vary between subjects. The spatial patterns of TMS-evoked potentials differ between subjects especially at long latencies, when probably more complex neuronal networks are active.

Requirements for Coregistration Accuracy in On-Scalp MEG.

Recent advances in magnetic sensing has made on-scalp magnetoencephalography (MEG) possible. In particular, optically-pumped magnetometers (OPMs) have reached sensitivity levels that enable their use in MEG. In contrast to the SQUID sensors used in current MEG systems, OPMs do not require cryogenic cooling and can thus be placed within millimetres from the head, enabling the construction of sensor arrays that conform to the shape of an individual's head. To properly estimate the location of neural sources within the brain, one must accurately know the position and orientation of sensors in relation to the head. With the adaptable on-scalp MEG sensor arrays, this coregistration becomes more challenging than in current SQUID-based MEG systems that use rigid sensor arrays. Here, we used simulations to quantify how accurately one needs to know the position and orientation of sensors in an on-scalp MEG system. The effects that different types of localisation errors have on forward modelling and source estimates obtained by minimum-norm estimation, dipole fitting, and beamforming are detailed. We found that sensor position errors generally have a larger effect than orientation errors and that these errors affect the localisation accuracy of superficial sources the most. To obtain similar or higher accuracy than with current SQUID-based MEG systems, RMS sensor position and orientation errors should be [Formula: see text] and [Formula: see text], respectively.

Truncated RAP-MUSIC (TRAP-MUSIC) for MEG and EEG source localization.

Electrically active brain regions can be located applying MUltiple SIgnal Classification (MUSIC) on magneto- or electroencephalographic (MEG; EEG) data. We introduce a new MUSIC method, called truncated recursively-applied-and-projected MUSIC (TRAP-MUSIC). It corrects a hidden deficiency of the conventional RAP-MUSIC algorithm, which prevents estimation of the true number of brain-signal sources accurately. The correction is done by applying a sequential dimension reduction to the signal-subspace projection. We show that TRAP-MUSIC significantly improves the performance of MUSIC-type localization; in particular, it successfully and robustly locates active brain regions and estimates their number. We compare TRAP-MUSIC and RAP-MUSIC in simulations with varying key parameters, e.g., signal-to-noise ratio, correlation between source time-courses, and initial estimate for the dimension of the signal space. In addition, we validate TRAP-MUSIC with measured MEG data. We suggest that with the proposed TRAP-MUSIC method, MUSIC-type localization could become more reliable and suitable for various online and offline MEG and EEG applications.