Distributed current estimates using cortical orientation constraints

Decoding kinematic information from beta-band motor rhythms of speech motor cortex: a methodological/analytic approach using concurrent speech movement tracking and magnetoencephalography

Introduction

Articulography and functional neuroimaging are two major tools for studying the neurobiology of speech production. Until now, however, it has generally not been feasible to use both in the same experimental setup because of technical incompatibilities between the two methodologies.

Methods

Here we describe results from a novel articulography system dubbed Magneto-articulography for the Assessment of Speech Kinematics (MASK), which is technically compatible with magnetoencephalography (MEG) brain scanning systems. In the present paper we describe our methodological and analytic approach for extracting brain motor activities related to key kinematic and coordination event parameters derived from time-registered MASK tracking measurements. Data were collected from 10 healthy adults with tracking coils on the tongue, lips, and jaw. Analyses targeted the gestural landmarks of reiterated utterances/ipa/ and /api/, produced at normal and faster rates.

Results

The results show that (1) Speech sensorimotor cortex can be reliably located in peri-rolandic regions of the left hemisphere; (2) mu (8-12 Hz) and beta band (13-30 Hz) neuromotor oscillations are present in the speech signals and contain information structures that are independent of those present in higher-frequency bands; and (3) hypotheses concerning the information content of speech motor rhythms can be systematically evaluated with multivariate pattern analytic techniques.

Discussion

These results show that MASK provides the capability, for deriving subject-specific articulatory parameters, based on well-established and robust motor control parameters, in the same experimental setup as the MEG brain recordings and in temporal and spatial co-register with the brain data. The analytic approach described here provides new capabilities for testing hypotheses concerning the types of kinematic information that are encoded and processed within specific components of the speech neuromotor system.

Characterizing directional dynamics of semantic prediction based on inter-regional temporal generalization

The event-related potential/field component N400(m) has been widely used as a neural index for semantic prediction. It has long been hypothesized that feedback information from inferior frontal areas plays a critical role in generating the N400. However, due to limitations in causal connectivity estimation, direct testing of this hypothesis has remained difficult. Here, magnetoencephalography (MEG) data was obtained during a classic N400 paradigm where the semantic predictability of a fixed target noun was manipulated in simple German sentences. To estimate causality, we implemented a novel approach based on machine learning and temporal generalization to estimate the effect of inferior frontal gyrus (IFG) on temporal areas. In this method, a support vector machine (SVM) classifier is trained on each time point of the neural activity in IFG to classify less predicted (LP) and highly predicted (HP) nouns and then tested on all time points of superior/middle temporal sub-regions activity (and vice versa, to establish spatio-temporal evidence for or against causality). The decoding accuracy was significantly above chance level when the classifier was trained on IFG activity and tested on future activity in superior and middle temporal gyrus (STG/MTG). The results present new evidence for a model predictive speech comprehension where predictive IFG activity is fed back to shape subsequent activity in STG/MTG, implying a feedback mechanism in N400 generation. In combination with the also observed strong feedforward effect from left STG/MTG to IFG, our findings provide evidence of dynamic feedback and feedforward influences between IFG and temporal areas during N400 generation.

Visual Stimuli Modulate Local Field Potentials But Drive No High-Frequency Activity in Human Auditory Cortex

Neuroimaging studies suggest cross-sensory visual influences in human auditory cortices (ACs). Whether these influences reflect active visual processing in human ACs, which drives neuronal firing and concurrent broadband high-frequency activity (BHFA; >70 Hz), or whether they merely modulate sound processing is still debatable. Here, we presented auditory, visual, and audiovisual stimuli to 16 participants (7 women, 9 men) with stereo-EEG depth electrodes implanted near ACs for presurgical monitoring. Anatomically normalized group analyses were facilitated by inverse modeling of intracranial source currents. Analyses of intracranial event-related potentials (iERPs) suggested cross-sensory responses to visual stimuli in ACs, which lagged the earliest auditory responses by several tens of milliseconds. Visual stimuli also modulated the phase of intrinsic low-frequency oscillations and triggered 15-30 Hz event-related desynchronization in ACs. However, BHFA, a putative correlate of neuronal firing, was not significantly increased in ACs after visual stimuli, not even when they coincided with auditory stimuli. Intracranial recordings demonstrate cross-sensory modulations, but no indication of active visual processing in human ACs.

Interictal magnetic signals in new-onset Rolandic epilepsy may help with timing of treatment selection

OBJECTIVE

With research progress on Rolandic epilepsy (RE), its "benign" nature has been phased out. Clinicians are exhibiting an increasing tendency toward a more assertive treatment approach for RE. Nonetheless, in clinical practice, delayed treatment remains common because of the "self-limiting" nature of RE. Therefore, this study aimed to identify an imaging marker to aid treatment decisions and select a more appropriate time for initiating therapy for RE.

METHODS

We followed up with children newly diagnosed with RE, classified them into medicated and non-medicated groups according to the follow-up results, and compared them with matched healthy controls. Before beginning follow-up visits, interictal magnetic data were collected using magnetoencephalography in treatment-naïve recently diagnosed patients. The spectral power of the whole brain during initial diagnosis was determined using minimum normative estimation combined with the Welch technique.

RESULTS

A difference was observed in the magnetic source intensity within the left caudal anterior cingulate and precentral and postcentral gyri in the delta band between the medicated and non-medicated groups. The results revealed good discriminatory ability within the receiver operator characteristic curve. In the medicated group, there was a specific change in the frontotemporal magnetic source intensity, which shifted from high to low frequencies, compared with the healthy control group.

SIGNIFICANCE

The intensity of the precentral gyrus magnetic source within the delta band showed good specificity. Considering the rigor of initial treatment, the intensity of the precentral gyrus magnetic source can provide some help as an imaging marker for initial RE treatment, particularly for the timing of treatment initiation.

A role for retro-splenial cortex in the task-related P3 network

OBJECTIVE

The P3 is an event-related response observed in relation to task-relevant sensory events. Despite its ubiquitous presence, the neural generators of the P3 are controversial and not well identified.

METHODS

We compared source analysis of combined magneto- and electroencephalography (M/EEG) data with functional magnetic resonance imaging (fMRI) and simulation studies to better understand the sources of the P3 in an auditory oddball paradigm.

RESULTS

Our results suggest that the dominant source of the classical, postero-central P3 lies in the retro-splenial cortex of the ventral cingulate gyrus. A second P3 source in the anterior insular cortex contributes little to the postero-central maximum. Multiple other sources in the auditory, somatosensory, and anterior midcingulate cortex are active in an overlapping time window but can be functionally dissociated based on their activation time courses.

CONCLUSIONS

The retro-splenial cortex is a dominant source of the parietal P3 maximum in EEG.

SIGNIFICANCE

These results provide a new perspective for the interpretation of the extensive research based on the P3 response.

Both stronger and weaker cerebro-cerebellar functional connectivity patterns during processing of spoken sentences in autism spectrum disorder

Cerebellar differences have long been documented in autism spectrum disorder (ASD), yet the extent to which such differences might impact language processing in ASD remains unknown. To investigate this, we recorded brain activity with magnetoencephalography (MEG) while ASD and age-matched typically developing (TD) children passively processed spoken meaningful English and meaningless Jabberwocky sentences. Using a novel source localization approach that allows higher resolution MEG source localization of cerebellar activity, we found that, unlike TD children, ASD children showed no difference between evoked responses to meaningful versus meaningless sentences in right cerebellar lobule VI. ASD children also had atypically weak functional connectivity in the meaningful versus meaningless speech condition between right cerebellar lobule VI and several left-hemisphere sensorimotor and language regions in later time windows. In contrast, ASD children had atypically strong functional connectivity for in the meaningful versus meaningless speech condition between right cerebellar lobule VI and primary auditory cortical areas in an earlier time window. The atypical functional connectivity patterns in ASD correlated with ASD severity and the ability to inhibit involuntary attention. These findings align with a model where cerebro-cerebellar speech processing mechanisms in ASD are impacted by aberrant stimulus-driven attention, which could result from atypical temporal information and predictions of auditory sensory events by right cerebellar lobule VI.

A role for retro-splenial cortex in the task-related P3 network

Objective

The P3 is an event-related response observed in relation to task-relevant sensory events. Despite its ubiquitous presence, the neural generators of the P3 are controversial and not well identified.

Methods

