Source Time Courses Research Articles

EEG dynamic source imaging using a regularized optimization with spatio-temporal constraints.

One of the most important needs in neuroimaging is brain dynamic source imaging with high spatial and temporal resolution. EEG source imaging estimates the underlying sources from EEG recordings, which provides enhanced spatial resolution with intrinsically high temporal resolution. To ensure identifiability in the underdetermined source reconstruction problem, constraints on EEG sources are essential. This paper introduces a novel method for estimating source activities based on spatio-temporal constraints and a dynamic source imaging algorithm. The method enhances time resolution by incorporating temporal evolution of neural activity into a regularization function. Additionally, two spatial regularization constraints based on and norms are applied in the transformed domain to address both focal and spread neural activities, achieved through spatial gradient and Laplacian transform. Performance evaluation, conducted quantitatively using synthetic datasets, discusses the influence of parameters such as source extent, number of sources, correlation level, and SNR level on temporal and spatial metrics. Results demonstrate that the proposed method provides superior spatial and temporal reconstructions compared to state-of-the-art inverse solutions including STRAPS, sLORETA, SBL, dSPM, and MxNE. This improvement is attributed to the simultaneous integration of transformed spatial and temporal constraints. When applied to a real auditory ERP dataset, our algorithm accurately reconstructs brain source time series and locations, effectively identifying the origins of auditory evoked potentials. In conclusion, our proposed method with spatio-temporal constraints outperforms the state-of-the-art algorithms in estimating source distribution and time courses.

Resting state connectivity biomarkers of seizure freedom after epilepsy surgery

Alterations in brain networks may cause the lowering of the seizure threshold and hypersynchronization that underlie the recurrence of unprovoked seizures in epilepsy. The aim of this work is to estimate functional network characteristics, which may help predicting outcome of epilepsy surgery.Twenty patients were studied (11 females, 9 males, mean age 33 years) with scalp-recorded HD-EEG in resting state (eyes closed, no interictal discharges) before intracranial evaluation, which allowed the precise determination of the epileptogenic zone. Dipole source time courses in the brain were estimated using Weighted Minimum Norm Estimate based on HD-EEG signals. Information inflow and outflow of atlas-based brain regions were computed using partial directed connectivity. A set of graph measures for pairwise connections in standard EEG frequency bands was calculated.After epilepsy surgery 10 patients were seizure-free (Engel 1a) and 10 patients continued suffering from seizures (Engel outcome worse than 1a). Inflow of the regions containing the epileptogenic zone in the beta and delta frequency bands was significantly lower in patients who achieved seizure-freedom after surgery, compared with patients who continued to have seizures (p = 0.012, and p = 0.026, respectively). Average path length in the beta frequency band was significantly higher in patients who achieved seizure freedom (p = 0.012). In the delta frequency band, local efficiency and clustering coefficient were significantly higher in patients who achieved seizure freedom (0.033, 0.046).In patients who achieved seizure freedom after surgery, the preoperative analysis of the epileptic network exhibited stronger separation of the region containing the seizure onset zone, with less inflow of information. In contrast, shorter paths within the epileptic network may facilitate hypersynchronous neuronal activity and thus the recurrence of seizures in non-seizure free patients. This study supports the hypothesis that epileptic network properties might help to define suitable candidates for epilepsy surgery.

FMRI-SI-STBF: An fMRI-informed Bayesian electromagnetic spatio-temporal extended source imaging

