Assessing and improving the spatial accuracy in MEG source localization by depth-weighted minimum-norm estimates

XDL-ESI: Electrophysiological Sources Imaging via explainable deep learning framework with validation on simultaneous EEG and iEEG

Electroencephalography (EEG) or Magnetoencephalography (MEG) source imaging aims to estimate the underlying activated brain sources to explain the observed EEG/MEG recordings. Solving the inverse problem of EEG/MEG Source Imaging (ESI) is challenging due to its ill-posed nature. To achieve a unique solution, it is essential to apply sophisticated regularization constraints to restrict the solution space. Traditionally, the design of regularization terms is based on assumptions about the spatiotemporal structure of the underlying source dynamics. In this paper, we propose a novel paradigm for ESI via an Explainable Deep Learning framework, termed as XDL-ESI, which connects the iterative optimization algorithm with deep learning architecture by unfolding the iterative updates with neural network modules. The proposed framework has the advantages of (1) establishing a data-driven approach to model the source solution structure instead of using hand-crafted regularization terms; (2) improving the robustness of source solutions by introducing a topological loss that leverages the geometric spatial information applying varying penalties on distinct localization errors; (3) improving the reconstruction efficiency and interpretability as it inherits the advantages from both the iterative optimization algorithms (interpretability) and deep learning approaches (function approximation). The proposed XDL-ESI framework provides an efficient, accurate, and interpretable paradigm to solve the ESI inverse problem with satisfactory performance in both simulated data and real clinical data. Specially, this approach is further validated using simultaneous EEG and intracranial EEG (iEEG).

Functional and structural maturation of auditory cortex from 2 months to 2 years old

BACKGROUND

In school-age children, the myelination of the auditory radiation thalamocortical pathway is associated with the latency of auditory evoked responses, with the myelination of thalamocortical axons facilitating the rapid propagation of acoustic information. Little is known regarding this auditory system function-structure association in infants and toddlers.

METHODS AND PARTICIPANTS

The present study tested the hypothesis that maturation of auditory radiation white-matter microstructure (e.g., fractional anisotropy (FA); measured using diffusion-weighted MRI) is associated with the latency of the infant auditory response (the P2m response, measured using magnetoencephalography, MEG) in a cross-sectional (N = 47, 2 to 24 months, 19 females) as well as longitudinal cohort (N = 18, 2 to 29 months, 8 females) of typically developing infants and toddlers. Of 18 longitudinal infants, 2 infants had data from 3 timepoints and 16 infants had data from 2 timepoints.

RESULTS

In the cross-sectional sample, non-linear maturation of P2m latency and auditory radiation diffusion measures were observed. Auditory radiation diffusion accounted for significant variance in P2m latency, even after removing the variance associated with age in both P2m latency and auditory radiation diffusion measures. In the longitudinal sample, latency and FA associations could be observed at the level of a single child.

CONCLUSIONS

Findings provide strong support for the hypothesis that an increase in thalamocortical neural conduction velocity, due to increased axon diameter and/or myelin maturation, contributes to a decrease in the infant P2m auditory evoked response latency.

SIGNIFICANCE

Infant multimodal brain imaging identifies brain mechanisms contributing to the rapid changes in neural circuit activity during the first two years of life.

Spelling principles matter: An ERP study investigating the processing of different types of pseudohomophones

Spelling in any writing system is governed by fundamental principles. We examined the processing of two types of pseudohomophones constructed from words whose spellings are based on different principles - on the traditional principle of writing, requiring memorization of their spelling, and on the morphological principle, allowing the determination of their spelling from another word with the same morpheme (root) to examine the dependence of the occurrence of orthography-phonology conflict on spelling principles. Event-related potentials were recorded from 22 volunteers during silent reading. Pseudohomophones based on the morphological principle increased the N400 amplitude, emphasizing semantic and morphological processing importance. The P600 component showed significant effects for differentiating words and pseudohomophones based on the traditional principle, predominantly indicating the involvement of memory and reanalysis processes. Source reconstruction demonstrates that both pseudohomophones activate the left inferior frontal gyrus. However, pseudohomophones based on the traditional principle additionally activate the right and left postcentral gyrus, indicating the involvement of additional areas in the differentiation process. The earlier differences for stimuli based on the morphological principle indicate access to smaller units (morphemes), whereas stimuli based on the traditional principle require whole word processing. Our findings underscore the significant role of spelling principles in orthographic processing.

Localized estimation of event-related neural source activity from simultaneous MEG-EEG with a recurrent neural network

Estimating intracranial current sources underlying the electromagnetic signals observed from extracranial sensors is a perennial challenge in non-invasive neuroimaging. Established solutions to this inverse problem treat time samples independently without considering the temporal dynamics of event-related brain processes. This paper describes current source estimation from simultaneously recorded magneto- and electro-encephalography (MEEG) using a recurrent neural network (RNN) that learns sequential relationships from neural data. The RNN was trained in two phases: (1) pre-training and (2) transfer learning with L1 regularization applied to the source estimation layer. Performance of using scaled labels derived from MEEG, magnetoencephalography (MEG), or electroencephalography (EEG) were compared, as were results from volumetric source space with free dipole orientation and surface source space with fixed dipole orientation. Exact low-resolution electromagnetic tomography (eLORETA) and mixed-norm L1/L2 (MxNE) source estimation methods were also applied to these data for comparison with the RNN method. The RNN approach outperformed other methods in terms of output signal-to-noise ratio, correlation and mean-squared error metrics evaluated against reference event-related field (ERF) and event-related potential (ERP) waveforms. Using MEEG labels with fixed-orientation surface sources produced the most consistent estimates. To estimate sources of ERF and ERP waveforms, the RNN generates temporal dynamics within its internal computational units, driven by sequential structure in neural data used as training labels. It thus provides a data-driven model of computational transformations from psychophysiological events into corresponding event-related neural signals, which is unique among MEEG source reconstruction solutions.

Changes in electrophysiological aperiodic activity during cognitive control in Parkinson's disease

Cognitive symptoms in Parkinson's disease are common and can significantly affect patients' quality of life. Therefore, there is an urgent clinical need to identify a signature derived from behavioural and/or neuroimaging indicators that could predict which patients are at increased risk for early and rapid cognitive decline. Recently, converging evidence identified that aperiodic activity of the EEG reflects meaningful physiological information associated with age, development, cognitive and perceptual states or pathologies. In this study, we aimed to investigate aperiodic activity in Parkinson's disease during cognitive control and characterize its possible association with behaviour. Here, we recorded high-density EEG in 30 healthy controls and 30 Parkinson's disease patients during a Simon task. We analysed task-related behavioural data in the context of the activation-suppression model and extracted aperiodic parameters (offset, exponent) at both scalp and source levels. Our results showed lower behavioural performances in cognitive control as well as higher offsets in patients in the parieto-occipital areas, suggesting increased excitability in Parkinson's disease. A small congruence effect on aperiodic parameters in pre- and post-central brain areas was also found, possibly associated with task execution. Significant differences in aperiodic parameters between the resting-state, pre- and post-stimulus phases were seen across the whole brain, which confirmed that the observed changes in aperiodic activity are linked to task execution. No correlation was found between aperiodic activity and behaviour or clinical features. Our findings provide evidence that EEG aperiodic activity in Parkinson's disease is characterized by greater offsets, and that aperiodic parameters differ depending on arousal state. However, our results do not support the hypothesis that the behaviour-related differences observed in Parkinson's disease are related to aperiodic changes. Overall, this study highlights the importance of considering aperiodic activity contributions in brain disorders and further investigating the relationship between aperiodic activity and behaviour.

Altered neural synchronization in response to 2 Hz amplitude-modulated tones in the auditory cortex of children with Autism Spectrum Disorder: An MEG study

OBJECTIVE

Some studies have hypothesized that atypical neural synchronization at the delta frequency band in the auditory cortex is associated with phonological and language skills in children with Autism Spectrum Disorder (ASD), but it is still poorly understood. This study investigated this neural activity and addressed the relationships between auditory response and behavioral measures of children with ASD.

METHODS

We used magnetoencephalography and individual brain models to investigate 2 Hz Auditory Steady-State Response (ASSR) in 20 primary-school-aged children with ASD and 20 age-matched typically developing (TD) controls.

RESULTS

