Spatiotemporal whole-brain activity and functional connectivity of melodies recognition.

Music is a non-verbal human language, built on logical, hierarchical structures, that offers excellent opportunities to explore how the brain processes complex spatiotemporal auditory sequences. Using the high temporal resolution of magnetoencephalography, we investigated the unfolding brain dynamics of 70 participants during the recognition of previously memorized musical sequences compared to novel sequences matched in terms of entropy and information content. Measures of both whole-brain activity and functional connectivity revealed a widespread brain network underlying the recognition of the memorized auditory sequences, which comprised primary auditory cortex, superior temporal gyrus, insula, frontal operculum, cingulate gyrus, orbitofrontal cortex, basal ganglia, thalamus, and hippocampus. Furthermore, while the auditory cortex responded mainly to the first tones of the sequences, the activity of higher-order brain areas such as the cingulate gyrus, frontal operculum, hippocampus, and orbitofrontal cortex largely increased over time during the recognition of the memorized versus novel musical sequences. In conclusion, using a wide range of analytical techniques spanning from decoding to functional connectivity and building on previous works, our study provided new insights into the spatiotemporal whole-brain mechanisms for conscious recognition of auditory sequences.

CutFEM-based MEG forward modeling improves source separability and sensitivity to quasi-radial sources: A somatosensory group study.

Source analysis of magnetoencephalography (MEG) data requires the computation of the magnetic fields induced by current sources in the brain. This so-called MEG forward problem includes an accurate estimation of the volume conduction effects in the human head. Here, we introduce the Cut finite element method (CutFEM) for the MEG forward problem. CutFEM's meshing process imposes fewer restrictions on tissue anatomy than tetrahedral meshes while being able to mesh curved geometries contrary to hexahedral meshing. To evaluate the new approach, we compare CutFEM with a boundary element method (BEM) that distinguishes three tissue compartments and a 6-compartment hexahedral FEM in an n = 19 group study of somatosensory evoked fields (SEF). The neural generators of the 20 ms post-stimulus SEF components (M20) are reconstructed using both an unregularized and a regularized inversion approach. Changing the forward model resulted in reconstruction differences of about 1 centimeter in location and considerable differences in orientation. The tested 6-compartment FEM approaches significantly increase the goodness of fit to the measured data compared with the 3-compartment BEM. They also demonstrate higher quasi-radial contributions for sources below the gyral crowns. Furthermore, CutFEM improves source separability compared with both other approaches. We conclude that head models with 6 compartments rather than 3 and the new CutFEM approach are valuable additions to MEG source reconstruction, in particular for sources that are predominantly radial.

A dynamic generative model can extract interpretable oscillatory components from multichannel neurophysiological recordings.

Modern neurophysiological recordings are performed using multichannel sensor arrays that are able to record activity in an increasingly high number of channels numbering in the 100s to 1000s. Often, underlying lower-dimensional patterns of activity are responsible for the observed dynamics, but these representations are difficult to reliably identify using existing methods that attempt to summarize multivariate relationships in a post hoc manner from univariate analyses or using current blind source separation methods. While such methods can reveal appealing patterns of activity, determining the number of components to include, assessing their statistical significance, and interpreting them requires extensive manual intervention and subjective judgment in practice. These difficulties with component selection and interpretation occur in large part because these methods lack a generative model for the underlying spatio-temporal dynamics. Here, we describe a novel component analysis method anchored by a generative model where each source is described by a bio-physically inspired state-space representation. The parameters governing this representation readily capture the oscillatory temporal dynamics of the components, so we refer to it as oscillation component analysis. These parameters - the oscillatory properties, the component mixing weights at the sensors, and the number of oscillations - all are inferred in a data-driven fashion within a Bayesian framework employing an instance of the expectation maximization algorithm. We analyze high-dimensional electroencephalography and magnetoencephalography recordings from human studies to illustrate the potential utility of this method for neuroscience data.

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.

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.

