A new approach to neuroimaging with magnetoencephalography

Neural dynamics of verbal working memory processing in children and adolescents

Development of cognitive functions and the underlying neurophysiology is evident throughout childhood and adolescence, with higher order processes such as working memory (WM) being some of the last cognitive faculties to fully mature. Previous functional neuroimaging studies of the neurodevelopment of WM have largely focused on overall regional activity levels rather than the temporal dynamics of neural component recruitment. In this study, we used magnetoencephalography (MEG) to examine the neural dynamics of WM in a large cohort of children and adolescents who were performing a high-load, modified verbal Sternberg WM task. Consistent with previous studies in adults, our findings indicated left-lateralized activity throughout the task period, beginning in the occipital cortices and spreading anterior to include temporal and prefrontal cortices during later encoding and into maintenance. During maintenance, the occipital alpha increase that has been widely reported in adults was found to be relatively weak in this developmental sample, suggesting continuing development of this component of neural processing, which was supported by correlational analyses. Intriguingly, we also found sex-specific developmental effects in alpha responses in the right inferior frontal region during encoding and in parietal and occipital cortices during maintenance. These findings suggested a developmental divergence between males and females in the maturation of neural circuitry serving WM during the transition from childhood to adolescence.

Spatial frequency supports the emergence of categorical representations in visual cortex during natural scene perception

In navigating our environment, we rapidly process and extract meaning from visual cues. However, the relationship between visual features and categorical representations in natural scene perception is still not well understood. Here, we used natural scene stimuli from different categories and filtered at different spatial frequencies to address this question in a passive viewing paradigm. Using representational similarity analysis (RSA) and cross-decoding of magnetoencephalography (MEG) data, we show that categorical representations emerge in human visual cortex at ∼180 ms and are linked to spatial frequency processing. Furthermore, dorsal and ventral stream areas reveal temporally and spatially overlapping representations of low and high-level layer activations extracted from a feedforward neural network. Our results suggest that neural patterns from extrastriate visual cortex switch from low-level to categorical representations within 200 ms, highlighting the rapid cascade of processing stages essential in human visual perception.

Spatiotemporal dynamics in human visual cortex rapidly encode the emotional content of faces

Recognizing emotion in faces is important in human interaction and survival, yet existing studies do not paint a consistent picture of the neural representation supporting this task. To address this, we collected magnetoencephalography (MEG) data while participants passively viewed happy, angry and neutral faces. Using time-resolved decoding of sensor-level data, we show that responses to angry faces can be discriminated from happy and neutral faces as early as 90 ms after stimulus onset and only 10 ms later than faces can be discriminated from scrambled stimuli, even in the absence of differences in evoked responses. Time-resolved relevance patterns in source space track expression-related information from the visual cortex (100 ms) to higher-level temporal and frontal areas (200-500 ms). Together, our results point to a system optimised for rapid processing of emotional faces and preferentially tuned to threat, consistent with the important evolutionary role that such a system must have played in the development of human social interactions.

Cognitive loading via mental arithmetic modulates effects of blink-related oscillations on precuneus and ventral attention network regions

Blink-related oscillations (BROs) have been linked with environmental monitoring processes associated with blinking, with cortical activations in the bilateral precuneus. Although BROs have been described under resting and passive fixation conditions, little is known about their characteristics under cognitive loading. To address this, we investigated BRO effects during both mental arithmetic (MA) and passive fixation (PF) tasks using magnetoencephalography (n =20), while maintaining the same sensory environment in both tasks. Our results confirmed the presence of BRO effects in both MA and PF tasks, with similar characteristics including blink-related increase in global field power and blink-related activation of the bilateral precuneus. In addition, cognitive loading due to MA also modulated BRO effects by decreasing BRO-induced cortical activations in key brain regions including the bilateral anterior precuneus. Interestingly, blinking during MA-but not PF-activated regions of the ventral attention network (i.e., right supramarginal gyrus and inferior frontal gyrus), suggesting possible recruitment of these areas for blink processing under cognitive loading conditions. Time-frequency analysis revealed a consistent pattern of BRO-related effects in the precuneus in both tasks, but with task-related functional segregation within the anterior and posterior subregions. Based on these findings, we postulate a potential neurocognitive mechanism for blink processing in the precuneus. This study is the first investigation of BRO effects under cognitive loading, and our results provide compelling new evidence for the important cognitive implications of blink-related processing in the human brain.

Haptic exploration attenuates and alters somatosensory cortical oscillations

KEY POINTS

Several behavioural studies have shown the sensory perceptions are reduced during movement; yet the neurophysiological reason for this is not clear. Participants underwent stimulation of the median nerve when either sitting quietly (i.e. passive stimulation condition) or performing haptic exploration of a ball with the left hand. Magnetoencephalographic brain imaging and advanced beamforming methods were used to identify the differences in somatosensory cortical responses. We show that the neural populations active during the passive stimulation condition were strongly gated during the haptic exploration task. These results imply that the reduced haptic perceptions might be governed by gating of certain somatosensory neural populations.

ABSTRACT

Several behavioural studies have shown that children have reduced sensory perceptions during movement; however, the neurophysiological nexus for these altered perceptions remains unknown. We used magnetoencephalographic brain imaging and advanced beamforming methods to address this knowledge gap. In our experiment, a cohort of children (aged 10-18 years) underwent stimulation of the median nerve when either sitting quietly (i.e. passive stimulation condition) or performing haptic exploration of a ball with the left hand. Our results revealed two novel observations. First, there was a relationship between the child's age and the strength of the beta (18-26 Hz) response seen within the somatosensory cortices during the passive stimulation condition. This suggests that there may be an age-dependent change in the processing of peripheral feedback by the somatosensory cortices. Second, all of the cortical regions that were active during the passive stimulation condition were almost completely gated during the haptic task. Instead, the haptic task involved neural oscillations within Brodmann area 2, which is known to convey less spatially precise tactile information but is involved in the processing of more complex somatosensations across the respective digits. These results imply that the reduced somatosensory perceptions seen during movements in healthy children may be related to the gating of certain neural generators, as well as activation of haptic-specific neural generators within the somatosensory cortices. The utilization of such haptic-specific circuits during development may lead to the enhanced somatosensory processing during haptic exploration seen in healthy adults.

The developmental trajectory of sensorimotor cortical oscillations

Numerous studies of motor control have confirmed beta and gamma oscillations in the primary motor cortices during basic movements. These responses include a robust beta decrease that precedes and extends through movement onset, a transient gamma response that coincides with the movement, and a post-movement beta rebound (PMBR) response that occurs after movement offset. While the existence of these responses has been confirmed by many studies, very few studies have examined their developmental trajectory. In the current study, we utilized magnetoencephalography (MEG) to investigate age-related changes in sensorimotor cortical oscillations in a large cross-section of children and adolescents (n = 94; age range = 9 -15 years-old). All participants performed a stimulus detection task with their right finger and the resulting MEG data were examined using oscillatory analysis methods and imaged using a beamformer. Consistent with adult studies, these youth participants exhibited characteristic beta (16-24 Hz) decreases prior to and during movement, as well as PMBR responses following movement offset, and a transient gamma (74-84 Hz) response during movement execution. Our primary findings were that the strength of the PMBR increased with age, while the strength of the gamma synchronization decreased with chronological age. In addition, the strength of each motor-related oscillatory response was significantly correlated with the power of spontaneous activity in the same frequency range and same voxel. This was the case for all three oscillatory responses. In conclusion, we investigated motor-related oscillatory activity in the largest cohort of children and adolescents reported to date, and our results indicated that beta and gamma cortical oscillations continue to develop as children transition into adolescents, and that these responses may not be fully matured until young to middle adulthood.

Load modulates the alpha and beta oscillatory dynamics serving verbal working memory

A network of predominantly left-lateralized brain regions has been linked to verbal working memory (VWM) performance. However, the impact of memory load on the oscillatory dynamics serving VWM is far less understood. To further investigate this, we had 26 healthy adults perform a high-load (6 letter) and low-load (4 letter) variant of a VWM task while undergoing magnetoencephalography (MEG). MEG data were evaluated in the time-frequency domain and significant oscillatory responses spanning the encoding and maintenance phases were reconstructed using a beamformer. To determine the impact of load on the neural dynamics, the resulting images were examined using paired-samples t-tests and virtual sensor analyses. Our results indicated stronger increases in frontal theta activity in the high- relative to low-load condition during early encoding. Stronger decreases in alpha/beta activity were also observed during encoding in bilateral posterior cortices during the high-load condition, and the strength of these load effects increased as encoding progressed. During maintenance, stronger decreases in alpha activity in the left inferior frontal gyrus, middle temporal gyrus, supramarginal gyrus, and inferior parietal cortices were detected during high- relative to low-load performance, with the strength of these load effects remaining largely static throughout maintenance. Finally, stronger increases in occipital alpha activity were observed during maintenance in the high-load condition, and the strength of these effects grew stronger with time during the first half of maintenance, before dissipating during the latter half of maintenance. Notably, this was the first study to utilize a whole-brain approach to statistically evaluate the temporal dynamics of load-related oscillatory differences during encoding and maintenance processes, and our results highlight the importance of spatial, temporal, and spectral specificity in this regard.