We compared source analysis of combined magneto- and electroencephalography (M/EEG) data with functional magnetic resonance imaging (fMRI) and simulation studies to better understand the sources of the P3 in an auditory oddball paradigm.

Results

Our results suggest that the dominant source of the classical, postero-central P3 lies in the retro-splenial cortex of the ventral cingulate gyrus. A second P3 source in the anterior insular cortex contributes little to the postero-central maximum. Multiple other sources in the auditory, somatosensory, and anterior midcingulate cortex are active in an overlapping time window but can be functionally dissociated based on their activation time courses.

Conclusion

The retro-splenial cortex is a dominant source of the parietal P3 maximum in EEG.

Significance

These results provide a new perspective for the interpretation of the extensive research based on the P3 response.

Targeted Modulation of Human Brain Interregional Effective Connectivity With Spike-Timing Dependent Plasticity

OBJECTIVE

The ability to selectively up- or downregulate interregional brain connectivity would be useful for research and clinical purposes. Toward this aim, cortico-cortical paired associative stimulation (ccPAS) protocols have been developed in which two areas are repeatedly stimulated with a millisecond-level asynchrony. However, ccPAS results in humans using bifocal transcranial magnetic stimulation (TMS) have been variable, and the mechanisms remain unproven. In this study, our goal was to test whether ccPAS mechanism is spike-timing-dependent plasticity (STDP).

MATERIALS AND METHODS

Eleven healthy participants received ccPAS to the left primary motor cortex (M1) → right M1 with three different asynchronies (5 milliseconds shorter, equal to, or 5 milliseconds longer than the 9-millisecond transcallosal conduction delay) in separate sessions. To observe the neurophysiological effects, single-pulse TMS was delivered to the left M1 before and after ccPAS while cortico-cortical evoked responses were extracted from the contralateral M1 using source-resolved electroencephalography.

RESULTS

Consistent with STDP mechanisms, the effects on synaptic strengths flipped depending on the asynchrony. Further implicating STDP, control experiments suggested that the effects were unidirectional and selective to the targeted connection.

CONCLUSION

The results support the idea that ccPAS induces STDP and may selectively up- or downregulate effective connectivity between targeted regions in the human brain.

Language-related brain areas in childhood epilepsy with centrotemporal spikes studied with MEG

OBJECTIVE

Children with self-limited epilepsy with centrotemporal spikes (SeLECTS) typically indicate cognitive impairment with widespread speech impairment. We explored how epilepsy affects language-related brain areas and areas in their vicinity.

METHODS

Twenty-two children with SeLECTS and declined verbal comprehension (DVC), 21 with SeLECTS and normal verbal comprehension (NVC), and 23 healthy controls (HCs) underwent high-sampling magnetoencephalography recordings. According to a previous study, 24 language-related regions of interest were selected bilaterally, and the relative spectral power was estimated using a minimum norm estimate.

RESULTS

The highest mean power spectral density was observed in the delta band for the DVC group, in the theta band for the NVC group, and in the alpha band for HCs within language-specific brain regions. The distinctions between the DVC and NVC groups in the delta and theta frequency bands were primarily concentrated in the right linguistic brain area.

CONCLUSIONS

Children with SeLECTS may have developmental problems in language-related brain areas, with different developmental levels observed in the DVC, NVC, and HC groups. The DVC group could have inferior speech comprehension due to a more significant number of seizures and more left-sided spike locations.

SIGNIFICANCE

Children having SeLECTS showed impaired brain maturation, leading to associated language impairment.

Yttrium-Iron Garnet Magnetometer in MEG: Advance towards Multi-Channel Arrays

Recently, a new kind of sensor applicable in magnetoencephalography (MEG) has been presented: a solid-state yttrium-iron garnet magnetometer (YIGM). The feasibility of yttrium-iron garnet magnetometers (YIGMs) was demonstrated in an alpha-rhythm registration experiment. In this paper, we propose the analysis of lead-field matrices for different possible multi-channel on-scalp sensor layouts using YIGMs with respect to information theory. Real noise levels of the new sensor were used to compute signal-to-noise ratio (SNR) and total information capacity (TiC), and compared with corresponding metrics that can be obtained with well-established MEG systems based on superconducting quantum interference devices (SQUIDs) and optically pumped magnetometers (OPMs). The results showed that due to YIGMs' proximity to the subject's scalp, they outperform SQUIDs and OPMs at their respective noise levels in terms of SNR and TiC. However, the current noise levels of YIGM sensors are unfortunately insufficient for constructing a multichannel YIG-MEG system. This simulation study provides insight into the direction for further development of YIGM sensors to create a multi-channel MEG system, namely, by decreasing the noise levels of sensors.

The IAS-MEEG Package: A Flexible Inverse Source Reconstruction Platform for Reconstruction and Visualization of Brain Activity from M/EEG Data

We present a standalone Matlab software platform complete with visualization for the reconstruction of the neural activity in the brain from MEG or EEG data. The underlying inversion combines hierarchical Bayesian models and Krylov subspace iterative least squares solvers. The Bayesian framework of the underlying inversion algorithm allows to account for anatomical information and possible a priori belief about the focality of the reconstruction. The computational efficiency makes the software suitable for the reconstruction of lengthy time series on standard computing equipment. The algorithm requires minimal user provided input parameters, although the user can express the desired focality and accuracy of the solution. The code has been designed so as to favor the parallelization performed automatically by Matlab, according to the resources of the host computer. We demonstrate the flexibility of the platform by reconstructing activity patterns with supports of different sizes from MEG and EEG data. Moreover, we show that the software reconstructs well activity patches located either in the subcortical brain structures or on the cortex. The inverse solver and visualization modules can be used either individually or in combination. We also provide a version of the inverse solver that can be used within Brainstorm toolbox. All the software is available online by Github, including the Brainstorm plugin, with accompanying documentation and test data.

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.

Behavioral Studies Using Large-Scale Brain Networks – Methods and Validations

Mapping human behaviors to brain activity has become a key focus in modern cognitive neuroscience. As methods such as functional MRI (fMRI) advance cognitive scientists show an increasing interest in investigating neural activity in terms of functional connectivity and brain networks, rather than activation in a single brain region. Due to the noisy nature of neural activity, determining how behaviors are associated with specific neural signals is not well-established. Previous research has suggested graph theory techniques as a solution. Graph theory provides an opportunity to interpret human behaviors in terms of the topological organization of brain network architecture. Graph theory-based approaches, however, only scratch the surface of what neural connections relate to human behavior. Recently, the development of data-driven methods, e.g., machine learning and deep learning approaches, provide a new perspective to study the relationship between brain networks and human behaviors across the whole brain, expanding upon past literatures. In this review, we sought to revisit these data-driven approaches to facilitate our understanding of neural mechanisms and build models of human behaviors. We start with the popular graph theory approach and then discuss other data-driven approaches such as connectome-based predictive modeling, multivariate pattern analysis, network dynamic modeling, and deep learning techniques that quantify meaningful networks and connectivity related to cognition and behaviors. Importantly, for each topic, we discuss the pros and cons of the methods in addition to providing examples using our own data for each technique to describe how these methods can be applied to real-world neuroimaging data.

Ballistocardiogram suppression in concurrent EEG-MRI by dynamic modeling of heartbeats

The ballistocardiogram (BCG), the induced electric potentials by the head motion originating from heartbeats, is a prominent source of noise in electroencephalography (EEG) data during magnetic resonance imaging (MRI). Although methods have been proposed to suppress the BCG artifact, more work considering the variability of cardiac cycles and head motion across time and subjects is needed to provide highly robust correction. Here, a method called "dynamic modeling of heartbeats" (DMH) is proposed to reduce BCG artifacts in EEG data recorded inside an MRI system. The DMH method models BCG artifacts by combining EEG points at time instants with similar dynamics. The modeled BCG artifact is then subtracted from the EEG recording to suppress the BCG artifact. Performance of DMH was tested and specifically compared with the Optimal Basis Set (OBS) method on EEG data recorded inside a 3T MRI system with either no MRI acquisition (Inside-MRI), echo-planar imaging (EPI-EEG), or fast MRI acquisition using simultaneous multi-slice and inverse imaging methods (SMS-InI-EEG). In a steady-state visual evoked response (SSVEP) paradigm, the 15-Hz oscillatory neuronal activity at the visual cortex after DMH processing was about 130% of that achieved by OBS processing for Inside-MRI, SMS-InI-EEG, and EPI-EEG conditions. The DMH method is computationally efficient for suppressing BCG artifacts and in the future may help to improve the quality of EEG data recorded in high-field MRI systems for neuroscientific and clinical applications.