Homogeneous B0 coil design method for open-access ultra-low field magnetic resonance imaging: A simulation study.

A multimodal brain function measurement system integrating functional magnetic resonance imaging (fMRI) and magnetoencephalography (MEG) is expected to be a tool that will provide new insights into neuroscience. To integrate fMRI and MEG, an ultra-low-field MRI (ULF-MRI) scanner that can generate a static magnetic field (B0) with an electromagnetic coil and turn off the B0 during MEG measurements is desirable. While electromagnetic B0 coil has the above advantages, it also has a trade-off between size and the broadness of the magnetic field homogeneity. In this study, we proposed a method for designing a B0 multi-stage circular coil arrangement that determines the number of coils required to maximize magnetic field homogeneity and minimize the total wiring length of the coils. The optimized multi-stage coil arrangement had an external shape of 600 mm in diameter and a maximum height of 600 mm, with an aperture of 600 mm in diameter and 300 mm in height. The magnetic field homogeneity was <100 ppm over a 210 mm diameter spherical volume (DSV). Compared to a previous two coil pairs arrangement with the same magnetic field homogeneity, the diameter was 1/1.9 times smaller, indicating that the newly designed B0 coil arrangement realized a smaller size and wider magnetic field homogeneity.

Reliability of dynamic causal modelling of resting-state magnetoencephalography.

This study assesses the reliability of resting-state dynamic causal modelling (DCM) of magnetoencephalography (MEG) under conductance-based canonical microcircuit models, in terms of both posterior parameter estimates and model evidence. We use resting-state MEG data from two sessions, acquired 2 weeks apart, from a cohort with high between-subject variance arising from Alzheimer's disease. Our focus is not on the effect of disease, but on the reliability of the methods (as within-subject between-session agreement), which is crucial for future studies of disease progression and drug intervention. To assess the reliability of first-level DCMs, we compare model evidence associated with the covariance among subject-specific free energies (i.e., the 'quality' of the models) with versus without interclass correlations. We then used parametric empirical Bayes (PEB) to investigate the differences between the inferred DCM parameter probability distributions at the between subject level. Specifically, we examined the evidence for or against parameter differences (i) within-subject, within-session, and between-epochs; (ii) within-subject between-session; and (iii) within-site between-subjects, accommodating the conditional dependency among parameter estimates. We show that for data acquired close in time, and under similar circumstances, more than 95% of inferred DCM parameters are unlikely to differ, speaking to mutual predictability over sessions. Using PEB, we show a reciprocal relationship between a conventional definition of 'reliability' and the conditional dependency among inferred model parameters. Our analyses confirm the reliability and reproducibility of the conductance-based DCMs for resting-state neurophysiological data. In this respect, the implicit generative modelling is suitable for interventional and longitudinal studies of neurological and psychiatric disorders.

Multi-Modal Electrophysiological Source Imaging With Attention Neural Networks Based on Deep Fusion of EEG and MEG.

The process of reconstructing underlying cortical and subcortical electrical activities from Electroencephalography (EEG) or Magnetoencephalography (MEG) recordings is called Electrophysiological Source Imaging (ESI). Given the complementarity between EEG and MEG in measuring radial and tangential cortical sources, combined EEG/MEG is considered beneficial in improving the reconstruction performance of ESI algorithms. Traditional algorithms mainly emphasize incorporating predesigned neurophysiological priors to solve the ESI problem. Deep learning frameworks aim to directly learn the mapping from scalp EEG/MEG measurements to the underlying brain source activities in a data-driven manner, demonstrating superior performance compared to traditional methods. However, most of the existing deep learning approaches for the ESI problem are performed on a single modality of EEG or MEG, meaning the complementarity of these two modalities has not been fully utilized. How to fuse the EEG and MEG in a more principled manner under the deep learning paradigm remains a challenging question. This study develops a Multi-Modal Deep Fusion (MMDF) framework using Attention Neural Networks (ANN) to fully leverage the complementary information between EEG and MEG for solving the ESI inverse problem, which is termed as MMDF-ANN. Specifically, our proposed brain source imaging approach consists of four phases, including feature extraction, weight generation, deep feature fusion, and source mapping. Our experimental results on both synthetic dataset and real dataset demonstrated that using a fusion of EEG and MEG can significantly improve the source localization accuracy compared to using a single-modality of EEG or MEG. Compared to the benchmark algorithms, MMDF-ANN demonstrated good stability when reconstructing sources with extended activation areas and situations of EEG/MEG measurements with a low signal-to-noise ratio.