Reliability of Static and Dynamic Network Metrics in the Resting-State: A MEG-Beamformed Connectivity Analysis

The resting activity of the brain can be described by so-called intrinsic connectivity networks (ICNs), which consist of spatially and temporally distributed, but functionally connected, nodes. The coordinated activity of the resting state can be explored via magnetoencephalography (MEG) by studying frequency-dependent functional brain networks at the source level. Although many algorithms for the analysis of brain connectivity have been proposed, the reliability of network metrics derived from both static and dynamic functional connectivity is still unknown. This is a particular problem for studies of associations between ICN metrics and personality variables or other traits, and for studies of differences between patient and control groups, which both depend critically on the reliability of the metrics used. A detailed investigation of the reliability of metrics derived from resting-state MEG repeat scans is therefore a prerequisite for the development of connectomic biomarkers. Here, we first estimated both static (SFC) and dynamic functional connectivity (DFC) after beamforming source reconstruction using the imaginary part of the phase locking index (iPLV) and the correlation of the amplitude envelope (CorEnv). Using our approach, functional network microstates (FCμstates) were derived from the DFC and chronnectomics were computed from the evolution of FCμstates across experimental time. In both temporal scales, the reliability of network metrics (SFC), the FCμstates and the related chronnectomics were evaluated for every frequency band. Chronnectomic statistics and FCμstates were generally more reliable than node-wise static network metrics. CorEnv-based network metrics were more reproducible at the static approach. The reliability of chronnectomics have been evaluated also in a second dataset. This study encourages the analysis of MEG resting-state via DFC.

IFCN-endorsed practical guidelines for clinical magnetoencephalography (MEG)

Magnetoencephalography (MEG) records weak magnetic fields outside the human head and thereby provides millisecond-accurate information about neuronal currents supporting human brain function. MEG and electroencephalography (EEG) are closely related complementary methods and should be interpreted together whenever possible. This manuscript covers the basic physical and physiological principles of MEG and discusses the main aspects of state-of-the-art MEG data analysis. We provide guidelines for best practices of patient preparation, stimulus presentation, MEG data collection and analysis, as well as for MEG interpretation in routine clinical examinations. In 2017, about 200 whole-scalp MEG devices were in operation worldwide, many of them located in clinical environments. Yet, the established clinical indications for MEG examinations remain few, mainly restricted to the diagnostics of epilepsy and to preoperative functional evaluation of neurosurgical patients. We are confident that the extensive ongoing basic MEG research indicates potential for the evaluation of neurological and psychiatric syndromes, developmental disorders, and the integrity of cortical brain networks after stroke. Basic and clinical research is, thus, paving way for new clinical applications to be identified by an increasing number of practitioners of MEG.

Aberrant oscillatory dynamics during somatosensory processing in HIV-infected adults

While the arrival of combination antiretroviral therapy significantly decreased the prevalence of HIV-associated dementia, between 35 and 70% of all infected adults continue to develop some form of cognitive impairment. These deficits appears to affect multiple neural subsystems, but the mechanisms and extent of damage are not fully understood. In the current study, we utilized magnetoencephalography (MEG), advanced oscillatory analysis methods, and a paired-pulse somatosensory stimulation paradigm to interrogate pre-attentive inhibitory processing in 43 HIV-infected adults and 28 demographically-matched uninfected controls. MEG responses were imaged using a beamformer, and time series data were extracted from the peak voxel in grand-averaged functional brain images to quantify the dynamics of sensory gating, oscillatory power, spontaneous power, and other neural indices. We found a significantly weakened response to the second stimulation compared to the first across groups, indicating significant sensory gating irrespective of HIV-infection. Interestingly, HIV-infected participants exhibited reduced neural responses in the 20-75 Hz gamma range to each somatosensory stimulation compared to uninfected controls, and exhibited significant alterations in peak gamma frequency in response to the second stimulation. Finally, HIV-infected participants also had significantly stronger spontaneous activity in the gamma range (i.e., 20-75 Hz) during the baseline period before stimulation onset. In conclusion, while HIV-infected participants had the capacity to efficiently gate somatosensory input, their overall oscillatory responses were weaker, spontaneous baseline activity was stronger, and their response to the second stimulation had an altered peak gamma frequency. We propose that this pattern of deficits suggests dysfunction in the somatosensory cortices, which is potentially secondary to accelerated aging.

Neurophysiological changes in the visuomotor network after practicing a motor task

Although it is well appreciated that practicing a motor task updates the associated internal model, it is still unknown how the cortical oscillations linked with the motor action change with practice. The present study investigates the short-term changes (e.g., fast motor learning) in the α- and β-event-related desynchronizations (ERD) associated with the production of a motor action. To this end, we used magnetoencephalography to identify changes in the α- and β-ERD in healthy adults after participants practiced a novel isometric ankle plantarflexion target-matching task. After practicing, the participants matched the targets faster and had improved accuracy, faster force production, and a reduced amount of variability in the force output when trying to match the target. Parallel with the behavioral results, the strength of the β-ERD across the motor-planning and execution stages was reduced after practice in the sensorimotor and occipital cortexes. No pre/postpractice changes were found in the α-ERD during motor planning or execution. Together, these outcomes suggest that fast motor learning is associated with a decrease in β-ERD power. The decreased strength likely reflects a more refined motor plan, a reduction in neural resources needed to perform the task, and/or an enhancement of the processes that are involved in the visuomotor transformations that occur before the onset of the motor action. These results may augment the development of neurologically based practice strategies and/or lead to new practice strategies that increase motor learning. NEW & NOTEWORTHY We aimed to determine the effects of practice on the movement-related cortical oscillatory activity. Following practice, we found that the performance of the ankle plantarflexion target-matching task improved and the power of the β-oscillations decreased in the sensorimotor and occipital cortexes. These novel findings capture the β-oscillatory activity changes in the sensorimotor and occipital cortexes that are coupled with behavioral changes to demonstrate the effects of motor learning.

Efficacy and mechanism of sub-sensory sacral (optimised) neuromodulation in adults with faecal incontinence: study protocol for a randomised controlled trial

Background

Faecal incontinence (FI) is a substantial health problem with a prevalence of approximately 8% in community-dwelling populations. Sacral neuromodulation (SNM) is considered the first-line surgical treatment option in adults with FI in whom conservative therapies have failed. The clinical efficacy of SNM has never been rigorously determined in a trial setting and the underlying mechanism of action remains unclear.

Methods/design

The design encompasses a multicentre, randomised, double-blind crossover trial and cohort follow-up study. Ninety participants will be randomised to one of two groups (SNM/SHAM or SHAM/SNM) in an allocation ratio of 1:1. The main inclusion criteria will be adults aged 18–75 years meeting Rome III and ICI definitions of FI, who have failed non-surgical treatments to the UK standard, who have a minimum of eight FI episodes in a 4-week screening period, and who are clinically suitable for SNM. The primary objective is to estimate the clinical efficacy of sub-sensory SNM vs. SHAM at 32 weeks based on the primary outcome of frequency of FI episodes using a 4-week paper diary, using mixed Poisson regression analysis on the intention-to-treat principle. The study is powered (0.9) to detect a 30% reduction in frequency of FI episodes between sub-sensory SNM and SHAM stimulation over a 32-week crossover period.

