MEG/EEG source reconstruction, statistical evaluation, and visualization with NUTMEG

Neural correlates of motor learning: Network communication versus local oscillations

Learning new motor skills through training, also termed motor learning, is central for everyday life. Current training strategies recommend intensive task-repetitions aimed at inducing local activation of motor areas, associated with changes in oscillation amplitudes ("event-related power") during training. More recently, another neural mechanism was suggested to influence motor learning: modulation of functional connectivity (FC), that is, how much spatially separated brain regions communicate with each other before and during training. The goal of the present study was to compare the impact of these two neural processing types on motor learning. We measured EEG before, during, and after a finger-tapping task (FTT) in 20 healthy subjects. The results showed that training gain, long-term expertise (i.e., average motor performance), and consolidation were all predicted by whole-brain alpha- and beta-band FC at motor areas, striatum, and mediotemporal lobe (MTL). Local power changes during training did not predict any dependent variable. Thus, network dynamics seem more crucial than local activity for motor sequence learning, and training techniques should attempt to facilitate network interactions rather than local cortical activation.

Neural mechanisms underlying improved new-word learning with high-density transcranial direct current stimulation

Neurobehavioral studies have provided evidence for the effectiveness of anodal tDCS on language production, by stimulation of the left Inferior Frontal Gyrus (IFG) or of left Temporo-Parietal Junction (TPJ). However, tDCS is currently not used in clinical practice outside of trials, because behavioral effects have been inconsistent and underlying neural effects unclear. Here, we propose to elucidate the neural correlates of verb and noun learning and to determine if they can be modulated with anodal high-definition (HD) tDCS stimulation. Thirty-six neurotypical participants were randomly allocated to anodal HD-tDCS over either the left IFG, the left TPJ, or sham stimulation. On day one, participants performed a naming task (pre-test). On day two, participants underwent a new-word learning task with rare nouns and verbs concurrently to HD-tDCS for 20 min. The third day consisted of a post-test of naming performance. EEG was recorded at rest and during naming on each day. Verb learning was significantly facilitated by left IFG stimulation. HD-tDCS over the left IFG enhanced functional connectivity between the left IFG and TPJ and this correlated with improved learning. HD-tDCS over the left TPJ enabled stronger local activation of the stimulated area (as indexed by greater alpha and beta-band power decrease) during naming, but this did not translate into better learning. Thus, tDCS can induce local activation or modulation of network interactions. Only the enhancement of network interactions, but not the increase in local activation, leads to robust improvement of word learning. This emphasizes the need to develop new neuromodulation methods influencing network interactions. Our study suggests that this may be achieved through behavioral activation of one area and concomitant activation of another area with HD-tDCS.

Abnormal gamma phase-amplitude coupling in the parahippocampal cortex is associated with network hyperexcitability in Alzheimer's disease

While animal models of Alzheimer's disease (AD) have shown altered gamma oscillations (∼40 Hz) in local neural circuits, the low signal-to-noise ratio of gamma in the resting human brain precludes its quantification via conventional spectral estimates. Phase-amplitude coupling (PAC) indicating the dynamic integration between the gamma amplitude and the phase of low-frequency (4-12 Hz) oscillations is a useful alternative to capture local gamma activity. In addition, PAC is also an index of neuronal excitability as the phase of low-frequency oscillations that modulate gamma amplitude, effectively regulates the excitability of local neuronal firing. In this study, we sought to examine the local neuronal activity and excitability using gamma PAC, within brain regions vulnerable to early AD pathophysiology-entorhinal cortex and parahippocampus, in a clinical population of patients with AD and age-matched controls. Our clinical cohorts consisted of a well-characterized cohort of AD patients (n = 50; age, 60 ± 8 years) with positive AD biomarkers, and age-matched, cognitively unimpaired controls (n = 35; age, 63 ± 5.8 years). We identified the presence or the absence of epileptiform activity in AD patients (AD patients with epileptiform activity, AD-EPI+, n = 20; AD patients without epileptiform activity, AD-EPI-, n = 30) using long-term electroencephalography (LTM-EEG) and 1-hour long magnetoencephalography (MEG) with simultaneous EEG. Using the source reconstructed MEG data, we computed gamma PAC as the coupling between amplitude of the gamma frequency (30-40 Hz) with phase of the theta (4-8 Hz) and alpha (8-12 Hz) frequency oscillations, within entorhinal and parahippocampal cortices. We found that patients with AD have reduced gamma PAC in the left parahippocampal cortex, compared to age-matched controls. Furthermore, AD-EPI+ patients showed greater reductions in gamma PAC than AD-EPI- in bilateral parahippocampal cortices. In contrast, entorhinal cortices did not show gamma PAC abnormalities in patients with AD. Our findings demonstrate the spatial patterns of altered gamma oscillations indicating possible region-specific manifestations of network hyperexcitability within medial temporal lobe regions vulnerable to AD pathophysiology. Greater deficits in AD-EPI+ suggests that reduced gamma PAC is a sensitive index of network hyperexcitability in AD patients. Collectively, the current results emphasize the importance of investigating the role of neural circuit hyperexcitability in early AD pathophysiology and explore its potential as a modifiable contributor to AD pathobiology.

Impaired long-range excitatory time scale predicts abnormal neural oscillations and cognitive deficits in Alzheimer's disease

BACKGROUND

Alzheimer's disease (AD) is the most common form of dementia, progressively impairing cognitive abilities. While neuroimaging studies have revealed functional abnormalities in AD, how these relate to aberrant neuronal circuit mechanisms remains unclear. Using magnetoencephalography imaging we documented abnormal local neural synchrony patterns in patients with AD. To identify global abnormal biophysical mechanisms underlying the spatial and spectral electrophysiological patterns in AD, we estimated the parameters of a biophysical spectral graph model (SGM).

METHODS

SGM is an analytic neural mass model that describes how long-range fiber projections in the brain mediate the excitatory and inhibitory activity of local neuronal subpopulations. Unlike other coupled neuronal mass models, the SGM is linear, available in closed-form, and parameterized by a small set of biophysical interpretable global parameters. This facilitates their rapid and unambiguous inference which we performed here on a well-characterized clinical population of patients with AD (N = 88, age = 62.73 +/- 8.64 years) and a cohort of age-matched controls (N = 88, age = 65.07 +/- 9.92 years).

RESULTS

Patients with AD showed significantly elevated long-range excitatory neuronal time scales, local excitatory neuronal time scales and local inhibitory neural synaptic strength. The long-range excitatory time scale had a larger effect size, compared to local excitatory time scale and inhibitory synaptic strength and contributed highest for the accurate classification of patients with AD from controls. Furthermore, increased long-range time scale was associated with greater deficits in global cognition.

CONCLUSIONS

These results demonstrate that long-range excitatory time scale of neuronal activity, despite being a global measure, is a key determinant in the local spectral signatures and cognition in the human brain, and how it might be a parsimonious factor underlying altered neuronal activity in AD. Our findings provide new insights into mechanistic links between abnormal local spectral signatures and global connectivity measures in AD.

Global MEG Resting State Functional Connectivity in Children with Autism and Sensory Processing Dysfunction

Sensory processing dysfunction not only affects most individuals with autism spectrum disorder (ASD), but at least 5% of children without ASD also experience dysfunctional sensory processing. Our understanding of the relationship between sensory dysfunction and resting state brain activity is still emerging. This study compared long-range resting state functional connectivity of neural oscillatory behavior in children aged 8-12 years with autism spectrum disorder (ASD; N=18), those with sensory processing dysfunction (SPD; N=18) who do not meet ASD criteria, and typically developing control participants (TDC; N=24) using magnetoencephalography (MEG). Functional connectivity analyses were performed in the alpha and beta frequency bands, which are known to be implicated in sensory information processing. Group differences in functional connectivity and associations between sensory abilities and functional connectivity were examined. Distinct patterns of functional connectivity differences between ASD and SPD groups were found only in the beta band, but not in the alpha band. In both alpha and beta bands, ASD and SPD cohorts differed from the TDC cohort. Somatosensory cortical beta-band functional connectivity was associated with tactile processing abilities, while higher-order auditory cortical alpha-band functional connectivity was associated with auditory processing abilities. These findings demonstrate distinct long-range neural synchrony alterations in SPD and ASD that are associated with sensory processing abilities. Neural synchrony measures could serve as potential sensitive biomarkers for ASD and SPD.

