How to assess the accuracy of volume conduction models? A validation study with stereotactic EEG data.

Introduction: Volume conduction models of the human head are used in various neuroscience fields, such as for source reconstruction in EEG and MEG, and for modeling the effects of brain stimulation. Numerous studies have quantified the accuracy and sensitivity of volume conduction models by analyzing the effects of the geometrical and electrical features of the head model, the sensor model, the source model, and the numerical method. Most studies are based on simulations as it is hard to obtain sufficiently detailed measurements to compare to models. The recording of stereotactic EEG during electric stimulation mapping provides an opportunity for such empirical validation. Methods: In the study presented here, we used the potential distribution of volume-conducted artifacts that are due to cortical stimulation to evaluate the accuracy of finite element method (FEM) volume conduction models. We adopted a widely used strategy for numerical comparison, i.e., we fixed the geometrical description of the head model and the mathematical method to perform simulations, and we gradually altered the head models, by increasing the level of detail of the conductivity profile. We compared the simulated potentials at different levels of refinement with the measured potentials in three epilepsy patients. Results: Our results show that increasing the level of detail of the volume conduction head model only marginally improves the accuracy of the simulated potentials when compared to in-vivo sEEG measurements. The mismatch between measured and simulated potentials is, throughout all patients and models, maximally 40 microvolts (i.e., 10% relative error) in 80% of the stimulation-recording combination pairs and it is modulated by the distance between recording and stimulating electrodes. Discussion: Our study suggests that commonly used strategies used to validate volume conduction models based solely on simulations might give an overly optimistic idea about volume conduction model accuracy. We recommend more empirical validations to be performed to identify those factors in volume conduction models that have the highest impact on the accuracy of simulated potentials. We share the dataset to allow researchers to further investigate the mismatch between measurements and FEM models and to contribute to improving volume conduction models.

The Swedish National Facility for Magnetoencephalography Parkinson's disease dataset.

Parkinson's disease (PD) is characterised by a loss of dopamine and dopaminergic cells. The consequences hereof are widespread network disturbances in brain function. It is an ongoing topic of investigation how the disease-related changes in brain function manifest in PD relate to clinical symptoms. We present The Swedish National Facility for Magnetoencephalography Parkinson's Disease Dataset (NatMEG-PD) as an Open Science contribution to identify the functional neural signatures of Parkinson's disease and contribute to diagnosis and treatment. The dataset contains whole-head magnetoencephalographic (MEG) recordings from 66 well-characterised PD patients on their regular dose of dopamine replacement therapy and 68 age- and sex-matched healthy controls. NatMEG-PD contains three-minute eyes-closed resting-state MEG, MEG during an active movement task, and MEG during passive movements. The data includes anonymised MRI for source analysis and clinical scores. MEG data is rich in nature and can be used to explore numerous functional features. By sharing these data, we hope other researchers will contribute to advancing our understanding of the relationship between brain activity and disease state or symptoms.

Parallel or sequential? Decoding conceptual and phonological/phonetic information from MEG signals during language production.

Speaking requires the temporally coordinated planning of core linguistic information, from conceptual meaning to articulation. Recent neurophysiological results suggested that these operations involve a cascade of neural events with subsequent onset times, whilst competing evidence suggests early parallel neural activation. To test these hypotheses, we examined the sources of neuromagnetic activity recorded from 34 participants overtly naming 134 images from 4 object categories (animals, tools, foods and clothes). Within each category, word length and phonological neighbourhood density were co-varied to target phonological/phonetic processes. Multivariate pattern analyses (MVPA) searchlights in source space decoded object categories in occipitotemporal and middle temporal cortex, and phonological/phonetic variables in left inferior frontal (BA 44) and motor cortex early on. The findings suggest early activation of multiple variables due to intercorrelated properties and interactivity of processing, thus raising important questions about the representational properties of target words during the preparatory time enabling overt speaking.

Heterogeneous Cortical Effects of Spinal Cord Stimulation.