Secondary objectives include: measurement of established and new clinical outcomes after 1 year of therapy using new (2017 published) optimised therapy (with standardised SNM-lead placement); validation of new electronic outcome measures (events) and a device to record them, and identification of potential biological effects of SNM on underlying anorectal afferent neuronal pathophysiology (hypothesis: SNM leads to increased frequency of perceived transient anal sphincter relaxations; improved conscious sensation of defaecatory urge and cortical/subcortical changes in afferent responses to anorectal electrical stimulation (main techniques: high-resolution anorectal manometry and magnetoencephalography).

Discussion

This trial will determine clinical effect size for sub-sensory chronic electrical stimulation of the sacral innervation. It will provide experimental evidence of modifiable afferent neurophysiology that may aid future patient selection as well as a basic understanding of the pathophysiology of FI.

Trial registration

International Standard Randomised Controlled Trial Number: ISRCTN98760715. Registered on 15 September 2017.

Electronic supplementary material

The online version of this article (10.1186/s13063-018-2689-1) contains supplementary material, which is available to authorized users.

Improving the Nulling Beamformer Using Subspace Suppression

Magnetoencephalography (MEG) captures the magnetic fields generated by neuronal current sources with sensors outside the head. In MEG analysis these current sources are estimated from the measured data to identify the locations and time courses of neural activity. Since there is no unique solution to this so-called inverse problem, multiple source estimation techniques have been developed. The nulling beamformer (NB), a modified form of the linearly constrained minimum variance (LCMV) beamformer, is specifically used in the process of inferring interregional interactions and is designed to eliminate shared signal contributions, or cross-talk, between regions of interest (ROIs) that would otherwise interfere with the connectivity analyses. The nulling beamformer applies the truncated singular value decomposition (TSVD) to remove small signal contributions from a ROI to the sensor signals. However, ROIs with strong crosstalk will have high separating power in the weaker components, which may be removed by the TSVD operation. To address this issue we propose a new method, the nulling beamformer with subspace suppression (NBSS). This method, controlled by a tuning parameter, reweights the singular values of the gain matrix mapping from source to sensor space such that components with high overlap are reduced. By doing so, we are able to measure signals between nearby source locations with limited cross-talk interference, allowing for reliable cortical connectivity analysis between them. In two simulations, we demonstrated that NBSS reduces cross-talk while retaining ROIs' signal power, and has higher separating power than both the minimum norm estimate (MNE) and the nulling beamformer without subspace suppression. We also showed that NBSS successfully localized the auditory M100 event-related field in primary auditory cortex, measured from a subject undergoing an auditory localizer task, and suppressed cross-talk in a nearby region in the superior temporal sulcus.

Virtual localization of the seizure onset zone: Using non-invasive MEG virtual electrodes at stereo-EEG electrode locations in refractory epilepsy patients

In some patients with medically refractory epilepsy, EEG with intracerebrally placed electrodes (stereo-electroencephalography, SEEG) is needed to locate the seizure onset zone (SOZ) for successful epilepsy surgery. SEEG has limitations and entails risk of complications because of its invasive character. Non-invasive magnetoencephalography virtual electrodes (MEG-VEs) may overcome SEEG limitations and optimize electrode placement making SOZ localization safer. Our purpose was to assess whether interictal activity measured by MEG-VEs and SEEG at identical anatomical locations were comparable, and whether MEG-VEs activity properties could determine the location of a later resected brain area (RA) as an approximation of the SOZ. We analyzed data from nine patients who underwent MEG and SEEG evaluation, and surgery for medically refractory epilepsy. MEG activity was retrospectively reconstructed using beamforming to obtain VEs at the anatomical locations corresponding to those of SEEG electrodes. Spectral, functional connectivity and functional network properties were obtained for both, MEG-VEs and SEEG time series, and their correlation and reliability were established. Based on these properties, the approximation of the SOZ was characterized by the differences between RA and non-RA (NRA). We found significant positive correlation and reliability between MEG-VEs and SEEG spectral measures (particularly in delta [0.5-4 Hz], alpha2 [10-13 Hz], and beta [13-30 Hz] bands) and broadband functional connectivity. Both modalities showed significantly slower activity and a tendency towards increased broadband functional connectivity in the RA compared to the NRA. Our findings show that spectral and functional connectivity properties of non-invasively obtained MEG-VEs match those of invasive SEEG recordings, and can characterize the SOZ. This suggests that MEG-VEs might be used for optimal SEEG planning and fewer depth electrode implantations, making the localization of the SOZ safer and more successful.

Altered Brain Dynamics in Patients With Type 1 Diabetes During Working Memory Processing

It is now generally accepted that diabetes increases the risk for cognitive impairment, but the precise mechanisms are poorly understood. A critical problem in linking diabetes to cognitive impairment is that patients often have multiple comorbidities (e.g., obesity, hypertension) that have been independently linked to cognitive deficits. In the study reported here we focused on young adults with and without type 1 diabetes who were virtually free of such comorbidities. The two groups were matched on major health and demographic factors, and all participants completed a verbal working memory task during magnetoencephalographic brain imaging. We hypothesized that patients would have altered neural dynamics in verbal working memory processing and that these differences would directly relate to clinical disease measures. Accordingly, we found that patients had significantly stronger neural responses in the superior parietal cortices during memory encoding and significantly weaker activity in parietal-occipital regions during maintenance compared with control subjects. Moreover, disease duration and glycemic control were both significantly correlated with neural responses in various brain regions. In conclusion, young healthy adults with type 1 diabetes already have aberrant neural processing relative to their peers without diabetes, using compensatory responses to perform the task, and glucose management and duration may play a central role.

Beta Oscillatory Dynamics in the Prefrontal and Superior Temporal Cortices Predict Spatial Working Memory Performance

The oscillatory dynamics serving spatial working memory (SWM), and how such dynamics relate to performance, are poorly understood. To address these topics, the present study recruited 22 healthy adults to perform a SWM task during magnetoencephalography (MEG). The resulting MEG data were transformed into the time-frequency domain, and significant oscillatory responses were imaged using a beamformer. Voxel time series data were extracted from the cluster peaks to quantify the dynamics, while whole-brain partial correlation maps were computed to identify regions where oscillatory strength varied with accuracy on the SWM task. The results indicated transient theta oscillations in spatially distinct subregions of the prefrontal cortices at the onset of encoding and maintenance, which may underlie selection of goal-relevant information. Additionally, strong and persistent decreases in alpha and beta oscillations were observed throughout encoding and maintenance in parietal, temporal, and occipital regions, which could serve sustained attention and maintenance processes during SWM performance. The neuro-behavioral correlations revealed that beta activity within left dorsolateral prefrontal control regions and bilateral superior temporal integration regions was negatively correlated with SWM accuracy. Notably, this is the first study to employ a whole-brain approach to significantly link neural oscillations to behavioral performance in the context of SWM.

Group Analysis in FieldTrip of Time-Frequency Responses: A Pipeline for Reproducibility at Every Step of Processing, Going From Individual Sensor Space Representations to an Across-Group Source Space Representation

An important aim of an analysis pipeline for magnetoencephalographic (MEG) data is that it allows for the researcher spending maximal effort on making the statistical comparisons that will answer his or her questions. The example question being answered here is whether the so-called beta rebound differs between novel and repeated stimulations. Two analyses are presented: going from individual sensor space representations to, respectively, an across-group sensor space representation and an across-group source space representation. The data analyzed are neural responses to tactile stimulations of the right index finger in a group of 20 healthy participants acquired from an Elekta Neuromag System. The processing steps covered for the first analysis are MaxFiltering the raw data, defining, preprocessing and epoching the data, cleaning the data, finding and removing independent components related to eye blinks, eye movements and heart beats, calculating participants' individual evoked responses by averaging over epoched data and subsequently removing the average response from single epochs, calculating a time-frequency representation and baselining it with non-stimulation trials and finally calculating a grand average, an across-group sensor space representation. The second analysis starts from the grand average sensor space representation and after identification of the beta rebound the neural origin is imaged using beamformer source reconstruction. This analysis covers reading in co-registered magnetic resonance images, segmenting the data, creating a volume conductor, creating a forward model, cutting out MEG data of interest in the time and frequency domains, getting Fourier transforms and estimating source activity with a beamformer model where power is expressed relative to MEG data measured during periods of non-stimulation. Finally, morphing the source estimates onto a common template and performing group-level statistics on the data are covered. Functions for saving relevant figures in an automated and structured manner are also included. The protocol presented here can be applied to any research protocol where the emphasis is on source reconstruction of induced responses where the underlying sources are not coherent.

tDCS Modulates Visual Gamma Oscillations and Basal Alpha Activity in Occipital Cortices: Evidence from MEG

Transcranial direct-current stimulation (tDCS) is now a widely used method for modulating the human brain, but the resulting physiological effects are not understood. Recent studies have combined magnetoencephalography (MEG) with simultaneous tDCS to evaluate online changes in occipital alpha and gamma oscillations, but no study to date has quantified the offline (i.e., after tDCS) alterations in these responses. Thirty-five healthy adults received active or sham anodal tDCS to the occipital cortices, and then completed a visual stimulation paradigm during MEG that is known to elicit robust gamma and alpha oscillations. The resulting MEG data were imaged and peak voxel time series were extracted to evaluate tDCS effects. We found that tDCS to the occipital increased the amplitude of local gamma oscillations, and basal alpha levels during the baseline. tDCS was also associated with network-level effects, including increased gamma oscillations in the prefrontal cortex, parietal, and other visual attention regions. Finally, although tDCS did not modulate peak gamma frequency, this variable was inversely correlated with gamma amplitude, which is consistent with a GABA-gamma link. In conclusion, tDCS alters gamma oscillations and basal alpha levels. The net offline effects on gamma activity are consistent with the view that anodal tDCS decreases local GABA.

Non-invasive Investigation of Human Hippocampal Rhythms Using Magnetoencephalography: A Review

Hippocampal rhythms are believed to support crucial cognitive processes including memory, navigation, and language. Due to the location of the hippocampus deep in the brain, studying hippocampal rhythms using non-invasive magnetoencephalography (MEG) recordings has generally been assumed to be methodologically challenging. However, with the advent of whole-head MEG systems in the 1990s and development of advanced source localization techniques, simulation and empirical studies have provided evidence that human hippocampal signals can be sensed by MEG and reliably reconstructed by source localization algorithms. This paper systematically reviews simulation studies and empirical evidence of the current capacities and limitations of MEG “deep source imaging” of the human hippocampus. Overall, these studies confirm that MEG provides a unique avenue to investigate human hippocampal rhythms in cognition, and can bridge the gap between animal studies and human hippocampal research, as well as elucidate the functional role and the behavioral correlates of human hippocampal oscillations.

Increased Small-World Network Topology Following Deployment-Acquired Traumatic Brain Injury Associated with the Development of Post-Traumatic Stress Disorder

Cross-sectional and longitudinal studies in active duty and veteran cohorts have both demonstrated that deployment-acquired traumatic brain injury (TBI) is an independent risk factor for developing post-traumatic stress disorder (PTSD), beyond confounds such as combat exposure, physical injury, predeployment TBI, and pre-deployment psychiatric symptoms. This study investigated how resting-state brain networks differ between individuals who developed PTSD and those who did not following deployment-acquired TBI. Participants included postdeployment veterans with deployment-acquired TBI history both with and without current PTSD diagnosis. Graph metrics, including small-worldness, clustering coefficient, and modularity, were calculated from individually constructed whole-brain networks based on 5-min eyes-open resting-state magnetoencephalography (MEG) recordings. Analyses were adjusted for age and premorbid IQ. Results demonstrated that participants with current PTSD displayed higher levels of small-worldness, F(1,12) = 5.364, p < 0.039, partial eta squared = 0.309, and Cohen's d = 0.972, and clustering coefficient, F(1, 12) = 12.204, p < 0.004, partial eta squared = 0.504, and Cohen's d = 0.905, than participants without current PTSD. There were no between-group differences in modularity or the number of modules present. These findings are consistent with a hyperconnectivity hypothesis of the effect of TBI history on functional networks rather than a disconnection hypothesis, demonstrating increased levels of clustering coefficient rather than a decrease as might be expected; however, these results do not account for potential changes in brain structure. These results demonstrate the potential pathological sequelae of changes in functional brain networks following deployment-acquired TBI and represent potential neurobiological changes associated with deployment-acquired TBI that may increase the risk of subsequently developing PTSD.

The neural effects of positively and negatively re-experiencing mental fatigue sensation: a magnetoencephalography study

Fatigue sensation is an essential biological alarm that urges us to take rest to avoid disrupting homeostasis and thus plays an important role in maintaining well-being. However, there are situations in which the anticipation of unpleasant fatigue sensation undesirably reduces motivation for activity. The aim of this study was to examine whether thinking positively about the fatigue sensation would increase motivation to accomplish the workload. Fourteen healthy male volunteers participated in this study and performed a two-back test for 30 min to induce mental fatigue sensation. After their subjective level of fatigue had recovered to the baseline level, they re-experienced the fatigue sensation experienced in the two-back test positively, negatively, and without any modification (i.e., re-experienced the fatigue sensation as it was). The level of motivation to perform another two-back test they felt during the re-experiencing was assessed. The neural activity related to the re-experiencing was recorded using magnetoencephalography. The level of the motivation to perform another two-back test was increased by positively re-experiencing the fatigue sensation. The increase in delta band power in Brodmann area 7 was positively associated with the increase in motivation. These results show that positive thinking about fatigue sensation can enhance motivation and suggest that this enhanced motivation may have some effects on visual attention system.

The Hierarchy of Brain Networks Is Related to Insulin Growth Factor-1 in a Large, Middle-Aged, Healthy Cohort: An Exploratory Magnetoencephalography Study

Recently, a large study demonstrated that lower serum levels of insulin growth factor-1 (IGF-1) relate to brain atrophy and to a greater risk for developing Alzheimer's disease in a healthy elderly population. We set out to test if functional brain networks relate to IGF-1 levels in the middle aged. Hence, we studied the association between IGF-1 and magnetoencephalography-based functional network characteristics in a middle-aged population. The functional connections between brain areas were estimated for six frequency bands (delta, theta, alpha1, alpha2, beta, gamma) using the phase lag index. Subsequently, the topology of the frequency-specific functional networks was characterized using the minimum spanning tree. Our results showed that lower levels of serum IGF-1 relate to a globally less integrated functional network in the beta and theta band. The associations remained significant when correcting for gender and systemic effects of IGF-1 that might indirectly affect the brain. The value of this exploratory study is the demonstration that lower levels of IGF-1 are associated with brain network topology in the middle aged.

The lifespan trajectory of neural oscillatory activity in the motor system

•

The lifespan trajectory of resting and motor-related beta oscillations is unknown.

•

These beta dynamics were examined in participants aged 9–75 years using MEG imaging.

•

Resting beta levels and motor-related beta oscillations follow unique trajectories.

•

The dynamic relationship between these two measures predicts motor performance.

Numerous studies connect beta oscillations in the motor cortices to volitional movement, and beta is known to be aberrant in multiple movement disorders. However, the dynamic interplay between these beta oscillations, motor performance, and spontaneous beta power (e.g., during rest) in the motor cortices remains unknown. This study utilized magnetoencephalography (MEG) to investigate these three parameters and their lifespan trajectory in 57 healthy participants aged 9–75 years old. Movement-related beta activity was imaged using a beamforming approach, and voxel time series data were extracted from the peak voxels in the primary motor cortices. Our results indicated that spontaneous beta power during rest followed a quadratic lifespan trajectory, while movement-related beta oscillations linearly increased with age. Follow-on analyses showed that spontaneous beta power and the beta minima during movement, together, significantly predicted task performance above and beyond the effects of age. These data are the first to show lifespan trajectories among measures of beta activity in the motor cortices, and suggest that the healthy brain compensates for age-related increases in spontaneous beta activity by increasing the strength of beta oscillations within the motor cortices which, when successful, enables normal motor performance into later life.

Abnormal language-related oscillatory responses in primary progressive aphasia

Patients with Primary Progressive Aphasia (PPA) may react to linguistic stimuli differently than healthy controls, reflecting degeneration of language networks and engagement of compensatory mechanisms. We used magnetoencephalography (MEG) to evaluate oscillatory neural responses in sentence comprehension, in patients with PPA and age-matched controls. Participants viewed sentences containing semantically and syntactically anomalous words that evoke distinct oscillatory responses. For age-matched controls, semantic anomalies elicited left-lateralized 8–30 Hz power decreases distributed along ventral brain regions, whereas syntactic anomalies elicited bilateral power decreases in both ventral and dorsal regions. In comparison to controls, patients with PPA showed altered patterns of induced oscillations, characterized by delayed latencies and attenuated amplitude, which were correlated with linguistic impairment measured offline. The recruitment of right hemisphere temporo-parietal areas (also found in controls) was correlated with preserved semantic processing abilities, indicating that preserved neural activity in these regions was able to support successful semantic processing. In contrast, syntactic processing was more consistently impaired in PPA, regardless of neural activity patterns, suggesting that this domain of language is particularly vulnerable to the neuronal loss. In addition, we found that delayed peak latencies of oscillatory responses were associated with lower accuracy for detecting semantic anomalies, suggesting that language deficits observed in PPA may be linked to delayed or slowed information processing.

•

Evaluated induced oscillations in patients with PPA using MEG.

•

PPA patients showed delayed latencies and attenuated amplitude of responses.

•

Preserved right hemisphere regions support semantic processing.

•

Delayed latencies of oscillatory responses associated with impaired performance.

•

Language deficits in PPA linked to delayed or slowed information processing.

Neural activity induced by visual food stimuli presented out of awareness: a preliminary magnetoencephalography study

Obesity is a major public health problem in modern society. Appetitive behavior has been proposed to be partially driven by unconscious decision-making processes and thus, targeting the unconscious cognitive processes related to eating behavior is essential to develop strategies for overweight individuals and obese patients. Here, we presented food pictures below the threshold of awareness to healthy male volunteers and examined neural activity related to appetitive behavior using magnetoencephalography. We found that, among participants who did not recognize food pictures during the experiment, an index of heart rate variability assessed by electrocardiography (low-frequency component power/high-frequency component power ratio, LF/HF) just after picture presentation was increased compared with that just before presentation, and the increase in LF/HF was negatively associated with the score for cognitive restraint of food intake. In addition, increased LF/HF was negatively associated with increased alpha band power in Brodmann area (BA) 47 caused by food pictures presented below the threshold of awareness, and level of cognitive restraint was positively associated with increased alpha band power in BA13. Our findings may provide valuable clues to the development of methods assessing unconscious regulation of appetite and offer avenues for further study of the neural mechanisms related to eating behavior.

Oscillatory dynamics in the dorsal and ventral attention networks during the reorienting of attention

The ability to reorient attention within the visual field is central to daily functioning, and numerous fMRI studies have shown that the dorsal and ventral attention networks (DAN, VAN) are critical to such processes. However, despite the instantaneous nature of attentional shifts, the dynamics of oscillatory activity serving attentional reorientation remain poorly characterized. In this study, we utilized magnetoencephalography (MEG) and a Posner task to probe the dynamics of attentional reorienting in 29 healthy adults. MEG data were transformed into the time-frequency domain and significant oscillatory responses were imaged using a beamformer. Voxel time series were then extracted from peak voxels in the functional beamformer images. These time series were used to quantify the dynamics of attentional reorienting, and to compute dynamic functional connectivity. Our results indicated strong increases in theta and decreases in alpha and beta activity across many nodes in the DAN and VAN. Interestingly, theta responses were generally stronger during trials that required attentional reorienting relative to those that did not, while alpha and beta oscillations were more dynamic, with many regions exhibiting significantly stronger responses during non-reorienting trials initially, and the opposite pattern during later processing. Finally, stronger functional connectivity was found following target presentation (575-700 ms) between bilateral superior parietal lobules during attentional reorienting. In sum, these data show that visual attention is served by multiple cortical regions within the DAN and VAN, and that attentional reorienting processes are often associated with spectrally-specific oscillations that have largely distinct spatiotemporal dynamics.

Consistency of magnetoencephalographic functional connectivity and network reconstruction using a template versus native MRI for co-registration

INTRODUCTION

Studies using functional connectivity and network analyses based on magnetoencephalography (MEG) with source localization are rapidly emerging in neuroscientific literature. However, these analyses currently depend on the availability of costly and sometimes burdensome individual MR scans for co-registration. We evaluated the consistency of these measures when using a template MRI, instead of native MRI, for the analysis of functional connectivity and network topology.

METHODS

Seventeen healthy participants underwent resting-state eyes-closed MEG and anatomical MRI. These data were projected into source space using an atlas-based peak voxel and a centroid beamforming approach either using (1) participants' native MRIs or (2) the Montreal Neurological Institute's template. For both methods, time series were reconstructed from 78 cortical atlas regions. Relative power was determined in six classical frequency bands per region and globally averaged. Functional connectivity (phase lag index) between each pair of regions was calculated. The adjacency matrices were then used to reconstruct functional networks, of which regional and global metrics were determined. Intraclass correlation coefficients were calculated and Bland-Altman plots were made to quantify the consistency and potential bias of the use of template versus native MRI.

RESULTS

Co-registration with the template yielded largely consistent relative power, connectivity, and network estimates compared to native MRI.

DISCUSSION

These findings indicate that there is no (systematic) bias or inconsistency between template and native MRI co-registration of MEG. They open up possibilities for retrospective and prospective analyses to MEG datasets in the general population that have no native MRIs available. Hum Brain Mapp, 2017. © 2017 The Authors Human Brain Mapping Published by Wiley Periodicals, Inc. Hum Brain Mapp 39:104-119, 2018. © 2017 Wiley Periodicals, Inc.

Contrasting functional imaging parametric maps: The mislocation problem and alternative solutions

In the field of neuroimaging, researchers often resort to contrasting parametric maps to identify differences between conditions or populations. Unfortunately, contrast patterns mix effects related to amplitude and location differences and tend to peak away from sources of genuine brain activity to an extent that scales with the smoothness of the maps. Here, we illustrate this mislocation problem on source maps reconstructed from magnetoencephalographic recordings and propose a novel, dedicated location-comparison method. In realistic simulations, contrast mislocation was on average ∼10 mm when genuine sources were placed at the same location, and was still above 5 mm when sources were 20 mm apart. The dedicated location-comparison method achieved a sensitivity of ∼90% when inter-source distance was 12 mm. Its benefit is also illustrated on real brain-speech entrainment data. In conclusion, contrasts of parametric maps provide precarious information for source location. To specifically address the question of location difference, one should turn to dedicated methods as the one proposed here.

Neurofeedback for Tinnitus Treatment - Review and Current Concepts

An effective treatment to completely alleviate chronic tinnitus symptoms has not yet been discovered. However, recent developments suggest that neurofeedback (NFB), a method already popular in the treatment of other psychological and neurological disorders, may provide a suitable alternative. NFB is a non-invasive method generally based on electrophysiological recordings and visualizing of certain aspects of brain activity as positive or negative feedback that enables patients to voluntarily control their brain activity and thus triggers them to unlearn typical neural activity patterns related to tinnitus. The purpose of this review is to summarize and discuss previous findings of neurofeedback treatment studies in the field of chronic tinnitus. In doing so, also an overview about the underlying theories of tinnitus emergence is presented and results of resting-state EEG and MEG studies summarized and critically discussed. To date, neurofeedback as well as electrophysiological tinnitus studies lack general guidelines that are crucial to produce more comparable and consistent results. Even though neurofeedback has already shown promising results for chronic tinnitus treatment, further research is needed in order to develop more sophisticated protocols that are able to tackle the individual needs of tinnitus patients more specifically.

Evidence for unconscious regulation of performance in fatigue

Since fatigue is prevalent in modern societies, it is necessary to clarify the neural mechanisms of fatigue. The regulation of performance through fatigue sensation is one of the mechanisms that decreases performance in fatigue. However, it is unknown whether subjective feeling of fatigue is necessary for the regulation of performance. Here, we examined whether decreased performance occurs without increased fatigue sensation through the experiment which was designed to test if fatigue can be learned unconsciously. Healthy male volunteers performed a fatigue-inducing hand-grip task for 10 min while viewing a target image presented without awareness. On the next day, they viewed a control and the target images presented with awareness and the neural activity caused by viewing the images was measured using magnetoencephalography. Results showed the level of fatigue sensation was not altered but grip-strength was decreased by viewing the target image on the second day. The level of beta band power in Brodmann's area 31 was increased by viewing the target image and this increase was negatively associated with the decrease of grip-strength caused in the hand-grip task. These findings demonstrated that fatigue can be learned unconsciously and that there is a mechanism to decrease performance without fatigue sensation.

Modulation of Neural Oscillatory Activity during Dynamic Face Processing

Various neuroimaging and neurophysiological methods have been used to examine neural activation patterns in response to faces. However, much of previous research has relied on static images of faces, which do not allow a complete description of the temporal structure of face-specific neural activities to be made. More recently, insights are emerging from fMRI studies about the neural substrates that underpin our perception of naturalistic dynamic face stimuli, but the temporal and spectral oscillatory activity associated with processing dynamic faces has yet to be fully characterized. Here, we used MEG and beamformer source localization to examine the spatiotemporal profile of neurophysiological oscillatory activity in response to dynamic faces. Source analysis revealed a number of regions showing enhanced activation in response to dynamic relative to static faces in the distributed face network, which were spatially coincident with regions that were previously identified with fMRI. Furthermore, our results demonstrate that perception of realistic dynamic facial stimuli activates a distributed neural network at varying time points facilitated by modulations in low-frequency power within alpha and beta frequency ranges (8-30 Hz). Naturalistic dynamic face stimuli may provide a better means of representing the complex nature of perceiving facial expressions in the real world, and neural oscillatory activity can provide additional insights into the associated neural processes.

The Effect of Head Model Simplification on Beamformer Source Localization

Beamformers are a widely-used tool in brain analysis with magnetoencephalography (MEG) and electroencephalography (EEG). For the construction of the beamformer filters realistic head volume conductor modeling is necessary for accurately computing the EEG and MEG leadfields, i.e., for solving the EEG and MEG forward problem. In this work, we investigate the influence of including realistic head tissue compartments into a finite element method (FEM) model on the beamformer's localization ability. Specifically, we investigate the effect of including cerebrospinal fluid, gray matter, and white matter distinction, as well as segmenting the skull bone into compacta and spongiosa, and modeling white matter anisotropy. We simulate an interictal epileptic measurement with white sensor noise. Beamformer filters are constructed with unit gain, unit array gain, and unit noise gain constraint. Beamformer source positions are determined by evaluating power and excess sample kurtosis (g2) of the source-waveforms at all source space nodes. For both modalities, we see a strong effect of modeling the cerebrospinal fluid and white and gray matter. Depending on the source position, both effects can each be in the magnitude of centimeters, rendering their modeling necessary for successful localization. Precise skull modeling mainly effected the EEG up to a few millimeters, while both modalities could profit from modeling white matter anisotropy to a smaller extent of 5-10 mm. The unit noise gain or neural activity index beamformer behaves similarly to the array gain beamformer when noise strength is sufficiently high. Variance localization seems more robust against modeling errors than kurtosis.

Neural basis of individual differences in the response to mental stress: a magnetoencephalography study

Stress is a risk factor for the onset of mental disorders. Although stress response varies across individuals, the mechanism of individual differences remains unclear. Here, we investigated the neural basis of individual differences in response to mental stress using magnetoencephalography (MEG). Twenty healthy male volunteers completed the Temperament and Character Inventory (TCI). The experiment included two types of tasks: a non-stress-inducing task and a stress-inducing task. During these tasks, participants passively viewed non-stress-inducing images and stress-inducing images, respectively, and MEG was recorded. Before and after each task, MEG and electrocardiography were recorded and subjective ratings were obtained. We grouped participants according to Novelty seeking (NS) - tendency to be exploratory, and Harm avoidance (HA) - tendency to be cautious. Participants with high NS and low HA (n = 10) assessed by TCI had a different neural response to stress than those with low NS and high HA (n = 10). Event-related desynchronization (ERD) in the beta frequency band was observed only in participants with high NS and low HA in the brain region extending from Brodmann's area 31 (including the posterior cingulate cortex and precuneus) from 200 to 350 ms after the onset of picture presentation in the stress-inducing task. Individual variation in personality traits (NS and HA) was associated with the neural response to mental stress. These findings increase our understanding of the psychological and neural basis of individual differences in the stress response, and will contribute to development of the psychotherapeutic approaches to stress-related disorders.

Spontaneous Blinks Activate the Precuneus: Characterizing Blink-Related Oscillations Using Magnetoencephalography

Spontaneous blinking occurs 15–20 times per minute. Although blinking has often been associated with its physiological role of corneal lubrication, there is now increasing behavioral evidence suggesting that blinks are also modulated by cognitive processes such as attention and information processing. Recent low-density electroencephalography (EEG) studies have reported so-called blink-related oscillations (BROs) associated with spontaneous blinking at rest. Delta-band (0.5–4 Hz) BROs are thought to originate from the precuneus region involved in environmental monitoring and awareness, with potential clinical utility in evaluation of disorders of consciousness. However, the neural mechanisms of BROs have not been elucidated. Using magnetoencephalography (MEG), we characterized delta-band BROs in 36 healthy individuals while controlling for background brain activity. Results showed that, compared to pre-blink baseline, delta-band BROs resulted in increased global field power (p < 0.001) and time-frequency spectral power (p < 0.05) at the sensor level, peaking at ~250 ms post-blink maximum. Source localization showed that spontaneous blinks activated the bilateral precuneus (p < 0.05 FWE), and source activity within the precuneus was also consistent with sensor-space results. Crucially, these effects were only observed in the blink condition and were absent in the control condition, demonstrating that results were due to spontaneous blinks rather than as part of the inherent brain activity. The current study represents the first MEG examination of BROs. Our findings suggest that spontaneous blinks activate the precuneus regions consistent with environmental monitoring and awareness, and provide important neuroimaging support for the cognitive role of spontaneous blinks.

A mean field model for movement induced changes in the beta rhythm

In electrophysiological recordings of the brain, the transition from high amplitude to low amplitude signals are most likely caused by a change in the synchrony of underlying neuronal population firing patterns. Classic examples of such modulations are the strong stimulus-related oscillatory phenomena known as the movement related beta decrease (MRBD) and post-movement beta rebound (PMBR). A sharp decrease in neural oscillatory power is observed during movement (MRBD) followed by an increase above baseline on movement cessation (PMBR). MRBD and PMBR represent important neuroscientific phenomena which have been shown to have clinical relevance. Here, we present a parsimonious model for the dynamics of synchrony within a synaptically coupled spiking network that is able to replicate a human MEG power spectrogram showing the evolution from MRBD to PMBR. Importantly, the high-dimensional spiking model has an exact mean field description in terms of four ordinary differential equations that allows considerable insight to be obtained into the cause of the experimentally observed time-lag from movement termination to the onset of PMBR (∼ 0.5 s), as well as the subsequent long duration of PMBR (∼ 1 - 10 s). Our model represents the first to predict these commonly observed and robust phenomena and represents a key step in their understanding, in health and disease.

The peak frequency of motor-related gamma oscillations is modulated by response competition

Movement execution generally occurs in an environment with numerous distractors, and requires the selection of a motor plan from multiple possible alternatives. However, the impact of such distractors on cortical motor function during movement remains largely unknown. Previous studies have identified two movement-related oscillatory responses that are critical to motor planning and execution, and these responses include the peri-movement beta event-related desynchronization (ERD) and the movement-related gamma synchronization (MRGS). In the current study, we investigate how visual distractors cuing alternative movements modulate the beta ERD and MRGS responses. To this end, we recorded magnetoencephalography (MEG) during an arrow-based version of the Eriksen flanker task in 42 healthy adults. All MEG data were transformed in to the time-frequency domain and the beta ERD and MRGS responses were imaged using a beamformer. Virtual sensors (voxel time series) were then extracted from the peak voxels of each response for the congruent and incongruent flanker conditions separately, and these data were examined for conditional differences during the movement. Our results indicated that participants exhibited the classic "flanker effect," as they responded significantly slower during incongruent relative to congruent trials. Our most important MEG finding was a significant increase in the peak frequency of the MRGS in the incongruent compared to the congruent condition, with no conditional effect on response amplitude. In addition, we found significantly stronger peri-movement beta ERD responses in the ipsilateral motor cortex during incongruent compared to congruent trials, but no conditional effect on frequency. These data are the first to show that the peak frequency of the MRGS response is linked to the task parameters, and varies from trial to trial in individual participants. More globally, these data suggest that beta and gamma oscillations are modulated by visual distractors causing response competition.

Attention training improves aberrant neural dynamics during working memory processing in veterans with PTSD

Posttraumatic stress disorder (PTSD) is associated with executive functioning deficits, including disruptions in working memory (WM). Recent studies suggest that attention training reduces PTSD symptomatology, but the underlying neural mechanisms are unknown. We used high-density magnetoencephalography (MEG) to evaluate whether attention training modulates brain regions serving WM processing in PTSD. Fourteen veterans with PTSD completed a WM task during a 306-sensor MEG recording before and after 8 sessions of attention training treatment. A matched comparison sample of 12 combat-exposed veterans without PTSD completed the same WM task during a single MEG session. To identify the spatiotemporal dynamics, each group's data were transformed into the time-frequency domain, and significant oscillatory brain responses were imaged using a beamforming approach. All participants exhibited activity in left hemispheric language areas consistent with a verbal WM task. Additionally, veterans with PTSD and combat-exposed healthy controls each exhibited oscillatory responses in right hemispheric homologue regions (e.g., right Broca's area); however, these responses were in opposite directions. Group differences in oscillatory activity emerged in the theta band (4-8 Hz) during encoding and in the alpha band (9-12 Hz) during maintenance and were significant in right prefrontal and right supramarginal and inferior parietal regions. Importantly, following attention training, these significant group differences were reduced or eliminated. This study provides initial evidence that attention training improves aberrant neural activity in brain networks serving WM processing.