Setting up EEG Source Imaging in Practice

SUMMARY

Adding EEG source imaging to a clinical practice has clear advantages over visual inspection of EEG. This article offers insight on incorporating EEG source imaging into an EEG laboratory and the best practices for producing optimal source analysis results.

Impaired long-range excitatory time scale predicts abnormal neural oscillations and cognitive deficits in Alzheimer's disease

Alzheimer's disease (AD) is the most common form of dementia, progressively impairing memory and cognition. While neuroimaging studies have revealed functional abnormalities in AD, how these relate to aberrant neuronal circuit mechanisms remains unclear. Using magnetoencephalography imaging we documented abnormal local neural synchrony patterns in patients with AD. To identify abnormal biophysical mechanisms underlying these abnormal electrophysiological patterns, we estimated the parameters of a spectral graph-theory model (SGM). SGM is an analytic model that describes how long-range fiber projections in the brain mediate the excitatory and inhibitory activity of local neuronal subpopulations. The long-range excitatory time scale was associated with greater deficits in global cognition and was able to distinguish AD patients from controls with high accuracy. These results demonstrate that long-range excitatory time scale of neuronal activity, despite being a global measure, is a key determinant in the spatiospectral signatures and cognition in AD.

Dynamic functional connectivity MEG features of Alzheimer's disease

Dynamic resting state functional connectivity (RSFC) characterizes time-varying fluctuations of functional brain network activity. While many studies have investigated static functional connectivity, it has been unclear whether features of dynamic functional connectivity are associated with neurodegenerative diseases. Popular sliding-window and clustering methods for extracting dynamic RSFC have various limitations that prevent extracting reliable features to address this question. Here, we use a novel and robust time-varying dynamic network (TVDN) approach to extract the dynamic RSFC features from high resolution magnetoencephalography (MEG) data of participants with Alzheimer's disease (AD) and matched controls. The TVDN algorithm automatically and adaptively learns the low-dimensional spatiotemporal manifold of dynamic RSFC and detects dynamic state transitions in data. We show that amongst all the functional features we investigated, the dynamic manifold features are the most predictive of AD. These include: the temporal complexity of the brain network, given by the number of state transitions and their dwell times, and the spatial complexity of the brain network, given by the number of eigenmodes. These dynamic features have higher sensitivity and specificity in distinguishing AD from healthy subjects than the existing benchmarks do. Intriguingly, we found that AD patients generally have higher spatial complexity but lower temporal complexity compared with healthy controls. We also show that graph theoretic metrics of dynamic component of TVDN are significantly different in AD versus controls, while static graph metrics are not statistically different. These results indicate that dynamic RSFC features are impacted in neurodegenerative disease like Alzheimer's disease, and may be crucial to understanding the pathophysiological trajectory of these diseases.

Impaired long-range excitatory time scale predicts abnormal neural oscillations and cognitive deficits in Alzheimer's disease

Alzheimer's disease (AD) is the most common form of dementia, progressively impairing memory and cognition. While neuroimaging studies have revealed functional abnormalities in AD, how these relate to aberrant neuronal circuit mechanisms remains unclear. Using magnetoencephalography imaging we documented abnormal local neural synchrony patterns in patients with AD. To identify abnormal biophysical mechanisms underlying these abnormal electrophysiological patterns, we estimated the parameters of a spectral graph-theory model (SGM). SGM is an analytic model that describes how long-range fiber projections in the brain mediate the excitatory and inhibitory activity of local neuronal subpopulations. The long-range excitatory time scale was associated with greater deficits in global cognition and was able to distinguish AD patients from controls with high accuracy. These results demonstrate that long-range excitatory time scale of neuronal activity, despite being a global measure, is a key determinant in the spatiospectral signatures and cognition in AD.

Measurement of Frontal Midline Theta Oscillations using OPM-MEG

Optically pumped magnetometers (OPMs) are an emerging lightweight and compact sensor that can measure magnetic fields generated by the human brain. OPMs enable construction of wearable magnetoencephalography (MEG) systems, which offer advantages over conventional instrumentation. However, when trying to measure signals at low frequency, higher levels of inherent sensor noise, magnetic interference and movement artefact introduce a significant challenge. Accurate characterisation of low frequency brain signals is important for neuroscientific, clinical, and paediatric MEG applications and consequently, demonstrating the viability of OPMs in this area is critical. Here, we undertake measurement of theta band (4-8 Hz) neural oscillations and contrast a newly developed 174 channel triaxial wearable OPM-MEG system with conventional (cryogenic-MEG) instrumentation. Our results show that visual steady state responses at 4 Hz, 6 Hz and 8 Hz can be recorded using OPM-MEG with a signal-to-noise ratio (SNR) that is not significantly different to conventional MEG. Moreover, we measure frontal midline theta oscillations during a 2-back working memory task, again demonstrating comparable SNR for both systems. We show that individual differences in both the amplitude and spatial signature of induced frontal-midline theta responses are maintained across systems. Finally, we show that our OPM-MEG results could not have been achieved without a triaxial sensor array, or the use of postprocessing techniques. Our results demonstrate the viability of OPMs for characterising theta oscillations and add weight to the argument that OPMs can replace cryogenic sensors as the fundamental building block of MEG systems.

EEG Beta functional connectivity decrease in the left amygdala correlates with the affective pain in fibromyalgia: A pilot study

Fibromyalgia (FM) is a major chronic pain disease with prominent affective disturbances, and pain-associated changes in neurotransmitters activity and in brain connectivity. However, correlates of affective pain dimension lack. The primary goal of this correlational cross-sectional case-control pilot study was to find electrophysiological correlates of the affective pain component in FM. We examined the resting-state EEG spectral power and imaginary coherence in the beta (β) band (supposedly indexing the GABAergic neurotransmission) in 16 female patients with FM and 11 age-adjusted female controls. FM patients displayed lower functional connectivity in the High β (Hβ, 20-30 Hz) sub-band than controls (p = 0.039) in the left basolateral complex of the amygdala (p = 0.039) within the left mesiotemporal area, in particular, in correlation with a higher affective pain component level (r = 0.50, p = 0.049). Patients showed higher Low β (Lβ, 13-20 Hz) relative power than controls in the left prefrontal cortex (p = 0.001), correlated with ongoing pain intensity (r = 0.54, p = 0.032). For the first time, GABA-related connectivity changes correlated with the affective pain component are shown in the amygdala, a region highly involved in the affective regulation of pain. The β power increase in the prefrontal cortex could be compensatory to pain-related GABAergic dysfunction.

Brainstorm-DUNEuro: An integrated and user-friendly Finite Element Method for modeling electromagnetic brain activity

Human brain activity generates scalp potentials (electroencephalography - EEG), intracranial potentials (iEEG), and external magnetic fields (magnetoencephalography - MEG). These electrophysiology (e-phys) signals can often be measured simultaneously for research and clinical applications. The forward problem involves modeling these signals at their sensors for a given equivalent current dipole configuration within the brain. While earlier researchers modeled the head as a simple set of isotropic spheres, today's magnetic resonance imaging (MRI) data allow for a detailed anatomic description of brain structures and anisotropic characterization of tissue conductivities. We present a complete pipeline, integrated into the Brainstorm software, that allows users to automatically generate an individual and accurate head model based on the subject's MRI and calculate the electromagnetic forward solution using the finite element method (FEM). The head model generation is performed by integrating the latest tools for MRI segmentation and FEM mesh generation. The final head model comprises the five main compartments: white-matter, gray-matter, CSF, skull, and scalp. The anisotropic brain conductivity model is based on the effective medium approach (EMA), which estimates anisotropic conductivity tensors from diffusion-weighted imaging (DWI) data. The FEM electromagnetic forward solution is obtained through the DUNEuro library, integrated into Brainstorm, and accessible with either a user-friendly graphical interface or scripting. With tutorials and example data sets available in an open-source format on the Brainstorm website, this integrated pipeline provides access to advanced FEM tools for electromagnetic modeling to a broader neuroscience community.

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.