FLUX: A pipeline for MEG analysis

Magnetoencephalography (MEG) allows for quantifying modulations of human neuronal activity on a millisecond time scale while also making it possible to estimate the location of the underlying neuronal sources. The technique relies heavily on signal processing and source modelling. To this end, there are several open-source toolboxes developed by the community. While these toolboxes are powerful as they provide a wealth of options for analyses, the many options also pose a challenge for reproducible research as well as for researchers new to the field. The FLUX pipeline aims to make the analyses steps and setting explicit for standard analysis done in cognitive neuroscience. It focuses on quantifying and source localization of oscillatory brain activity, but it can also be used for event-related fields and multivariate pattern analysis. The pipeline is derived from the Cogitate consortium addressing a set of concrete cognitive neuroscience questions. Specifically, the pipeline including documented code is defined for MNE Python (a Python toolbox) and FieldTrip (a Matlab toolbox), and a data set on visuospatial attention is used to illustrate the steps. The scripts are provided as notebooks implemented in Jupyter Notebook and MATLAB Live Editor providing explanations, justifications and graphical outputs for the essential steps. Furthermore, we also provide suggestions for text and parameter settings to be used in registrations and publications to improve replicability and facilitate pre-registrations. The FLUX can be used for education either in self-studies or guided workshops. We expect that the FLUX pipeline will strengthen the field of MEG by providing some standardization on the basic analysis steps and by aligning approaches across toolboxes. Furthermore, we also aim to support new researchers entering the field by providing education and training. The FLUX pipeline is not meant to be static; it will evolve with the development of the toolboxes and with new insights. Furthermore, with the anticipated increase in MEG systems based on the Optically Pumped Magnetometers, the pipeline will also evolve to embrace these developments.

Multiple Source Detection based on Spatial Clustering and Its Applications on Wearable OPM-MEG

OBJECTIVE

Magnetoencephalography (MEG) is a non-invasive technique that measures the magnetic fields of brain activity. In particular, a new type of optically pumped magnetometer (OPM)-based wearable MEG system has been developed in recent years. Source localization in MEG can provide signals and locations of brain activity. However, conventional source localization methods face the difficulty of accurately estimating multiple sources. The present study presented a new parametric method to estimate the number of sources and localize multiple sources. In addition, we applied the proposed method to a constructed wearable OPM-MEG system.

METHODS

We used spatial clustering of the dipole spatial distribution to detect sources. The spatial distribution of dipoles was obtained by segmenting the MEG data temporally into slices and then estimating the parameters of the dipoles on each data slice using the particle swarm optimization algorithm. Spatial clustering was performed using the spatial-temporal density-based spatial clustering of applications with a noise algorithm. The performance of our approach for detecting multiple sources was compared with that of four typical benchmark algorithms using the OPM-MEG sensor configuration.

RESULTS

The simulation results showed that the proposed method had the best performance for detecting multiple sources. Moreover, the effectiveness of the method was verified by a multimodel sensory stimuli experiment on a real constructed 31-channel OPM-MEG.

CONCLUSION

Our study provides an effective method for the detection of multiple sources.

SIGNIFICANCE

With the improvement of the source localization methods, MEG may have a wider range of applications in neuroscience and clinical research.

Aging-Related Changes in Cortical Sources of Sleep Oscillatory Neural Activity Following Motor Learning Reflect Contributions of Cortical Thickness and Pre-sleep Functional Activity

Oscillatory neural activity during sleep, such as that in the delta and sigma bands, is important for motor learning consolidation. This activity is reduced with typical aging, and this reduction may contribute to aging-related declines in motor learning consolidation. Evidence suggests that brain regions involved in motor learning contribute to oscillatory neural activity during subsequent sleep. However, aging-related differences in regional contributions to sleep oscillatory activity following motor learning are unclear. To characterize these differences, we estimated the cortical sources of consolidation-related oscillatory activity using individual anatomical information in young and older adults during non-rapid eye movement sleep after motor learning and analyzed them in light of cortical thickness and pre-sleep functional brain activation. High-density electroencephalogram was recorded from young and older adults during a midday nap, following completion of a functional magnetic resonance imaged serial reaction time task as part of a larger experimental protocol. Sleep delta activity was reduced with age in a left-weighted motor cortical network, including premotor cortex, primary motor cortex, supplementary motor area, and pre-supplementary motor area, as well as non-motor regions in parietal, temporal, occipital, and cingulate cortices. Sleep theta activity was reduced with age in a similar left-weighted motor network, and in non-motor prefrontal and middle cingulate regions. Sleep sigma activity was reduced with age in left primary motor cortex, in a non-motor right-weighted prefrontal-temporal network, and in cingulate regions. Cortical thinning mediated aging-related sigma reductions in lateral orbitofrontal cortex and frontal pole, and partially mediated delta reductions in parahippocampal, fusiform, and lingual gyri. Putamen, caudate, and inferior parietal cortex activation prior to sleep predicted frontal and motor cortical contributions to sleep delta and theta activity in an age-moderated fashion, reflecting negative relationships in young adults and positive or absent relationships in older adults. Overall, these results support the local sleep hypothesis that brain regions active during learning contribute to consolidation-related neural activity during subsequent sleep and demonstrate that sleep oscillatory activity in these regions is reduced with aging.

MEG correlates of temporal regularity relevant to pitch perception in human auditory cortex

We recorded neural responses in human participants to three types of pitch-evoking regular stimuli at rates below and above the lower limit of pitch using magnetoencephalography (MEG). These bandpass filtered (1–4 kHz) stimuli were harmonic complex tones (HC), click trains (CT), and regular interval noise (RIN). Trials consisted of noise-regular-noise (NRN) or regular-noise-regular (RNR) segments in which the repetition rate (or fundamental frequency F0) was either above (250 Hz) or below (20 Hz) the lower limit of pitch. Neural activation was estimated and compared at the senor and source levels.

The pitch-relevant regular stimuli (F0 = 250 Hz) were all associated with marked evoked responses at around 140 ms after noise-to-regular transitions at both sensor and source levels. In particular, greater evoked responses to pitch-relevant stimuli than pitch-irrelevant stimuli (F0 = 20 Hz) were localized along the Heschl's sulcus around 140 ms. The regularity-onset responses for RIN were much weaker than for the other types of regular stimuli (HC, CT). This effect was localized over planum temporale, planum polare, and lateral Heschl's gyrus. Importantly, the effect of pitch did not interact with the stimulus type. That is, we did not find evidence to support different responses for different types of regular stimuli from the spatiotemporal cluster of the pitch effect (∼140 ms).

The current data demonstrate cortical sensitivity to temporal regularity relevant to pitch that is consistently present across different pitch-relevant stimuli in the Heschl's sulcus between Heschl's gyrus and planum temporale, both of which have been identified as a “pitch center” based on different modalities.

Laminar dynamics of high amplitude beta bursts in human motor cortex

Motor cortical activity in the beta frequency range is one of the strongest and most studied movement-related neural signals. At the single trial level, beta band activity is often characterized by transient, high amplitude, bursting events rather than slowly modulating oscillations. The timing of these bursting events is tightly linked to behavior, suggesting a more dynamic functional role for beta activity than previously believed. However, the neural mechanisms underlying beta bursts in sensorimotor circuits are poorly understood. To address this, we here leverage and extend recent developments in high precision MEG for temporally resolved laminar analysis of burst activity, combined with a neocortical circuit model that simulates the biophysical generators of the electrical currents which drive beta bursts. This approach pinpoints the generation of beta bursts in human motor cortex to distinct excitatory synaptic inputs to deep and superficial cortical layers, which drive current flow in opposite directions. These laminar dynamics of beta bursts in motor cortex align with prior invasive animal recordings within the somatosensory cortex, and suggest a conserved mechanism for somatosensory and motor cortical beta bursts. More generally, we demonstrate the ability for uncovering the laminar dynamics of event-related neural signals in human non-invasive recordings. This provides important constraints to theories about the functional role of burst activity for movement control in health and disease, and crucial links between macro-scale phenomena measured in humans and micro-circuit activity recorded from animal models.

Bilateral Representation of Sensorimotor Responses in Benign Adult Familial Myoclonus Epilepsy: An MEG Study