First, we found a between-group difference in the localization of the auditory response, so as the topology of 2 Hz ASSR was more superior and posterior in TD children when comparing to children with ASD. Second, the power of 2 Hz ASSR was reduced in the ASD group. Finally, we observed a significant association between the amplitude of neural response and language skills in children with ASD.

CONCLUSIONS

The study provided the evidence of reduced neural response in children with ASD and its relation to language skills.

SIGNIFICANCE

These findings may inform future interventions targeting auditory and language impairments in ASD population.

Convolutional neural networks can identify brain interactions involved in decoding spatial auditory attention

Human listeners have the ability to direct their attention to a single speaker in a multi-talker environment. The neural correlates of selective attention can be decoded from a single trial of electroencephalography (EEG) data. In this study, leveraging the source-reconstructed and anatomically-resolved EEG data as inputs, we sought to employ CNN as an interpretable model to uncover task-specific interactions between brain regions, rather than simply to utilize it as a black box decoder. To this end, our CNN model was specifically designed to learn pairwise interaction representations for 10 cortical regions from five-second inputs. By exclusively utilizing these features for decoding, our model was able to attain a median accuracy of 77.56% for within-participant and 65.14% for cross-participant classification. Through ablation analysis together with dissecting the features of the models and applying cluster analysis, we were able to discern the presence of alpha-band-dominated inter-hemisphere interactions, as well as alpha- and beta-band dominant interactions that were either hemisphere-specific or were characterized by a contrasting pattern between the right and left hemispheres. These interactions were more pronounced in parietal and central regions for within-participant decoding, but in parietal, central, and partly frontal regions for cross-participant decoding. These findings demonstrate that our CNN model can effectively utilize features known to be important in auditory attention tasks and suggest that the application of domain knowledge inspired CNNs on source-reconstructed EEG data can offer a novel computational framework for studying task-relevant brain interactions.

Emotion brain network topology in healthy subjects following passive listening to different auditory stimuli

A large body of research establishes the efficacy of musical intervention in many aspects of physical, cognitive, communication, social, and emotional rehabilitation. However, the underlying neural mechanisms for musical therapy remain elusive. This study aimed to investigate the potential neural correlates of musical therapy, focusing on the changes in the topology of emotion brain network. To this end, a Bayesian statistical approach and a cross-over experimental design were employed together with two resting-state magnetoencephalography (MEG) as controls. MEG recordings of 30 healthy subjects were acquired while listening to five auditory stimuli in random order. Two resting-state MEG recordings of each subject were obtained, one prior to the first stimulus (pre) and one after the final stimulus (post). Time series at the level of brain regions were estimated using depth-weighted minimum norm estimation (wMNE) source reconstruction method and the functional connectivity between these regions were computed. The resultant connectivity matrices were used to derive two topological network measures: transitivity and global efficiency which are important in gauging the functional segregation and integration of brain network respectively. The differences in these measures between pre- and post-stimuli resting MEG were set as the equivalence regions. We found that the network measures under all auditory stimuli were equivalent to the resting state network measures in all frequency bands, indicating that the topology of the functional brain network associated with emotional regulation in healthy subjects remains unchanged following these auditory stimuli. This suggests that changes in the emotion network topology may not be the underlying neural mechanism of musical therapy. Nonetheless, further studies are required to explore the neural mechanisms of musical interventions especially in the populations with neuropsychiatric disorders.

EEG/MEG source imaging of deep brain activity within the maximum entropy on the mean framework: Simulations and validation in epilepsy

Electro/Magneto-EncephaloGraphy (EEG/MEG) source imaging (EMSI) of epileptic activity from deep generators is often challenging due to the higher sensitivity of EEG/MEG to superficial regions and to the spatial configuration of subcortical structures. We previously demonstrated the ability of the coherent Maximum Entropy on the Mean (cMEM) method to accurately localize the superficial cortical generators and their spatial extent. Here, we propose a depth-weighted adaptation of cMEM to localize deep generators more accurately. These methods were evaluated using realistic MEG/high-density EEG (HD-EEG) simulations of epileptic activity and actual MEG/HD-EEG recordings from patients with focal epilepsy. We incorporated depth-weighting within the MEM framework to compensate for its preference for superficial generators. We also included a mesh of both hippocampi, as an additional deep structure in the source model. We generated 5400 realistic simulations of interictal epileptic discharges for MEG and HD-EEG involving a wide range of spatial extents and signal-to-noise ratio (SNR) levels, before investigating EMSI on clinical HD-EEG in 16 patients and MEG in 14 patients. Clinical interictal epileptic discharges were marked by visual inspection. We applied three EMSI methods: cMEM, depth-weighted cMEM and depth-weighted minimum norm estimate (MNE). The ground truth was defined as the true simulated generator or as a drawn region based on clinical information available for patients. For deep sources, depth-weighted cMEM improved the localization when compared to cMEM and depth-weighted MNE, whereas depth-weighted cMEM did not deteriorate localization accuracy for superficial regions. For patients' data, we observed improvement in localization for deep sources, especially for the patients with mesial temporal epilepsy, for which cMEM failed to reconstruct the initial generator in the hippocampus. Depth weighting was more crucial for MEG (gradiometers) than for HD-EEG. Similar findings were found when considering depth weighting for the wavelet extension of MEM. In conclusion, depth-weighted cMEM improved the localization of deep sources without or with minimal deterioration of the localization of the superficial sources. This was demonstrated using extensive simulations with MEG and HD-EEG and clinical MEG and HD-EEG for epilepsy patients.

Test-retest reliability of corticokinematic coherence in young children with cerebral palsy: An observational longitudinal study

OBJECTIVES

To assess the test-retest reliability of the corticokinematic coherence (CKC), an electrophysiological marker of proprioception, in children with cerebral palsy (CP).

METHODS

Electroencephalography (EEG) signals from 15 children with unilateral or bilateral CP aged 23 to 53 months were recorded in two sessions 3 months apart using 128-channel EEG caps. During each session, children's fingers were moved at 2 Hz by an experimenter, in separate recordings for the more-affected (MA) and less-affected (LA) hands. The CKC was computed at the electrode and source levels, at movement frequency F0 (2 Hz) and its first harmonic F1 (4 Hz). A two-way mixed-effects model intraclass-correlation coefficient (ICC) was computed for the maximum CKC strength across electrodes at F0 and F1 obtained during the two sessions.

RESULTS

ICC of the CKC strength acquired from LA and MA hands pooled together were respectively 0.51 (95% CI: 0.30-0.68) at F0 and 0.96 (95% CI: 0.93-0.98) at F1. The mean distances separating the CKC peaks in the source space at the two evaluation times were in the order of a centimeter.

CONCLUSION

CKC is a robust electrophysiologic marker to study the longitudinal changes in cortical processing of proprioceptive afferences in young children with CP.

The maturation of infant and toddler visual cortex neural activity and associations with fine motor performance

Our understanding of how visual cortex neural processes mature during infancy and toddlerhood is limited. Using magnetoencephalography (MEG), the present study investigated the development of visual evoked responses (VERs) in both cross-sectional and longitudinal samples of infants and toddlers 2 months to 3 years. Brain space analyses focused on N1m and P1m latency, as well as the N1m-to-P1m amplitude. Associations between VER measures and developmental quotient (DQ) scores in the cognitive/visual and fine motor domains were also examined. Results showed a nonlinear decrease in N1m and P1m latency as a function of age, characterized by rapid changes followed by slower progression, with the N1m latency plateauing at 6-7 months and the P1m latency plateauing at 8-9 months. The N1m-to-P1m amplitude also exhibited a non-linear decrease, with strong responses observed in younger infants (∼2-3 months) and then a gradual decline. Associations between N1m and P1m latency and fine motor DQ scores were observed, suggesting that infants with faster visual processing may be better equipped to perform fine motor tasks. The present findings advance our understanding of the maturation of the infant visual system and highlight the relationship between the maturation of visual system and fine motor skills.

Highlights

The infant N1m and P1m latency shows a nonlinear decrease.N1m latency decreases precede P1m latency decreases.N1m-to-P1m amplitude shows a nonlinear decrease, with stronger responses in younger than older infants.N1m and P1m latency are associated with fine motor DQ.

Functional and structural maturation of auditory cortex from 2 months to 2 years old