Flexible adjustment of anticipations in human outcome processing

To sense whether thoughts refer to current reality or not, a capacity called orbitofrontal reality filtering, depends on an orbitofrontal signal when anticipated outcomes fail to occur. Here, we explored the flexibility and precision of outcome processing in a deterministic reversal learning task. Healthy subjects decided which one of two colored squares hid a target stimulus. Brain activity was measured with high-density electroencephalography. Stimuli resembling, but not identical with, the target stimuli were initially processed like different stimuli from 210 to 250 ms, irrespective of behavioral relevance. From 250 ms on, they were processed according to behavioral relevance: If they required a subsequent switch, they were processed like different stimuli; if they had been declared potential targets, they were treated like true targets. Stimuli requiring a behavioral switch induced strong theta activity in orbitofrontal, ventromedial, and medial temporal regions. The study indicates flexible adaptation of anticipations but precise processing of outcomes, mainly determined by behavioral relevance.

Altered excitatory and inhibitory neuronal subpopulation parameters are distinctly associated with tau and amyloid in Alzheimer's disease

Background

Neuronal- and circuit-level abnormalities of excitation and inhibition are shown to be associated with tau and amyloid-beta (Aβ) in preclinical models of Alzheimer's disease (AD). These relationships remain poorly understood in patients with AD.

Methods

Using empirical spectra from magnetoencephalography and computational modeling (neural mass model), we examined excitatory and inhibitory parameters of neuronal subpopulations and investigated their specific associations to regional tau and Aβ, measured by positron emission tomography, in patients with AD.

Results

Patients with AD showed abnormal excitatory and inhibitory time-constants and neural gains compared to age-matched controls. Increased excitatory time-constants distinctly correlated with higher tau depositions while increased inhibitory time-constants distinctly correlated with higher Aβ depositions.

Conclusions

Our results provide critical insights about potential mechanistic links between abnormal neural oscillations and cellular correlates of impaired excitatory and inhibitory synaptic functions associated with tau and Aβ in patients with AD.

Funding

This study was supported by the National Institutes of Health grants: K08AG058749 (KGR), F32AG050434-01A1 (KGR), K23 AG038357 (KAV), P50 AG023501, P01 AG19724 (BLM), P50-AG023501 (BLM and GDR), R01 AG045611 (GDR); AG034570, AG062542 (WJ); NS100440 (SSN), DC176960 (SSN), DC017091 (SSN), AG062196 (SSN); a grant from John Douglas French Alzheimer's Foundation (KAV); grants from Larry L. Hillblom Foundation: 2015-A-034-FEL (KGR), 2019-A-013-SUP (KGR); grants from the Alzheimer's Association: AARG-21-849773 (KGR); PCTRB-13-288476 (KAV), and made possible by Part the CloudTM (ETAC-09-133596); a grant from Tau Consortium (GDR and WJJ), and a gift from the S. D. Bechtel Jr. Foundation.

Neuronal synchrony abnormalities associated with subclinical epileptiform activity in early-onset Alzheimer's disease

Since the first demonstrations of network hyperexcitability in scientific models of Alzheimer's disease, a growing body of clinical studies have identified subclinical epileptiform activity and associated cognitive decline in patients with Alzheimer's disease. An obvious problem presented in these studies is lack of sensitive measures to detect and quantify network hyperexcitability in human subjects. In this study we examined whether altered neuronal synchrony can be a surrogate marker to quantify network hyperexcitability in patients with Alzheimer's disease. Using magnetoencephalography (MEG) at rest, we studied 30 Alzheimer's disease patients without subclinical epileptiform activity, 20 Alzheimer's disease patients with subclinical epileptiform activity and 35 age-matched controls. Presence of subclinical epileptiform activity was assessed in patients with Alzheimer's disease by long-term video-EEG and a 1-h resting MEG with simultaneous EEG. Using the resting-state source-space reconstructed MEG signal, in patients and controls we computed the global imaginary coherence in alpha (8-12 Hz) and delta-theta (2-8 Hz) oscillatory frequencies. We found that Alzheimer's disease patients with subclinical epileptiform activity have greater reductions in alpha imaginary coherence and greater enhancements in delta-theta imaginary coherence than Alzheimer's disease patients without subclinical epileptiform activity, and that these changes can distinguish between Alzheimer's disease patients with subclinical epileptiform activity and Alzheimer's disease patients without subclinical epileptiform activity with high accuracy. Finally, a principal component regression analysis showed that the variance of frequency-specific neuronal synchrony predicts longitudinal changes in Mini-Mental State Examination in patients and controls. Our results demonstrate that quantitative neurophysiological measures are sensitive biomarkers of network hyperexcitability and can be used to improve diagnosis and to select appropriate patients for the right therapy in the next-generation clinical trials. The current results provide an integrative framework for investigating network hyperexcitability and network dysfunction together with cognitive and clinical correlates in patients with Alzheimer's disease.

Temporal specificity of abnormal neural oscillations during phonatory events in laryngeal dystonia

Laryngeal dystonia is a debilitating disorder of voicing in which the laryngeal muscles are intermittently in spasm resulting in involuntary interruptions during speech. The central pathophysiology of laryngeal dystonia, underlying computational impairments in vocal motor control, remains poorly understood. Although prior imaging studies have found aberrant activity in the CNS during phonation in patients with laryngeal dystonia, it is not known at what timepoints during phonation these abnormalities emerge and what function may be impaired. To investigate this question, we recruited 22 adductor laryngeal dystonia patients (15 female, age range = 28.83-72.46 years) and 18 controls (eight female, age range = 27.40-71.34 years). We leveraged the fine temporal resolution of magnetoencephalography to monitor neural activity around glottal movement onset, subsequent voice onset and after the onset of pitch feedback perturbations. We examined event-related beta-band (12-30 Hz) and high-gamma-band (65-150 Hz) neural oscillations. Prior to glottal movement onset, we observed abnormal frontoparietal motor preparatory activity. After glottal movement onset, we observed abnormal activity in the somatosensory cortex persisting through voice onset. Prior to voice onset and continuing after, we also observed abnormal activity in the auditory cortex and the cerebellum. After pitch feedback perturbation onset, we observed no differences between controls and patients in their behavioural responses to the perturbation. But in patients, we did find abnormal activity in brain regions thought to be involved in the auditory feedback control of vocal pitch (premotor, motor, somatosensory and auditory cortices). Our study results confirm the abnormal processing of somatosensory feedback that has been seen in other studies. However, there were several remarkable findings in our study. First, patients have impaired vocal motor activity even before glottal movement onset, suggesting abnormal movement preparation. These results are significant because (i) they occur before movement onset, abnormalities in patients cannot be ascribed to deficits in vocal performance and (ii) they show that neural abnormalities in laryngeal dystonia are more than just abnormal responses to sensory feedback during phonation as has been hypothesized in some previous studies. Second, abnormal auditory cortical activity in patients begins even before voice onset, suggesting abnormalities in setting up auditory predictions before the arrival of auditory feedback at voice onset. Generally, activation abnormalities identified in key brain regions within the speech motor network around various phonation events not only provide temporal specificity to neuroimaging phenotypes in laryngeal dystonia but also may serve as potential therapeutic targets for neuromodulation.

A unified view on beamformers for M/EEG source reconstruction

Beamforming is a popular method for functional source reconstruction using magnetoencephalography (MEG) and electroencephalography (EEG) data. Beamformers, which were first proposed for MEG more than two decades ago, have since been applied in hundreds of studies, demonstrating that they are a versatile and robust tool for neuroscience. However, certain characteristics of beamformers remain somewhat elusive and there currently does not exist a unified documentation of the mathematical underpinnings and computational subtleties of beamformers as implemented in the most widely used academic open source software packages for MEG analysis (Brainstorm, FieldTrip, MNE, and SPM). Here, we provide such documentation that aims at providing the mathematical background of beamforming and unifying the terminology. Beamformer implementations are compared across toolboxes and pitfalls of beamforming analyses are discussed. Specifically, we provide details on handling rank deficient covariance matrices, prewhitening, the rank reduction of forward fields, and on the combination of heterogeneous sensor types, such as magnetometers and gradiometers. The overall aim of this paper is to contribute to contemporary efforts towards higher levels of computational transparency in functional neuroimaging.

