The influence of wakefulness fluctuations on brain networks involved in centrotemporal spike occurrence.

OBJECTIVE: Drowsiness has been implicated in the modulation of centro-temporal spikes (CTS) in Self-limited epilepsy with Centro-Temporal Spikes (SeLECTS). Here, we explore this relationship and whether fluctuations in wakefulness influence the brain networks involved in CTS generation. METHODS: Functional MRI (fMRI) and electroencephalography (EEG) was simultaneously acquired in 25 SeLECTS. A multispectral EEG index quantified drowsiness ('EWI': EEG Wakefulness Index). EEG (Pearson Correlation, Cross Correlation, Trend Estimation, Granger Causality) and fMRI (PPI: psychophysiological interactions) analytic approaches were adopted to explore respectively: (a) the relationship between EWI and changes in CTS frequency and (b) the functional connectivity of the networks involved in CTS generation and wakefulness oscillations. EEG analyses were repeated on a sample of routine EEG from the same patient's cohort. RESULTS: No correlation was found between EWI fluctuations and CTS density during the EEG-fMRI recordings, while they showed an anticorrelated trend when drowsiness was followed by proper sleep in routine EEG traces. According to PPI findings, EWI fluctuations modulate the connectivity between the brain networks engaged by CTS and the left frontal operculum. CONCLUSIONS: While CTS frequency per se seems unrelated to drowsiness, wakefulness oscillations modulate the connectivity between CTS generators and key regions of the language circuitry, a cognitive function often impaired in SeLECTS. SIGNIFICANCE: This work advances our understanding of (a) interaction between CTS occurrence and vigilance fluctuations and (b) possible mechanisms responsible for language disruption in SeLECTS.

Propofol Reversibly Attenuates Short-Range Microstate Ordering and 20 Hz Microstate Oscillations.

Microstate sequences summarize the changing voltage patterns measured by electroencephalography, using a clustering approach to reduce the high dimensionality of the underlying data. A common approach is to restrict the pattern matching step to local maxima of the global field power (GFP) and to interpolate the microstate fit in between. In this study, we investigate how the anesthetic propofol affects microstate sequence periodicity and predictability, and how these metrics are changed by interpolation. We performed two frequency analyses on microstate sequences, one based on time-lagged mutual information, the other based on Fourier transform methodology, and quantified the effects of interpolation. Resting-state microstate sequences had a 20 Hz frequency peak related to dominant 10 Hz (alpha) rhythms, and the Fourier approach demonstrated that all five microstate classes followed this frequency. The 20 Hz periodicity was reversibly attenuated under moderate propofol sedation, as shown by mutual information and Fourier analysis. Characteristic microstate frequencies could only be observed in non-interpolated microstate sequences and were masked by smoothing effects of interpolation. Information-theoretic analysis revealed faster microstate dynamics and larger entropy rates under propofol, whereas Shannon entropy did not change significantly. In moderate sedation, active information storage decreased for non-interpolated sequences. Signatures of non-equilibrium dynamics were observed in non-interpolated sequences, but no changes were observed between sedation levels. All changes occurred while subjects were able to perform an auditory perception task. In summary, we show that low dose propofol reversibly increases the randomness of microstate sequences and attenuates microstate oscillations without correlation to cognitive task performance. Microstate dynamics between GFP peaks reflect physiological processes that are not accessible in interpolated sequences.

Frequency Analysis of EEG Microstate Sequences in Wakefulness and NREM Sleep.

The majority of EEG microstate analyses concern wakefulness, and the existing sleep studies have focused on changes in spatial microstate properties and on microstate transitions between adjacent time points, the shortest available time scale. We present a more extensive time series analysis of unsmoothed EEG microstate sequences in wakefulness and non-REM sleep stages across many time scales. Very short time scales are assessed with Markov tests, intermediate time scales by the entropy rate and long time scales by a spectral analysis which identifies characteristic microstate frequencies. During the descent from wakefulness to sleep stage N3, we find that the increasing mean microstate duration is a gradual phenomenon explained by a continuous slowing of microstate dynamics as described by the relaxation time of the transition probability matrix. The finite entropy rate, which considers longer microstate histories, shows that microstate sequences become more predictable (less random) with decreasing vigilance level. Accordingly, the Markov property is absent in wakefulness but in sleep stage N3, 10/19 subjects have microstate sequences compatible with a second-order Markov process. A spectral microstate analysis is performed by comparing the time-lagged mutual information coefficients of microstate sequences with the autocorrelation function of the underlying EEG. We find periodic microstate behavior in all vigilance states, linked to alpha frequencies in wakefulness, theta activity in N1, sleep spindle frequencies in N2, and in the delta frequency band in N3. In summary, we show that EEG microstates are a dynamic phenomenon with oscillatory properties that slow down in sleep and are coupled to specific EEG frequencies across several sleep stages.