In school-age children, the myelination of the auditory radiation thalamocortical pathway is associated with the latency of auditory evoked responses, with the myelination of thalamocortical axons facilitating the rapid propagation of acoustic information. Little is known regarding this auditory system function-structure association in infants and toddlers. The present study tested the hypothesis that maturation of auditory radiation white-matter microstructure (e.g., fractional anisotropy (FA); measured using diffusion-weighted MRI) is associated with the latency of the infant auditory response (P2m measured using magnetoencephalography, MEG) in a cross-sectional (2 to 24 months) as well as longitudinal cohort (2 to 29 months) of typically developing infants and toddlers. In the cross-sectional sample, non-linear maturation of P2m latency and auditory radiation diffusion measures were observed. After removing the variance associated with age in both P2m latency and auditory radiation diffusion measures, auditory radiation still accounted for significant variance in P2m latency. In the longitudinal sample, latency and FA associations could be observed at the level of a single child. Findings provide strong support for a contribution of auditory radiation white matter to rapid cortical auditory encoding processes in infants.

Abnormalities in both stimulus-induced and baseline MEG alpha oscillations in the auditory cortex of children with Autism Spectrum Disorder

The neurobiology of Autism Spectrum Disorder (ASD) is hypothetically related to the imbalance between neural excitation (E) and inhibition (I). Different studies have revealed that alpha-band (8-12 Hz) activity in magneto- and electroencephalography (MEG and EEG) may reflect E and I processes and, thus, can be of particular interest in ASD research. Previous findings indicated alterations in event-related and baseline alpha activity in different cortical systems in individuals with ASD, and these abnormalities were associated with core and co-occurring conditions of ASD. However, the knowledge on auditory alpha oscillations in this population is limited. This MEG study investigated stimulus-induced (Event-Related Desynchronization, ERD) and baseline alpha-band activity (both periodic and aperiodic) in the auditory cortex and also the relationships between these neural activities and behavioral measures of children with ASD. Ninety amplitude-modulated tones were presented to two groups of children: 20 children with ASD (5 girls, Mage = 10.03, SD = 1.7) and 20 typically developing controls (9 girls, Mage = 9.11, SD = 1.3). Children with ASD had a bilateral reduction of alpha-band ERD, reduced baseline aperiodic-adjusted alpha power, and flattened aperiodic exponent in comparison to TD children. Moreover, lower raw baseline alpha power and aperiodic offset in the language-dominant left auditory cortex were associated with better language skills of children with ASD measured in formal assessment. The findings highlighted the alterations of E / I balance metrics in response to basic auditory stimuli in children with ASD and also provided evidence for the contribution of low-level processing to language difficulties in ASD.

Changes in visually and auditory attended audiovisual speech processing in cochlear implant users: A longitudinal ERP study

Limited auditory input, whether caused by hearing loss or by electrical stimulation through a cochlear implant (CI), can be compensated by the remaining senses. Specifically for CI users, previous studies reported not only improved visual skills, but also altered cortical processing of unisensory visual and auditory stimuli. However, in multisensory scenarios, it is still unclear how auditory deprivation (before implantation) and electrical hearing experience (after implantation) affect cortical audiovisual speech processing. Here, we present a prospective longitudinal electroencephalography (EEG) study which systematically examined the deprivation- and CI-induced alterations of cortical processing of audiovisual words by comparing event-related potentials (ERPs) in postlingually deafened CI users before and after implantation (five weeks and six months of CI use). A group of matched normal-hearing (NH) listeners served as controls. The participants performed a word-identification task with congruent and incongruent audiovisual words, focusing their attention on either the visual (lip movement) or the auditory speech signal. This allowed us to study the (top-down) attention effect on the (bottom-up) sensory cortical processing of audiovisual speech. When compared to the NH listeners, the CI candidates (before implantation) and the CI users (after implantation) exhibited enhanced lipreading abilities and an altered cortical response at the N1 latency range (90-150 ms) that was characterized by a decreased theta oscillation power (4-8 Hz) and a smaller amplitude in the auditory cortex. After implantation, however, the auditory-cortex response gradually increased and developed a stronger intra-modal connectivity. Nevertheless, task efficiency and activation in the visual cortex was significantly modulated in both groups by focusing attention on the visual as compared to the auditory speech signal, with the NH listeners additionally showing an attention-dependent decrease in beta oscillation power (13-30 Hz). In sum, these results suggest remarkable deprivation effects on audiovisual speech processing in the auditory cortex, which partially reverse after implantation. Although even experienced CI users still show distinct audiovisual speech processing compared to NH listeners, pronounced effects of (top-down) direction of attention on (bottom-up) audiovisual processing can be observed in both groups. However, NH listeners but not CI users appear to show enhanced allocation of cognitive resources in visually as compared to auditory attended audiovisual speech conditions, which supports our behavioural observations of poorer lipreading abilities and reduced visual influence on audition in NH listeners as compared to CI users.

Impaired Cortical Tracking of Speech in Children with Developmental Language Disorder

In developmental language disorder (DLD), learning to comprehend and express oneself with spoken language is impaired, but the reason for this remains unknown. Using millisecond-scale magnetoencephalography recordings combined with machine learning models, we investigated whether the possible neural basis of this disruption lies in poor cortical tracking of speech. The stimuli were common spoken Finnish words (e.g., dog, car, hammer) and sounds with corresponding meanings (e.g., dog bark, car engine, hammering). In both children with DLD (10 boys and 7 girls) and typically developing (TD) control children (14 boys and 3 girls), aged 10-15 years, the cortical activation to spoken words was best modeled as time-locked to the unfolding speech input at ∼100 ms latency between sound and cortical activation. Amplitude envelope (amplitude changes) and spectrogram (detailed time-varying spectral content) of the spoken words, but not other sounds, were very successfully decoded based on time-locked brain responses in bilateral temporal areas; based on the cortical responses, the models could tell at ∼75-85% accuracy which of the two sounds had been presented to the participant. However, the cortical representation of the amplitude envelope information was poorer in children with DLD compared with TD children at longer latencies (at ∼200-300 ms lag). We interpret this effect as reflecting poorer retention of acoustic-phonetic information in short-term memory. This impaired tracking could potentially affect the processing and learning of words as well as continuous speech. The present results offer an explanation for the problems in language comprehension and acquisition in DLD.

Binocularly suppressed stimuli induce brain activities related to aesthetic emotions

Introduction

Aesthetic emotions are a class of emotions aroused by evaluating aesthetically appealing objects or events. While evolutionary aesthetics suggests the adaptive roles of these emotions, empirical assessments are lacking. Previous neuroscientific studies have demonstrated that visual stimuli carrying evolutionarily important information induce neural responses even when presented non-consciously. To examine the evolutionary importance of aesthetic emotions, we conducted a neuroscientific study using magnetoencephalography (MEG) to measure induced neural responses to non-consciously presented portrait paintings categorised as biological and non-biological and examined associations between the induced responses and aesthetic ratings.

Methods

MEG and pre-rating data were collected from 23 participants. The pre-rating included visual analogue scales for object saliency, facial saliency, liking, and beauty scores, in addition to 'biologi-ness,' which was used for subcategorising stimuli into biological and non-biological. The stimuli were presented non-consciously using a continuous flash suppression paradigm or consciously using binocular presentation without flashing masks, while dichotomic behavioural responses were obtained (beauty or non-beauty). Time-frequency decomposed MEG data were used for correlation analysis with pre-rating scores for each category.

Results

Behavioural data revealed that saliency scores of non-consciously presented stimuli influenced dichotomic responses (beauty or non-beauty). MEG data showed that non-consciously presented portrait paintings induced spatiotemporally distributed low-frequency brain activities associated with aesthetic ratings, which were distinct between the biological and non-biological categories and conscious and non-conscious conditions.

Conclusion

Aesthetic emotion holds evolutionary significance for humans. Neural pathways are sensitive to visual images that arouse aesthetic emotion in distinct ways for biological and non-biological categories, which are further influenced by consciousness. These differences likely reflect the diversity in mechanisms of aesthetic processing, such as processing fluency, active elaboration, and predictive processing. The aesthetic processing of non-conscious stimuli appears to be characterised by fluency-driven affective processing, while top-down regulatory processes are suppressed. This study provides the first empirical evidence supporting the evolutionary significance of aesthetic processing.