Macroscopic brain dynamics beyond contralateral primary motor cortex for movement prediction.

This study investigates the complex relationship between upper limb movement direction and macroscopic neural signals in the brain, which is critical for understanding brain-computer interfaces (BCI). Conventional BCI research has primarily focused on a local area, such as the contralateral primary motor cortex (M1), relying on the population-based decoding method with microelectrode arrays. In contrast, macroscopic approaches such as electroencephalography (EEG) and magnetoencephalography (MEG) utilize numerous electrodes to cover broader brain regions. This study probes the potential differences in the mechanisms of microscopic and macroscopic methods. It is important to determine which neural activities effectively predict movements. To investigate this, we analyzed MEG data from nine right-handed participants while performing arm-reaching tasks. We employed dynamic statistical parametric mapping (dSPM) to estimate source activity and built a decoding model composed of long short-term memory (LSTM) and a multilayer perceptron to predict movement trajectories. This model achieved a high correlation coefficient of 0.79 between actual and predicted trajectories. Subsequently, we identified brain regions sensitive to predicting movement direction using the integrated gradients (IG) method, which assesses the predictive contribution of each source activity. The resulting salience map demonstrated a distribution without significant differences across motor-related regions, including M1. Predictions based solely on M1 activity yielded a correlation coefficient of 0.42, nearly half as effective as predictions incorporating all source activities. This suggests that upper limb movements are influenced by various factors such as movement coordination, planning, body and target position recognition, and control, beyond simple muscle activity. All of the activities are needed in the decoding model using macroscopic signals. Our findings also revealed that contralateral and ipsilateral hemispheres contribute equally to movement prediction, implying that BCIs could potentially benefit patients with brain damage in the contralateral hemisphere by utilizing brain signals from the ipsilateral hemisphere. In conclusion, this study demonstrates that macroscopic activity from large brain regions significantly contributes to predicting upper limb movement. Non-invasive BCI systems would require a comprehensive collection of neural signals from multiple brain regions.

Multitask learning of a biophysically-detailed neuron model.

The human brain operates at multiple levels, from molecules to circuits, and understanding these complex processes requires integrated research efforts. Simulating biophysically-detailed neuron models is a computationally expensive but effective method for studying local neural circuits. Recent innovations have shown that artificial neural networks (ANNs) can accurately predict the behavior of these detailed models in terms of spikes, electrical potentials, and optical readouts. While these methods have the potential to accelerate large network simulations by several orders of magnitude compared to conventional differential equation based modelling, they currently only predict voltage outputs for the soma or a select few neuron compartments. Our novel approach, based on enhanced state-of-the-art architectures for multitask learning (MTL), allows for the simultaneous prediction of membrane potentials in each compartment of a neuron model, at a speed of up to two orders of magnitude faster than classical simulation methods. By predicting all membrane potentials together, our approach not only allows for comparison of model output with a wider range of experimental recordings (patch-electrode, voltage-sensitive dye imaging), it also provides the first stepping stone towards predicting local field potentials (LFPs), electroencephalogram (EEG) signals, and magnetoencephalography (MEG) signals from ANN-based simulations. While LFP and EEG are an important downstream application, the main focus of this paper lies in predicting dendritic voltages within each compartment to capture the entire electrophysiology of a biophysically-detailed neuron model. It further presents a challenging benchmark for MTL architectures due to the large amount of data involved, the presence of correlations between neighbouring compartments, and the non-Gaussian distribution of membrane potentials.