OBJECTIVES: The understanding of the cortical effects of spinal cord stimulation (SCS) remains limited. Multiple studies have investigated the effects of SCS in resting-state electroencephalography. However, owing to the large variation in reported outcomes, we aimed to describe the differential cortical responses between two types of SCS and between responders and nonresponders using magnetoencephalography (MEG). MATERIALS AND METHODS: We conducted 5-minute resting-state MEG recordings in 25 patients with chronic pain with active SCS in three sessions, each after a one-week exposure to tonic, burst, or sham SCS. We extracted six spectral features from the measured neurophysiological signals: the alpha peak frequency; alpha power ratio (power 7-9 Hz/power 9-11 Hz); and average power in the theta (4-7.5 Hz), alpha (8-12.5 Hz), beta (13-30 Hz), and low-gamma (30.5-60 Hz) frequency bands. We compared these features (using nonparametric permutation t-tests) for MEG sensor and cortical map effects across stimulation paradigms, between participants who reported low (< 5, responders) vs high (≥ 5, nonresponders) pain scores, and in three representative participants. RESULTS: We found statistically significant (p < 0.05, false discovery rate corrected) increased MEG sensor signal power below 3 Hz in response to burst SCS compared with tonic and sham SCS. We did not find statistically significant differences (all p > 0.05) between the power spectra of responders and nonresponders. Our data did not show statistically significant differences in the spectral features of interest among the three stimulation paradigms or between responders and nonresponders. These results were confirmed by the MEG cortical maps. However, we did identify certain trends in the MEG source maps for all comparisons and several features, with substantial variation across participants. CONCLUSIONS: The considerable variation in cortical responses to the various SCS treatment options necessitates studies with sample sizes larger than commonly reported in the field and more personalized treatment plans. Studies with a finer stratification between responders and nonresponders are required to advance the knowledge on SCS treatment effects.

fNIRS is sensitive to leg activity in the primary motor cortex after systemic artifact correction.