Simultaneous MEG and EEG source imaging of electrophysiological activity in response to acute transcranial photobiomodulation

Introduction

Transcranial photobiomodulation (tPBM) is a non-invasive neuromodulation technique that improves human cognition. The effects of tPBM of the right forehead on neurophysiological activity have been previously investigated using EEG in sensor space. However, the spatial resolution of these studies is limited. Magnetoencephalography (MEG) is known to facilitate a higher spatial resolution of brain source images. This study aimed to image post-tPBM effects in brain space based on both MEG and EEG measurements across the entire human brain.

Methods

MEG and EEG scans were concurrently acquired for 6 min before and after 8-min of tPBM delivered using a 1,064-nm laser on the right forehead of 25 healthy participants. Group-level changes in both the MEG and EEG power spectral density with respect to the baseline (pre-tPBM) were quantified and averaged within each frequency band in the sensor space. Constrained modeling was used to generate MEG and EEG source images of post-tPBM, followed by cluster-based permutation analysis for family wise error correction (p < 0.05).

Results

The 8-min tPBM enabled significant increases in alpha (8-12 Hz) and beta (13-30 Hz) powers across multiple cortical regions, as confirmed by MEG and EEG source images. Moreover, tPBM-enhanced oscillations in the beta band were located not only near the stimulation site but also in remote cerebral regions, including the frontal, parietal, and occipital regions, particularly on the ipsilateral side.

Discussion

MEG and EEG results shown in this study demonstrated that tPBM modulates neurophysiological activity locally and in distant cortical areas. The EEG topographies reported in this study were consistent with previous observations. This study is the first to present MEG and EEG evidence of the electrophysiological effects of tPBM in the brain space, supporting the potential utility of tPBM in treating neurological diseases through the modulation of brain oscillations.

Neural correlates of retrospective memory confidence during face-name associative learning

The ability to accurately assess one's own memory performance during learning is essential for adaptive behavior, but the brain mechanisms underlying this metamemory function are not well understood. We investigated the neural correlates of memory accuracy and retrospective memory confidence in a face-name associative learning task using magnetoencephalography in healthy young adults (n = 32). We found that high retrospective confidence was associated with stronger occipital event-related fields during encoding and widespread event-related fields during retrieval compared to low confidence. On the other hand, memory accuracy was linked to medial temporal activities during both encoding and retrieval, but only in low-confidence trials. A decrease in oscillatory power at alpha/beta bands in the parietal regions during retrieval was associated with higher memory confidence. In addition, representational similarity analysis at the single-trial level revealed distributed but differentiable neural activities associated with memory accuracy and confidence during both encoding and retrieval. In summary, our study unveiled distinct neural activity patterns related to memory confidence and accuracy during associative learning and underscored the crucial role of parietal regions in metamemory.

No Differences in Auditory Steady-State Responses in Children with Autism Spectrum Disorder and Typically Developing Children

Auditory steady-state response (ASSR) has been studied as a potential biomarker for abnormal auditory sensory processing in autism spectrum disorder (ASD), with mixed results. Motivated by prior somatosensory findings of group differences in inter-trial coherence (ITC) between ASD and typically developing (TD) individuals at twice the steady-state stimulation frequency, we examined ASSR at 25 and 50 as well as 43 and 86 Hz in response to 25-Hz and 43-Hz auditory stimuli, respectively, using magnetoencephalography. Data were recorded from 22 ASD and 31 TD children, ages 6-17 years. ITC measures showed prominent ASSRs at the stimulation and double frequencies, without significant group differences. These results do not support ASSR as a robust ASD biomarker of abnormal auditory processing in ASD. Furthermore, the previously observed atypical double-frequency somatosensory response in ASD did not generalize to the auditory modality. Thus, the hypothesis about modality-independent abnormal local connectivity in ASD was not supported.

Dorsal brain activity reflects the severity of menopausal symptoms

OBJECTIVE

The severity of menopausal symptoms, despite being triggered by hormonal imbalance, does not directly correspond to hormone levels in the blood; thus, the level of unpleasantness is assessed using subjective questionnaires in clinical practice. To provide better treatments, alternative objective assessments have been anticipated to support medical interviews and subjective assessments. This study aimed to develop a new objective measurement for assessing unpleasantness.

METHODS

Fourteen participants with menopausal symptoms and two age-matched participants who visited our outpatient section were enrolled. Resting-state brain activity was measured using magnetoencephalography. The level of unpleasantness of menopausal symptoms was measured using the Kupperman Kohnenki Shogai Index. The blood level of follicle-stimulating hormone and luteinizing hormone were also measured. Correlation analyses were performed between the oscillatory power of brain activity, index score, and hormone levels.

RESULTS

The level of unpleasantness of menopausal symptoms was positively correlated with high-frequency oscillatory powers in the parietal and bordering cortices (alpha; P = 0.016, beta; P = 0.015, low gamma; P = 0.010). The follicle-stimulating hormone blood level was correlated with high-frequency oscillatory powers in the dorsal part of the cortex (beta; P = 0.008, beta; P = 0.005, low gamma; P = 0.017), whereas luteinizing hormone blood level was not correlated.

CONCLUSION

Resting-state brain activity can serve as an objective measurement of unpleasantness associated with menopausal symptoms, which aids the selection of appropriate treatment and monitors its outcome.

CiftiStorm pipeline: facilitating reproducible EEG/MEG source connectomics

We present CiftiStorm, an electrophysiological source imaging (ESI) pipeline incorporating recently developed methods to improve forward and inverse solutions. The CiftiStorm pipeline produces Human Connectome Project (HCP) and megconnectome-compliant outputs from dataset inputs with varying degrees of spatial resolution. The input data can range from low-sensor-density electroencephalogram (EEG) or magnetoencephalogram (MEG) recordings without structural magnetic resonance imaging (sMRI) to high-density EEG/MEG recordings with an HCP multimodal sMRI compliant protocol. CiftiStorm introduces a numerical quality control of the lead field and geometrical corrections to the head and source models for forward modeling. For the inverse modeling, we present a Bayesian estimation of the cross-spectrum of sources based on multiple priors. We facilitate ESI in the T1w/FSAverage32k high-resolution space obtained from individual sMRI. We validate this feature by comparing CiftiStorm outputs for EEG and MRI data from the Cuban Human Brain Mapping Project (CHBMP) acquired with technologies a decade before the HCP MEG and MRI standardized dataset.

Connectome spectrum electromagnetic tomography: A method to reconstruct electrical brain source networks at high-spatial resolution

Connectome spectrum electromagnetic tomography (CSET) combines diffusion MRI-derived structural connectivity data with well-established graph signal processing tools to solve the M/EEG inverse problem. Using simulated EEG signals from fMRI responses, and two EEG datasets on visual-evoked potentials, we provide evidence supporting that (i) CSET captures realistic neurophysiological patterns with better accuracy than state-of-the-art methods, (ii) CSET can reconstruct brain responses more accurately and with more robustness to intrinsic noise in the EEG signal. These results demonstrate that CSET offers high spatio-temporal accuracy, enabling neuroscientists to extend their research beyond the current limitations of low sampling frequency in functional MRI and the poor spatial resolution of M/EEG.

OPM-MEG Measuring Phase Synchronization on Source Time Series: Application in Rhythmic Median Nerve Stimulation

The magnetoencephalogram (MEG) based on array optically pumped magnetometers (OPMs) has the potential of replacing conventional cryogenic superconducting quantum interference device. Phase synchronization is a common method for measuring brain oscillations and functional connectivity. Verifying the feasibility and fidelity of OPM-MEG in measuring phase synchronization will help its widespread application in the study of aforementioned neural mechanisms. The analysis method on source-level time series can weaken the influence of instantaneous field spread effect. In this paper, the OPM-MEG was used for measuring the evoked responses of 20Hz rhythmic and arrhythmic median nerve stimulation, and the inter-trial phase synchronization (ITPS) and inter-reginal phase synchronization (IRPS) of primary somatosensory cortex (SI) and secondary somatosensory cortex (SII) were analysed. The results find that under rhythmic condition, the evoked responses of SI and SII show continuous oscillations and the effect of resetting phase. The values of ITPS and IRPS significantly increase at the stimulation frequency of 20Hz and its harmonic of 40Hz, whereas the arrhythmic stimulation does not exhibit this phenomenon. Moreover, in the initial stage of stimulation, the ITPS and IRPS values are significantly higher at Mu rhythm in the rhythmic condition compared to arrhythmic. In conclusion, the results demonstrate the ability of OPM-MEG in measuring phase pattern and functional connectivity on source-level, and may also prove beneficial for the study on the mechanism of rhythmic stimulation therapy for rehabilitation.