Multimodal functional neuroimaging by integrating functional magnetic resonance imaging (fMRI) and electroencephalography (EEG) has the promise of recovering brain activities with high spatiotemporal resolution, which is crucial for neuroscience research and clinical diagnosis. However, the misalignment of the localizations between fMRI and EEG activities may degrade the accuracy of the fMRI-constrained EEG source imaging (ESI) technique. To leverage the complementary spatiotemporal resolution of fMRI and EEG in a data-driven fashion, we propose an asymmetric approach for EEG/fMRI fusion, termed fMRI-informed source imaging based on spatiotemporal basis functions (fMRI-SI-STBF). fMRI-SI-STBF employs the covariance components (CCs) derived from clusters defined by fMRI and EEG signals as spatial priors within the empirical Bayesian framework. Additionally, fMRI-SI-STBF represents the current source matrix as a linear combination of several unknown temporal basis functions (TBFs) by matrix decomposition. The relative contribution of each of the fMRI-informed and EEG-informed CCs, as well as the number and profiles of the TBFs, are all automatically determined based on the EEG data using variational Bayesian inference. Our results demonstrate that fMRI-SI-STBF can effectively utilize valid fMRI information for ESI and is robust to invalid fMRI priors. This robustness is essential for practical ESI since the validity of fMRI priors is often unclear considering that fMRI is an indirect measure of neural activity. Moreover, fMRI-SI-STBF can achieve performance improvement by incorporating temporal constraints compared to methods that use spatial constraints only. For the numerical simulations, fMRI-SI-STBF reconstructs the source extents, locations and time courses more accurately than existing EEG-fMRI ESI methods (i.e., fwMNE, fMRI-SI-SBF) and ESI methods without fMRI priors (i.e., wMNE, LORETA, SBL, SI-STBF, SI-SBF), indicated by the smaller spatial dispersion (average SD <5 mm), distance of localization error (average DLE <2 mm), shape error (average SE <0.9) and larger model evidence values.

Clinical Validation of the Champagne Algorithm for Evoked Response Source Localization in Magnetoencephalography

Magnetoencephalography (MEG) is a robust method for non-invasive functional brain mapping of sensory cortices due to its exceptional spatial and temporal resolution. The clinical standard for MEG source localization of functional landmarks from sensory evoked responses is the equivalent current dipole (ECD) localization algorithm, known to be sensitive to initialization, noise, and manual choice of the number of dipoles. Recently many automated and robust algorithms have been developed, including the Champagne algorithm, an empirical Bayesian algorithm, with powerful abilities for MEG source reconstruction and time course estimation (Wipf et al. 2010; Owen et al. 2012). Here, we evaluate automated Champagne performance in a clinical population of tumor patients where there was minimal failure in localizing sensory evoked responses using the clinical standard, ECD localization algorithm. MEG data of auditory evoked potentials and somatosensory evoked potentials from 21 brain tumor patients were analyzed using Champagne, and these results were compared with equivalent current dipole (ECD) fit. Across both somatosensory and auditory evoked field localization, we found there was a strong agreement between Champagne and ECD localizations in all cases. Given resolution of 8mm voxel size, peak source localizations from Champagne were below 10mm of ECD peak source localization. The Champagne algorithm provides a robust and automated alternative to manual ECD fits for clinical localization of sensory evoked potentials and can contribute to improved clinical MEG data processing workflows.

Unisensory and Multisensory Responses in Fetal Alcohol Spectrum Disorders (FASD): Effects of Spatial Congruence

While it is generally accepted that structural and functional brain deficits underlie the behavioral deficits associated with Fetal Alcohol Spectrum Disorders (FASD), the degree to which these problems are expressed in sensory pathology is unknown. Electrophysiological measures indicate that neural processing is delayed in visual and auditory domains. Furthermore, multiple reports of white matter deficits due to prenatal alcohol exposure indicate altered cortical connectivity in individuals with FASD. Multisensory integration requires close coordination between disparate cortical areas leading us to hypothesize that individuals with FASD will have impaired multisensory integration relative to healthy control (HC) participants. Participants’ neurophysiological responses were recorded using magnetoencephalography (MEG) during passive unisensory or simultaneous, spatially congruent or incongruent multisensory auditory and somatosensory stimuli. Source timecourses from evoked responses were estimated using multi-dipole spatiotemporal modeling. Auditory M100 response latency was faster for the multisensory relative to the unisensory condition but no group differences were observed. M200 auditory latency to congruent stimuli was earlier and congruent amplitude was larger in participants with FASD relative to controls. Somatosensory M100 response latency was faster in right hemisphere for multisensory relative to unisensory stimulation in both groups. FASD participants’ somatosensory M200 responses were delayed by 13 ms, but only for the unisensory presentation of the somatosensory stimulus. M200 results indicate that unisensory and multisensory processing is altered in FASD; it remains to be seen if the multisensory response represents a normalization of the unisensory deficits.