BACKGROUND: functional near-infrared spectroscopy (fNIRS) is an increasingly popular tool to study cortical activity during movement and gait that requires further validation. This study aimed to assess (1) whether fNIRS can detect the difficult-to-measure leg area of the primary motor cortex (M1) and distinguish it from the hand area; and (2) whether fNIRS can differentiate between automatic (i.e., not requiring one's attention) and non-automatic movement processes. Special attention was attributed to systemic artifacts (i.e., changes in blood pressure, heart rate, breathing) which were assessed and corrected by short channels, i.e., fNIRS channels which are mainly sensitive to superficial scalp hemodynamics. METHODS: Twenty-three seated, healthy participants tapped four fingers on a keyboard or tapped the right foot on four squares on the floor in a specific order given by a 12-digit sequence (e.g., 434141243212). Two different sequences were executed: a beforehand learned (i.e., automatic) version and a newly learned (i.e., non-automatic) version. A 36-channel fNIRS device including 12 short channels covered multiple motor-related cortical areas including M1. The fNIRS data were analyzed with a general linear model (GLM). Correlation between the expected functional hemodynamic responses (i.e. task regressor) and the short channels (i.e. nuisance regressors), necessitated performing a separate short channel regression instead of integrating them in the GLM. RESULTS: Consistent with the M1 somatotopy, we found significant HbO increases of very large effect size in the lateral M1 channels during finger tapping (Cohen's d = 1.35, p<0.001) and significant HbO increases of moderate effect size in the medial M1 channels during foot tapping (Cohen's d = 0.8, p<0.05). The cortical activity differences between automatic and non-automatic tasks were not significantly different. Importantly, leg movements produced large systemic fluctuations, which were adequately removed by the use of all available short channels. DISCUSSION: Our results indicate that fNIRS is sensitive to leg activity in M1, though the sensitivity is lower than for finger activity and requires rigorous correction for systemic fluctuations. We furthermore highlight that systemic artifacts may result in an unreliable GLM analysis when short channels show signals that are similar to the expected hemodynamic responses.

Introduction to the shared near infrared spectroscopy format.

Significance: Functional near-infrared spectroscopy (fNIRS) is a popular neuroimaging technique with proliferating hardware platforms, analysis approaches, and software tools. There has not been a standardized file format for storing fNIRS data, which has hindered the sharing of data as well as the adoption and development of software tools. Aim: We endeavored to design a file format to facilitate the analysis and sharing of fNIRS data that is flexible enough to meet the community's needs and sufficiently defined to be implemented consistently across various hardware and software platforms. Approach: The shared NIRS format (SNIRF) specification was developed in consultation with the academic and commercial fNIRS community and the Society for functional Near Infrared Spectroscopy. Results: The SNIRF specification defines a format for fNIRS data acquired using continuous wave, frequency domain, time domain, and diffuse correlation spectroscopy devices. Conclusions: We present the SNIRF along with validation software and example datasets. Support for reading and writing SNIRF data has been implemented by major hardware and software platforms, and the format has found widespread use in the fNIRS community.

Estimating the influence of stroke lesions on MEG source reconstruction.

Source reconstruction of magnetoencephalography (MEG) has been used to assess brain reorganization after brain damage, such as stroke. Lesions result in parts of the brain having an electrical conductivity that differs from the normal values. The effect this has on the forward solutions (i.e., the propagation of electric currents and magnetic fields generated by cortical activity) is well predictable. However, their influence on source localization results is not well characterized and understood. This is specifically a concern for patient studies with asymmetric (i.e., within one hemisphere) lesions focusing on asymmetric and lateralized brain activity, such as language. In particular, it is good practice to consider the level of geometrical detail that is necessary to compute and interpret reliable source reconstruction results. To understand the effect of lesions on source estimates and propose recommendations to researchers working with clinical data, in this study we consider the trade off between improved accuracy and the additional effort to compute more realistic head models, with the aim to answer the question whether the additional effort is worth it. We simulated and analyzed the effects of a stroke lesion (i.e., an asymmetrically distributed CSF-filled cavity) in the head model with three different sizes and locations when performing MEG source reconstruction using a finite element method (FEM). We compared the effect of the lesion with a homogeneous head model that neglects the lesion. We computed displacement and attenuation/amplification maps to quantify the localization errors and signal magnitude modulation. We conclude that brain lesions leading to asymmetrically distributed CSF-filled cavities should be modeled when performing MEG source reconstruction, especially when investigating deep sources or post-stroke hemispheric lateralization of functions. The strongest effects are not only visible in perilesional areas, but can extend up to 20 mm from the lesion. Bigger lesions lead to stronger effects impacting larger areas, independently from the lesion location. Lastly, we conclude that more priority should be given to usability and accessibility of the required computational tools, to allow researchers with less technical expertise to use the improved methods that are available but currently not widely adopted yet.

Advances in human intracranial electroencephalography research, guidelines and good practices.

Since the second-half of the twentieth century, intracranial electroencephalography (iEEG), including both electrocorticography (ECoG) and stereo-electroencephalography (sEEG), has provided an intimate view into the human brain. At the interface between fundamental research and the clinic, iEEG provides both high temporal resolution and high spatial specificity but comes with constraints, such as the individual's tailored sparsity of electrode sampling. Over the years, researchers in neuroscience developed their practices to make the most of the iEEG approach. Here we offer a critical review of iEEG research practices in a didactic framework for newcomers, as well addressing issues encountered by proficient researchers. The scope is threefold: (i) review common practices in iEEG research, (ii) suggest potential guidelines for working with iEEG data and answer frequently asked questions based on the most widespread practices, and (iii) based on current neurophysiological knowledge and methodologies, pave the way to good practice standards in iEEG research. The organization of this paper follows the steps of iEEG data processing. The first section contextualizes iEEG data collection. The second section focuses on localization of intracranial electrodes. The third section highlights the main pre-processing steps. The fourth section presents iEEG signal analysis methods. The fifth section discusses statistical approaches. The sixth section draws some unique perspectives on iEEG research. Finally, to ensure a consistent nomenclature throughout the manuscript and to align with other guidelines, e.g., Brain Imaging Data Structure (BIDS) and the OHBM Committee on Best Practices in Data Analysis and Sharing (COBIDAS), we provide a glossary to disambiguate terms related to iEEG research.

The Time Course of Language Production as Revealed by Pattern Classification of MEG Sensor Data.

Language production involves a complex set of computations, from conceptualization to articulation, which are thought to engage cascading neural events in the language network. However, recent neuromagnetic evidence suggests simultaneous meaning-to-speech mapping in picture naming tasks, as indexed by early parallel activation of frontotemporal regions to lexical semantic, phonological, and articulatory information. Here we investigate the time course of word production, asking to what extent such "earliness" is a distinctive property of the associated spatiotemporal dynamics. Using MEG, we recorded the neural signals of 34 human subjects (26 males) overtly naming 134 images from four semantic object categories (animals, foods, tools, clothes). Within each category, we covaried word length, as quantified by the number of syllables contained in a word, and phonological neighborhood density to target lexical and post-lexical phonological/phonetic processes. Multivariate pattern analyses searchlights in sensor space distinguished the stimulus-locked spatiotemporal responses to object categories early on, from 150 to 250 ms after picture onset, whereas word length was decoded in left frontotemporal sensors at 250-350 ms, followed by the latency of phonological neighborhood density (350-450 ms). Our results suggest a progression of neural activity from posterior to anterior language regions for the semantic and phonological/phonetic computations preparing overt speech, thus supporting serial cascading models of word production.SIGNIFICANCE STATEMENT Current psycholinguistic models make divergent predictions on how a preverbal message is mapped onto articulatory output during the language planning. Serial models predict a cascading sequence of hierarchically organized neural computations from conceptualization to articulation. In contrast, parallel models posit early simultaneous activation of multiple conceptual, phonological, and articulatory information in the language system. Here we asked whether such earliness is a distinctive property of the neural dynamics of word production. The combination of the millisecond precision of MEG with multivariate pattern analyses revealed subsequent onset times for the neural events supporting semantic and phonological/phonetic operations, progressing from posterior occipitotemporal to frontal sensor areas. The findings bring new insights for refining current theories of language production.

Microscopy-BIDS: An Extension to the Brain Imaging Data Structure for Microscopy Data.

The Brain Imaging Data Structure (BIDS) is a specification for organizing, sharing, and archiving neuroimaging data and metadata in a reusable way. First developed for magnetic resonance imaging (MRI) datasets, the community-led specification evolved rapidly to include other modalities such as magnetoencephalography, positron emission tomography, and quantitative MRI (qMRI). In this work, we present an extension to BIDS for microscopy imaging data, along with example datasets. Microscopy-BIDS supports common imaging methods, including 2D/3D, ex/in vivo, micro-CT, and optical and electron microscopy. Microscopy-BIDS also includes comprehensible metadata definitions for hardware, image acquisition, and sample properties. This extension will facilitate future harmonization efforts in the context of multi-modal, multi-scale imaging such as the characterization of tissue microstructure with qMRI.

Sharing individualised template MRI data for MEG source reconstruction: A solution for open data while keeping subject confidentiality.

The increasing requirements for adoption of FAIR data management and sharing original research data from neuroimaging studies can be at odds with protecting the anonymity of the research participants due to the person-identifiable anatomical features in the data. We propose a solution to this dilemma for anatomical MRIs used in MEG source analysis. In MEG analysis, the channel-level data is reconstructed to the source-level using models derived from anatomical MRIs. Sharing data, therefore, requires sharing the anatomical MRI to replicate the analysis. The suggested solution is to replace the individual anatomical MRIs with individualised warped templates that can be used to carry out the MEG source analysis and that provide sufficient geometrical similarity to the original participants' MRIs. First, we demonstrate how the individualised template warping can be implemented with one of the leading open-source neuroimaging analysis toolboxes. Second, we compare results from four different MEG source reconstruction methods performed with an individualised warped template to those using the participant's original MRI. While the source reconstruction results are not numerically identical, there is a high similarity between the results for single dipole fits, dynamic imaging of coherent sources beamforming, and atlas-based virtual channel beamforming. There is a moderate similarity between minimum-norm estimates, as anticipated due to this method being anatomically constrained and dependent on the exact morphological features of the cortical sheet. We also compared the morphological features of the warped template to those of the original MRI. These showed a high similarity in grey matter volume and surface area, but a low similarity in the average cortical thickness and the mean folding index within cortical parcels. Taken together, this demonstrates that the results obtained by MEG source reconstruction can be preserved with the warped templates, whereas the anatomical and morphological fingerprint is sufficiently altered to protect the anonymity of research participants. In cases where participants consent to sharing anatomical MRI data, it remains preferable to share the original defaced data with an appropriate data use agreement. In cases where participants did not consent to share their MRIs, the individualised warped MRI template offers a good compromise in sharing data for reuse while retaining anonymity for research participants.

Deep-Breathing Biofeedback Trainability in a Virtual-Reality Action Game: A Single-Case Design Study With Police Trainers.

It is widely recognized that police performance may be hindered by psychophysiological state changes during acute stress. To address the need for awareness and control of these physiological changes, police academies in many countries have implemented Heart-Rate Variability (HRV) biofeedback training. Despite these trainings now being widely delivered in classroom setups, they typically lack the arousing action context needed for successful transfer to the operational field, where officers must apply learned skills, particularly when stress levels rise. The study presented here aimed to address this gap by training physiological control skills in an arousing decision-making context. We developed a Virtual-Reality (VR) breathing-based biofeedback training in which police officers perform deep and slow diaphragmatic breathing in an engaging game-like action context. This VR game consisted of a selective shoot/don't shoot game designed to assess response inhibition, an impaired capacity in high arousal situations. Biofeedback was provided based on adherence to a slow breathing pace: the slower and deeper the breathing, the less constrained peripheral vision became, facilitating accurate responses to the in-game demands. A total of nine male police trainers completed 10 sessions over a 4-week period as part of a single-case experimental ABAB study-design (i.e., alternating sessions with and without biofeedback). Results showed that eight out of nine participants showed improved breathing control in action, with a positive effect on breathing-induced low frequency HRV, while also improving their in-game behavioral performance. Critically, the breathing-based skill learning transferred to subsequent sessions in which biofeedback was not presented. Importantly, all participants remained highly engaged throughout the training. Altogether, our study showed that our VR environment can be used to train breathing regulation in an arousing and active decision-making context.

Brain structural and functional correlates to defense-related inhibition of muscle sympathetic nerve activity in man.

An individual's blood pressure (BP) reactivity to stress is linked to increased risk of hypertension and cardiovascular disease. However, inter- and intra-individual BP variability makes understanding the coupling between stress, BP reactivity, and long-term outcomes challenging. Previous microneurographic studies of sympathetic signaling to muscle vasculature (i.e. muscle sympathetic nerve activity, MSNA) have established a neural predictor for an individual's BP reactivity during short-lasting stress. Unfortunately, this method is invasive, technically demanding, and time-consuming and thus not optimal for widespread use. Potential central nervous system correlates have not been investigated. We used MagnetoEncephaloGraphy and Magnetic Resonance Imaging to search for neural correlates to sympathetic response profiles within the central autonomic network and sensorimotor (Rolandic) regions in 20 healthy young males. The main correlates include (a) Rolandic beta rebound and an anterior cingulate cortex (ACC) response elicited by sudden stimulation and (b) cortical thickness in the ACC. Our findings highlight the involvement of the ACC in reactions to stress entailing peripheral sympathetic responses to environmental stimuli. The Rolandic response furthermore indicates a surprisingly strong link between somatosensory and autonomic processes. Our results thus demonstrate the potential in using non-invasive neuroimaging-based measures of stress-related MSNA reactions, previously assessed only using invasive microneurography.

Are alpha and beta oscillations spatially dissociated over the cortex in context-driven spoken-word production?

Decreases in oscillatory alpha- and beta-band power have been consistently found in spoken-word production. These have been linked to both motor preparation and conceptual-lexical retrieval processes. However, the observed power decreases have a broad frequency range that spans two "classic" (sensorimotor) bands: alpha and beta. It remains unclear whether alpha- and beta-band power decreases contribute independently when a spoken word is planned. Using a re-analysis of existing magnetoencephalography data, we probed whether the effects in alpha and beta bands are spatially distinct. Participants read a sentence that was either constraining or non-constraining toward the final word, which was presented as a picture. In separate blocks participants had to name the picture or score its predictability via button press. Irregular-resampling auto-spectral analysis (IRASA) was used to isolate the oscillatory activity in the alpha and beta bands from the background 1-over-f spectrum. The sources of alpha- and beta-band oscillations were localized based on the participants' individualized peak frequencies. For both tasks, alpha- and beta-power decreases overlapped in left posterior temporal and inferior parietal cortex, regions that have previously been associated with conceptual and lexical processes. The spatial distributions of the alpha and beta power effects were spatially similar in these regions to the extent we could assess it. By contrast, for left frontal regions, the spatial distributions differed between alpha and beta effects. Our results suggest that for conceptual-lexical retrieval, alpha and beta oscillations do not dissociate spatially and, thus, are distinct from the classical sensorimotor alpha and beta oscillations.

A unified view on beamformers for M/EEG source reconstruction.

Beamforming is a popular method for functional source reconstruction using magnetoencephalography (MEG) and electroencephalography (EEG) data. Beamformers, which were first proposed for MEG more than two decades ago, have since been applied in hundreds of studies, demonstrating that they are a versatile and robust tool for neuroscience. However, certain characteristics of beamformers remain somewhat elusive and there currently does not exist a unified documentation of the mathematical underpinnings and computational subtleties of beamformers as implemented in the most widely used academic open source software packages for MEG analysis (Brainstorm, FieldTrip, MNE, and SPM). Here, we provide such documentation that aims at providing the mathematical background of beamforming and unifying the terminology. Beamformer implementations are compared across toolboxes and pitfalls of beamforming analyses are discussed. Specifically, we provide details on handling rank deficient covariance matrices, prewhitening, the rank reduction of forward fields, and on the combination of heterogeneous sensor types, such as magnetometers and gradiometers. The overall aim of this paper is to contribute to contemporary efforts towards higher levels of computational transparency in functional neuroimaging.

Enhancing reproducibility in developmental EEG research: BIDS, cluster-based permutation tests, and effect sizes.

Developmental research using electroencephalography (EEG) offers valuable insights in brain processes early in life, but at the same time, applying this sensitive technique to young children who are often non-compliant and have short attention spans comes with practical limitations. It is thus of particular importance to optimally use the limited resources to advance our understanding of development through reproducible and replicable research practices. Here, we describe methodological approaches that help maximize the reproducibility of developmental EEG research. We discuss how to transform EEG data into the standardized Brain Imaging Data Structure (BIDS) which organizes data according to the FAIR data sharing principles. We provide a tutorial on how to use cluster-based permutation testing to analyze developmental EEG data. This versatile test statistic solves the multiple comparison problem omnipresent in EEG analysis and thereby substantially decreases the risk of reporting false discoveries. Finally, we describe how to quantify effect sizes, in particular of cluster-based permutation results. Reporting effect sizes conveys a finding's impact and robustness which in turn informs future research. To demonstrate these methodological approaches to data organization, analysis and report, we use a publicly accessible infant EEG dataset and provide a complete copy of the analysis code.

The Human Connectome Project: A retrospective.

The Human Connectome Project (HCP) was launched in 2010 as an ambitious effort to accelerate advances in human neuroimaging, particularly for measures of brain connectivity; apply these advances to study a large number of healthy young adults; and freely share the data and tools with the scientific community. NIH awarded grants to two consortia; this retrospective focuses on the "WU-Minn-Ox" HCP consortium centered at Washington University, the University of Minnesota, and University of Oxford. In just over 6 years, the WU-Minn-Ox consortium succeeded in its core objectives by: 1) improving MR scanner hardware, pulse sequence design, and image reconstruction methods, 2) acquiring and analyzing multimodal MRI and MEG data of unprecedented quality together with behavioral measures from more than 1100 HCP participants, and 3) freely sharing the data (via the ConnectomeDB database) and associated analysis and visualization tools. To date, more than 27 Petabytes of data have been shared, and 1538 papers acknowledging HCP data use have been published. The "HCP-style" neuroimaging paradigm has emerged as a set of best-practice strategies for optimizing data acquisition and analysis. This article reviews the history of the HCP, including comments on key events and decisions associated with major project components. We discuss several scientific advances using HCP data, including improved cortical parcellations, analyses of connectivity based on functional and diffusion MRI, and analyses of brain-behavior relationships. We also touch upon our efforts to develop and share a variety of associated data processing and analysis tools along with detailed documentation, tutorials, and an educational course to train the next generation of neuroimagers. We conclude with a look forward at opportunities and challenges facing the human neuroimaging field from the perspective of the HCP consortium.