Is source-resolved magnetoencephalographic mismatch negativity a viable biomarker for early psychosis?

Mismatch negativity (MMN) is an auditory event-related response reflecting the pre-attentive detection of novel stimuli and is a biomarker of cortical dysfunction in schizophrenia (SZ). MMN to pitch (pMMN) and to duration (dMMN) deviant stimuli are impaired in chronic SZ, but it is less clear if MMN is reduced in first-episode SZ, with inconsistent findings in scalp-level EEG studies. Here, we investigated the neural generators of pMMN and dMMN with MEG recordings in 26 first-episode schizophrenia spectrum (FEsz) and 26 matched healthy controls (C). We projected MEG inverse solutions into precise functionally meaningful auditory cortex areas. MEG-derived MMN sources were in bilateral primary auditory cortex (A1) and belt areas. In A1, pMMN FEsz reduction showed a trend towards statistical significance (F(1,50) = 3.31; p = .07), and dMMN was reduced in FEsz (F(1,50) = 4.11; p = .04). Hypothesis-driven comparisons at each hemisphere revealed dMMN reduction in FEsz occurred in the left (t(56) = 2.23; p = .03; d = .61) but not right (t(56) = 1.02; p = .31; d = .28) hemisphere, with a moderate effect size. The added precision of MEG source solution with high-resolution MRI and parcellation of A1 may be requisite to detect the emerging pathophysiology and indicates a critical role for left hemisphere pathology at psychosis onset. However, the moderate effect size in left A1, albeit larger than reported in scalp MMN meta-analyses, casts doubt on the clinical utility of MMN for differential diagnosis, as a majority of patients will overlap with the healthy individual's distribution.

Verb generation for presurgical mapping: Gaining specificity

Verb generation is among the most frequently used tasks in presurgical mapping. Because this task involves many processes, the overall brain effects are not specific. While it is necessary to identify the whole network involving noun comprehension or semantic retrieval and lexical selection to produce the verb, isolation of those components is also crucial. Here, we present data from four patients undergoing presurgical brain mapping. The study implied a reanalysis of magnetoencephalography data with a recategorization of the used items. It aimed to extract the task component that relies on the inferior frontal gyrus (IFG). The task could be applied with higher specificity when targeting frontal areas. For that, we based item classification on the selection demands imposed by the noun. It is a robust finding that the IFG carries out this selection and that a quantitative index can be calculated for each noun, which depends on the selection effort (Proceedings of the National Academy of Sciences of the United States of America, 1997; 94(26):14792-14797, Proceedings of the National Academy of Sciences of the United States of America, 1998; 95(26):15855-15860). Data showed focality and specificity, with a correlation between this derived index and source activations in the inferior frontal gyrus for all patients. Strikingly, we detected when the right-hemisphere homologue area was involved in the selection process in two patients showing reorganization or language right lateralization. The present data are a step towards a dissection of broad specific tasks frequently used in presurgical protocols.

Deep EEG source localization via EMD-based fMRI high spatial frequency

Brain imaging with a high-spatiotemporal resolution is crucial for accurate brain-function mapping. Electroencephalography (EEG) and functional Magnetic Resonance Imaging (fMRI) are two popular neuroimaging modalities with complementary features that record brain function with high temporal and spatial resolution, respectively. One popular non-invasive way to obtain data with both high spatial and temporal resolutions is to combine the fMRI activation map and EEG data to improve the spatial resolution of the EEG source localization. However, using the whole fMRI map may cause spurious results for the EEG source localization, especially for deep brain regions. Considering the head's conductivity, deep regions' sources with low activity are unlikely to be detected by the EEG electrodes at the scalp. In this study, we use fMRI's high spatial-frequency component to identify the local high-intensity activations that are most likely to be captured by the EEG. The 3D Empirical Mode Decomposition (3D-EMD), a data-driven method, is used to decompose the fMRI map into its spatial-frequency components. Different validation measurements for EEG source localization show improved performance for the EEG inverse-modeling informed by the fMRI's high-frequency spatial component compared to the fMRI-informed EEG source-localization methods. The level of improvement varies depending on the voxels' intensity and their distribution. Our experimental results also support this conclusion.

Perceptual Grouping in Autism Spectrum Disorder: An Exploratory Magnetoencephalography Study

Visual information is organised according to visual grouping principles. In visual grouping tasks individuals with ASD have shown equivocal performance. We explored neural correlates of Gestalt grouping in individuals with and without ASD. Neuromagnetic activity of individuals with (15) and without (18) ASD was compared during a visual grouping task testing grouping by proximity versus similarity. Individuals without ASD showed stronger evoked responses with earlier peaks in response to both grouping types indicating an earlier neuronal differentiation between grouping principles in individuals without ASD. In contrast, individuals with ASD showed particularly prolonged processing of grouping by similarity suggesting a high demand of neural resources. The neuronal processing differences found could explain less efficient grouping performance observed behaviourally in ASD.

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.

Oscillatory characteristics of resting-state magnetoencephalography reflect pathological and symptomatic conditions of cognitive impairment

Background

Dementia and mild cognitive impairment are characterised by symptoms of cognitive decline, which are typically assessed using neuropsychological assessments (NPAs), such as the Mini-Mental State Examination (MMSE) and Frontal Assessment Battery (FAB). Magnetoencephalography (MEG) is a novel clinical assessment technique that measures brain activities (summarised as oscillatory parameters), which are associated with symptoms of cognitive impairment. However, the relevance of MEG and regional cerebral blood flow (rCBF) data obtained using single-photon emission computed tomography (SPECT) has not been examined using clinical datasets. Therefore, this study aimed to investigate the relationships among MEG oscillatory parameters, clinically validated biomarkers computed from rCBF, and NPAs using outpatient data retrieved from hospital records.

Methods

Clinical data from 64 individuals with mixed pathological backgrounds were retrieved and analysed. MEG oscillatory parameters, including relative power (RP) from delta to high gamma bands, mean frequency, individual alpha frequency, and Shannon's spectral entropy, were computed for each cortical region. For SPECT data, three pathological parameters-'severity', 'extent', and 'ratio'-were computed using an easy z-score imaging system (eZIS). As for NPAs, the MMSE and FAB scores were retrieved.

Results

MEG oscillatory parameters were correlated with eZIS parameters. The eZIS parameters associated with Alzheimer's disease pathology were reflected in theta power augmentation and slower shift of the alpha peak. Moreover, MEG oscillatory parameters were found to reflect NPAs. Global slowing and loss of diversity in neural oscillatory components correlated with MMSE and FAB scores, whereas the associations between eZIS parameters and NPAs were sparse.

Conclusion

MEG oscillatory parameters correlated with both SPECT (i.e. eZIS) parameters and NPAs, supporting the clinical validity of MEG oscillatory parameters as pathological and symptomatic indicators. The findings indicate that various components of MEG oscillatory characteristics can provide valuable pathological and symptomatic information, making MEG data a rich resource for clinical examinations of patients with cognitive impairments. SPECT (i.e. eZIS) parameters showed no correlations with NPAs. The results contributed to a better understanding of the characteristics of electrophysiological and pathological examinations for patients with cognitive impairments, which will help to facilitate their co-use in clinical application, thereby improving patient care.

Measuring neuronal avalanches to inform brain-computer interfaces

Large-scale interactions among multiple brain regions manifest as bursts of activations called neuronal avalanches, which reconfigure according to the task at hand and, hence, might constitute natural candidates to design brain-computer interfaces (BCIs). To test this hypothesis, we used source-reconstructed magneto/electroencephalography during resting state and a motor imagery task performed within a BCI protocol. To track the probability that an avalanche would spread across any two regions, we built an avalanche transition matrix (ATM) and demonstrated that the edges whose transition probabilities significantly differed between conditions hinged selectively on premotor regions in all subjects. Furthermore, we showed that the topology of the ATMs allows task-decoding above the current gold standard. Hence, our results suggest that neuronal avalanches might capture interpretable differences between tasks that can be used to inform brain-computer interfaces.

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.