Complexity Measures for EEG Microstate Sequences: Concepts and Algorithms.

EEG microstate sequence analysis quantifies properties of ongoing brain electrical activity which is known to exhibit complex dynamics across many time scales. In this report we review recent developments in quantifying microstate sequence complexity, we classify these approaches with regard to different complexity concepts, and we evaluate excess entropy as a yet unexplored quantity in microstate research. We determined the quantities entropy rate, excess entropy, Lempel-Ziv complexity (LZC), and Hurst exponents on Potts model data, a discrete statistical mechanics model with a temperature-controlled phase transition. We then applied the same techniques to EEG microstate sequences from wakefulness and non-REM sleep stages and used first-order Markov surrogate data to determine which time scales contributed to the different complexity measures. We demonstrate that entropy rate and LZC measure the Kolmogorov complexity (randomness) of microstate sequences, whereas excess entropy and Hurst exponents describe statistical complexity which attains its maximum at intermediate levels of randomness. We confirmed the equivalence of entropy rate and LZC when the LZ-76 algorithm is used, a result previously reported for neural spike train analysis (Amigó et al., Neural Comput 16:717-736, https://doi.org/10.1162/089976604322860677 , 2004). Surrogate data analyses prove that entropy-based quantities and LZC focus on short-range temporal correlations, whereas Hurst exponents include short and long time scales. Sleep data analysis reveals that deeper sleep stages are accompanied by a decrease in Kolmogorov complexity and an increase in statistical complexity. Microstate jump sequences, where duplicate states have been removed, show higher randomness, lower statistical complexity, and no long-range correlations. Regarding the practical use of these methods, we suggest that LZC can be used as an efficient entropy rate estimator that avoids the estimation of joint entropies, whereas entropy rate estimation via joint entropies has the advantage of providing excess entropy as the second parameter of the same linear fit. We conclude that metrics of statistical complexity are a useful addition to microstate analysis and address a complexity concept that is not yet covered by existing microstate algorithms while being actively explored in other areas of brain research.

Low-dimensional organization of global brain states of reduced consciousness.

Brain states are frequently represented using a unidimensional scale measuring the richness of subjective experience (level of consciousness). This description assumes a mapping between the high-dimensional space of whole-brain configurations and the trajectories of brain states associated with changes in consciousness, yet this mapping and its properties remain unclear. We combine whole-brain modeling, data augmentation, and deep learning for dimensionality reduction to determine a mapping representing states of consciousness in a low-dimensional space, where distances parallel similarities between states. An orderly trajectory from wakefulness to patients with brain injury is revealed in a latent space whose coordinates represent metrics related to functional modularity and structure-function coupling, increasing alongside loss of consciousness. Finally, we investigate the effects of model perturbations, providing geometrical interpretation for the stability and reversibility of states. We conclude that conscious awareness depends on functional patterns encoded as a low-dimensional trajectory within the vast space of brain configurations.

Deep Brain Stimulation of the Anterior Nucleus of the Thalamus in Drug-Resistant Epilepsy in the MORE Multicenter Patient Registry.

BACKGROUND AND OBJECTIVES: The efficacy of deep brain stimulation of the anterior nucleus of the thalamus (ANT DBS) in patients with drug-resistant epilepsy (DRE) was demonstrated in the double-blind Stimulation of the Anterior Nucleus of the Thalamus for Epilepsy randomized controlled trial. The Medtronic Registry for Epilepsy (MORE) aims to understand the safety and longer-term effectiveness of ANT DBS therapy in routine clinical practice. METHODS: MORE is an observational registry collecting prospective and retrospective clinical data. Participants were at least 18 years old, with focal DRE recruited across 25 centers from 13 countries. They were followed for at least 2 years in terms of seizure frequency (SF), responder rate (RR), health-related quality of life (Quality of Life in Epilepsy Inventory 31), depression, and safety outcomes. RESULTS: Of the 191 patients recruited, 170 (mean [SD] age of 35.6 [10.7] years, 43% female) were implanted with DBS therapy and met all eligibility criteria. At baseline, 38% of patients reported cognitive impairment. The median monthly SF decreased by 33.1% from 15.8 at baseline to 8.8 at 2 years (p < 0.0001) with 32.3% RR. In the subgroup of 47 patients who completed 5 years of follow-up, the median monthly SF decreased by 55.1% from 16 at baseline to 7.9 at 5 years (p < 0.0001) with 53.2% RR. High-volume centers (>10 implantations) had 42.8% reduction in median monthly SF by 2 years in comparison with 25.8% in low-volume center. In patients with cognitive impairment, the reduction in median monthly SF was 26.0% by 2 years compared with 36.1% in patients without cognitive impairment. The most frequently reported adverse events were changes (e.g., increased frequency/severity) in seizure (16%), memory impairment (patient-reported complaint, 15%), depressive mood (patient-reported complaint, 13%), and epilepsy (12%). One definite sudden unexpected death in epilepsy case was reported. DISCUSSION: The MORE registry supports the effectiveness and safety of ANT DBS therapy in a real-world setting in the 2 years following implantation. CLASSIFICATION OF EVIDENCE: This study provides Class IV evidence that ANT DBS reduces the frequency of seizures in patients with drug-resistant focal epilepsy. TRIAL REGISTRATION INFORMATION: MORE ClinicalTrials.gov Identifier: NCT01521754, first posted on January 31, 2012.

Data-driven discovery of canonical large-scale brain dynamics.

Human behavior and cognitive function correlate with complex patterns of spatio-temporal brain dynamics, which can be simulated using computational models with different degrees of biophysical realism. We used a data-driven optimization algorithm to determine and classify the types of local dynamics that enable the reproduction of different observables derived from functional magnetic resonance recordings. The phase space analysis of the resulting equations revealed a predominance of stable spiral attractors, which optimized the similarity to the empirical data in terms of the synchronization, metastability, and functional connectivity dynamics. For stable limit cycles, departures from harmonic oscillations improved the fit in terms of functional connectivity dynamics. Eigenvalue analyses showed that proximity to a bifurcation improved the accuracy of the simulation for wakefulness, whereas deep sleep was associated with increased stability. Our results provide testable predictions that constrain the landscape of suitable biophysical models, while supporting noise-driven dynamics close to a bifurcation as a canonical mechanism underlying the complex fluctuations that characterize endogenous brain activity.

fMRI spectral signatures of sleep.

Sleep can be distinguished from wake by changes in brain electrical activity, typically assessed using electroencephalography (EEG). The hallmark of nonrapid-eye-movement (NREM) sleep is the shift from high-frequency, low-amplitude wake EEG to low-frequency, high-amplitude sleep EEG dominated by spindles and slow waves. Here we identified signatures of sleep in brain hemodynamic activity, using simultaneous functional MRI (fMRI) and EEG. We found that, at the transition from wake to sleep, fMRI blood oxygen level-dependent (BOLD) activity evolved from a mixed-frequency pattern to one dominated by two distinct oscillations: a low-frequency (<0.1 Hz) oscillation prominent in light sleep and correlated with the occurrence of spindles, and a high-frequency oscillation (>0.1 Hz) prominent in deep sleep and correlated with the occurrence of slow waves. The two oscillations were both detectable across the brain but exhibited distinct spatiotemporal patterns. During the falling-asleep process, the low-frequency oscillation first appeared in the thalamus, then the posterior cortex, and lastly the frontal cortex, while the high-frequency oscillation first appeared in the midbrain, then the frontal cortex, and lastly the posterior cortex. During the waking-up process, both oscillations disappeared first from the thalamus, then the frontal cortex, and lastly the posterior cortex. The BOLD oscillations provide local signatures of spindle and slow wave activity. They may be employed to monitor the regional occurrence of sleep or wakefulness, track which regions are the first to fall asleep or wake up at the wake-sleep transitions, and investigate local homeostatic sleep processes.

Unifying turbulent dynamics framework distinguishes different brain states.

Significant advances have been made by identifying the levels of synchrony of the underlying dynamics of a given brain state. This research has demonstrated that non-conscious dynamics tend to be more synchronous than in conscious states, which are more asynchronous. Here we go beyond this dichotomy to demonstrate that different brain states are underpinned by dissociable spatiotemporal dynamics. We investigated human neuroimaging data from different brain states (resting state, meditation, deep sleep and disorders of consciousness after coma). The model-free approach was based on Kuramoto's turbulence framework using coupled oscillators. This was extended by a measure of the information cascade across spatial scales. Complementarily, the model-based approach used exhaustive in silico perturbations of whole-brain models fitted to these measures. This allowed studying of the information encoding capabilities in given brain states. Overall, this framework demonstrates that elements from turbulence theory provide excellent tools for describing and differentiating between brain states.

New-onset refractory status epilepticus (NORSE) and febrile infection-related epilepsy syndrome (FIRES) of unknown aetiology: A comparison of the incomparable?

PURPOSE: The discussion is ongoing whether new-onset refractory status epilepticus (NORSE) in adults and febrile infection-related epilepsy syndrome (FIRES) in children are one syndrome if the aetiology is unknown. In this study we will compare an adult cohort with NORSE and a paediatric cohort with FIRES in order to determine if they are similar or different. METHODS: We retrospectively compared 18 adults with NORSE versus 48 children with FIRES, both cohorts without identifiable cause despite extensive investigations. We analyzed demographic and clinical data using Mann-Whitney-U and χ2- tests. RESULTS: NORSE affected more women (78% vs. 42%; P = 0.009) than in FIRES. Median acute hospital stay was longer in FIRES (35 days [interquartile range, IQR=36] vs. 20 days [IQR=19]; P<0.001). FIRES was treated more frequently with coma therapy (82% vs. 28%; P<0.001) and with a higher median number of antiseizure medicines (7 [IQR=5] vs. 4 [IQR=2]; P<0.001). Children with FIRES showed a higher cerebrospinal fluid (CSF) cell count (10 cells/µl; P = 0.002) but a lower CSF protein level than adults with NORSE (48 mg/dl; P = 0.028). Immunotherapy was administered more frequently in FIRES (73% vs. 22%; P<0.001) than in NORSE. Group differences in number of antiseizure medicines after hospital stay (P = 0.229) and in overall mortality (P = 0.327) were not significant. CONCLUSION: In our explorative comparison, differences prevailed. NORSE and FIRES should be compared prospectively in age-matched cohorts.

Unconsciousness reconfigures modular brain network dynamics.

The dynamic core hypothesis posits that consciousness is correlated with simultaneously integrated and differentiated assemblies of transiently synchronized brain regions. We represented time-dependent functional interactions using dynamic brain networks and assessed the integrity of the dynamic core by means of the size and flexibility of the largest multilayer module. As a first step, we constrained parameter selection using a newly developed benchmark for module detection in heterogeneous temporal networks. Next, we applied a multilayer modularity maximization algorithm to dynamic brain networks computed from functional magnetic resonance imaging (fMRI) data acquired during deep sleep and under propofol anesthesia. We found that unconsciousness reconfigured network flexibility and reduced the size of the largest spatiotemporal module, which we identified with the dynamic core. Our results represent a first characterization of modular brain network dynamics during states of unconsciousness measured with fMRI, adding support to the dynamic core hypothesis of human consciousness.

Nonequilibrium brain dynamics as a signature of consciousness.

The cognitive functions of human and nonhuman primates rely on the dynamic interplay of distributed neural assemblies. As such, it seems unlikely that cognition can be supported by macroscopic brain dynamics at the proximity of equilibrium. We confirmed this hypothesis by investigating electrocorticography data from nonhuman primates undergoing different states of unconsciousness (sleep, and anesthesia with propofol, ketamine, and ketamine plus medetomidine), and functional magnetic resonance imaging data from humans, both during deep sleep and under propofol anesthesia. Systematically, all states of reduced consciousness unfolded at higher proximity to equilibrium compared to conscious wakefulness, as demonstrated by the computation of entropy production and the curl of probability flux in phase space. Our results establish nonequilibrium macroscopic brain dynamics as a robust signature of consciousness, opening the way for the characterization of cognition and awareness using tools from statistical mechanics.

Decoding brain states on the intrinsic manifold of human brain dynamics across wakefulness and sleep.

Current state-of-the-art functional magnetic resonance imaging (fMRI) offers remarkable imaging quality and resolution, yet, the intrinsic dimensionality of brain dynamics in different states (wakefulness, light and deep sleep) remains unknown. Here we present a method to reveal the low dimensional intrinsic manifold underlying human brain dynamics, which is invariant of the high dimensional spatio-temporal representation of the neuroimaging technology. By applying this intrinsic manifold framework to fMRI data acquired in wakefulness and sleep, we reveal the nonlinear differences between wakefulness and three different sleep stages, and successfully decode these different brain states with a mean accuracy across participants of 96%. Remarkably, a further group analysis shows that the intrinsic manifolds of all participants share a common topology. Overall, our results reveal the intrinsic manifold underlying the spatiotemporal dynamics of brain activity and demonstrate how this manifold enables the decoding of different brain states such as wakefulness and various sleep stages.

Perturbations in dynamical models of whole-brain activity dissociate between the level and stability of consciousness.

Consciousness transiently fades away during deep sleep, more stably under anesthesia, and sometimes permanently due to brain injury. The development of an index to quantify the level of consciousness across these different states is regarded as a key problem both in basic and clinical neuroscience. We argue that this problem is ill-defined since such an index would not exhaust all the relevant information about a given state of consciousness. While the level of consciousness can be taken to describe the actual brain state, a complete characterization should also include its potential behavior against external perturbations. We developed and analyzed whole-brain computational models to show that the stability of conscious states provides information complementary to their similarity to conscious wakefulness. Our work leads to a novel methodological framework to sort out different brain states by their stability and reversibility, and illustrates its usefulness to dissociate between physiological (sleep), pathological (brain-injured patients), and pharmacologically-induced (anesthesia) loss of consciousness.

Noise-driven multistability vs deterministic chaos in phenomenological semi-empirical models of whole-brain activity.

An outstanding open problem in neuroscience is to understand how neural systems are capable of producing and sustaining complex spatiotemporal dynamics. Computational models that combine local dynamics with in vivo measurements of anatomical and functional connectivity can be used to test potential mechanisms underlying this complexity. We compared two conceptually different mechanisms: noise-driven switching between equilibrium solutions (modeled by coupled Stuart-Landau oscillators) and deterministic chaos (modeled by coupled Rossler oscillators). We found that both models struggled to simultaneously reproduce multiple observables computed from the empirical data. This issue was especially manifested in the case of noise-driven dynamics close to a bifurcation, which imposed overly strong constraints on the optimal model parameters. In contrast, the chaotic model could produce complex behavior over a range of parameters, thus being capable of capturing multiple observables at the same time with good performance. Our observations support the view of the brain as a non-equilibrium system able to produce endogenous variability. We presented a simple model capable of jointly reproducing functional connectivity computed at different temporal scales. Besides adding to our conceptual understanding of brain complexity, our results inform and constrain the future development of biophysically realistic large-scale models.

Signature of consciousness in brain-wide synchronization patterns of monkey and human fMRI signals.

During the sleep-wake cycle, the brain undergoes profound dynamical changes, which manifest subjectively as transitions between conscious experience and unconsciousness. Yet, neurophysiological signatures that can objectively distinguish different consciousness states based are scarce. Here, we show that differences in the level of brain-wide signals can reliably distinguish different stages of sleep and anesthesia from the awake state in human and monkey fMRI resting state data. Moreover, a whole-brain computational model can faithfully reproduce changes in global synchronization and other metrics such as functional connectivity, structure-function relationship, integration and segregation across vigilance states. We demonstrate that the awake brain is close to a Hopf bifurcation, which naturally coincides with the emergence of globally correlated fMRI signals. Furthermore, simulating lesions of individual brain areas highlights the importance of connectivity hubs in the posterior brain and subcortical nuclei for maintaining the model in the awake state, as predicted by graph-theoretical analyses of structural data.

Generative Embeddings of Brain Collective Dynamics Using Variational Autoencoders.

We consider the problem of encoding pairwise correlations between coupled dynamical systems in a low-dimensional latent space based on few distinct observations. We use variational autoencoders (VAEs) to embed temporal correlations between coupled nonlinear oscillators that model brain states in the wake-sleep cycle into a two-dimensional manifold. Training a VAE with samples generated using two different parameter combinations results in an embedding that encodes the repertoire of collective dynamics, as well as the topology of the underlying connectivity network. We first follow this approach to infer the trajectory of brain states measured from wakefulness to deep sleep from the two end points of this trajectory; then, we show that the same architecture was capable of representing the pairwise correlations of generic Landau-Stuart oscillators coupled by complex network topology.

Connectivity dynamics from wakefulness to sleep.

Interest in time-resolved connectivity in fMRI has grown rapidly in recent years. The most widely used technique for studying connectivity changes over time utilizes a sliding windows approach. There has been some debate about the utility of shorter versus longer windows, the use of fixed versus adaptive windows, as well as whether observed resting state dynamics during wakefulness may be predominantly due to changes in sleep state and subject head motion. In this work we use an independent component analysis (ICA)-based pipeline applied to concurrent EEG/fMRI data collected during wakefulness and various sleep stages and show: 1) connectivity states obtained from clustering sliding windowed correlations of resting state functional network time courses well classify the sleep states obtained from EEG data, 2) using shorter sliding windows instead of longer non-overlapping windows improves the ability to capture transition dynamics even at windows as short as 30 â€‹s, 3) motion appears to be mostly associated with one of the states rather than spread across all of them 4) a fixed tapered sliding window approach outperforms an adaptive dynamic conditional correlation approach, and 5) consistent with prior EEG/fMRI work, we identify evidence of multiple states within the wakeful condition which are able to be classified with high accuracy. Classification of wakeful only states suggest the presence of time-varying changes in connectivity in fMRI data beyond sleep state or motion. Results also inform about advantageous technical choices, and the identification of different clusters within wakefulness that are separable suggest further studies in this direction.