Beyond Stereo-EEG: Is It Worth Combining Stereo-EEG With Other Diagnostic Methods?

SUMMARY: Stereo-EEG is a widely used method to improve the diagnostic precision of presurgical workup in patients with refractory epilepsy. Its ability to detect epileptic activity and identify epileptic networks largely depends on the chosen implantation strategy. Even in an ideal situation, electrodes record activity generated in <10% of the brain and contacts only record from brain tissue in their immediate proximity. In this article, the authors discuss how recording stereo-EEG simultaneously with other diagnostic methods can improve its diagnostic value in clinical and research settings. It can help overcome the limited spatial coverage of intracranial recording and better understand the sources of epileptic activity. Simultaneous scalp EEG is the most widely available method, often used to understand large epileptic networks, seizure propagation, and EEG activity occurring on the contralateral hemisphere. Simultaneous magnetoencephalography allows for more precise source localization and identification of deep sources outside the stereo-EEG coverage. Finally, simultaneous functional MRI can highlight metabolic changes following epileptic activity and help understand the widespread network changes associated with interictal activity. This overview highlights advantages and methodological challenges for all these methods. Clinical use and research applications are presented for each approach.

Multi-scale entropy assessment of magnetoencephalography signals in schizophrenia.

Schizophrenia is a severe disruption in cognition and emotion, affecting fundamental human functions. In this study, we applied Multi-Scale Entropy analysis to resting-state Magnetoencephalography data from 54 schizophrenia patients and 98 healthy controls. This method quantifies the temporal complexity of the signal across different time scales using the concept of sample entropy. Results show significantly higher sample entropy in schizophrenia patients, primarily in central, parietal, and occipital lobes, peaking at time scales equivalent to frequencies between 15 and 24 Hz. To disentangle the contributions of the amplitude and phase components, we applied the same analysis to a phase-shuffled surrogate signal. The analysis revealed that most differences originate from the amplitude component in the δ, α, and ß power bands. While the phase component had a smaller magnitude, closer examination reveals clear spatial patterns and significant differences across specific brain regions. We assessed the potential of multi-scale entropy as a schizophrenia biomarker by comparing its classification performance to conventional spectral analysis and a cognitive task (the n-back paradigm). The discriminative power of multi-scale entropy and spectral features was similar, with a slight advantage for multi-scale entropy features. The results of the n-back test were slightly below those obtained from multi-scale entropy and spectral features.

Contingent magnetic variation and beta-band oscillations in sensorimotor temporal decision-making.

The ability to accurately encode the temporal information of sensory events and hence to make prompt action is fundamental to humans' prompt behavioral decision-making. Here we examined the ability of ensemble coding (averaging multiple inter-intervals in a sound sequence) and subsequent immediate reproduction of target duration at half, equal, or double that of the perceived mean interval in a sensorimotor loop. With magnetoencephalography (MEG), we found that the contingent magnetic variation (CMV) in the central scalp varied as a function of the averaging tasks, with a faster rate for buildup amplitudes and shorter peak latencies in the "half" condition as compared to the "double" condition. ERD (event-related desynchronization) -to-ERS (event-related synchronization) latency was shorter in the "half" condition. A robust beta band (15-23 Hz) power suppression and recovery between the final tone and the action of key pressing was found for time reproduction. The beta modulation depth (i.e., the ERD-to-ERS power difference) was larger in motor areas than in primary auditory areas. Moreover, results of phase slope index (PSI) indicated that beta oscillations in the left supplementary motor area (SMA) led those in the right superior temporal gyrus (STG), showing SMA to STG directionality for the processing of sequential (temporal) auditory interval information. Our findings provide the first evidence to show that CMV and beta oscillations predict the coupling between perception and action in time averaging.