Feasibility of Reconstructing Source Functional Connectivity with Low-Density EEG

OBJECTIVES

Functional connectivity (FC) is increasingly used as target for neuromodulation and enhancement of performance. A reliable assessment of FC with electroencephalography (EEG) currently requires a laboratory environment with high-density montages and a long preparation time. This study investigated the feasibility of reconstructing source FC with a low-density EEG montage towards a usage in real life applications.

METHODS

Source FC was reconstructed with inverse solutions and quantified as node degree of absolute imaginary coherence in alpha frequencies. We used simulated coherent point sources as well as two real datasets to investigate the impact of electrode density (19 vs. 128 electrodes) and usage of template vs. individual MRI-based head models on localization accuracy. In addition, we checked whether low-density EEG is able to capture inter-individual variations in coherence strength.

RESULTS

In numerical simulations as well as real data, a reduction of the number of electrodes led to less reliable reconstructions of coherent sources and of coupling strength. Yet, when comparing different approaches to reconstructing FC from 19 electrodes, source FC obtained with beamformers outperformed sensor FC, FC computed after independent component analysis, and source FC obtained with sLORETA. In particular, only source FC based on beamformers was able to capture neural correlates of motor behavior.

CONCLUSION

Reconstructions of FC from low-density EEG is challenging, but may be feasible when using source reconstructions with beamformers.

Disrupted resting-state EEG alpha-band interactions as a novel marker for the severity of visual field deficits after brain lesion

OBJECTIVE

Homonymous visual field deficits (HFVDs) are frequent following brain lesions. Current restoration treatments aim at activating areas of residual vision through numerous stimuli, but show limited effect. Recent findings suggest that spontaneous neural α-band coupling is more efficient for enabling visual perception in healthy humans than task-induced activations. Here, we evaluated whether it is also associated with the severity of HFVD.

METHODS

Ten patients with HFVDs after brain damage in the subacute to chronic stage and ten matched healthy controls underwent visual stimulation with alternating checkerboards and electroencephalography recordings of stimulation-induced power changes and of spontaneous neural interactions during rest.

RESULTS

Visual areas of the affected hemisphere showed reduced event-related power decrease in α and β frequency bands, but also reduced spontaneous α-band interactions during rest, as compared to contralesional areas and healthy controls. A multivariate stepwise regression retained the degree of disruption of spontaneous interactions, but not the reduced task-induced power changes as predictor for the severity of the visual deficit.

CONCLUSIONS

Spontaneous α-band interactions of visual areas appear as a better marker for the severity of HFVDs than task-induced activations.

SIGNIFICANCE

Treatment attempts of HFVDs should try to enhance spontaneous α-band coupling of structurally intact ipsilesional areas.

Unification of sparse Bayesian learning algorithms for electromagnetic brain imaging with the majorization minimization framework

Methods for electro- or magnetoencephalography (EEG/MEG) based brain source imaging (BSI) using sparse Bayesian learning (SBL) have been demonstrated to achieve excellent performance in situations with low numbers of distinct active sources, such as event-related designs. This paper extends the theory and practice of SBL in three important ways. First, we reformulate three existing SBL algorithms under the majorization-minimization (MM) framework. This unification perspective not only provides a useful theoretical framework for comparing different algorithms in terms of their convergence behavior, but also provides a principled recipe for constructing novel algorithms with specific properties by designing appropriate bounds of the Bayesian marginal likelihood function. Second, building on the MM principle, we propose a novel method called LowSNR-BSI that achieves favorable source reconstruction performance in low signal-to-noise-ratio (SNR) settings. Third, precise knowledge of the noise level is a crucial requirement for accurate source reconstruction. Here we present a novel principled technique to accurately learn the noise variance from the data either jointly within the source reconstruction procedure or using one of two proposed cross-validation strategies. Empirically, we could show that the monotonous convergence behavior predicted from MM theory is confirmed in numerical experiments. Using simulations, we further demonstrate the advantage of LowSNR-BSI over conventional SBL in low-SNR regimes, and the advantage of learned noise levels over estimates derived from baseline data. To demonstrate the usefulness of our novel approach, we show neurophysiologically plausible source reconstructions on averaged auditory evoked potential data.

Neural dynamics of semantic categorization in semantic variant of primary progressive aphasia

Semantic representations are processed along a posterior-to-anterior gradient reflecting a shift from perceptual (e.g., it has eight legs) to conceptual (e.g., venomous spiders are rare) information. One critical region is the anterior temporal lobe (ATL): patients with semantic variant primary progressive aphasia (svPPA), a clinical syndrome associated with ATL neurodegeneration, manifest a deep loss of semantic knowledge. We test the hypothesis that svPPA patients perform semantic tasks by over-recruiting areas implicated in perceptual processing. We compared MEG recordings of svPPA patients and healthy controls during a categorization task. While behavioral performance did not differ, svPPA patients showed indications of greater activation over bilateral occipital cortices and superior temporal gyrus, and inconsistent engagement of frontal regions. These findings suggest a pervasive reorganization of brain networks in response to ATL neurodegeneration: the loss of this critical hub leads to a dysregulated (semantic) control system, and defective semantic representations are seemingly compensated via enhanced perceptual processing.

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.

Interhemispheric Auditory Cortical Synchronization in Asymmetric Hearing Loss

Objectives:

Auditory cortical activation of the two hemispheres to monaurally presented tonal stimuli has been shown to be asynchronous in normal hearing (NH) but synchronous in the extreme case of adult-onset asymmetric hearing loss (AHL) with single-sided deafness. We addressed the wide knowledge gap between these two anchoring states of interhemispheric temporal organization. The objectives of this study were as follows: (1) to map the trajectory of interhemispheric temporal reorganization from asynchrony to synchrony using magnitude of interaural threshold difference as the independent variable in a cross-sectional study and (2) to evaluate reversibility of interhemispheric synchrony in association with hearing in noise performance by amplifying the aidable poorer ear in a repeated measures, longitudinal study.

Design:

The cross-sectional and longitudinal cohorts were comprised of 49 subjects (AHL; N = 21; 11 male, 10 female; mean age = 48 years) and NH (N = 28; 16 male, 12 female; mean age = 45 years). The maximum interaural threshold difference of the two cohorts spanned from 0 to 65 dB. Magnetoencephalography analyses focused on latency of the M100 peak response from auditory cortex in both hemispheres between 50 msec and 150 msec following monaural tonal stimulation at the frequency (0.5, 1, 2, 3, or 4 kHz) corresponding to the maximum and minimum interaural threshold difference for better and poorer ears separately. The longitudinal AHL cohort was drawn from three subjects in the cross-sectional AHL cohort (all male; ages 49 to 60 years; varied AHL etiologies; no amplification for at least 2 years). All longitudinal study subjects were treated by monaural amplification of the poorer ear and underwent repeated measures examination of the M100 response latency and quick speech in noise hearing in noise performance at baseline, and postamplification months 3, 6, and 12.

Results:

The M100 response peak latency values in the ipsilateral hemisphere lagged those in the contralateral hemisphere for all stimulation conditions. The mean (SD) interhemispheric latency difference values (ipsilateral less contralateral) to better ear stimulation for three categories of maximum interaural threshold difference were as follows: NH (≤ 10 dB)—8.6 (3.0) msec; AHL (15 to 40 dB)—3.0 (1.2) msec; AHL (≥ 45 dB)—1.4 (1.3) msec. In turn, the magnitude of difference values were used to define interhemispheric temporal organization states of asynchrony, mixed asynchrony and synchrony, and synchrony, respectively. Amplification of the poorer ear in longitudinal subjects drove interhemispheric organization change from baseline synchrony to postamplification asynchrony and hearing in noise performance improvement in those with baseline impairment over a 12-month period.

Conclusions:

Interhemispheric temporal organization in AHL was anchored between states of asynchrony in NH and synchrony in single-sided deafness. For asymmetry magnitudes between 15 and 40 dB, the intermediate mixed state of asynchrony and synchrony was continuous and reversible. Amplification of the poorer ear in AHL improved hearing in noise performance and restored normal temporal organization of auditory cortices in the two hemispheres. The return to normal interhemispheric asynchrony from baseline synchrony and improvement in hearing following monoaural amplification of the poorer ear evolved progressively over a 12-month period.