Deep Source Localization with Magnetoencephalography Based on Sensor Array Decomposition and Beamforming

In recent years, the source localization technique of magnetoencephalography (MEG) has played a prominent role in cognitive neuroscience and in the diagnosis and treatment of neurological and psychological disorders. However, locating deep brain activities such as in the mesial temporal structures, especially in preoperative evaluation of epilepsy patients, may be more challenging. In this work we have proposed a modified beamforming approach for finding deep sources. First, an iterative spatiotemporal signal decomposition was employed for reconstructing the sensor arrays, which could characterize the intrinsic discriminant features for interpreting sensor signals. Next, a sensor covariance matrix was estimated under the new reconstructed space. Then, a well-known vector beamforming approach, which was a linearly constraint minimum variance (LCMV) approach, was applied to compute the solution for the inverse problem. It can be shown that the proposed source localization approach can give better localization accuracy than two other commonly-used beamforming methods (LCMV, MUSIC) in simulated MEG measurements generated with deep sources. Further, we applied the proposed approach to real MEG data recorded from ten patients with medically-refractory mesial temporal lobe epilepsy (mTLE) for finding epileptogenic zone(s), and there was a good agreement between those findings by the proposed approach and the clinical comprehensive results.

Veterans with post-traumatic stress disorder exhibit altered emotional processing and attentional control during an emotional Stroop task