The neurophysiological brain-fingerprint of Parkinson's disease.

BACKGROUND: Research in healthy young adults shows that characteristic patterns of brain activity define individual "brain-fingerprints" that are unique to each person. However, variability in these brain-fingerprints increases in individuals with neurological conditions, challenging the clinical relevance and potential impact of the approach. Our study shows that brain-fingerprints derived from neurophysiological brain activity are associated with pathophysiological and clinical traits of individual patients with Parkinson's disease (PD). METHODS: We created brain-fingerprints from task-free brain activity recorded through magnetoencephalography in 79 PD patients and compared them with those from two independent samples of age-matched healthy controls (N = 424 total). We decomposed brain activity into arrhythmic and rhythmic components, defining distinct brain-fingerprints for each type from recording durations of up to 4 min and as short as 30 s. FINDINGS: The arrhythmic spectral components of cortical activity in patients with Parkinson's disease are more variable over short periods, challenging the definition of a reliable brain-fingerprint. However, by isolating the rhythmic components of cortical activity, we derived brain-fingerprints that distinguished between patients and healthy controls with about 90% accuracy. The most prominent cortical features of the resulting Parkinson's brain-fingerprint are mapped to polyrhythmic activity in unimodal sensorimotor regions. Leveraging these features, we also demonstrate that Parkinson's symptom laterality can be decoded directly from cortical neurophysiological activity. Furthermore, our study reveals that the cortical topography of the Parkinson's brain-fingerprint aligns with that of neurotransmitter systems affected by the disease's pathophysiology. INTERPRETATION: The increased moment-to-moment variability of arrhythmic brain-fingerprints challenges patient differentiation and explains previously published results. We outline patient-specific rhythmic brain signaling features that provide insights into both the neurophysiological signature and symptom laterality of Parkinson's disease. Thus, the proposed definition of a rhythmic brain-fingerprint of Parkinson's disease may contribute to novel, refined approaches to patient stratification. Symmetrically, we discuss how rhythmic brain-fingerprints may contribute to the improved identification and testing of therapeutic neurostimulation targets. FUNDING: Data collection and sharing for this project was provided by the Quebec Parkinson Network (QPN), the Pre-symptomatic Evaluation of Novel or Experimental Treatments for Alzheimer's Disease (PREVENT-AD; release 6.0) program, the Cambridge Centre for Aging Neuroscience (Cam-CAN), and the Open MEG Archives (OMEGA). The QPN is funded by a grant from Fonds de Recherche du Québec - Santé (FRQS). PREVENT-AD was launched in 2011 as a $13.5 million, 7-year public-private partnership using funds provided by McGill University, the FRQS, an unrestricted research grant from Pfizer Canada, the Levesque Foundation, the Douglas Hospital Research Centre and Foundation, the Government of Canada, and the Canada Fund for Innovation. The Brainstorm project is supported by funding to SB from the NIH (R01-EB026299-05). Further funding to SB for this study included a Discovery grant from the Natural Sciences and Engineering Research Council of Canada of Canada (436355-13), and the CIHR Canada research Chair in Neural Dynamics of Brain Systems (CRC-2017-00311).

Source localization comparison and combination of OPM-MEG and fMRI to detect sensorimotor cortex responses.