Spontaneous Network Coupling Enables Efficient Task Performance without Local Task-Induced Activations

Neurobehavioral studies in humans have long concentrated on changes in local activity levels during repetitive executions of a task. Spontaneous neural coupling within extended networks has latterly been found to also influence performance. Here, we intend to uncover the underlying mechanisms, the relative importance, and the interaction between spontaneous coupling and task-induced activations. To do so, we recorded two groups of healthy participants (male and female) during rest and while they performed either a visual perception or a motor sequence task. We demonstrate that, for both tasks, stronger activations during the task as well as greater network coupling through spontaneous α rhythms at rest predict performance. However, high performers present an absence of classical task-induced activations and, instead, stronger spontaneous network coupling. Activations were thus a compensation mechanism needed only in subjects with lower spontaneous network interactions. This challenges classical models of neural processing and calls for new strategies in attempts to train and enhance performance.SIGNIFICANCE STATEMENT Our findings challenge the widely accepted notion that task-induced activations are of paramount importance for behavior. This will have an important impact on interpretations of human neurobehavioral research. They further link the widely used techniques of quantifying network communication in the brain with classical neuroscience methods and demonstrate possible ways of how network communication influences human behavior. Traditional training methods attempt to enhance neural activations through task repetitions. Our findings suggest a more efficient neural target for learning: enhancing spontaneous neural interactions. This will be of major interest for a large variety of scientific fields with very broad applications in schools, work, and others.

NUTMEG: Open Source Software for M/EEG Source Reconstruction

Neurodynamic Utility Toolbox for Magnetoencephalo- and Electroencephalography (NUTMEG) is an open-source MATLAB-based toolbox for the analysis and reconstruction of magnetoencephalography/electroencephalography data in source space. NUTMEG includes a variety of options for the user in data import, preprocessing, source reconstruction, and functional connectivity. A group analysis toolbox allows the user to run a variety of inferential statistics on their data in an easy-to-use GUI-driven format. Importantly, NUTMEG features an interactive five-dimensional data visualization platform. A key feature of NUTMEG is the availability of a large menu of interference cancelation and source reconstruction algorithms. Each NUTMEG operation acts as a stand-alone MATLAB function, allowing the package to be easily adaptable and scripted for the more advanced user for interoperability with other software toolboxes. Therefore, NUTMEG enables a wide range of users access to a complete "sensor-to- source-statistics" analysis pipeline.

Resting-State SEEG May Help Localize Epileptogenic Brain Regions

BACKGROUND

Stereotactic electroencephalography (SEEG) is a minimally invasive neurosurgical method to localize epileptogenic brain regions in epilepsy but requires days in the hospital with interventions to trigger several seizures.

OBJECTIVE

To make initial progress in the development of network analysis methods to identify epileptogenic brain regions using brief, resting-state SEEG data segments, without requiring seizure recordings.

METHODS

In a cohort of 15 adult focal epilepsy patients undergoing SEEG, we evaluated functional connectivity (alpha-band imaginary coherence) across sampled regions using brief (2 min) resting-state data segments. Bootstrapped logistic regression was used to generate a model to predict epileptogenicity of individual regions.

RESULTS

Compared to nonepileptogenic structures, we found increased functional connectivity within epileptogenic regions (P < .05) and between epileptogenic areas and other structures (P < .01, paired t-tests, corrected). Epileptogenic areas also demonstrated higher clustering coefficient (P < .01) and betweenness centrality (P < .01), and greater decay of functional connectivity with distance (P < .05, paired t-tests, corrected). Our functional connectivity model to predict epileptogenicity of individual regions demonstrated an area under the curve of 0.78 and accuracy of 80.4%.

CONCLUSION

Our study represents a preliminary step towards defining resting-state SEEG functional connectivity patterns to help localize epileptogenic brain regions ahead of neurosurgical treatment without requiring seizure recordings.

The Visual Word Form Area compensates for auditory working memory dysfunction in schizophrenia

Auditory working memory impairments feature prominently in schizophrenia. However, the existence of altered and perhaps compensatory neural dynamics, sub-serving auditory working memory, remains largely unexplored. We compared the dynamics of induced high gamma power (iHGP) across cortex in humans during speech-sound working memory in individuals with schizophrenia (SZ) and healthy comparison subjects (HC) using magnetoencephalography (MEG). SZ showed similar task performance to HC while utilizing different brain regions. During encoding of speech sounds, SZ lacked the correlation of iHGP with task performance in posterior superior temporal gyrus (STGp) that was observed in healthy subjects. Instead, SZ recruited the visual word form area (VWFA) during both stimulus encoding and response preparation. Importantly, VWFA activity during encoding correlated with the magnitude of SZ hallucinations, task performance and an independent measure of verbal working memory. These findings suggest that VWFA plasticity is harnessed to compensate for STGp dysfunction in schizophrenia patients with hallucinations.

Global resting-state functional connectivity of neural oscillations in tinnitus with and without hearing loss

This study examined global resting-state functional connectivity of neural oscillations in individuals with chronic tinnitus and normal and impaired hearing. We tested the hypothesis that distinct neural oscillatory networks are engaged in tinnitus with and without hearing loss. In both tinnitus groups, with and without hearing loss, we identified multiple frequency band-dependent regions of increased and decreased global functional connectivity. We also found that the auditory domain of tinnitus severity, assayed by the Tinnitus Functional Index, was associated with global functional connectivity in both auditory and nonauditory regions. These findings provide candidate biomarkers to target and monitor treatments for tinnitus with and without hearing loss.

Neurophysiological signatures in Alzheimer's disease are distinctly associated with TAU, amyloid-β accumulation, and cognitive decline

Neural synchrony is intricately balanced in the normal resting brain but becomes altered in Alzheimer's disease (AD). To determine the neurophysiological manifestations associated with molecular biomarkers of AD neuropathology, in patients with AD, we used magnetoencephalographic imaging (MEGI) and positron emission tomography with amyloid-beta (Aβ) and TAU tracers. We found that alpha oscillations (8 to 12 Hz) were hyposynchronous in occipital and posterior temporoparietal cortices, whereas delta-theta oscillations (2 to 8 Hz) were hypersynchronous in frontal and anterior temporoparietal cortices, in patients with AD compared to age-matched controls. Regional patterns of alpha hyposynchrony were unique in each neurobehavioral phenotype of AD, whereas the regional patterns of delta-theta hypersynchrony were similar across the phenotypes. Alpha hyposynchrony strongly colocalized with TAU deposition and was modulated by the degree of TAU tracer uptake. In contrast, delta-theta hypersynchrony colocalized with both TAU and Aβ depositions and was modulated by both TAU and Aβ tracer uptake. Furthermore, alpha hyposynchrony but not delta-theta hypersynchrony was correlated with the degree of global cognitive dysfunction in patients with AD. The current study demonstrates frequency-specific neurophysiological signatures of AD pathophysiology and suggests that neurophysiological measures from MEGI are sensitive indices of network disruptions mediated by TAU and Aβ and associated cognitive decline. These findings facilitate the pursuit of novel therapeutic approaches toward normalizing network synchrony in AD.

Intervention-specific patterns of cortical function plasticity during auditory encoding in people with schizophrenia

Schizophrenia is a neurocognitive illness characterized by behavioral and neural impairments in both early auditory processing and higher order verbal working memory. Previously we have shown intervention-specific cognitive performance improvements with computerized, targeted training of auditory processing (AT) when compared to a computer games (CG) control intervention that emphasized visual processing. To investigate spatiotemporal changes in patterns of neural activity specific to the AT intervention, the current study used magnetoencephalography (MEG) imaging to derive induced high gamma band oscillations (HGO) during auditory encoding, before and after 50 h (∼10 weeks) of exposure to either the AT or CG intervention. During stimulus encoding, AT intervention-specific changes in high gamma activity occurred in left middle frontal and left middle-superior temporal cortices. In contrast, CG intervention-specific changes were observed in right medial frontal and supramarginal gyri during stimulus encoding, and in bilateral temporal cortices during response preparation. These data reveal that, in schizophrenia, intensive exposure to either training of auditory processing or exposure to visuospatial activities produces significant but complementary patterns of cortical function plasticity within a distributed fronto-temporal network. These results underscore the importance of delineating the specific neuroplastic effects of targeted behavioral interventions to ensure desired neurophysiological changes and avoid unintended consequences on neural system functioning.