BACKGROUND

Post-traumatic stress disorder (PTSD) is often associated with attention allocation and emotional regulation difficulties, but the brain dynamics underlying these deficits are unknown. The emotional Stroop task (EST) is an ideal means to monitor these difficulties, because participants are asked to attend to non-emotional aspects of the stimuli. In this study, we used magnetoencephalography (MEG) and the EST to monitor attention allocation and emotional regulation during the processing of emotionally charged stimuli in combat veterans with and without PTSD.

METHOD

A total of 31 veterans with PTSD and 20 without PTSD performed the EST during MEG. Three categories of stimuli were used, including combat-related, generally threatening and neutral words. MEG data were imaged in the time-frequency domain and the network dynamics were probed for differences in processing threatening and non-threatening words.

RESULTS

Behaviorally, veterans with PTSD were significantly slower in responding to combat-related relative to neutral and generally threatening words. Veterans without PTSD exhibited no significant differences in responding to the three different word types. Neurophysiologically, we found a significant three-way interaction between group, word type and time period across multiple brain regions. Follow-up testing indicated stronger theta-frequency (4-8 Hz) responses in the right ventral prefrontal (0.4-0.8 s) and superior temporal cortices (0.6-0.8 s) of veterans without PTSD compared with those with PTSD during the processing of combat-related words.