Patients with cortical reflex myoclonus manifest typical neurophysiologic characteristics due to primary sensorimotor cortex (S1/M1) hyperexcitability, namely, contralateral giant somatosensory-evoked potentials/fields and a C-reflex (CR) in the stimulated arm. Some patients show a CR in both arms in response to unilateral stimulation, with about 10-ms delay in the non-stimulated compared with the stimulated arm. This bilateral C-reflex (BCR) may reflect strong involvement of bilateral S1/M1. However, the significance and exact pathophysiology of BCR within 50 ms are yet to be established because it is difficult to identify a true ipsilateral response in the presence of the giant component in the contralateral hemisphere. We hypothesized that in patients with BCR, bilateral S1/M1 activity will be detected using MEG source localization and interhemispheric connectivity will be stronger than in healthy controls (HCs) between S1/M1 cortices. We recruited five patients with cortical reflex myoclonus with BCR and 15 HCs. All patients had benign adult familial myoclonus epilepsy. The median nerve was electrically stimulated unilaterally. Ipsilateral activity was investigated in functional regions of interest that were determined by the N20m response to contralateral stimulation. Functional connectivity was investigated using weighted phase-lag index (wPLI) in the time-frequency window of 30–50 ms and 30–100 Hz. Among seven of the 10 arms of the patients who showed BCR, the average onset-to-onset delay between the stimulated and the non-stimulated arm was 8.4 ms. Ipsilateral S1/M1 activity was prominent in patients. The average time difference between bilateral cortical activities was 9.4 ms. The average wPLI was significantly higher in the patients compared with HCs in specific cortico-cortical connections. These connections included precentral-precentral, postcentral-precentral, inferior parietal (IP)-precentral, and IP-postcentral cortices interhemispherically (contralateral region-ipsilateral region), and precentral-IP and postcentral-IP intrahemispherically (contralateral region-contralateral region). The ipsilateral response in patients with BCR may be a pathologically enhanced motor response homologous to the giant component, which was too weak to be reliably detected in HCs. Bilateral representation of sensorimotor responses is associated with disinhibition of the transcallosal inhibitory pathway within homologous motor cortices, which is mediated by the IP. IP may play a role in suppressing the inappropriate movements seen in cortical myoclonus.

Composition within and between Languages in the Bilingual Mind: MEG Evidence from Korean/English Bilinguals

The ability of the human brain to build complex expressions from simpler parts is fascinating, but the ability of the bilingual brain to do so is perhaps even more remarkable. When highly proficient bilinguals converse, they can fluidly switch from one language to another even inside sentences. Thus, they build expressions using words from more than one language. How are bilinguals able to compose words across different languages in real time? While robust evidence has implicated the left anterior temporal lobe (LATL) for the composition of words within one language, we do not know how the LATL, or other regions implicated for composition, operates when the language switches. We also do not know whether prefrontal regions associated with language control are recruited for language switching during composition. We addressed these questions with magnetoencephalography measurements in bilinguals who are fluent in two typologically distant languages, English and Korean. We observed early composition effects in the LATL at ∼200 ms that were unaffected by either language or orthography switching, which was also varied (Hangul vs Roman alphabet). Thus, the combinatory mechanism at 200 ms housed in the anterior temporal cortex appears blind to the language in which its input concepts are expressed. However, in later time windows, language and orthography switching interacted both in regions implicated for composition [LATL, ventromedial prefrontal cortex, left inferior frontal gyrus (LIFG)] as well as in regions associated with language control (ACC, LIFG). This establishes a starting point for understanding how bilingual brains code switch: words are initially combined without consideration of which language they come from, but language switching affects later processing.

How the brain negotiates divergent executive processing demands: Evidence of network reorganization in fleeting brain states

During performance in everyday contexts, multiple networks draw from shared executive resources to maintain attention, regulate arousal, and solve problems. At times, requirements for attention and self-regulation appear to be in competition. How does the brain attempt to resolve conflicts arising from such divergent processing demands? Here we demonstrate that the brain is capable of managing multiple processes via rapidly cycling between functional brain states over time, as it is typically regarded. Treating the brain as a complex system, comprising relationships within and between functional networks, we implemented Hidden Markov Modeling (HMM) on electroencephalographic (EEG) data to identify nonlinear brain states in both intra and internetwork synchrony that produced better performance for women subjects who were tasked with solving difficult problems under autobiographically-relevant, evaluative stress. Prior work often found that emotion-regulation and default-mode network (ERN and DMN) activity conflicted with the frontoparietal network's (FPN) ability to facilitate executive functioning necessary for problem solving. Contrastingly, we discovered that fleeting, nonlinear states dominated by FPN and ERN internetwork synchrony supported optimum performance generally, while during stress, states dominated by ERN and DMN intranetwork synchrony were more important for performance. These results imply that the brain may be capable of resolving competing processes through networks' cooperative dynamics. Further, data suggests a novel role for DMN as a mechanism for integrating external threats with internal, self-referent processing during evaluative stress within the observed population.

Co-registration Comparison of On-Scalp Magnetoencephalography and Magnetic Resonance Imaging

Magnetoencephalography (MEG) can non-invasively measure the electromagnetic activity of the brain. A new type of MEG, on-scalp MEG, has attracted the attention of researchers recently. Compared to the conventional SQUID-MEG, on-scalp MEG constructed with optically pumped magnetometers is wearable and has a high signal-to-noise ratio. While the co-registration between MEG and magnetic resonance imaging (MRI) significantly influences the source localization accuracy, co-registration error requires assessment, and quantification. Recent studies have evaluated the co-registration error of on-scalp MEG mainly based on the surface fit error or the repeatability error of different measurements, which do not reflect the true co-registration error. In this study, a three-dimensional-printed reference phantom was constructed to provide the ground truth of MEG sensor locations and orientations relative to MRI. The co-registration performances of commonly used three devices—electromagnetic digitization system, structured-light scanner, and laser scanner—were compared and quantified by the indices of final co-registration errors in the reference phantom and human experiments. Furthermore, the influence of the co-registration error on the performance of source localization was analyzed via simulations. The laser scanner had the best co-registration accuracy (rotation error of 0.23° and translation error of 0.76 mm based on the phantom experiment), whereas the structured-light scanner had the best cost performance. The results of this study provide recommendations and precautions for researchers regarding selecting and using an appropriate device for the co-registration of on-scalp MEG and MRI.

Synchronization patterns reveal neuronal coding of working memory content

Neuronal oscillations are suggested to play an important role in auditory working memory (WM), but their contribution to content-specific representations has remained unclear. Here, we measure magnetoencephalography during a retro-cueing task with parametric ripple-sound stimuli, which are spectrotemporally similar to speech but resist non-auditory memory strategies. Using machine learning analyses, with rigorous between-subject cross-validation and non-parametric permutation testing, we show that memorized sound content is strongly represented in phase-synchronization patterns between subregions of auditory and frontoparietal cortices. These phase-synchronization patterns predict the memorized sound content steadily across the studied maintenance period. In addition to connectivity-based representations, there are indices of more local, “activity silent” representations in auditory cortices, where the decoding accuracy of WM content significantly increases after task-irrelevant “impulse stimuli.” Our results demonstrate that synchronization patterns across auditory sensory and association areas orchestrate neuronal coding of auditory WM content. This connectivity-based coding scheme could also extend beyond the auditory domain.

Mamashli et al. use machine learning analyses of human magnetoencephalography (MEG) recordings to study “working memory,” maintenance of information in mind over brief periods of time. Their results show that the human brain maintains working memory content in transient functional connectivity patterns across sensory and association areas.

Mirror Visual Feedback Induces M1 Excitability by Disengaging Functional Connections of Perceptuo-Motor-Attentional Processes during Asynchronous Bimanual Movement: A Magnetoencephalographic Study

Mirror visual feedback (MVF) has been shown to increase the excitability of the primary motor cortex (M1) during asynchronous bimanual movement. However, the functional networks underlying this process remain unclear. We recruited 16 healthy volunteers to perform asynchronous bimanual movement, that is, their left hand performed partial range of movement while their right hand performed normal full range of movement. Their ongoing brain activities were recorded by whole-head magnetoencephalography during the movement. Participants were required to keep both hands stationary in the control condition. In the other two conditions, participants were required to perform asynchronous bimanual movement with MVF (Asy_M) and without MVF (Asy_w/oM). Greater M1 excitability was found under Asy_M than under Asy_w/oM. More importantly, when receiving MVF, the visual cortex reduced its functional connection to brain regions associated with perceptuo-motor-attentional process (i.e., M1, superior temporal gyrus, and dorsolateral prefrontal cortex). This is the first study to demonstrate a global functional network of MVF during asynchronous bimanual movement, providing a foundation for future research to examine the neural mechanisms of mirror illusion in motor control.