Cross-modal plasticity in adult single-sided deafness revealed by alpha band resting-state functional connectivity

Single-sided deafness (SSD) or profound unilateral hearing loss is the condition where the transfer of acoustic information to the brain is restricted to one ear. SSD impairment is most evident under adverse acoustic environments with overlapping interference, which burdens cognitive resources. It is known that bilateral deafness induces cross-modal brain plasticity within visual cortical areas. Here we investigate whether similar cross-modal plasticity is observed in adult-onset SSD. In SSD patients (n ​= ​29) and matched controls (n ​= ​29) we estimated voxel level resting-state power and functional connectivity in the alpha band (8-12 ​Hz) from magnetoencephalography (MEG) data. We examined both global functional connectivity (mean functional connectivity of each voxel with the rest of the brain), and seeded functional connectivity of primary auditory cortices (A1), primary visual cortices (V1) and posterior cingulate cortex (PCC) of the default mode network (DMN). Power reduction was observed in left auditory cortex. Global functional connectivity showed reduction in frontal cortices and enhancement in visual cortex. Seeded functional connectivity of auditory cortices showed reduction in temporal, frontal and occipital regions, and enhancement in parietal cortex. Interestingly, seeded functional connectivity of visual cortices showed enhancement in visual cortices, inferior parietal lobe, post-central gyrus, and the precuneus, and reduction in auditory cortex. Seeded functional connectivity of PCC showed reduction in frontal cortical regions that are part of the DMN, attention, and working memory networks. Adult-onset SSD exhibited widespread cross-modal brain plasticity involving alterations in auditory, visual, attention, working memory and default mode networks.

Sensorimotor Cortical Oscillations during Movement Preparation in 16p11.2 Deletion Carriers

Sensorimotor deficits are prevalent in many neurodevelopmental disorders like autism, including one of its common genetic etiologies, a 600 kb reciprocal deletion/duplication at 16p11.2. We have previously shown that copy number variations of 16p11.2 impact regional brain volume, white matter integrity, and early sensory responses in auditory cortex. Here, we test the hypothesis that abnormal cortical neurophysiology is present when genes in the 16p11.2 region are haploinsufficient, and in humans that this in turn may account for behavioral deficits specific to deletion carriers. We examine sensorimotor cortical network activity in males and females with 16p11.2 deletions compared with both typically developing individuals, and those with duplications of 16p11.2, using magnetoencephalographic imaging during preparation of overt speech or hand movements in tasks designed to be easy for all participants. In deletion carriers, modulation of beta oscillations (12-30 Hz) were increased during both movement types over effector-specific regions of motor cortices compared with typically developing individuals or duplication carriers, with no task-related performance differences between cohorts, even when corrected for their own cognitive and sensorimotor deficits. Reduced left hemispheric language specialization was observed in deletion carriers but not in duplication carriers. Neural activity over sensorimotor cortices in deletion carriers was linearly related to clinical measures of speech and motor impairment. These findings link insufficient copy number repeats at 16p11.2 to excessive neural activity (e.g., increased beta oscillations) in motor cortical networks for speech and hand motor control. These results have significant implications for understanding the neural basis of autism and related neurodevelopmental disorders.SIGNIFICANCE STATEMENT The recurrent ∼600 kb deletion at 16p11.2 (BP4-BP5) is one of the most common genetic etiologies of ASD and, more generally, of neurodevelopmental disorders. Here, we use high-resolution magnetoencephalographic imaging (MEG-I) to define with millisecond precision the underlying neurophysiological signature of motor impairments for individuals with 16p11.2 deletions. We identify significant increases in beta (12-30 Hz) suppression in sensorimotor cortices related to performance during speech and hand movement tasks. These findings not only provide a neurophysiological phenotype for the clinical presentation of this genetic deletion, but also guide our understanding of how genetic variation encodes for neural oscillatory dynamics.

Beta-band activity in medial prefrontal cortex predicts source memory encoding and retrieval accuracy

Reality monitoring is defined as the ability to distinguish internally self-generated information from externally-derived information. The medial prefrontal cortex (mPFC) is a key brain region subserving reality monitoring and has been shown to be activated specifically during the retrieval of self-generated information. However, it is unclear if mPFC is activated during the encoding of self-generated information into memory. If so, it is important to understand whether successful retrieval of self-generated information critically depends on enhanced neural activity within mPFC during initial encoding of this self-generated information. We used magnetoencephalographic imaging (MEGI) to determine the timing and location of cortical activity during a reality-monitoring task involving self generated contextual source memory encoding and retrieval. We found both during encoding and retrieval of self-generated information, when compared to externally-derived information, mPFC showed significant task induced oscillatory power modulation in the beta-band. During initial encoding of self-generated information, greater mPFC beta-band power reductions occurred within a time window of −700 ms to −500 ms prior to vocalization. This increased activity in mPFC was not observed during encoding of externally-derived information. Additionally, increased mPFC activity during encoding of self-generated information predicted subsequent retrieval accuracy of this self-generated information. Beta-band activity in mPFC was also observed during the initial retrieval of self-generated information within a time window of 300 to 500 ms following stimulus onset and correlated with accurate retrieval performance of self-generated information. Together, these results further highlight the importance of mPFC in mediating the initial generation and awareness of participants’ internal thoughts.

Neural correlates of abnormal auditory feedback processing during speech production in Alzheimer's disease

Accurate integration of sensory inputs and motor commands is essential to achieve successful behavioral goals. A robust model of sensorimotor integration is the pitch perturbation response, in which speakers respond rapidly to shifts of the pitch in their auditory feedback. In a previous study, we demonstrated abnormal sensorimotor integration in patients with Alzheimer's disease (AD) with an abnormally enhanced behavioral response to pitch perturbation. Here we examine the neural correlates of the abnormal pitch perturbation response in AD patients, using magnetoencephalographic imaging. The participants phonated the vowel /α/ while a real-time signal processor briefly perturbed the pitch (100 cents, 400 ms) of their auditory feedback. We examined the high-gamma band (65-150 Hz) responses during this task. AD patients showed significantly reduced left prefrontal activity during the early phase of perturbation and increased right middle temporal activity during the later phase of perturbation, compared to controls. Activity in these brain regions significantly correlated with the behavioral response. These results demonstrate that impaired prefrontal modulation of speech-motor-control network and additional recruitment of right temporal regions are significant mediators of aberrant sensorimotor integration in patients with AD. The abnormal neural integration mechanisms signify the contribution of cortical network dysfunction to cognitive and behavioral deficits in AD.

EEG Source Imaging: A Practical Review of the Analysis Steps

The electroencephalogram (EEG) is one of the oldest technologies to measure neuronal activity of the human brain. It has its undisputed value in clinical diagnosis, particularly (but not exclusively) in the identification of epilepsy and sleep disorders and in the evaluation of dysfunctions in sensory transmission pathways. With the advancement of digital technologies, the analysis of EEG has moved from pure visual inspection of amplitude and frequency modulations over time to a comprehensive exploration of the temporal and spatial characteristics of the recorded signals. Today, EEG is accepted as a powerful tool to capture brain function with the unique advantage of measuring neuronal processes in the time frame in which these processes occur, namely in the sub-second range. However, it is generally stated that EEG suffers from a poor spatial resolution that makes it difficult to infer to the location of the brain areas generating the neuronal activity measured on the scalp. This statement has challenged a whole community of biomedical engineers to offer solutions to localize more precisely and more reliably the generators of the EEG activity. High-density EEG systems combined with precise information of the head anatomy and sophisticated source localization algorithms now exist that convert the EEG to a true neuroimaging modality. With these tools in hand and with the fact that EEG still remains versatile, inexpensive and portable, electrical neuroimaging has become a widely used technology to study the functions of the pathological and healthy human brain. However, several steps are needed to pass from the recording of the EEG to 3-dimensional images of neuronal activity. This review explains these different steps and illustrates them in a comprehensive analysis pipeline integrated in a stand-alone freely available academic software: Cartool. The information about how the different steps are performed in Cartool is only meant as a suggestion. Other EEG source imaging software may apply similar or different approaches to the different steps.