CONCLUSIONS

Our data indicated that veterans with PTSD exhibited deficits in attention allocation and emotional regulation when processing trauma cues, while those without PTSD were able to regulate emotion by directing attention away from threat.

The heritability of multi-modal connectivity in human brain activity

Patterns of intrinsic human brain activity exhibit a profile of functional connectivity that is associated with behaviour and cognitive performance, and deteriorates with disease. This paper investigates the relative importance of genetic factors and the common environment between twins in determining this functional connectivity profile. Using functional magnetic resonance imaging (fMRI) on 820 subjects from the Human Connectome Project, and magnetoencephalographic (MEG) recordings from a subset, the heritability of connectivity among 39 cortical regions was estimated. On average over all connections, genes account for about 15% of the observed variance in fMRI connectivity (and about 10% in alpha-band and 20% in beta-band oscillatory power synchronisation), which substantially exceeds the contribution from the environment shared between twins. Therefore, insofar as twins share a common upbringing, it appears that genes, rather than the developmental environment, have the dominant role in determining the coupling of neuronal activity.

Directional information flow in patients with Alzheimer's disease. A source-space resting-state MEG study

In a recent magnetoencephalography (MEG) study, we found posterior-to-anterior information flow over the cortex in higher frequency bands in healthy subjects, with a reversed pattern in the theta band. A disruption of information flow may underlie clinical symptoms in Alzheimer's disease (AD). In AD, highly connected regions (hubs) in posterior areas are mostly disrupted. We therefore hypothesized that in AD the information flow from these hub regions would be disturbed. We used resting-state MEG recordings from 27 early-onset AD patients and 26 healthy controls. Using beamformer-based virtual electrodes, we estimated neuronal oscillatory activity for 78 cortical regions of interest (ROIs) and 12 subcortical ROIs of the AAL atlas, and calculated the directed phase transfer entropy (dPTE) as a measure of information flow between these ROIs. Group differences were evaluated using permutation tests and, for the AD group, associations between dPTE and general cognition or CSF biomarkers were determined using Spearman correlation coefficients. We confirmed the previously reported posterior-to-anterior information flow in the higher frequency bands in the healthy controls, and found it to be disturbed in the beta band in AD. Most prominently, the information flow from the precuneus and the visual cortex, towards frontal and subcortical structures, was decreased in AD. These disruptions did not correlate with cognitive impairment or CSF biomarkers. We conclude that AD pathology may affect the flow of information between brain regions, particularly from posterior hub regions, and that changes in the information flow in the beta band indicate an aspect of the pathophysiological process in AD.