From letters to composed concepts: A magnetoencephalography study of reading

Language comprehension requires the recognition of individual words and the combination of their meanings to yield complex concepts or interpretations. This combinatory process often requires the insertion of unstated semantic material between words, based on thematic or feature knowledge. For example, the phrase horse barn is not interpreted as a blend of a horse and a barn, but specifically a barn where horses are kept. Previous neuroscientific evidence suggests that left posterior and anterior temporal cortex underpin thematic and feature‐based concept knowledge, respectively, but much remains unclear about how these areas contribute to combinatory language processing. Using magnetoencephalography, we contrasted source‐localized responses to modifier‐noun phrases involving thematic relations versus feature modifications, while also examining how lower‐level orthographic processing fed composition. Participants completed three procedures examining responses to letter‐strings, adjective‐noun phrases, and noun–noun combinations that varied the semantic relations between words. We found that sections of the left anterior temporal lobe, posterior temporal lobe, and cortex surrounding the angular gyrus were all engaged in the minimal composition of adjective‐noun phrases, a more distributed network than in most prior studies of minimal composition. Of these regions, only the left posterior temporal lobe was additionally sensitive to implicit thematic relations between composing words, suggesting that it houses a specialized relational processing component in a wider composition network. We additionally identified a left occipitotemporal progression from orthographic to lexical processing, feeding ventral anterior areas engaged in the combination of word meanings. Finally, by examining source signal leakage, we characterized the degree to which these responses could be distinguished from one another using source estimation.

One Way or Another: Cortical Language Areas Flexibly Adapt Processing Strategies to Perceptual And Contextual Properties of Speech

Cortical circuits rely on the temporal regularities of speech to optimize signal parsing for sound-to-meaning mapping. Bottom-up speech analysis is accelerated by top-down predictions about upcoming words. In everyday communications, however, listeners are regularly presented with challenging input-fluctuations of speech rate or semantic content. In this study, we asked how reducing speech temporal regularity affects its processing-parsing, phonological analysis, and ability to generate context-based predictions. To ensure that spoken sentences were natural and approximated semantic constraints of spontaneous speech we built a neural network to select stimuli from large corpora. We analyzed brain activity recorded with magnetoencephalography during sentence listening using evoked responses, speech-to-brain synchronization and representational similarity analysis. For normal speech theta band (6.5-8 Hz) speech-to-brain synchronization was increased and the left fronto-temporal areas generated stronger contextual predictions. The reverse was true for temporally irregular speech-weaker theta synchronization and reduced top-down effects. Interestingly, delta-band (0.5 Hz) speech tracking was greater when contextual/semantic predictions were lower or if speech was temporally jittered. We conclude that speech temporal regularity is relevant for (theta) syllabic tracking and robust semantic predictions while the joint support of temporal and contextual predictability reduces word and phrase-level cortical tracking (delta).

Identifying Mild Cognitive Impairment in Parkinson's Disease With Electroencephalogram Functional Connectivity

Objective

Cognitive impairment occurs frequently in Parkinson’s disease (PD) and negatively impacts the patient’s quality of life. However, its pathophysiological mechanism remains unclear, hindering the development of new therapies. Changes in brain connectivity are related to cognitive impairment in patients with PD, with the dorsolateral prefrontal cortex (DLPFC) being considered the essential region related to PD cognitive impairment. Nevertheless, few studies have focused on the global connectivity responsible for communication with the DLPFC node, the posterior division of the middle frontal gyrus (PMFG) in patients with PD; this was the focus of this study.

Methods

We applied resting-state electroencephalography (EEG) and calculated a reliable functional connectivity measurement, the debiased weighted phase lag index (dWPLI), to examine inter-regional functional connectivity in 68 patients with PD who were classified into two groups according to their cognitive condition.

Results

We observed that altered left and right PMFG-based functional connectivity associated with cognitive impairment in patients with PD in the theta frequency bands under the eyes closed condition (r = −0.426, p < 0.001 and r = −0.437, p < 0.001, respectively). Exploratory results based on the MoCA subdomains indicated that poorer visuospatial function was associated with higher right PMFG-based functional connectivity (r = −0.335, p = 0.005), and poorer attention function was associated with higher left and right PMFG-based functional connectivity (r = −0.380, p = 0.001 and r = −0.256, p = 0.035, respectively). Further analysis using logistic regression and receiver operating characteristic (ROC) curves found that this abnormal functional connectivity was an independent risk factor for cognitive impairment [odds ratio (OR): 2.949, 95% confidence interval (CI): 1.294–6.725, p = 0.01 for left PMFG; OR: 11.278, 95% CI: 2.578–49.335, p = 0.001 for right PMFG, per 0.1 U], and provided moderate classification power to discriminate between cognitive abilities in patients with PD [area under the ROC curve (AUC) = 0.770 for left PMFG; AUC = 0.809 for right PMFG].

Conclusion

These preliminary findings indicate that abnormal PMFG-based functional connectivity patterns associated with cognitive impairment in the theta frequency bands under the eyes closed condition and altered functional connectivity patterns have the potential to act as reliable biomarkers for identifying cognitive impairment in patients with PD.

Enhancing visual motion discrimination by desynchronizing bifocal oscillatory activity

Visual motion discrimination involves reciprocal interactions in the alpha band between the primary visual cortex (V1) and mediotemporal areas (V5/MT). We investigated whether modulating alpha phase synchronization using individualized multisite transcranial alternating current stimulation (tACS) over V5 and V1 regions would improve motion discrimination. We tested 3 groups of healthy subjects with the following conditions: (1) individualized In-Phase V1_alpha-V5_alpha tACS (0° lag), (2) individualized Anti-Phase V1_alpha-V5_alpha tACS (180° lag) and (3) sham tACS. Motion discrimination and EEG activity were recorded before, during and after tACS. Performance significantly improved in the Anti-Phase group compared to the In-Phase group 10 and 30 min after stimulation. This result was explained by decreases in bottom-up alpha-V1 gamma-V5 phase-amplitude coupling. One possible explanation of these results is that Anti-Phase V1_alpha-V5_alpha tACS might impose an optimal phase lag between stimulation sites due to the inherent speed of wave propagation, hereby supporting optimized neuronal communication.

Altered maturation and atypical cortical processing of spoken sentences in autism spectrum disorder

Autism spectrum disorder (ASD) is associated with widespread receptive language impairments, yet the neural mechanisms underlying these deficits are poorly understood. Neuroimaging has shown that processing of socially-relevant sounds, including speech and non-speech, is atypical in ASD. However, it is unclear how the presence of lexical-semantic meaning affects speech processing in ASD. Here, we recorded magnetoencephalography data from individuals with ASD (N = 22, ages 7-17, 4 females) and typically developing (TD) peers (N = 30, ages 7-17, 5 females) during unattended listening to meaningful auditory speech sentences and meaningless jabberwocky sentences. After adjusting for age, ASD individuals showed stronger responses to meaningless jabberwocky sentences than to meaningful speech sentences in the same left temporal and parietal language regions where TD individuals exhibited stronger responses to meaningful speech. Maturational trajectories of meaningful speech responses were atypical in temporal, but not parietal, regions in ASD. Temporal responses were associated with ASD severity, while parietal responses were associated with aberrant involuntary attentional shifting in ASD. Our findings suggest a receptive speech processing dysfunction in ASD, wherein unattended meaningful speech elicits abnormal engagement of the language system, while unattended meaningless speech, filtered out in TD individuals, engages the language system through involuntary attention capture.

Auditory cues facilitate object movement processing in human extrastriate visual cortex during simulated self-motion: A pilot study

Visual segregation of moving objects is a considerable computational challenge when the observer moves through space. Recent psychophysical studies suggest that directionally congruent, moving auditory cues can substantially improve parsing object motion in such settings, but the exact brain mechanisms and visual processing stages that mediate these effects are still incompletely known. Here, we utilized multivariate pattern analyses (MVPA) of MRI-informed magnetoencephalography (MEG) source estimates to examine how crossmodal auditory cues facilitate motion detection during the observer's self-motion. During MEG recordings, participants identified a target object that moved either forward or backward within a visual scene that included nine identically textured objects simulating forward observer translation. Auditory motion cues 1) improved the behavioral accuracy of target localization, 2) significantly modulated the MEG source activity in the areas V2 and human middle temporal complex (hMT+), and 3) increased the accuracy at which the target movement direction could be decoded from hMT+ activity using MVPA. The increase of decoding accuracy by auditory cues in hMT+ was significant also when superior temporal activations in or near auditory cortices were regressed out from the hMT+ source activity to control for source estimation biases caused by point spread. Taken together, these results suggest that parsing object motion from self-motion-induced optic flow in the human extrastriate visual cortex can be facilitated by crossmodal influences from auditory system.