Cerebellar EEG source localization reveals age-related compensatory activity moderated by genetic risk for Alzheimer's disease

The apolipoprotein-E (APOE) ε4 allele is the greatest genetic risk factor for late-onset Alzheimer's disease (AD), but alone it is not sufficiently predictive. Because neuropathological changes associated with AD begin decades before cognitive symptoms, neuroimaging of healthy, cognitively intact ε4 carriers (ε4+) may enable early characterization of patterns associated with risk for future decline. Research in the cerebral cortex highlights a period of compensatory recruitment in elders and ε4+, which serves to maintain cognitive functioning. Yet, AD-related changes may occur even earlier in the cerebellum. Advances in electroencephalography (EEG) source localization now allow effective modeling of cerebellar activity. Importantly, healthy aging and AD are associated with declines in both cerebellar functions and executive functioning (EF). However, it is not known whether cerebellar activity can detect pre-symptomatic AD risk. Thus, the current study analyzed cerebellar EEG source localization during an EF-dependent stop-signal task (i.e., inhibitory control) in healthy, intact older adults (Mage  = 80 years; 20 ε4+, 25 ε4-). Task performance was comparable between groups. Older age predicted greater activity in left crus II and lobule VIIb during the P300 window (i.e., performance evaluation), consistent with age-related compensation. Age*ε4 moderations specifically showed that compensatory patterns were evident only in ε4-, suggesting that cerebellar compensatory resources may already be depleted in healthy ε4+ elders. Thus, the posterolateral cerebellum is sensitive to AD-related neural deficits in healthy elders. Characterization of these patterns may be essential for the earliest possible detection of AD risk, which would enable critical early intervention prior to symptom onset.

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.

The Characteristics of Power Spectral Density in Bipolar Disorder at the Resting State

Bipolar disorder (BD) is a common psychiatric disorder, but its pathophysiology is not fully elucidated. The current study focused on its electrophysiological characteristics, especially power spectral density (PSD). Resting state with eyes opened magnetoencephalography data were collected from 21 patients with BD and 22 healthy controls. The whole brain's PSD was calculated from source reconstructed waveforms at each frequency band (delta: 1-3 Hz, theta: 4-7 Hz, alpha: 8-12 Hz, low beta: 13-19 Hz, high beta: 20-29 Hz, and gamma: 30-80 Hz). We compared PSD values on the marked vertices at each frequency band between healthy and patient groups using a Mann-Whitney rank test to examine the relationship between significantly different PSD and clinical measures. The PSD in patients with BD was significantly decreased in lower frequency bands, mainly in the default mode network (DMN) areas (bilateral medial prefrontal cortex, bilateral precuneus, left inferior parietal lobe, and right temporal cortex in the alpha band) and salience network areas (SAL; left anterior insula [AI] at the delta band, anterior cingulate cortex at the theta band, and right AI at the alpha band). No significant differences in PSD were observed at low beta and high beta. PSD was not correlated with age or other clinical scales. Altered PSDs of the DMN and SAL were observed in the delta, theta, and alpha bands. These alterations contribute to the vulnerability of BD through the disturbance of self-referential mental activity and switching between the default mode and frontoparietal networks.

Connectivity-based Meta-Bands: A new approach for automatic frequency band identification in connectivity analyses

The majority of electroencephalographic (EEG) and magnetoencephalographic (MEG) studies filter and analyse neural signals in specific frequency ranges, known as "canonical" frequency bands. However, this segmentation, is not exempt from limitations, mainly due to the lack of adaptation to the neural idiosyncrasies of each individual. In this study, we introduce a new data-driven method to automatically identify frequency ranges based on the topological similarity of the frequency-dependent functional neural network. The resting-state neural activity of 195 cognitively healthy subjects from three different databases (MEG: 123 subjects; EEG1: 27 subjects; EEG2: 45 subjects) was analysed. In a first step, MEG and EEG signals were filtered with a narrow-band filter bank (1 Hz bandwidth) from 1 to 70 Hz with a 0.5 Hz step. Next, the connectivity in each of these filtered signals was estimated using the orthogonalized version of the amplitude envelope correlation to obtain the frequency-dependent functional neural network. Finally, a community detection algorithm was used to identify communities in the frequency domain showing a similar network topology. We have called this approach the "Connectivity-based Meta-Bands" (CMB) algorithm. Additionally, two types of synthetic signals were used to configure the hyper-parameters of the CMB algorithm. We observed that the classical approaches to band segmentation are partially aligned with the underlying network topologies at group level for the MEG signals, but they are missing individual idiosyncrasies that may be biasing previous studies, as revealed by our methodology. On the other hand, the sensitivity of EEG signals to reflect this underlying frequency-dependent network structure is limited, revealing a simpler frequency parcellation, not aligned with that defined by the "canonical" frequency bands. To the best of our knowledge, this is the first study that proposes an unsupervised band segmentation method based on the topological similarity of functional neural network across frequencies. This methodology fully accounts for subject-specific patterns, providing more robust and personalized analyses, and paving the way for new studies focused on exploring the frequency-dependent structure of brain connectivity.

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.

WRA-MTSI: A Robust Extended Source Imaging Algorithm Based on Multi-Trial EEG

OBJECTIVE

Reconstructing brain activities from electroencephalography (EEG) signals is crucial for studying brain functions and their abnormalities. However, since EEG signals are nonstationary and vulnerable to noise, brain activities reconstructed from single-trial EEG data are often unstable, and significant variability may occur across different EEG trials even for the same cognitive task.

METHODS

In an effort to leverage the shared information across the EEG data of multiple trials, this paper proposes a multi-trial EEG source imaging method based on Wasserstein regularization, termed WRA-MTSI. In WRA-MTSI, Wasserstein regularization is employed to perform multi-trial source distribution similarity learning, and the structured sparsity constraint is enforced to enable accurate estimation of the source extents, locations and time series. The resulting optimization problem is solved by a computationally efficient algorithm based on the alternating direction method of multipliers (ADMM).

RESULTS

Both numerical simulations and real EEG data analysis demonstrate that WRA-MTSI outperforms existing single-trial ESI methods (e.g., wMNE, LORETA, SISSY, and SBL) in mitigating the influence of artifacts in EEG data. Moreover, WRA-MTSI yields superior performance compared to other state-of-the-art multi-trial ESI methods (e.g., group lasso, the dirty model, and MTW) in estimating source extents.

CONCLUSION AND SIGNIFICANCE

WRA-MTSI may serve as an effective robust EEG source imaging method in the presence of multi-trial noisy EEG data.

Effects of noise and noise reduction on audiovisual speech perception in cochlear implant users: An ERP study

OBJECTIVE

Hearing with a cochlear implant (CI) is difficult in noisy environments, but the use of noise reduction algorithms, specifically ForwardFocus, can improve speech intelligibility. The current event-related potentials (ERP) study examined the electrophysiological correlates of this perceptual improvement.

METHODS

Ten bimodal CI users performed a syllable-identification task in auditory and audiovisual conditions, with syllables presented from the front and stationary noise presented from the sides. Brainstorm was used for spatio-temporal evaluation of ERPs.

RESULTS

CI users revealed an audiovisual benefit as reflected by shorter response times and greater activation in temporal and occipital regions at P2 latency. However, in auditory and audiovisual conditions, background noise hampered speech processing, leading to longer response times and delayed auditory-cortex-activation at N1 latency. Nevertheless, activating ForwardFocus resulted in shorter response times, reduced listening effort and enhanced superior-frontal-cortex-activation at P2 latency, particularly in audiovisual conditions.

CONCLUSIONS

ForwardFocus enhances speech intelligibility in audiovisual speech conditions by potentially allowing the reallocation of attentional resources to relevant auditory speech cues.

SIGNIFICANCE

This study shows for CI users that background noise and ForwardFocus differentially affect spatio-temporal cortical response patterns, both in auditory and audiovisual speech conditions.

Event-Related Desynchronization of MEG Alpha-Band Oscillations during Simultaneous Presentation of Audio and Visual Stimuli in Children with Autism Spectrum Disorder