Hierarchical multiscale Bayesian algorithm for robust MEG/EEG source reconstruction

In this paper, we present a novel hierarchical multiscale Bayesian algorithm for electromagnetic brain imaging using magnetoencephalography (MEG) and electroencephalography (EEG). In particular, we present a solution to the source reconstruction problem for sources that vary in spatial extent. We define sensor data measurements using a generative probabilistic graphical model that is hierarchical across spatial scales of brain regions and voxels. We then derive a novel Bayesian algorithm for probabilistic inference with this graphical model. This algorithm enables robust reconstruction of sources that have different spatial extent, from spatially contiguous clusters of dipoles to isolated dipolar sources. We compare the new algorithm with several representative benchmarks on both simulated and real brain activities. The source locations and the correct estimation of source time courses used for the simulated data are chosen to test the performance on challenging source configurations. In simulations, performance of the novel algorithm shows superiority to several existing benchmark algorithms. We also demonstrate that the new algorithm is more robust to correlated brain activity present in real MEG and EEG data and is able to resolve distinct and functionally relevant brain areas with real MEG and EEG datasets.

A Comparative Study of the Robustness of Frequency-Domain Connectivity Measures to Finite Data Length.

In this work we use numerical simulation to investigate how the temporal length of the data affects the reliability of the estimates of brain connectivity from EEG time-series. We assume that the neural sources follow a stable MultiVariate AutoRegressive model, and consider three connectivity metrics: imaginary part of coherency (IC), generalized partial directed coherence (gPDC) and frequency-domain granger causality (fGC). In order to assess the statistical significance of the estimated values, we use the surrogate data test by generating phase-randomized and autoregressive surrogate data. We first consider the ideal case where we know the source time courses exactly. Here we show how, expectedly, even exact knowledge of the source time courses is not sufficient to provide reliable estimates of the connectivity when the number of samples gets small; however, while gPDC and fGC tend to provide a larger number of false positives, the IC becomes less sensitive to the presence of connectivity. Then we proceed with more realistic simulations, where the source time courses are estimated using eLORETA, and the EEG signal is affected by biological noise of increasing intensity. Using the ideal case as a reference, we show that the impact of biological noise on IC estimates is qualitatively different from the impact on gPDC and fGC.

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.

Information-Theoretic Evidence for Predictive Coding in the Face-Processing System.

Predictive coding suggests that the brain infers the causes of its sensations by combining sensory evidence with internal predictions based on available prior knowledge. However, the neurophysiological correlates of (pre)activated prior knowledge serving these predictions are still unknown. Based on the idea that such preactivated prior knowledge must be maintained until needed, we measured the amount of maintained information in neural signals via the active information storage (AIS) measure. AIS was calculated on whole-brain beamformer-reconstructed source time courses from MEG recordings of 52 human subjects during the baseline of a Mooney face/house detection task. Preactivation of prior knowledge for faces showed as α-band-related and β-band-related AIS increases in content-specific areas; these AIS increases were behaviorally relevant in the brain's fusiform face area. Further, AIS allowed decoding of the cued category on a trial-by-trial basis. Our results support accounts indicating that activated prior knowledge and the corresponding predictions are signaled in low-frequency activity (<30 Hz).SIGNIFICANCE STATEMENT Our perception is not only determined by the information our eyes/retina and other sensory organs receive from the outside world, but strongly depends also on information already present in our brains, such as prior knowledge about specific situations or objects. A currently popular theory in neuroscience, predictive coding theory, suggests that this prior knowledge is used by the brain to form internal predictions about upcoming sensory information. However, neurophysiological evidence for this hypothesis is rare, mostly because this kind of evidence requires strong a priori assumptions about the specific predictions the brain makes and the brain areas involved. Using a novel, assumption-free approach, we find that face-related prior knowledge and the derived predictions are represented in low-frequency brain activity.