Modeling regional changes in dynamic stability during sleep and wakefulness.

Global brain states are frequently placed within a unidimensional continuum by correlational studies, ranging from states of deep unconsciousness to ordinary wakefulness. An alternative is their multidimensional and mechanistic characterization in terms of different cognitive capacities, using computational models to reproduce the underlying neural dynamics. We explore this alternative by introducing a semi-empirical model linking regional activation and long-range functional connectivity in the different brain states visited during the natural wake-sleep cycle. Our model combines functional magnetic resonance imaging (fMRI) data, in vivo estimates of structural connectivity, and anatomically-informed priors to constrain the independent variation of regional activation. The best fit to empirical data was achieved using priors based on functionally coherent networks, with the resulting model parameters dividing the cortex into regions presenting opposite dynamical behavior. Frontoparietal regions approached a bifurcation from dynamics at a fixed point governed by noise, while sensorimotor regions approached a bifurcation from oscillatory dynamics. In agreement with human electrophysiological experiments, sleep onset induced subcortical deactivation with low correlation, which was subsequently reversed for deeper stages. Finally, we introduced periodic forcing of variable intensity to simulate external perturbations, and identified the key regions relevant for the recovery of wakefulness from deep sleep. Our model represents sleep as a state with diminished perceptual gating and the latent capacity for global accessibility that is required for rapid arousals. To the extent that the qualitative characterization of local dynamics is exhausted by the dichotomy between unstable and stable behavior, our work highlights how expanding the model parameter space can describe states of consciousness in terms of multiple dimensions with interpretations given by the choice of anatomically-informed priors.