Alpha-band (8-12 Hz) event-related desynchronization (ERD) or a decrease in alpha power in electro- and magnetoencephalography (EEG and MEG) reflects the involvement of a neural tissue in information processing. It is known that most children with autism spectrum disorder (ASD) have difficulties in information processing, and, thus, investigation of alpha oscillations is of particular interest in this population. Previous studies have demonstrated alterations in this neural activity in individuals with ASD; however, little is known about alpha ERD during simultaneous presentation of auditory and visual stimuli in children with and without ASD. As alpha oscillations are intimately related to attention, and attention deficit is one of the common co-occurring conditions of ASD, we predict that children with ASD can have altered alpha ERD in one of the sensory domains. In the present study, we used MEG to investigate alpha ERD in groups of 20 children with ASD and 20 age-matched typically developing controls. Simple amplitude-modulated tones were presented together with a fixation cross appearing on the screen. The results showed that children with ASD had a bilateral reduction in alpha-band ERD in the auditory but not visual cortex. Moreover, alterations in the auditory cortex were associated with a higher presence of autistic traits measured in behavioral assessment.

An artificial intelligence-based pipeline for automated detection and localisation of epileptic sources from magnetoencephalography

Objective.Magnetoencephalography (MEG) is a powerful non-invasive diagnostic modality for presurgical epilepsy evaluation. However, the clinical utility of MEG mapping for localising epileptic foci is limited by its low efficiency, high labour requirements, and considerable interoperator variability. To address these obstacles, we proposed a novel artificial intelligence-based automated magnetic source imaging (AMSI) pipeline for automated detection and localisation of epileptic sources from MEG data.Approach.To expedite the analysis of clinical MEG data from patients with epilepsy and reduce human bias, we developed an autolabelling method, a deep-learning model based on convolutional neural networks and a hierarchical clustering method based on a perceptual hash algorithm, to enable the coregistration of MEG and magnetic resonance imaging, the detection and clustering of epileptic activity, and the localisation of epileptic sources in a highly automated manner. We tested the capability of the AMSI pipeline by assessing MEG data from 48 epilepsy patients.Main results.The AMSI pipeline was able to rapidly detect interictal epileptiform discharges with 93.31% ± 3.87% precision based on a 35-patient dataset (with sevenfold patientwise cross-validation) and robustly rendered accurate localisation of epileptic activity with a lobar concordance of 87.18% against interictal and ictal stereo-electroencephalography findings in a 13-patient dataset. We also showed that the AMSI pipeline accomplishes the necessary processes and delivers objective results within a much shorter time frame (∼12 min) than traditional manual processes (∼4 h).Significance.The AMSI pipeline promises to facilitate increased utilisation of MEG data in the clinical analysis of patients with epilepsy.

Maturation of auditory cortex neural responses during infancy and toddlerhood

The infant auditory system rapidly matures across the first years of life, with a primary goal of obtaining ever-more-accurate real-time representations of the external world. Our understanding of how left and right auditory cortex neural processes develop during infancy, however, is meager, with few studies having the statistical power to detect potential hemisphere and sex differences in primary/secondary auditory cortex maturation. Using infant magnetoencephalography (MEG) and a cross-sectional study design, left and right auditory cortex P2m responses to pure tones were examined in 114 typically developing infants and toddlers (66 males, 2 to 24 months). Non-linear maturation of P2m latency was observed, with P2m latencies decreasing rapidly as a function of age during the first year of life, followed by slower changes between 12 and 24 months. Whereas in younger infants auditory tones were encoded more slowly in the left than right hemisphere, similar left and right P2m latencies were observed by ∼21 months of age due to faster maturation rate in the left than right hemisphere. No sex differences in the maturation of the P2m responses were observed. Finally, an earlier left than right hemisphere P2m latency predicted better language performance in older infants (12 to 24 months). Findings indicate the need to consider hemisphere when examining the maturation of auditory cortex neural activity in infants and toddlers and show that the pattern of left-right hemisphere P2m maturation is associated with language performance.

An in-vivo validation of ESI methods with focal sources

Electrophysiological source imaging (ESI) aims at reconstructing the precise origin of brain activity from measurements of the electric field on the scalp. Across laboratories/research centers/hospitals, ESI is performed with different methods, partly due to the ill-posedness of the underlying mathematical problem. However, it is difficult to find systematic comparisons involving a wide variety of methods. Further, existing comparisons rarely take into account the variability of the results with respect to the input parameters. Finally, comparisons are typically performed using either synthetic data, or in-vivo data where the ground-truth is only roughly known. We use an in-vivo high-density EEG dataset recorded during intracranial single pulse electrical stimulation, in which the true sources are substantially dipolar and their locations are precisely known. We compare ten different ESI methods, using their implementation in the MNE-Python package: MNE, dSPM, LORETA, sLORETA, eLORETA, LCMV beamformers, irMxNE, Gamma Map, SESAME and dipole fitting. We perform comparisons under multiple choices of input parameters, to assess the accuracy of the best reconstruction, as well as the impact of such parameters on the localization performance. Best reconstructions often fall within 1 cm from the true source, with most accurate methods hitting an average localization error of 1.2 cm and outperforming least accurate ones erring by 2.5 cm. As expected, dipolar and sparsity-promoting methods tend to outperform distributed methods. For several distributed methods, the best regularization parameter turned out to be the one in principle associated with low SNR, despite the high SNR of the available dataset. Depth weighting played no role for two out of the six methods implementing it. Sensitivity to input parameters varied widely between methods. While one would expect high variability being associated with low localization error at the best solution, this is not always the case, with some methods producing highly variable results and high localization error, and other methods producing stable results with low localization error. In particular, recent dipolar and sparsity-promoting methods provide significantly better results than older distributed methods. As we repeated the tests with "conventional" (32 channels) and dense (64, 128, 256 channels) EEG recordings, we observed little impact of the number of channels on localization accuracy; however, for distributed methods denser montages provide smaller spatial dispersion. Overall findings confirm that EEG is a reliable technique for localization of point sources and therefore reinforce the importance that ESI may have in the clinical context, especially when applied to identify the surgical target in potential candidates for epilepsy surgery.

MEG language mapping using a novel automatic ECD algorithm in comparison with MNE, dSPM, and DICS beamformer

Introduction

The single equivalent current dipole (sECD) is the standard clinical procedure for presurgical language mapping in epilepsy using magnetoencephalography (MEG). However, the sECD approach has not been widely used in clinical assessments, mainly because it requires subjective judgements in selecting several critical parameters. To address this limitation, we developed an automatic sECD algorithm (AsECDa) for language mapping.

Methods

The localization accuracy of the AsECDa was evaluated using synthetic MEG data. Subsequently, the reliability and efficiency of AsECDa were compared to three other common source localization methods using MEG data recorded during two sessions of a receptive language task in 21 epilepsy patients. These methods include minimum norm estimation (MNE), dynamic statistical parametric mapping (dSPM), and dynamic imaging of coherent sources (DICS) beamformer.

Results

For the synthetic single dipole MEG data with a typical signal-to-noise ratio, the average localization error of AsECDa was less than 2 mm for simulated superficial and deep dipoles. For the patient data, AsECDa showed better test-retest reliability (TRR) of the language laterality index (LI) than MNE, dSPM, and DICS beamformer. Specifically, the LI calculated with AsECDa revealed excellent TRR between the two MEG sessions across all patients (Cor = 0.80), while the LI for MNE, dSPM, DICS-event-related desynchronization (ERD) in the alpha band, and DICS-ERD in the low beta band ranged lower (Cor = 0.71, 0.64, 0.54, and 0.48, respectively). Furthermore, AsECDa identified 38% of patients with atypical language lateralization (i.e., right lateralization or bilateral), compared to 73%, 68%, 55%, and 50% identified by DICS-ERD in the low beta band, DICS-ERD in the alpha band, MNE, and dSPM, respectively. Compared to other methods, AsECDa's results were more consistent with previous studies that reported atypical language lateralization in 20-30% of epilepsy patients.

Discussion

Our study suggests that AsECDa is a promising approach for presurgical language mapping, and its fully automated nature makes it easy to implement and reliable for clinical evaluations.

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.

Listening to two speakers: Capacity and tradeoffs in neural speech tracking during Selective and Distributed Attention