Resting-State Magnetoencephalography Reveals Different Patterns of Aberrant Functional Connectivity in Combat-Related Mild Traumatic Brain Injury.

Blast mild traumatic brain injury (mTBI) is a leading cause of sustained impairment in military service members and veterans. However, the mechanism of persistent disability is not fully understood. The present study investigated disturbances in brain functioning in mTBI participants using a source-imaging-based approach to analyze functional connectivity (FC) from resting-state magnetoencephalography (rs-MEG). Study participants included 26 active-duty service members or veterans who had blast mTBI with persistent post-concussive symptoms, and 22 healthy control active-duty service members or veterans. The source time courses from regions of interest (ROIs) were used to compute ROI to whole-brain (ROI-global) FC for different frequency bands using two different measures: 1) time-lagged cross-correlation and 2) phase-lock synchrony. Compared with the controls, blast mTBI participants showed increased ROI-global FC in beta, gamma, and low-frequency bands, but not in the alpha band. Sources of abnormally increased FC included the: 1) prefrontal cortex (right ventromedial prefrontal cortex [vmPFC], right rostral anterior cingulate cortex [rACC]), and left ventrolateral and dorsolateral prefrontal cortex; 2) medial temporal lobe (bilateral parahippocampus, hippocampus, and amygdala); and 3) right putamen and cerebellum. In contrast, the blast mTBI group also showed decreased FC of the right frontal pole. Group differences were highly consistent across the two different FC measures. FC of the left ventrolateral prefrontal cortex correlated with executive functioning and processing speed in mTBI participants. Altogether, our findings of increased and decreased regionalpatterns of FC suggest that disturbances in intrinsic brain connectivity may be the result of multiple mechanisms, and are associated with cognitive sequelae of the injury.

Noninvasive Electromagnetic Source Imaging and Granger Causality Analysis: An Electrophysiological Connectome (eConnectome) Approach.

Combined source-imaging techniques and directional connectivity analysis can provide useful information about the underlying brain networks in a noninvasive fashion. Source-imaging techniques have been used successfully to either determine the source of activity or to extract source time-courses for Granger causality analysis, previously. In this work, we utilize source-imaging algorithms to both find the network nodes [regions of interest (ROI)] and then extract the activation time series for further Granger causality analysis. The aim of this work is to find network nodes objectively from noninvasive electromagnetic signals, extract activation time-courses, and apply Granger analysis on the extracted series to study brain networks under realistic conditions. Source-imaging methods are used to identify network nodes and extract time-courses and then Granger causality analysis is applied to delineate the directional functional connectivity of underlying brain networks. Computer simulations studies where the underlying network (nodes and connectivity pattern) is known were performed; additionally, this approach has been evaluated in partial epilepsy patients to study epilepsy networks from interictal and ictal signals recorded by EEG and/or Magnetoencephalography (MEG). Localization errors of network nodes are less than 5 mm and normalized connectivity errors of ∼20% in estimating underlying brain networks in simulation studies. Additionally, two focal epilepsy patients were studied and the identified nodes driving the epileptic network were concordant with clinical findings from intracranial recordings or surgical resection. Our study indicates that combined source-imaging algorithms with Granger causality analysis can identify underlying networks precisely (both in terms of network nodes location and internodal connectivity). The combined source imaging and Granger analysis technique is an effective tool for studying normal or pathological brain conditions.

EEG and MEG: sensitivity to epileptic spike activity as function of source orientation and depth