Magnetoencephalography reveals increased slow-to-fast alpha power ratios in patients with chronic pain.

INTRODUCTION: Objective disease markers are a key for diagnosis and personalized interventions. In chronic pain, such markers are still not available, and therapy relies on individual patients' reports. However, several pain studies have reported group-based differences in functional magnetic resonance imaging, electroencephalography, and magnetoencephalography (MEG). OBJECTIVES: We aimed to explore spectral differences in resting-state MEG brain signals between patients with chronic pain and pain-free controls and to characterize the cortical and subcortical regions involved. METHODS: We estimated power spectral density over 5 minutes of resting-state MEG recordings in patients with chronic pain and controls and derived 7 spectral features at the sensor and source levels: alpha peak frequency, alpha power ratio (power 7-9 Hz divided by power 9-11 Hz), and average power in theta, alpha, beta, low-gamma, and high-gamma bands. We performed nonparametric permutation t tests (false discovery rate corrected) to assess between-group differences in these 7 spectral features. RESULTS: Twenty-one patients with chronic pain and 25 controls were included. No significant group differences were found in alpha peak frequency or average power in any frequency band. The alpha power ratio was significantly higher (P < 0.05) in patients with chronic pain at both the sensor and brain source levels. The brain regions showing significantly higher ratios included the occipital, parietal, temporal and frontal lobe areas, insular and cingulate cortex, and right thalamus. CONCLUSION: The alpha power ratio is a simple, promising signal marker of chronic pain, affecting an expansive range of cortical and subcortical regions, including known pain-processing areas.