MEG current source reconstruction using a meta-analysis fMRI prior

Magnetoencephalography (MEG) offers a unique way to noninvasively investigate millisecond-order cortical activities by mapping sensor signals (magnetic fields outside the head) to cortical current sources using current source reconstruction methods. Current source reconstruction is defined as an ill-posed inverse problem, since the number of sensors is less than the number of current sources. One powerful approach to solving this problem is to use functional MRI (fMRI) data as a spatial constraint, although it boosts the cost of measurement and the burden on subjects. Here, we show how to use the meta-analysis fMRI data synthesized from thousands of papers instead of the individually recorded fMRI data. To mitigate the differences between the meta-analysis and individual data, the former are imported as prior information of the hierarchical Bayesian estimation. Using realistic simulations, we found out the performance of current source reconstruction using meta-analysis fMRI data to be better than that using low-quality individual fMRI data and conventional methods. By applying experimental data of a face recognition task, we qualitatively confirmed that group analysis results using the meta-analysis fMRI data showed a tendency similar to the results using the individual fMRI data. Our results indicate that the use of meta-analysis fMRI data improves current source reconstruction without additional measurement costs. We assume the proposed method would have greater effect for modalities with lower measurement costs, such as optically pumped magnetometers.

Differential neurodynamics and connectivity in the dorsal and ventral visual pathways during perception of emotional crowds and individuals: a MEG study

Reading the prevailing emotion of groups of people ("crowd emotion") is critical to understanding their overall intention and disposition. It alerts us to potential dangers, such as angry mobs or panicked crowds, giving us time to escape. A critical aspect of processing crowd emotion is that it must occur rapidly, because delays often are costly. Although knowing the timing of neural events is crucial for understanding how the brain guides behaviors using coherent signals from a glimpse of multiple faces, this information is currently lacking in the literature on face ensemble coding. Therefore, we used magnetoencephalography to examine the neurodynamics in the dorsal and ventral visual streams and the periamygdaloid cortex to compare perception of groups of faces versus individual faces. Forty-six participants compared two groups of four faces or two individual faces with varying emotional expressions and chose which group or individual they would avoid. We found that the dorsal stream was activated as early as 68 msec after the onset of stimuli containing groups of faces. In contrast, the ventral stream was activated later and predominantly for individual face stimuli. The latencies of the dorsal stream activation peaks correlated with participants' response times for facial crowds. We also found enhanced connectivity earlier between the periamygdaloid cortex and the dorsal stream regions for crowd emotion perception. Our findings suggest that ensemble coding of facial crowds proceeds rapidly and in parallel by engaging the dorsal stream to mediate adaptive social behaviors, via a distinct route from single face perception.

Children with autism spectrum disorder show altered functional connectivity and abnormal maturation trajectories in response to inverted faces

The processing of information conveyed by faces is a critical component of social communication. While the neurophysiology of processing upright faces has been studied extensively in autism spectrum disorder (ASD), less is known about the neurophysiological abnormalities associated with processing inverted faces in ASD. We used magnetoencephalography (MEG) to study both long-range and local functional connectivity, with the latter assessed using local cross-frequency coupling, in response to inverted faces stimuli, in 7-18 years old individuals with ASD and age and IQ matched typically developing (TD) individuals. We found abnormally reduced coupling between the phase of the alpha rhythm and the amplitude of the gamma rhythm in the fusiform face area (FFA) in response to inverted faces, as well as reduced long-range functional connectivity between the FFA and the inferior frontal gyrus (IFG) in response to inverted faces in the ASD group. These group differences were absent in response to upright faces. The magnitude of functional connectivity between the FFA and the IFG was significantly correlated with the severity of ASD, and FFA-IFG long-range functional connectivity increased with age in TD group, but not in the ASD group. Our findings suggest that both local and long-range functional connectivity are abnormally reduced in children with ASD when processing inverted faces, and that the pattern of abnormalities associated with the processing of inverted faces differs from the pattern of upright faces in ASD, likely due to the presumed greater reliance on top-down regulations necessary for efficient processing of inverted faces. LAY SUMMARY: We found alterations in the neurophysiological responses to inverted faces in children with ASD, that were not reflected in the evoked responses, and were not observed in the responses to upright faces. These alterations included reduced local functional connectivity in the fusiform face area (FFA), and decreased long-range alpha-band modulated functional connectivity between the FFA and the left IFG. The magnitude of long-range functional connectivity between the FFA and the inferior frontal gyrus was correlated with the severity of ASD.

Distributed source modeling of intracranial stereoelectro-encephalographic measurements

Intracranial stereoelectroencephalography (sEEG) provides unsurpassed sensitivity and specificity for human neurophysiology. However, functional mapping of brain functions has been limited because the implantations have sparse coverage and differ greatly across individuals. Here, we developed a distributed, anatomically realistic sEEG source-modeling approach for within- and between-subject analyses. In addition to intracranial event-related potentials (iERP), we estimated the sources of high broadband gamma activity (HBBG), a putative correlate of local neural firing. Our novel approach accounted for a significant portion of the variance of the sEEG measurements in leave-one-out cross-validation. After logarithmic transformations, the sensitivity and signal-to-noise ratio were linearly inversely related to the minimal distance between the brain location and electrode contacts (slope≈-3.6). The signa-to-noise ratio and sensitivity in the thalamus and brain stem were comparable to those locations at the vicinity of electrode contact implantation. The HGGB source estimates were remarkably consistent with analyses of intracranial-contact data. In conclusion, distributed sEEG source modeling provides a powerful neuroimaging tool, which facilitates anatomically-normalized functional mapping of human brain using both iERP and HBBG data.

Context Matters: Situational Stress Impedes Functional Reorganization of Intrinsic Brain Connectivity during Problem-Solving

Extensive research has established a relationship between individual differences in brain activity in a resting state and individual differences in behavior. Conversely, when individuals are engaged in various tasks, certain task-evoked reorganization occurs in brain functional connectivity, which can consequently influence individuals' performance as well. Here, we show that resting state and task-dependent state brain patterns interact as a function of contexts engendering stress. Findings revealed that when the resting state connectome was examined during performance, the relationship between connectome strength and performance only remained for participants under stress (who also performed worse than all other groups on the math task), suggesting that stress preserved brain patterns indicative of underperformance whereas non-stressed individuals spontaneously transitioned out of these patterns. Results imply that stress may impede the reorganization of a functional network in task-evoked brain states. This hypothesis was subsequently verified using graph theory measurements on a functional network, independent of behavior. For participants under stress, the functional network showed less topological alterations compared to non-stressed individuals during the transition from resting state to task-evoked state. Implications are discussed for network dynamics as a function of context.

Classification of evoked responses to inverted faces reveals both spatial and temporal cortical response abnormalities in Autism spectrum disorder

The neurophysiology of face processing has been studied extensively in the context of social impairments associated with autism spectrum disorder (ASD), but the existing studies have concentrated mainly on univariate analyses of responses to upright faces, and, less frequently, inverted faces. The small number of existing studies on neurophysiological responses to inverted face in ASD have used univariate approaches, with divergent results. Here, we used a data-driven, classification-based, multivariate machine learning decoding approach to investigate the temporal and spatial properties of the neurophysiological evoked response for upright and inverted faces, relative to the neurophysiological evoked response for houses, a neutral stimulus. 21 (2 females) ASD and 29 (4 females) TD participants ages 7 to 19 took part in this study. Group level classification accuracies were obtained for each condition, using first the temporal domain of the evoked responses, and then the spatial distribution of the evoked responses on the cortical surface, each separately. We found that classification of responses to inverted neutral faces vs. houses was less accurate in ASD compared to TD, in both the temporal and spatial domains. In contrast, there were no group differences in the classification of evoked responses to upright neutral faces relative to houses. Using the classification in the temporal domain, lower decoding accuracies in ASD were found around 120 ms and 170 ms, corresponding the known components of the evoked responses to faces. Using the classification in the spatial domain, lower decoding accuracies in ASD were found in the right superior marginal gyrus (SMG), intra-parietal sulcus (IPS) and posterior superior temporal sulcus (pSTS), but not in core face processing areas. Importantly, individual classification accuracies from both the temporal and spatial classifiers correlated with ASD severity, confirming the relevance of the results to the ASD phenotype.