BACKGROUND AND OBJECTIVES: The exploration of various neuroimaging techniques have become focal points within the field of neuroscience research. Magnetoencephalography based on optically pumped magnetometers (OPM-MEG) has shown significant potential to be the next generation of functional neuroimaging with the advantages of high signal intensity and flexible sensor arrangement. In this study, we constructed a 31-channel OPM-MEG system and performed a preliminary comparison of the temporal and spatial relationship between magnetic responses measured by OPM-MEG and blood-oxygen-level-dependent signals detected by functional magnetic resonance imaging (fMRI) during a grasping task. METHODS: For OPM-MEG, the ß-band (15-30 Hz) oscillatory activities can be reliably detected across multiple subjects and multiple session runs. To effectively localize the inhibitory oscillatory activities, a source power-spectrum ratio-based imaging method was proposed. This approach was compared with conventional source imaging methods, such as minimum norm-type and beamformer methods, and was applied in OPM-MEG source analysis. Subsequently, the spatial and temporal responses at the source-level between OPM-MEG and fMRI were analyzed. RESULTS: The effectiveness of the proposed method was confirmed through simulations compared to benchmark methods. Our demonstration revealed an average spatial separation of 10.57 ± 4.41 mm between the localization results of OPM-MEG and fMRI across four subjects. Furthermore, the fMRI-constrained OPM-MEG localization results indicated a more focused imaging extent. CONCLUSIONS: Taken together, the performance exhibited by OPM-MEG positions it as a potential instrument for functional surgery assessment.

Potential association between suicide risk, aggression, impulsivity, and the somatosensory system.

Aggression and impulsivity are linked to suicidal behaviors, but their relationship to the suicidal crisis remains unclear. This magnetoencephalography (MEG) study investigated the link between aggression, impulsivity, and resting-state MEG power and connectivity. Four risk groups were enrolled: high-risk (HR; n = 14), who had a recent suicidal crisis; lower-risk (LR; n = 41), who had a history of suicide attempts but no suicide attempt or ideation in the past year; clinical control (CC; n = 38), who had anxiety/mood disorders but no suicidal history; and minimal risk (MR; n = 28), who had no psychiatric/suicidal history. No difference in resting-state MEG power was observed between the groups. Individuals in the HR group with high self-reported aggression and impulsivity scores had reduced MEG power in regions responsible for sensory/emotion regulation vs. those in the HR group with low scores. The HR group also showed downregulated bidirectional glutamatergic feedback between the precuneus (PRE) and insula (INS) compared to the LR, CC, and MR groups. High self-reported impulsivity was linked to reduced PRE to INS feedback, whereas high risk-taking impulsivity was linked to upregulated INS to postcentral gyrus (PCG) and PCG to INS feedback. These preliminary findings suggest that glutamatergic-mediated sensory and emotion-regulation processes may function as potential suicide risk markers.

Expanding the clinical application of OPM-MEG using an effective automatic suppression method for the dental brace metal artifact.

Optically pumped magnetometer magnetoencephalography (OPM-MEG) holds significant promise for clinical functional brain imaging due to its superior spatiotemporal resolution. However, effectively suppressing metallic artifacts, particularly from devices such as orthodontic braces and vagal nerve stimulators remains a major challenge, hindering the wider clinical application of wearable OPM-MEG devices. A comprehensive analysis of metal artifact characteristics from time, frequency, and time-frequency perspectives was conducted for the first time using an OPM-MEG device in clinical medicine. This study focused on patients with metal orthodontics, examining the modulation of metal artifacts by breath and head movement, the incomplete regular sub-Gaussian distribution, and the high absolute power ratio in the 0.5-8 Hz band. The existing metal artifact suppression algorithms applied to SQUID-MEG, such as fast independent component analysis (FastICA), information maximization (Infomax), and algorithms for multiple unknown signal extraction (AMUSE), exhibit limited efficacy. Consequently, this study introduced the second-order blind identification (SOBI) algorithm, which utilized multiple time delays for the component separation of OPM-MEG measurement signals. We modified the time delays of the SOBI method to improve its efficacy in separating artifact components, particularly those in the ultralow frequency range. This approach employs the frequency-domain absolute power ratio, root mean square (RMS) value, and mutual information methods to automate the artifact component screening process. The effectiveness of this method was validated through simulation experiments involving four subjects in both resting and evoked experiments. In addition, the proposed method was also validated by the actual OPM-MEG evoked experiments of three subjects. Comparative analyses were conducted against the FastICA, Infomax, and AMUSE algorithms. Evaluation metrics included normalized mean square error, normalized delta band power error, RMS error, and signal-to-noise ratio, demonstrating that the proposed method provides optimal suppression of metal artifacts. This advancement holds promise for enhancing data quality and expanding the clinical applications of OPM-MEG.