Simultaneous electroencephalography (EEG) and magnetoencephalography (MEG) recordings of neuronal activity from epileptic patients reveal situations in which either EEG or MEG or both modalities show visible interictal spikes. While different signal-to-noise ratios (SNRs) of the spikes in EEG and MEG have been reported, a quantitative relation of spike source orientation and depth as well as the background brain activity to the SNR has not been established.We investigated this quantitative relationship for both dipole and patch sources in an anatomically realistic cortex model.Altogether, 5600 dipole and 3300 patch sources were distributed on the segmented cortical surfaces of two volunteers. The sources were classified according to their quantified depths and orientations, ranging from 20 mm to 60 mm below the skin surface and radial and tangential, respectively. The source time-courses mimicked an interictal spike, and the simulated background activity emulated resting activity. Simulations were conducted with individual three-compartment boundary element models. The SNR was evaluated for 128 EEG, 102 MEG magnetometer, and 204 MEG gradiometer channels.For superficial dipole and superficial patch sources, EEG showed higher SNRs for dominantly radial orientations, and MEG showed higher values for dominantly tangential orientations. Gradiometers provided higher SNR than magnetometers for superficial sources, particularly for those with dominantly tangential orientations. The orientation dependent difference in SNR in EEG and MEG gradually changed as the sources were located deeper, where the interictal spikes generated higher SNRs in EEG compared to those in MEG for all source orientations. With deep sources, the SNRs in gradiometers and magnetometers were of the same order.To better detect spikes, both EEG and MEG should be used.

Closely Spaced MEG Source Localization and Functional Connectivity Analysis Using a New Prewhitening Invariance of Noise Space Algorithm.

This paper proposed a prewhitening invariance of noise space (PW-INN) as a new magnetoencephalography (MEG) source analysis method, which is particularly suitable for localizing closely spaced and highly correlated cortical sources under real MEG noise. Conventional source localization methods, such as sLORETA and beamformer, cannot distinguish closely spaced cortical sources, especially under strong intersource correlation. Our previous work proposed an invariance of noise space (INN) method to resolve closely spaced sources, but its performance is seriously degraded under correlated noise between MEG sensors. The proposed PW-INN method largely mitigates the adverse influence of correlated MEG noise by projecting MEG data to a new space defined by the orthogonal complement of dominant eigenvectors of correlated MEG noise. Simulation results showed that PW-INN is superior to INN, sLORETA, and beamformer in terms of localization accuracy for closely spaced and highly correlated sources. Lastly, source connectivity between closely spaced sources can be satisfactorily constructed from source time courses estimated by PW-INN but not from results of other conventional methods. Therefore, the proposed PW-INN method is a promising MEG source analysis to provide a high spatial-temporal characterization of cortical activity and connectivity, which is crucial for basic and clinical research of neural plasticity.

Minimum variance beamformer weights revisited

Adaptive minimum variance beamformers are widely used analysis tools in MEG and EEG. When the target brain activity presents in the form of spatially localized responses, the procedure usually involves two steps. First, positions and orientations of the sources of interest are determined. Second, the filter weights are calculated and source time courses reconstructed. This last step is the object of the current study. Despite different approaches utilized at the source localization stage, basic expressions for the weights have the same form, dictated by the minimum variance condition. These classic expressions involve covariance matrix of the measured field, which includes contributions from both the sources of interest and the noise background. We show analytically that the same weights can alternatively be obtained, if the full field covariance is replaced with that of the noise, provided the beamformer points to the true sources precisely. In practice, however, a certain mismatch is always inevitable. We show that such mismatch results in partial suppression of the true sources if the traditional weights are used. To avoid this effect, the “alternative” weights based on properly estimated noise covariance should be applied at the second, source time course reconstruction step. We demonstrate mathematically and using simulated and real data that in many situations the alternative weights provide significantly better time course reconstruction quality than the traditional ones. In particular, they a) improve source-level SNR and yield more accurately reconstructed waveforms; b) provide more accurate estimates of inter-source correlations; and c) reduce the adverse influence of the source correlations on the performance of single-source beamformers, which are used most often. Importantly, the alternative weights come at no additional computational cost, as the structure of the expressions remains the same.