Lateralized evoked responses in parietal cortex demonstrate visual short-term memory deficits in first-episode schizophrenia

Working memory dysfunction may be central to neurocognitive deficits in schizophrenia. Maintenance of visual information in working memory, or visual short-term memory (vSTM), is linked to general cognitive dysfunction and predicts functional outcome. Lateralized change-detection tasks afford investigation of the contralateral delay activity (CDA), a useful tool for investigating vSTM dysfunction. Previous work suggests "hyperfocusing" of attention in schizophrenia, such that CDA is increased when a single item is maintained in vSTM but reduced for multiple items. If observed early in the disease, vSTM dysfunction may be a key feature of schizophrenia or target for intervention. We investigated CDA during lateralized vSTM of one versus three items using sensor-level electroencephalography and source-level magnetoencephalography in 26 individuals at their first episode of schizophrenia-spectrum psychosis (FESz) and 26 matched healthy controls. FESz were unable to modulate CDA with increased memory load - high-load CDA was reduced and low-load CDA was increased compared to controls. Further, sources of CDA in posterior parietal cortex were reduced in FESz and indices of working memory were correlated with neurocognitive deficits and symptom severity. These results support working memory maintenance dysfunction as a central and early component to the disorder. Targeted intervention focusing on vSTM deficits may be warranted to alleviate downstream effects of this disability.

Distinct Patterns of Functional Connectivity During the Comprehension of Natural, Narrative Speech

Recent continuous task studies, such as narrative speech comprehension, show that fluctuations in brain functional connectivity (FC) are altered and enhanced compared to the resting state. Here, we characterized the fluctuations in FC during comprehension of speech and time-reversed speech conditions. The correlations of Hilbert envelope of source-level EEG data were used to quantify FC between spatially separate brain regions. A symmetric multivariate leakage correction was applied to address the signal leakage issue before calculating FC. The dynamic FC was estimated based on a sliding time window. Then, principal component analysis (PCA) was performed on individually concatenated and temporally concatenated FC matrices to identify FC patterns. We observed that the mode of FC induced by speech comprehension can be characterized with a single principal component. The condition-specific FC demonstrated decreased correlations between frontal and parietal brain regions and increased correlations between frontal and temporal brain regions. The fluctuations of the condition-specific FC characterized by a shorter time demonstrated that dynamic FC also exhibited condition specificity over time. The FC is dynamically reorganized and FC dynamic pattern varies along a single mode of variation during speech comprehension. The proposed analysis framework seems valuable for studying the reorganization of brain networks during continuous task experiments.

Building the Arrow of Time… Over Time: A Sequence of Brain Activity Mapping Imagined Events in Time and Space

When moving, the spatiotemporal unfolding of events is bound to our physical trajectory, and time and space become entangled in episodic memory. When imagining past or future events, or being in different geographical locations, the temporal and spatial dimensions of mental events can be independently accessed and manipulated. Using time-resolved neuroimaging, we characterized brain activity while participants ordered historical events from different mental perspectives in time (e.g., when imagining being 9 years in the future) or in space (e.g., when imagining being in Cayenne). We describe 2 neural signatures of temporal ordinality: an early brain response distinguishing whether participants were mentally in the past, the present or the future (self-projection in time), and a graded activity at event retrieval, indexing the mental distance between the representation of the self in time and the event. Neural signatures of ordinality and symbolic distances in time were distinct from those observed in the homologous spatial task: activity indicating spatial order and distances overlapped in latency in distinct brain regions. We interpret our findings as evidence that the conscious representation of time and space share algorithms (egocentric mapping, distance, and ordinality computations) but different implementations with a distinctive status for the psychological "time arrow."

A methodological framework for inverse-modeling of propagating cortical activity using MEG/EEG

The prevailing view on the dynamics of large-scale electrical activity in the human cortex is that it constitutes a functional network of discrete and localized circuits. Within this view, a natural way to analyse magnetoencephalographic (MEG) and electroencephalographic (EEG) data is by adopting methods from network theory. Invasive recordings, however, demonstrate that cortical activity is spatially continuous, rather than discrete, and exhibits propagation behavior. Furthermore, human cortical activity is known to propagate under a variety of conditions such as non-REM sleep, general anesthesia, and coma. Although several MEG/EEG studies have investigated propagating cortical activity, not much is known about the conditions under which such activity can be successfully reconstructed from MEG/EEG sensor-data. This study provides a methodological framework for inverse-modeling of propagating cortical activity. Within this framework, cortical activity is represented in the spatial frequency domain, which is more natural than the dipole domain when dealing with spatially continuous activity. We define angular power spectra, which show how the power of cortical activity is distributed across spatial frequencies, angular gain/phase spectra, which characterize the spatial filtering properties of linear inverse operators, and angular resolution matrices, which summarize how linear inverse operators leak signal within and across spatial frequencies. We adopt the framework to provide insight into the performance of several linear inverse operators in reconstructing propagating cortical activity from MEG/EEG sensor-data. We also describe how prior spatial frequency information can be incorporated into the inverse-modeling to obtain better reconstructions.

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

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

Reduced Dorsal Visual Oscillatory Activity During Working Memory Maintenance in the First-Episode Schizophrenia Spectrum

Cognitive deficits in people with schizophrenia are among the hardest to treat and strongly predict functional outcome. The ability to maintain sensory precepts in memory over a short delay is impacted early in the progression of schizophrenia and has been linked to reliable neurophysiological markers. Yet, little is known about the mechanisms of these deficits. Here, we investigated possible neurophysiological mechanisms of impaired visual short-term memory (vSTM, aka working memory maintenance) in the first-episode schizophrenia spectrum (FESz) using magnetoencephalography (MEG). Twenty-eight FESz and 25 matched controls performed a lateralized change detection task where they were cued to selectively attend and remember colors of circles presented in either the left or right peripheral visual field over a 1 s delay. Contralateral alpha suppression (CAS) during the delay period was used to assess selective attention to cued visual hemifields held in vSTM. Delay-period CAS was compared between FESz and controls and between trials presenting one vs three items per visual hemifield. CAS in dorsal visual cortex was reduced in FESz compared to controls in high-load trials, but not low-load trials. Group differences in CAS were found beginning 100 ms after the disappearance of the memory set, suggesting deficits were not due to the initial deployment of attention to the cued visual hemifield prior to stimulus presentation. CAS was not greater for high-load vs low-load trials in FESz subjects, although this effect was prominent in controls. Further, lateralized gamma (34-40 Hz) power emerged in dorsal visual cortex prior to the onset of CAS in controls but not FESz. Gamma power in this cluster differed between groups at both high and low load. CAS deficits observed in FESz were correlated with change detection accuracy, working memory function, estimated IQ, and negative symptoms. Our results implicate deficits in CAS in trials requiring broad, but not narrow, focus of attention to spatially distributed objects maintained in vSTM in FESz, possibly due to reduced ability to broadly distribute visuospatial attention (alpha) or disruption of object-location binding (gamma) during encoding/consolidation. This early pathophysiology may shed light upon mechanisms of emerging working memory deficits that are intrinsic to schizophrenia.

Prior Knowledge Guides Speech Segregation in Human Auditory Cortex

Segregating concurrent sound streams is a computationally challenging task that requires integrating bottom-up acoustic cues (e.g. pitch) and top-down prior knowledge about sound streams. In a multi-talker environment, the brain can segregate different speakers in about 100 ms in auditory cortex. Here, we used magnetoencephalographic (MEG) recordings to investigate the temporal and spatial signature of how the brain utilizes prior knowledge to segregate 2 speech streams from the same speaker, which can hardly be separated based on bottom-up acoustic cues. In a primed condition, the participants know the target speech stream in advance while in an unprimed condition no such prior knowledge is available. Neural encoding of each speech stream is characterized by the MEG responses tracking the speech envelope. We demonstrate that an effect in bilateral superior temporal gyrus and superior temporal sulcus is much stronger in the primed condition than in the unprimed condition. Priming effects are observed at about 100 ms latency and last more than 600 ms. Interestingly, prior knowledge about the target stream facilitates speech segregation by mainly suppressing the neural tracking of the non-target speech stream. In sum, prior knowledge leads to reliable speech segregation in auditory cortex, even in the absence of reliable bottom-up speech segregation cue.