Neural correlates of reality filtering in schizophrenia spectrum disorder

BACKGROUND

A false sense of reality is a characteristic of schizophrenia spectrum disorders (SSD). Reality confusion may also emanate from posterior orbitofrontal cortex (OFC) lesions, as evident in confabulations that patients act upon and disorientation. This confusion can be measured by repeated runs of a continuous recognition task (CRT): patients increase their false positive rate from the second run on, failing to realize that an item is not a repetition within the current run. Correct handling of these stimuli, a faculty called orbitofrontal reality filtering (ORFi), induces a distinct frontal potential at 200-300 ms, the "ORFi potential". Patients with schizophrenia have been reported to fail in this task, too. Here, we explored the electrophysiology of ORFi in SSD.

METHODS

Evoked potentials, source, and connectivity analyses derived from high-density electroencephalograms of 17 patients with SSD and 15 age-matched healthy controls performing two runs of a CRT.

RESULTS

Although the patients obtained normal performance, they did not normally express the frontal potential typical of ORFi between 200 and 300 ms. Coherence analysis demonstrated virtually absent functional connectivity in the theta band within the memory network in this period. Source analysis showed increased activity in left medial temporal and prefrontal regions in patients.

CONCLUSIONS

SSD patients appear to invoke compensatory resources to handle the challenges of reality filtering. An abnormal ORFi potential may be an early biomarker of SSD.

Functional and Structural Brain Plasticity in Adult Onset Single-Sided Deafness

Single-sided deafness (SSD) or profound unilateral hearing loss obligates the only serviceable ear to capture all acoustic information. This loss of binaural function taxes cognitive resources for accurate listening performance, especially under adverse environments or challenging tasks. We hypothesized that adults with SSD would manifest both functional and structural brain plasticity compared to controls with normal binaural hearing. We evaluated functional alterations using magnetoencephalographic imaging (MEGI) of brain activation during performance of a moderately difficult auditory syllable sequence reproduction task and assessed structural integrity using diffusion tensor imaging (DTI). MEGI showed the SSD cohort to have increased induced oscillations in the theta band over the left superior temporal cortex and decreased induced gamma band oscillations over the frontal and parietal cortices between 175 and 475 ms following stimulus onset. DTI showed the SSD cohort to have extensive fractional anisotropy (FA) reduction in both auditory and non-auditory tracts and regions. Overlaying functional and structural changes revealed by the two imaging techniques demonstrated close registration of cortical areas and white matter tracts that expressed brain plasticity. Hence, complete loss of input from one ear in adulthood triggers both functional and structural alterations to dorsal temporal and frontal-parietal areas.

Empirical Comparison of Distributed Source Localization Methods for Single-Trial Detection of Movement Preparation

The development of technologies for the treatment of movement disorders, like stroke, is still of particular interest in brain-computer interface (BCI) research. In this context, source localization methods (SLMs), that reconstruct the cerebral origin of brain activity measured outside the head, e.g., via electroencephalography (EEG), can add a valuable insight into the current state and progress of the treatment. However, in BCIs SLMs were often solely considered as advanced signal processing methods that are compared against other methods based on the classification performance alone. Though, this approach does not guarantee physiological meaningful results. We present an empirical comparison of three established distributed SLMs with the aim to use one for single-trial movement prediction. The SLMs wMNE, sLORETA, and dSPM were applied on data acquired from eight subjects performing voluntary arm movements. Besides the classification performance as quality measure, a distance metric was used to asses the physiological plausibility of the methods. For the distance metric, which is usually measured to the source position of maximum activity, we further propose a variant based on clusters that is better suited for the single-trial case in which several sources are likely and the actual maximum is unknown. The two metrics showed different results. The classification performance revealed no significant differences across subjects, indicating that all three methods are equally well-suited for single-trial movement prediction. On the other hand, we obtained significant differences in the distance measure, favoring wMNE even after correcting the distance with the number of reconstructed clusters. Further, distance results were inconsistent with the traditional method using the maximum, indicating that for wMNE the point of maximum source activity often did not coincide with the nearest activation cluster. In summary, the presented comparison might help users to select an appropriate SLM and to understand the implications of the selection. The proposed methodology pays attention to the particular properties of distributed SLMs and can serve as a framework for further comparisons.

The neural dynamics of deficient memory control in heavily traumatized refugees

Victims of war, torture and natural catastrophes are prone to develop posttraumatic stress disorder (PTSD). These individuals experience the recurrent, involuntary intrusion of traumatic memories. What neurocognitive mechanisms are driving this memory disorder? Here we show that PTSD symptoms in heavily traumatized refugees are related to deficits in the effective control of memory retrieval. In a think/no-think task, PTSD patients were unable to forget memories that they had previously tried to suppress when compared to control participants with the same trauma history but without PTSD. Deficits in voluntary forgetting were clinically relevant since they correlated with memory intrusions in everyday life. Magnetoencephalography (MEG) recorded during suppression attempts revealed that PTSD patients were unable to downregulate signatures of sensory long-term memory traces in the gamma frequency band (70-120 Hz). Thus, our data suggest that the inability to suppress unwanted memories through modulation of gamma activity is related to PTSD symptom severity.

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

Cognitive neuroscience questions are commonly tested with experiments that involve a cohort of subjects. The cohort can consist of a handful of subjects for small studies to hundreds or thousands of subjects in open datasets. While there exist various online resources to get started with the analysis of magnetoencephalography (MEG) or electroencephalography (EEG) data, such educational materials are usually restricted to the analysis of a single subject. This is in part because data from larger group studies are harder to share, but also analyses of such data often require subject-specific decisions which are hard to document. This work presents the results obtained by the reanalysis of an open dataset from Wakeman and Henson (2015) using the MNE software package. The analysis covers preprocessing steps, quality assurance steps, sensor space analysis of evoked responses, source localization, and statistics in both sensor and source space. Results with possible alternative strategies are presented and discussed at different stages such as the use of high-pass filtering versus baseline correction, tSSS vs. SSS, the use of a minimum norm inverse vs. LCMV beamformer, and the use of univariate or multivariate statistics. This aims to provide a comparative study of different stages of M/EEG analysis pipeline on the same dataset, with open access to all of the scripts necessary to reproduce this analysis.

Modulating functional connectivity after stroke with neurofeedback: Effect on motor deficits in a controlled cross-over study

Synchronization of neural activity as measured with functional connectivity (FC) is increasingly used to study the neural basis of brain disease and to develop new treatment targets. However, solid evidence for a causal role of FC in disease and therapy is lacking. Here, we manipulated FC of the ipsilesional primary motor cortex in ten chronic human stroke patients through brain-computer interface technology with visual neurofeedback. We conducted a double-blind controlled crossover study to test whether manipulation of FC through neurofeedback had a behavioral effect on motor performance. Patients succeeded in increasing FC in the motor cortex. This led to improvement in motor function that was significantly greater than during neurofeedback training of a control brain area and proportional to the degree of FC enhancement. This result provides evidence that FC has a causal role in neurological function and that it can be effectively targeted with therapy.

•

Stroke patients participated in clinical trial on neurofeedback of functional connectivity.

•

Patients learned to enhance synchrony of neural activity in their motor cortex.

•

This led to reduced motor impairment.

•

Evidence for a causal role of neural synchrony in neurological deficits and recovery.

Amplitude by Peak Interaction but No Evidence of Auditory Mismatch Response Deficits to Frequency Change in Preschool-Aged Children with Fetal Alcohol Spectrum Disorders

BACKGROUND

Prior studies indicate that the auditory mismatch response is sensitive to early alterations in brain development in multiple developmental disorders. Prenatal alcohol exposure is known to impact early auditory processing. The current study hypothesized alterations in the mismatch response in young children with fetal alcohol spectrum disorders (FASD).

METHODS