Tracking the neurodevelopmental trajectory of beta band oscillations with optically pumped magnetometer-based magnetoencephalography.

Neural oscillations mediate the coordination of activity within and between brain networks, supporting cognition and behaviour. How these processes develop throughout childhood is not only an important neuroscientific question but could also shed light on the mechanisms underlying neurological and psychiatric disorders. However, measuring the neurodevelopmental trajectory of oscillations has been hampered by confounds from instrumentation. In this paper, we investigate the suitability of a disruptive new imaging platform - optically pumped magnetometer-based magnetoencephalography (OPM-MEG) - to study oscillations during brain development. We show how a unique 192-channel OPM-MEG device, which is adaptable to head size and robust to participant movement, can be used to collect high-fidelity electrophysiological data in individuals aged between 2 and 34 years. Data were collected during a somatosensory task, and we measured both stimulus-induced modulation of beta oscillations in sensory cortex, and whole-brain connectivity, showing that both modulate significantly with age. Moreover, we show that pan-spectral bursts of electrophysiological activity drive task-induced beta modulation, and that their probability of occurrence and spectral content change with age. Our results offer new insights into the developmental trajectory of beta oscillations and provide clear evidence that OPM-MEG is an ideal platform for studying electrophysiology in neurodevelopment.

On the speech envelope in the cortical tracking of speech.

The synchronization between the speech envelope and neural activity in auditory regions, referred to as cortical tracking of speech (CTS), plays a key role in speech processing. The method selected for extracting the envelope is a crucial step in CTS measurement, and the absence of a consensus on best practices among the various methods can influence analysis outcomes and interpretation. Here, we systematically compare five standard envelope extraction methods the absolute value of Hilbert transform (absHilbert), gammatone filterbanks, heuristic approach, Bark scale, and vocalic energy), analyzing their impact on the CTS. We present performance metrics for each method based on the recording of brain activity from participants listening to speech in clear and noisy conditions, utilizing intracranial EEG, MEG and EEG data. As expected, we observed significant CTS in temporal brain regions below 10 Hz across all datasets, regardless of the extraction methods. In general, the gammatone filterbanks approach consistently demonstrated superior performance compared to other methods. Results from our study can guide scientists in the field to make informed decisions about the optimal analysis to extract the CTS, contributing to advancing the understanding of the neuronal mechanisms implicated in CTS.

Using a linear dynamic system to measure functional connectivity from M/EEG.

Objective.Measures of functional connectivity (FC) can elucidate which cortical regions work together in order to complete a variety of behavioral tasks. This study's primary objective was to expand a previously published model of measuring FC to include multiple subjects and several regions of interest. While FC has been more extensively investigated in vision and other sensorimotor tasks, it is not as well understood in audition. The secondary objective of this study was to investigate how auditory regions are functionally connected to other cortical regions when attention is directed to different distinct auditory stimuli.Approach.This study implements a linear dynamic system (LDS) to measure the structured time-lagged dependence across several cortical regions in order to estimate their FC during a dual-stream auditory attention task.Results.The model's output shows consistent functionally connected regions across different listening conditions, indicative of an auditory attention network that engages regardless of endogenous switching of attention or different auditory cues being attended.Significance.The LDS implemented in this study implements a multivariate autoregression to infer FC across cortical regions during an auditory attention task. This study shows how a first-order autoregressive function can reliably measure functional connectivity from M/EEG data. Additionally, the study shows how auditory regions engage with the supramodal attention network outlined in the visual attention literature.