Distinct Regional Oscillatory Connectivity Patterns During Auditory Target and Novelty Processing

Auditory attention allows us to focus on relevant target sounds in the acoustic environment while maintaining the capability to orient to unpredictable (novel) sound changes. An open question is whether orienting to expected vs. unexpected auditory events are governed by anatomically distinct attention pathways, respectively, or by differing communication patterns within a common system. To address this question, we applied a recently developed PeSCAR analysis method to evaluate spectrotemporal functional connectivity patterns across subregions of broader cortical regions of interest (ROIs) to analyze magnetoencephalography data obtained during a cued auditory attention task. Subjects were instructed to detect a predictable harmonic target sound embedded among standard tones in one ear and to ignore the standard tones and occasional unpredictable novel sounds presented in the opposite ear. Phase coherence of estimated source activity was calculated between subregions of superior temporal, frontal, inferior parietal, and superior parietal cortex ROIs. Functional connectivity was stronger in response to target than novel stimuli between left superior temporal and left parietal ROIs and between left frontal and right parietal ROIs, with the largest effects observed in the beta band (15-35 Hz). In contrast, functional connectivity was stronger in response to novel than target stimuli in inter-hemispheric connections between left and right frontal ROIs, observed in early time windows in the alpha band (8-12 Hz). Our findings suggest that auditory processing of expected target vs. unexpected novel sounds involves different spatially, temporally, and spectrally distributed oscillatory connectivity patterns across temporal, parietal, and frontal areas.

Discovering dynamic task-modulated functional networks with specific spectral modes using MEG

Efficient neuronal communication between brain regions through oscillatory synchronization at certain frequencies is necessary for cognition. Such synchronized networks are transient and dynamic, established on the timescale of milliseconds in order to support ongoing cognitive operations. However, few studies characterizing dynamic electrophysiological brain networks have simultaneously accounted for temporal non-stationarity, spectral structure, and spatial properties. Here, we propose an analysis framework for characterizing the large-scale phase-coupling network dynamics during task performance using magnetoencephalography (MEG). We exploit the high spatiotemporal resolution of MEG to measure time-frequency dynamics of connectivity between parcellated brain regions, yielding data in tensor format. We then use a tensor component analysis (TCA)-based procedure to identify the spatio-temporal-spectral modes of covariation among separate regions in the human brain. We validate our pipeline using MEG data recorded during a hand movement task, extracting a transient motor network with beta-dominant spectral mode, which is significantly modulated by the movement task. Next, we apply the proposed pipeline to explore brain networks that support cognitive operations during a working memory task. The derived results demonstrate the temporal formation and dissolution of multiple phase-coupled networks with specific spectral modes, which are associated with face recognition, vision, and movement. The proposed pipeline can characterize the spectro-temporal dynamics of functional connectivity in the brain on the subsecond timescale, commensurate with that of cognitive performance.

Predictive regression modeling with MEG/EEG: from source power to signals and cognitive states

Predicting biomedical outcomes from Magnetoencephalography and Electroencephalography (M/EEG) is central to applications like decoding, brain-computer-interfaces (BCI) or biomarker development and is facilitated by supervised machine learning. Yet, most of the literature is concerned with classification of outcomes defined at the event-level. Here, we focus on predicting continuous outcomes from M/EEG signal defined at the subject-level, and analyze about 600 MEG recordings from Cam-CAN dataset and about 1000 EEG recordings from TUH dataset. Considering different generative mechanisms for M/EEG signals and the biomedical outcome, we propose statistically-consistent predictive models that avoid source-reconstruction based on the covariance as representation. Our mathematical analysis and ground-truth simulations demonstrated that consistent function approximation can be obtained with supervised spatial filtering or by embedding with Riemannian geometry. Additional simulations revealed that Riemannian methods were more robust to model violations, in particular geometric distortions induced by individual anatomy. To estimate the relative contribution of brain dynamics and anatomy to prediction performance, we propose a novel model inspection procedure based on biophysical forward modeling. Applied to prediction of outcomes at the subject-level, the analysis revealed that the Riemannian model better exploited anatomical information while sensitivity to brain dynamics was similar across methods. We then probed the robustness of the models across different data cleaning options. Environmental denoising was globally important but Riemannian models were strikingly robust and continued performing well even without preprocessing. Our results suggest each method has its niche: supervised spatial filtering is practical for event-level prediction while the Riemannian model may enable simple end-to-end learning.

Estimates of cortical column orientation improve MEG source inversion

Determining the anatomical source of brain activity non-invasively measured from EEG or MEG sensors is challenging. In order to simplify the source localization problem, many techniques introduce the assumption that current sources lie on the cortical surface. Another common assumption is that this current flow is orthogonal to the cortical surface, thereby approximating the orientation of cortical columns. However, it is not clear which cortical surface to use to define the current source locations, and normal vectors computed from a single cortical surface may not be the best approximation to the orientation of cortical columns. We compared three different surface location priors and five different approaches for estimating dipole vector orientation, both in simulations and visual and motor evoked MEG responses. We show that models with source locations on the white matter surface and using methods based on establishing correspondences between white matter and pial cortical surfaces dramatically outperform models with source locations on the pial or combined pial/white surfaces and which use methods based on the geometry of a single cortical surface in fitting evoked visual and motor responses. These methods can be easily implemented and adopted in most M/EEG analysis pipelines, with the potential to significantly improve source localization of evoked responses.

Exploring Frequency-Dependent Brain Networks from Ongoing EEG Using Spatial ICA During Music Listening

Recently, exploring brain activity based on functional networks during naturalistic stimuli especially music and video represents an attractive challenge because of the low signal-to-noise ratio in collected brain data. Although most efforts focusing on exploring the listening brain have been made through functional magnetic resonance imaging (fMRI), sensor-level electro- or magnetoencephalography (EEG/MEG) technique, little is known about how neural rhythms are involved in the brain network activity under naturalistic stimuli. This study exploited cortical oscillations through analysis of ongoing EEG and musical feature during freely listening to music. We used a data-driven method that combined music information retrieval with spatial Fourier Independent Components Analysis (spatial Fourier-ICA) to probe the interplay between the spatial profiles and the spectral patterns of the brain network emerging from music listening. Correlation analysis was performed between time courses of brain networks extracted from EEG data and musical feature time series extracted from music stimuli to derive the musical feature related oscillatory patterns in the listening brain. We found brain networks of musical feature processing were frequency-dependent. Musical feature time series, especially fluctuation centroid and key feature, were associated with an increased beta activation in the bilateral superior temporal gyrus. An increased alpha oscillation in the bilateral occipital cortex emerged during music listening, which was consistent with alpha functional suppression hypothesis in task-irrelevant regions. We also observed an increased delta-beta oscillatory activity in the prefrontal cortex associated with musical feature processing. In addition to these findings, the proposed method seems valuable for characterizing the large-scale frequency-dependent brain activity engaged in musical feature processing.

Spatiotemporal dynamics of postoperative functional plasticity in patients with brain tumors in language areas

Postoperative functional neuroimaging provides a unique opportunity to investigate the neural mechanisms that facilitate language network reorganization. Previous studies in patients with low grade gliomas (LGGs) in language areas suggest that postoperative recovery is likely due to functional neuroplasticity in peritumoral and contra-tumoral healthy regions, but have attributed varying degrees of importance to specific regions. In this study, we used Magnetoencephalography (MEG) to investigate functional connectivity changes in peritumoral and contra-tumoral regions after brain tumor resection. MEG recordings of cortical activity during resting-state were obtained from 12 patients with LGGs in left-hemisphere language brain areas. MEG data were recorded before (Pre session), and 3 (Post_1 session) and 6 (Post_2 session) months after awake craniotomy. For each MEG session, we measured the functional connectivity of the peritumoral and contra-tumoral regions to the rest of the brain across the 1-100 Hz frequency band. We found that functional connectivity in the Post_1 and Post_2 sessions was higher than in the Pre session only in peritumoral regions and within the alpha frequency band. Functional connectivity in peritumoral regions did not differ between the Post_1 and Post_2 sessions. Alpha connectivity enhancement in peritumoral regions was observed in all patients regardless of the LGG location. Together, these results suggest that postoperative language functional reorganization occurs in peritumoral regions regardless of the location of the tumor and mostly develops within 3 months after surgery.