Speech comprehension is severely compromised when several people talk at once, due to limited perceptual and cognitive resources. In such circumstances, top-down attention mechanisms can actively prioritize processing of task-relevant speech. However, behavioral and neural evidence suggest that this selection is not exclusive, and the system may have sufficient capacity to process additional speech input as well. Here we used a data-driven approach to contrast two opposing hypotheses regarding the system's capacity to co-represent competing speech: Can the brain represent two speakers equally or is the system fundamentally limited, resulting in tradeoffs between them? Neural activity was measured using magnetoencephalography (MEG) as human participants heard concurrent speech narratives and engaged in two tasks: Selective Attention, where only one speaker was task-relevant and Distributed Attention, where both speakers were equally relevant. Analysis of neural speech-tracking revealed that both tasks engaged a similar network of brain regions involved in auditory processing, attentional control and speech processing. Interestingly, during both Selective and Distributed Attention the neural representation of competing speech showed a bias towards one speaker. This is in line with proposed 'bottlenecks' for co-representation of concurrent speech and suggests that good performance on distributed attention tasks may be achieved by toggling attention between speakers over time.

Effect of channel density, inverse solutions and connectivity measures on EEG resting-state networks reconstruction: A simulation study

Along with the study of brain activity evoked by external stimuli, the past two decades witnessed an increased interest in characterizing the spontaneous brain activity occurring during resting conditions. The identification of connectivity patterns in this so-called "resting-state" has been the subject of a great number of electrophysiology-based studies, using the Electro/Magneto-Encephalography (EEG/MEG) source connectivity method. However, no consensus has been reached yet regarding a unified (if possible) analysis pipeline, and several involved parameters and methods require cautious tuning. This is particularly challenging when different analytical choices induce significant discrepancies in results and drawn conclusions, thereby hindering the reproducibility of neuroimaging research. Hence, our objective in this study was to shed light on the effect of analytical variability on outcome consistency by evaluating the implications of parameters involved in the EEG source connectivity analysis on the accuracy of resting-state networks (RSNs) reconstruction. We simulated, using neural mass models, EEG data corresponding to two RSNs, namely the default mode network (DMN) and dorsal attentional network (DAN). We investigated the impact of five channel densities (19, 32, 64, 128, 256), three inverse solutions (weighted minimum norm estimate (wMNE), exact low-resolution brain electromagnetic tomography (eLORETA), and linearly constrained minimum variance (LCMV) beamforming) and four functional connectivity measures (phase-locking value (PLV), phase-lag index (PLI), and amplitude envelope correlation (AEC) with and without source leakage correction), on the correspondence between reconstructed and reference networks. We showed that, with different analytical choices related to the number of electrodes, source reconstruction algorithm, and functional connectivity measure, high variability is present in the results. More specifically, our results show that a higher number of EEG channels significantly increased the accuracy of the reconstructed networks. Additionally, our results showed significant variability in the performance of the tested inverse solutions and connectivity measures. Such methodological variability and absence of analysis standardization represent a critical issue for neuroimaging studies that should be prioritized. We believe that this work could be useful for the field of electrophysiology connectomics, by increasing awareness regarding the challenge of variability in methodological approaches and its implications on reported results.

Neuromagnetic 40 Hz Auditory Steady-State Response in the left auditory cortex is related to language comprehension in children with Autism Spectrum Disorder

Language impairment is comorbid in most children with Autism Spectrum Disorder (ASD), but its neural mechanisms are still poorly understood. Some studies hypothesize that the atypical low-level sensory perception in the auditory cortex accounts for the abnormal language development in these children. One of the potential non-invasive measures of such low-level perception can be the cortical gamma-band oscillations registered with magnetoencephalography (MEG), and 40 Hz Auditory Steady-State Response (40 Hz ASSR) is a reliable paradigm for eliciting auditory gamma response. Although there is research in children with and without ASD using 40 Hz ASSR, nothing is known about the relationship between this auditory response in children with ASD and their language abilities measured directly in formal assessment. In the present study, we used MEG and individual brain models to investigate 40 Hz ASSR in primary-school-aged children with and without ASD. It was also used to assess how the strength of the auditory response is related to language abilities of children with ASD, their non-verbal IQ, and social functioning. A total of 40 children were included in the study. The results demonstrated that 40 Hz ASSR was reduced in the right auditory cortex in children with ASD when comparing them to typically developing controls. Importantly, our study provides the first evidence of the association between 40 Hz ASSR in the language-dominant left auditory cortex and language comprehension in children with ASD. This link was domain-specific because the other brain-behavior correlations were non-significant.

Reliable evaluation of functional connectivity and graph theory measures in source-level EEG: How many electrodes are enough?

OBJECTIVE

Using EEG to characterise functional brain networks through graph theory has gained significant interest in clinical and basic research. However, the minimal requirements for reliable measures remain largely unaddressed. Here, we examined functional connectivity estimates and graph theory metrics obtained from EEG with varying electrode densities.

METHODS

EEG was recorded with 128 electrodes in 33 participants. The high-density EEG data were subsequently subsampled into three sparser montages (64, 32, and 19 electrodes). Four inverse solutions, four measures of functional connectivity, and five graph theory metrics were tested.

RESULTS

The correlation between the results obtained with 128-electrode and the subsampled montages decreased as a function of the number of electrodes. As a result of decreased electrode density, the network metrics became skewed: mean network strength and clustering coefficient were overestimated, while characteristic path length was underestimated.

CONCLUSIONS

Several graph theory metrics were altered when electrode density was reduced. Our results suggest that, for optimal balance between resource demand and result precision, a minimum of 64 electrodes should be utilised when graph theory metrics are used to characterise functional brain networks in source-reconstructed EEG data.

SIGNIFICANCE

Characterisation of functional brain networks derived from low-density EEG warrants careful consideration.

Central-positive complexes in ECT-induced seizures: Possible evidence for thalamocortical mechanisms

OBJECTIVE

Central-positive complexes (CPCs) are elicited during electroconvulsive therapy (ECT) as generalized high-amplitude waveforms with maximum positive voltage over the vertex. While these complexes have been qualitatively assessed in previous literature, quantitative analyses are lacking. This study aims to characterize CPCs across temporal, spatial, and spectral domains.

METHODS

High-density 64-electrode electroencephalogram (EEG) recordings during 50 seizures acquired from 11 patients undergoing right unilateral ECT allowed for evaluation of spatiotemporal characteristics of CPCs via source localization and spectral analysis.

RESULTS

Peak-amplitude CPC scalp topology was consistent across seizures, showing maximal positive polarity over the midline fronto-central region and maximal negative polarity over the suborbital regions. The sources of these peak potentials were localized to the bilateral medial thalamus and cingulate cortical regions. Delta, beta, and gamma oscillations were correlated with the peak amplitude of CPCs during seizures induced during ketamine, whereas delta and gamma oscillations were associated with CPC peaks during etomidate anesthesia (excluding the dose-charge titration).

CONCLUSIONS

Our findings demonstrate the consistency of CPC presence across participant, stimulus charge, time, and anesthetic agent, with peaks localized to bilateral medial thalamus and cingulate cortical regions and associated with delta, beta, and gamma band oscillations (depending on the anesthetic condition).

SIGNIFICANCE

The consistency and reproducibility of CPCs offers ECT as a new avenue for studying the dynamics of generalized seizure activity and thalamocortical networks.

Atypical cortical processing of bottom-up speech binding cues in children with autism spectrum disorders

Individuals with autism spectrum disorder (ASD) commonly display speech processing abnormalities. Binding of acoustic features of speech distributed across different frequencies into coherent speech objects is fundamental in speech perception. Here, we tested the hypothesis that the cortical processing of bottom-up acoustic cues for speech binding may be anomalous in ASD. We recorded magnetoencephalography while ASD children (ages 7-17) and typically developing peers heard sentences of sine-wave speech (SWS) and modulated SWS (MSS) where binding cues were restored through increased temporal coherence of the acoustic components and the introduction of harmonicity. The ASD group showed increased long-range feedforward functional connectivity from left auditory to parietal cortex with concurrent decreased local functional connectivity within the parietal region during MSS relative to SWS. As the parietal region has been implicated in auditory object binding, our findings support our hypothesis of atypical bottom-up speech binding in ASD. Furthermore, the long-range functional connectivity correlated with behaviorally measured auditory processing abnormalities, confirming the relevance of these atypical cortical signatures to the ASD phenotype. Lastly, the group difference in the local functional connectivity was driven by the youngest participants, suggesting that impaired speech binding in ASD might be ameliorated upon entering adolescence.