A symmetric multivariate leakage correction for MEG connectomes

Ambiguities in the source reconstruction of magnetoencephalographic (MEG) measurements can cause spurious correlations between estimated source time-courses. In this paper, we propose a symmetric orthogonalisation method to correct for these artificial correlations between a set of multiple regions of interest (ROIs). This process enables the straightforward application of network modelling methods, including partial correlation or multivariate autoregressive modelling, to infer connectomes, or functional networks, from the corrected ROIs. Here, we apply the correction to simulated MEG recordings of simple networks and to a resting-state dataset collected from eight subjects, before computing the partial correlations between power envelopes of the corrected ROItime-courses. We show accurate reconstruction of our simulated networks, and in the analysis of real MEGresting-state connectivity, we find dense bilateral connections within the motor and visual networks, together with longer-range direct fronto-parietal connections.

Influence of the head model on EEG and MEG source connectivity analyses

The results of brain connectivity analysis using reconstructed source time courses derived from EEG and MEG data depend on a number of algorithmic choices. While previous studies have investigated the influence of the choice of source estimation method or connectivity measure, the effects of the head modeling errors or simplifications have not been studied sufficiently.In the present simulation study, we investigated the influence of particular properties of the head model on the reconstructed source time courses as well as on source connectivity analysis in EEG and MEG. Therefore, we constructed a realistic head model and applied the finite element method to solve the EEG and MEG forward problems. We considered the distinction between white and gray matter, the distinction between compact and spongy bone, the inclusion of a cerebrospinal fluid (CSF) compartment, and the reduction to a simple 3-layer model comprising only the skin, skull, and brain. Source time courses were reconstructed using a beamforming approach and the source connectivity was estimated by the imaginary coherence (ICoh) and the generalized partial directed coherence (GPDC).Our results show that in both EEG and MEG, neglecting the white and gray matter distinction or the CSF causes considerable errors in reconstructed source time courses and connectivity analysis, while the distinction between spongy and compact bone is just of minor relevance, provided that an adequate skull conductivity value is used. Large inverse and connectivity errors are found in the same regions that show large topography errors in the forward solution. Moreover, we demonstrate that the very conservative ICoh is relatively safe from the crosstalk effects caused by imperfect head models, as opposed to the GPDC.

Eccentric and concentric loading of the triceps surae: an in vivo study of dynamic muscle and tendon biomechanical parameters.

Triceps surae eccentric exercise is more effective than concentric exercise for treating Achilles tendinopathy, however the mechanisms underpinning these effects are unclear. This study compared the biomechanical characteristics of eccentric and concentric exercises to identify differences in the tendon load response. Eleven healthy volunteers performed eccentric and concentric exercises on a force plate, with ultrasonography, motion tracking, and EMG applied to measure Achilles tendon force, lower limb movement, and leg muscle activation. Tendon length was ultrasonographically tracked and quantified using a novel algorithm. The Fourier transform of the ground reaction force was also calculated to investigate for tremor, or perturbations. Tendon stiffness and extension did not vary between exercise types (P = .43). However, tendon perturbations were significantly higher during eccentric than concentric exercises (25%-40% higher, P = .02). Furthermore, perturbations during eccentric exercises were found to be negatively correlated with the tendon stiffness (R2 = .59). The particular efficacy of eccentric exercise does not appear to result from variation in tendon stiffness or extension within a given session. However, varied perturbation magnitude may have a role in mediating the observed clinical effects. This property is subject-specific, with the source and clinical time-course of such perturbations requiring further research.

Voxel-ICA for reconstruction of source signal time-series and orientation in EEG and MEG