Automatic detection and visualisation of MEG ripple oscillations in epilepsy

High frequency oscillations (HFOs, 80–500 Hz) in invasive EEG are a biomarker for the epileptic focus. Ripples (80–250 Hz) have also been identified in non-invasive MEG, yet detection is impeded by noise, their low occurrence rates, and the workload of visual analysis. We propose a method that identifies ripples in MEG through noise reduction, beamforming and automatic detection with minimal user effort. We analysed 15 min of presurgical resting-state interictal MEG data of 25 patients with epilepsy. The MEG signal-to-noise was improved by using a cross-validation signal space separation method, and by calculating ~ 2400 beamformer-based virtual sensors in the grey matter. Ripples in these sensors were automatically detected by an algorithm optimized for MEG. A small subset of the identified ripples was visually checked. Ripple locations were compared with MEG spike dipole locations and the resection area if available. Running the automatic detection algorithm resulted in on average 905 ripples per patient, of which on average 148 ripples were visually reviewed. Reviewing took approximately 5 min per patient, and identified ripples in 16 out of 25 patients. In 14 patients the ripple locations showed good or moderate concordance with the MEG spikes. For six out of eight patients who had surgery, the ripple locations showed concordance with the resection area: 4/5 with good outcome and 2/3 with poor outcome. Automatic ripple detection in beamformer-based virtual sensors is a feasible non-invasive tool for the identification of ripples in MEG. Our method requires minimal user effort and is easily applicable in a clinical setting.

•

Cross-validation signal space separation and beamformer increase the SNR in MEG.

•

Automatic detection of MEG ripples in the time domain is feasible.

•

Our method identifies ripples with minimal user effort and is clinically applicable.

•

Automatically detected ripples are concordant with MEG spikes in 14/16 patients.

•

Automatically detected ripples are concordant with resection area in 6/8 patients.

Heightened amygdala responsiveness in s-carriers of 5-HTTLPR genetic polymorphism reflects enhanced cortical rather than subcortical inputs: An MEG study

Short allele carriers (S-carriers) of the serotonin transporter gene (5-HTTLPR) show an elevated amygdala response to emotional stimuli relative to long allele carriers (LL-homozygous). However, whether this reflects increased responsiveness of the amygdala generally or interactions between the amygdala and the specific input systems remains unknown. It is argued that the amygdala receives input via a quick subcortical and a slower cortical pathway. If the elevated amygdala response in S-carriers reflects generally increased amygdala responding, then group differences in amygdala should be seen across the amygdala response time course. However, if the difference is a secondary consequence of enhanced amygdala-cortical interactions, then group differences might only be present later in the amygdala response. Using magnetoencephalography (MEG), we found an enhanced amygdala response to fearful expressions starting 40-50 ms poststimulus. However, group differences in the amygdala were only seen 190-200 ms poststimulus, preceded by increased superior temporal sulcus (STS) responses in S-carriers from 130 to 140 ms poststimulus. An enhanced amygdala response to angry expressions started 260-270 ms poststimulus with group differences in the amygdala starting at 160-170 ms poststimulus onset, preceded by increased STS responses in S-carriers from 150 to 160 ms poststimulus. These suggest that enhanced amygdala responses in S-carriers might reflect enhanced STS-amygdala connectivity in S-carriers. Hum Brain Mapp 38:4313-4321, 2017. © 2017 Wiley Periodicals, Inc.

Large-scale network dynamics of beta-band oscillations underlie auditory perceptual decision-making