Magnetic Measurement of Electrically Evoked Muscle Responses With Optically Pumped Magnetometers.

OBJECTIVE: Electroneurography has been an essential method for assessing peripheral nerve disorders for decades. During this procedure, a nerve is briefly electrically excited, and nerve conduction properties are identified by indirect means from the behavior of the innervated muscle. The magnetic field of the resulting muscle response can also be recorded by novel, uncooled magnetometers, which have become very attractive for different medical applications over recent years. These highly sensitive magnetometers are called optically pumped magnetometers. METHODS: We performed unaveraged and averaged magnetic signal detection of electrically evoked muscle responses using optically pumped magnetometers. We then discussed the suitability of this procedure for clinical applications in the context of diagnostic value and in direct comparison with the current electrical gold standard. RESULTS: The magnetic detection of muscle responses is possible using optically pumped magnetometers. Our magnetic results (averaged and unaveraged) closely match those from electrical measurements. CONCLUSION: Optically pumped magnetometers provide an alternative, contactless technology for electrode-based motor studies, but they are currently not ready for routine clinical use. This costly technology requires additional earth magnetic shielding because this is a prerequisite for proper operation. Currently, there are no diagnostic advantages over electrical measurements. Additionally, the required measurement setup and procedure are much more complicated. SIGNIFICANCE: In contrast to already published proof-of-principle studies for magnetomyography, we report in detail the results of the magnetic measurements of electrically evoked muscle responses in a shielded environment by applying supramaximal stimulation and finally validate our findings with electroneurography data as a reference.