Participants in this study were 9 children with a FASD and 17 control children (Control) aged 3 to 6 years. Participants underwent magnetoencephalography and structural magnetic resonance imaging scans separately. We compared groups on neurophysiological mismatch negativity (MMN) responses to auditory stimuli measured using the auditory oddball paradigm. Frequent (1,000 Hz) and rare (1,200 Hz) tones were presented at 72 dB.

RESULTS

There was no significant group difference in MMN response latency or amplitude represented by the peak located ~200 ms after stimulus presentation in the difference time course between frequent and infrequent tones. Examining the time courses to the frequent and infrequent tones separately, repeated measures analysis of variance with condition (frequent vs. rare), peak (N100m and N200m), and hemisphere as within-subject factors and diagnosis and sex as the between-subject factors showed a significant interaction of peak by diagnosis (p = 0.001), with a pattern of decreased amplitude from N100m to N200m in Control children and the opposite pattern in children with FASD. However, no significant difference was found with the simple effects comparisons. No group differences were found in the response latencies of the rare auditory evoked fields.

CONCLUSIONS

The results indicate that there was no detectable effect of alcohol exposure on the amplitude or latency of the MMNm response to simple tones modulated by frequency change in preschool-aged children with FASD. However, while discrimination abilities to simple tones may be intact, early auditory sensory processing revealed by the interaction between N100m and N200m amplitude indicates that auditory sensory processing may be altered in children with FASD.

Choice of Magnetometers and Gradiometers after Signal Space Separation

BACKGROUND

Modern Elekta Neuromag MEG devices include 102 sensor triplets containing one magnetometer and two planar gradiometers. The first processing step is often a signal space separation (SSS), which provides a powerful noise reduction. A question commonly raised by researchers and reviewers relates to which data should be employed in analyses: (1) magnetometers only, (2) gradiometers only, (3) magnetometers and gradiometers together. The MEG community is currently divided with regard to the proper answer.

METHODS

First, we provide theoretical evidence that both gradiometers and magnetometers result from the backprojection of the same SSS components. Then, we compare resting state and task-related sensor and source estimations from magnetometers and gradiometers in real MEG recordings before and after SSS.

RESULTS

SSS introduced a strong increase in the similarity between source time series derived from magnetometers and gradiometers (r² = 0.3-0.8 before SSS and r² > 0.80 after SSS). After SSS, resting state power spectrum and functional connectivity, as well as visual evoked responses, derived from both magnetometers and gradiometers were highly similar (Intraclass Correlation Coefficient > 0.8, r² > 0.8).

CONCLUSIONS

After SSS, magnetometer and gradiometer data are estimated from a single set of SSS components (usually ≤ 80). Equivalent results can be obtained with both sensor types in typical MEG experiments.

Simultaneous Reality Filtering and Encoding of Thoughts: The Substrate for Distinguishing between Memories of Real Events and Imaginations?

Any thought, whether it refers to the present moment or reflects an imagination, is again encoded as a new memory trace. Orbitofrontal reality filtering (ORFi) denotes an on-line mechanism which verifies whether upcoming thoughts relate to ongoing reality or not. Its failure induces reality confusion with confabulations and disorientation. If the result of this process were simultaneously encoded, it would easily explain later distinction between memories relating to a past reality and memories relating to imagination, a faculty called reality monitoring. How the brain makes this distinction is unknown but much research suggests that it depends on processes active when information is encoded. Here we explored the precise timing between ORFi and encoding as well as interactions between the involved brain structures. We used high-density evoked potentials and two runs of a continuous recognition task (CRT) combining the challenges of ORFi and encoding. ORFi was measured by the ability to realize that stimuli appearing in the second run had not appeared in this run yet. Encoding was measured with immediately repeated stimuli, which has been previously shown to induce a signal emanating from the medial temporal lobe (MTL), which has a protective effect on the memory trace. We found that encoding, as measured with this task, sets in at about 210 ms after stimulus presentation, 35 ms before ORFi. Both processes end at about 330 ms. Both were characterized by increased coherence in the theta band in the MTL during encoding and in the orbitofrontal cortex (OFC) during ORFi. The study suggests a complex interaction between OFC and MTL allowing for thoughts to be re-encoded while they undergo ORFi. The combined influence of these two processes at 200-300 ms may leave a memory trace that allows for later effortless reality monitoring in most everyday situations.

Longitudinal Structural and Functional Differences Between Proportional and Poor Motor Recovery After Stroke

BACKGROUND

Evolution of motor function during the first months after stroke is stereotypically bifurcated, consisting of either recovery to about 70% of maximum possible improvement ("proportional recovery, PROP") or in little to no improvement ("poor recovery, POOR"). There is currently no evidence that any rehabilitation treatment will prevent POOR and favor PROP.

OBJECTIVE

To perform a longitudinal and multimodal assessment of functional and structural changes in brain organization associated with PROP.

METHODS

Fugl-Meyer Assessments of the upper extremity and high-density electroencephalography (EEG) were obtained from 63 patients, diffusion tensor imaging from 46 patients, at 2 and 4 weeks (T0) and at 3 months (T1) after stroke onset.

RESULTS

We confirmed the presence of 2 distinct recovery patterns (PROP and POOR) in our sample. At T0, PROP patients had greater integrity of the corticospinal tract (CST) and greater EEG functional connectivity (FC) between the affected hemisphere and rest of the brain, in particular between the ventral premotor and the primary motor cortex. POOR patients suffered from degradation of corticocortical and corticofugal fiber tracts in the affected hemisphere between T0 and T1, which was not observed in PROP patients. Better initial CST integrity correlated with greater initial global FC, which was in turn associated with less white matter degradation between T0 and T1.

CONCLUSIONS

These findings suggest links between initial CST integrity, systems-level cortical network plasticity, reduction of white matter atrophy, and clinical motor recovery after stroke. This identifies candidate treatment targets.

Distinct spatiotemporal patterns of neuronal functional connectivity in primary progressive aphasia variants

Primary progressive aphasia is a syndrome characterized by progressive loss of language abilities with three main phenotypic clinical presentations, including logopenic, non-fluent/agrammatic, and semantic variants. Previous imaging studies have shown unique anatomic impacts within language networks in each variant. However, direct measures of spontaneous neuronal activity and functional integrity of these impacted neural networks in primary progressive aphasia are lacking. The aim of this study was to characterize the spatial and temporal patterns of resting state neuronal synchronizations in primary progressive aphasia syndromes. We hypothesized that resting state brain oscillations will show unique deficits within language network in each variant of primary progressive aphasia. We examined 39 patients with primary progressive aphasia including logopenic variant (n = 14, age = 61 ± 9 years), non-fluent/agrammatic variant (n = 12, age = 71 ± 8 years) and semantic variant (n = 13, age = 65 ± 7 years) using magnetoencephalographic imaging, compared to a control group that was matched in age and gender to each primary progressive aphasia subgroup (n = 20, age = 65 ± 5 years). Each patient underwent a complete clinical evaluation including a comprehensive battery of language tests. We examined the whole-brain resting state functional connectivity as measured by imaginary coherence in each patient group compared to the control cohort, in three frequency oscillation bands-delta-theta (2-8 Hz); alpha (8-12 Hz); beta (12-30 Hz). Each variant showed a distinct spatiotemporal pattern of altered functional connectivity compared to age-matched controls. Specifically, we found significant hyposynchrony of alpha and beta frequency within the left posterior temporal and occipital cortices in patients with the logopenic variant, within the left inferior frontal cortex in patients with the non-fluent/agrammatic variant, and within the left temporo-parietal junction in patients with the semantic variant. Patients with logopenic variant primary progressive aphasia also showed significant hypersynchrony of delta-theta frequency within bilateral medial frontal and posterior parietal cortices. Furthermore, region of interest-based analyses comparing the spatiotemporal patterns of variant-specific regions of interest identified in comparison to age-matched controls showed significant differences between primary progressive aphasia variants themselves. We also found distinct patterns of regional spectral power changes in each primary progressive aphasia variant, compared to age-matched controls. Our results demonstrate neurophysiological signatures of network-specific neuronal dysfunction in primary progressive aphasia variants. The unique spatiotemporal patterns of neuronal synchrony signify diverse neurophysiological disruptions and pathological underpinnings of the language network in each variant.