In electroencephalography (EEG) and magnetoencephalography signal processing, scalar beamformers are a popular technique for reconstruction of the time-course of a brain source in a single time-series. A prerequisite for scalar beamformers, however, is that the orientation of the source must be known or estimated, whereas in reality the orientation of a brain source is often not known in advance and current techniques for estimation of brain source orientation are effective only for high signal-to-noise ratio (SNR) brain sources. As a result, vector beamformers are applied which do not need the orientation of the source and reconstruct the source time-course in three orthogonal (x, y, and z) directions. To obtain a single time-course, the vector magnitude of the three orthogonal outputs of the beamformer can be calculated at each time point (often called neural activity index, NAI). The NAI, however, is different from the actual time-course of a source since it contains only positive values. Moreover, in estimating the magnitude of the desired source, the background activity (undesired signals) in the beamformer outputs also become all positive values, which, when added to each other, leads to a drop in the SNR. This becomes a serious problem when the desired source is weak. We propose applying independent component analysis (ICA) to the orthogonal time-courses of a brain voxel, as reconstructed by a vector beamformer, to reconstruct the time-course of a desired source in a single time-series. This approach also provides a good estimation of dipole orientation. Simulated and real EEG data were used to demonstrate the performance of voxel-ICA and were compared with a scalar beamformer and the magnitude time-series of a vector beamformer. This approach is especially helpful when the desired source is weak and the orientation of the source cannot be estimated by other means.

Head movement compensation in real-time magnetoencephalographic recordings

Neurofeedback- and brain-computer interface (BCI)-based interventions can be implemented using real-time analysis of magnetoencephalographic (MEG) recordings. Head movement during MEG recordings, however, can lead to inaccurate estimates of brain activity, reducing the efficacy of the intervention. Most real-time applications in MEG have utilized analyses that do not correct for head movement. Effective means of correcting for head movement are needed to optimize the use of MEG in such applications. Here we provide preliminary validation of a novel analysis technique, real-time source estimation (rtSE), that measures head movement and generates corrected current source time course estimates in real-time. rtSE was applied while recording a calibrated phantom to determine phantom position localization accuracy and source amplitude estimation accuracy under stationary and moving conditions. Results were compared to off-line analysis methods to assess validity of the rtSE technique. The rtSE method allowed for accurate estimation of current source activity at the source-level in real-time, and accounted for movement of the source due to changes in phantom position. The rtSE technique requires modifications and specialized analysis of the following MEG work flow steps.•Data acquisition•Head position estimation•Source localization•Real-time source estimationThis work explains the technical details and validates each of these steps.

MEG source imaging method using fast L1 minimum-norm and its applications to signals with brain noise and human resting-state source amplitude images

The present study developed a fast MEG source imaging technique based on Fast Vector-based Spatio-Temporal Analysis using a L1-minimum-norm (Fast-VESTAL) and then used the method to obtain the source amplitude images of resting-state magnetoencephalography (MEG) signals for different frequency bands. The Fast-VESTAL technique consists of two steps. First, L1-minimum-norm MEG source images were obtained for the dominant spatial modes of sensor-waveform covariance matrix. Next, accurate source time-courses with millisecond temporal resolution were obtained using an inverse operator constructed from the spatial source images of Step 1. Using simulations, Fast-VESTAL's performance was assessed for its 1) ability to localize multiple correlated sources; 2) ability to faithfully recover source time-courses; 3) robustness to different SNR conditions including SNR with negative dB levels; 4) capability to handle correlated brain noise; and 5) statistical maps of MEG source images. An objective pre-whitening method was also developed and integrated with Fast-VESTAL to remove correlated brain noise. Fast-VESTAL's performance was then examined in the analysis of human median-nerve MEG responses. The results demonstrated that this method easily distinguished sources in the entire somatosensory network. Next, Fast-VESTAL was applied to obtain the first whole-head MEG source-amplitude images from resting-state signals in 41 healthy control subjects, for all standard frequency bands. Comparisons between resting-state MEG sources images and known neurophysiology were provided. Additionally, in simulations and cases with MEG human responses, the results obtained from using conventional beamformer technique were compared with those from Fast-VESTAL, which highlighted the beamformer's problems of signal leaking and distorted source time-courses.