Perceptual decisions vary in the speed at which we make them. Evidence suggests that translating sensory information into perceptual decisions relies on distributed interacting neural populations, with decision speed hinging on power modulations of the neural oscillations. Yet the dependence of perceptual decisions on the large-scale network organization of coupled neural oscillations has remained elusive. We measured magnetoencephalographic signals in human listeners who judged acoustic stimuli composed of carefully titrated clouds of tone sweeps. These stimuli were used in two task contexts, in which the participants judged the overall pitch or direction of the tone sweeps. We traced the large-scale network dynamics of the source-projected neural oscillations on a trial-by-trial basis using power-envelope correlations and graph-theoretical network discovery. In both tasks, faster decisions were predicted by higher segregation and lower integration of coupled beta-band (∼16-28 Hz) oscillations. We also uncovered the brain network states that promoted faster decisions in either lower-order auditory or higher-order control brain areas. Specifically, decision speed in judging the tone sweep direction critically relied on the nodal network configurations of anterior temporal, cingulate, and middle frontal cortices. Our findings suggest that global network communication during perceptual decision-making is implemented in the human brain by large-scale couplings between beta-band neural oscillations.

Children with cerebral palsy have altered oscillatory activity in the motor and visual cortices during a knee motor task

The neuroimaging literature on cerebral palsy (CP) has predominantly focused on identifying structural aberrations within the white matter (e.g., fiber track integrity), with very few studies examining neural activity within the key networks that serve the production of motor actions. The current investigation used high-density magnetoencephalography to begin to fill this knowledge gap by quantifying the temporal dynamics of the alpha and beta cortical oscillations in children with CP (age = 15.5 ± 3 years; GMFCS levels II–III) and typically developing (TD) children (age = 14.1 ± 3 years) during a goal-directed isometric target-matching task using the knee joint. Advanced beamforming methods were used to image the cortical oscillations during the movement planning and execution stages. Compared with the TD children, our results showed that the children with CP had stronger alpha and beta event-related desynchronization (ERD) within the primary motor cortices, premotor area, inferior parietal lobule, and inferior frontal gyrus during the motor planning stage. Differences in beta ERD amplitude extended through the motor execution stage within the supplementary motor area and premotor cortices, and a stronger alpha ERD was detected in the anterior cingulate. Interestingly, our results also indicated that alpha and beta oscillations were weaker in the children with CP within the occipital cortices and visual MT area during movement execution. These altered alpha and beta oscillations were accompanied by slower reaction times and substantial target matching errors in the children with CP. We also identified that the strength of the alpha and beta ERDs during the motor planning and execution stages were correlated with the motor performance. Lastly, our regression analyses suggested that the beta ERD within visual areas during motor execution primarily predicted the amount of motor errors. Overall, these data suggest that uncharacteristic alpha and beta oscillations within visuomotor cortical networks play a prominent role in the atypical motor actions exhibited by children with CP.

•

Children with CP performed an isometric task with the knee joint.

•

Children with CP had stronger alpha and beta ERD during motor planning.

•

These ERD differences extended through the motor execution period.

•

Occipital cortices and visual MT area alpha and beta ERD were weaker.

•

Altered alpha and beta ERD were accompanied by impaired motor actions.

Changes in resting-state directed connectivity in cortico-subcortical networks correlate with cognitive function in Parkinson's disease

OBJECTIVE

The pathophysiological mechanisms underlying Parkinson's disease (PD)-related cognitive decline and conversion to PD dementia are poorly understood. In the healthy human brain, stable patterns of posterior-to-anterior cortical information flow have recently been demonstrated in the higher frequency bands using magnetoencephalography (MEG). In this study we estimated PD-related changes in information flow patterns, as well as the contribution of subcortical regions.

METHODS

Resting-state MEG recordings were acquired in moderately advanced PD patients (n=34; mean Hoehn and Yahr-stage 2.5) and healthy controls (n=12). MEG signals were projected to both cortical and subcortical brain regions, following which we estimated the balance between incoming and outgoing information flow per region.

RESULTS

In PD patients, compared to controls, preferential beta band information outflow was significantly higher for the basal ganglia and frontotemporal cortical regions, and significantly lower for parieto-occipital regions. In addition, in patients, low preferential information outflow from occipital regions correlated with poor global cognitive performance.

CONCLUSION

In the PD brain, a shift in balance towards more anterior-to-posterior beta band information flow takes place and is associated with poorer cognitive performance.

SIGNIFICANCE

Our results indicate that a reversal of the physiological posterior-to-anterior information flow may be an important mechanism in PD-related cognitive decline.

Spatiotemporal oscillatory dynamics of visual selective attention during a flanker task

The flanker task is a test of visual selective attention that has been widely used to probe error monitoring, response conflict, and related constructs. However, to date, few studies have focused on the selective attention component of this task and imaged the underlying oscillatory dynamics serving task performance. In this study, 21 healthy adults successfully completed an arrow-based version of the Eriksen flanker task during magnetoencephalography (MEG). All MEG data were pre-processed and transformed into the time-frequency domain. Significant oscillatory brain responses were imaged using a beamforming approach, and voxel time series were extracted from the peak responses to identify the temporal dynamics. Across both congruent and incongruent flanker conditions, our results indicated robust decreases in alpha (9-12Hz) activity in medial and lateral occipital regions, bilateral parietal cortices, and cerebellar areas during task performance. In parallel, increases in theta (3-7Hz) oscillatory activity were detected in dorsal and ventral frontal regions, and the anterior cingulate. As per conditional effects, stronger alpha responses (i.e., greater desynchronization) were observed in parietal, occipital, and cerebellar cortices during incongruent relative to congruent trials, whereas the opposite pattern emerged for theta responses (i.e., synchronization) in the anterior cingulate, left dorsolateral prefrontal, and ventral prefrontal cortices. Interestingly, the peak latency of theta responses in these latter brain regions was significantly correlated with reaction time, and may partially explain the amplitude difference observed between congruent and incongruent trials. Lastly, whole-brain exploratory analyses implicated the frontal eye fields, right temporoparietal junction, and premotor cortices. These findings suggest that regions of both the dorsal and ventral attention networks contribute to visual selective attention processes during incongruent trials, and that such differential processes are transient and fully completed shortly after the behavioral response in most trials.

Altered sensorimotor cortical oscillations in individuals with multiple sclerosis suggests a faulty internal model

Multiple sclerosis (MS) is a demyelinating disease that results in a broad array of symptoms, including impaired motor performance. How such demyelination of fibers affects the inherent neurophysiological activity in motor circuits, however, remains largely unknown. Potentially, the movement errors associated with MS may be due to imperfections in the internal model used to make predictions of the motor output that will meet the task demands. Prior magnetoencephalographic (MEG) and electroencephalographic brain imaging experiments have established that the beta (15-30 Hz) oscillatory activity in the sensorimotor cortices is related to the control of movement. Specifically, it has been suggested that the strength of the post-movement beta rebound may indicate the certainty of the internal model. In this study, we used MEG to evaluate the neural oscillatory activity in the sensorimotor cortices of individuals with MS and healthy individuals during a goal-directed isometric knee force task. Our results showed no difference between the individuals with MS and healthy individuals in the beta activity during the planning and execution stages of movement. However, we did find that individuals with MS exhibited a weaker post-movement beta rebound in the pre/postcentral gyri relative to healthy controls. Additionally, we found that the behavioral performance of individuals with MS was aberrant, and related to the strength of the post-movement beta rebound. These results suggest that the internal model may be faulty in individuals with MS. Hum Brain Mapp 38:4009-4018, 2017. © 2017 Wiley Periodicals, Inc.

Shared Neural Mechanisms for the Prediction of Own and Partner Musical Sequences after Short-term Piano Duet Training

Predictive mechanisms in the human brain can be investigated using markers for prediction violations like the mismatch negativity (MMN). Short-term piano training increases the MMN for melodic and rhythmic deviations in the training material. This increase occurs only when the material is actually played, not when it is only perceived through listening, suggesting that learning predictions about upcoming musical events are derived from motor involvement. However, music is often performed in concert with others. In this case, predictions about upcoming actions from a partner are a crucial part of the performance. In the present experiment, we use magnetoencephalography (MEG) to measure MMNs to deviations in one's own and a partner's musical material after both engaged in musical duet training. Event-related field (ERF) results revealed that the MMN increased significantly for own and partner material suggesting a neural representation of the partner's part in a duet situation. Source analysis using beamforming revealed common activations in auditory, inferior frontal, and parietal areas, similar to previous results for single players, but also a pronounced contribution from the cerebellum. In addition, activation of the precuneus and the medial frontal cortex was observed, presumably related to the need to distinguish between own and partner material.