Assessing the balance between excitation and inhibition in chronic pain through the aperiodic component of EEG
eLife Assessment
Gil Ávila et al. evaluated the aperiodic component in the medial prefrontal cortex using resting-state EEG recordings from 149 individuals with chronic pain and 115 healthy participants. The authors present compelling evidence that the aperiodic component of the EEG does not differentiate between those with chronic pain and healthy individuals. The study was well-designed and rigorously conducted, and the clear and conclusive results provide important insights that can guide future research in the field of pain neuroscience.
https://doi.org/10.7554/eLife.101727.3.sa0Important: Findings that have theoretical or practical implications beyond a single subfield
- Landmark
- Fundamental
- Important
- Valuable
- Useful
Compelling: Evidence that features methods, data and analyses more rigorous than the current state-of-the-art
- Exceptional
- Compelling
- Convincing
- Solid
- Incomplete
- Inadequate
During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments
Abstract
Chronic pain is a prevalent and debilitating condition whose neural mechanisms are incompletely understood. An imbalance of cerebral excitation and inhibition (E/I), particularly in the medial prefrontal cortex (mPFC), is believed to represent a crucial mechanism in the development and maintenance of chronic pain. Thus, identifying a non-invasive, scalable marker of E/I could provide valuable insights into the neural mechanisms of chronic pain and aid in developing clinically useful biomarkers. Recently, the aperiodic component of the electroencephalography (EEG) power spectrum has been proposed to represent a non-invasive proxy for E/I. We, therefore, assessed the aperiodic component in the mPFC of resting-state EEG recordings in 149 people with chronic pain and 115 healthy participants. We found robust evidence against differences in the aperiodic component in the mPFC between people with chronic pain and healthy participants, and no correlation between the aperiodic component and pain intensity. These findings were consistent across different subtypes of chronic pain and were similarly found in a whole-brain analysis. Their robustness was supported by preregistration and multiverse analyses across many different methodological choices. Together, our results suggest that the EEG aperiodic component does not differentiate between people with chronic pain and healthy individuals. These findings and the rigorous methodological approach can guide future studies investigating non-invasive, scalable markers of cerebral dysfunction in people with chronic pain and beyond.
Introduction
Chronic pain is a highly prevalent disease (Kennedy et al., 2014) that severely decreases the quality of life of those who live with it, and imposes a substantial burden on society and healthcare systems (Rice et al., 2016). The neural mechanisms of chronic pain are incompletely understood, and effective therapies are lacking (Cohen et al., 2021). Therefore, advancing our understanding, assessment, and treatment of chronic pain is urgently needed.
Human and animal studies have revealed structural and functional changes in an extended network of brain areas in chronic pain, including the medial prefrontal cortex (mPFC; Baliki and Apkarian, 2015; Kuner and Flor, 2017). Cellular and microcircuit studies in rodent models of chronic pain have indicated imbalances between neural excitation and inhibition in this network (Bliss et al., 2016; Kummer et al., 2020; Shiers and Price, 2020; Tan and Kuner, 2021). Correspondingly, positron emission tomography (PET) and magnetic resonance spectroscopy (MRS) studies in humans have suggested alterations in the excitation/inhibition (E/I) balance in experimental and chronic pain (Pasanta et al., 2023; Peek et al., 2020; Zhao et al., 2017). Furthermore, changes in E/I likely relate to central sensitization (Woolf, 2011), which has been observed in many chronic pain states, including nociplastic pain conditions (Kaplan et al., 2024). Thus, assessing cerebral excitation and inhibition promises insights into the neural mechanisms of chronic pain. Moreover, such insights might aid in developing clinically useful biomarkers of chronic pain. However, non-invasive, scalable, whole-brain measurements of the E/I balance in people with chronic pain have not been performed so far.
Recently, novel electroencephalography (EEG) measures have been proposed as non-invasive proxies for E/I with high translational potential (Ahmad et al., 2022). Specifically, the aperiodic exponent of the EEG power spectrum has been highlighted as a measure of the E/I ratio (Gao et al., 2017; Lombardi et al., 2017). The power spectrum of EEG signals (Figure 1) typically features narrow band peaks over a power distribution following a 1/fχ shape, with power decreasing linearly as frequency increases on a log-log scale (Buzsáki and Draguhn, 2004; Miller et al., 2009). Traditional EEG analyses have focused on oscillatory activity (peaks in the power spectrum) and activity at specific frequency bands. However, recent studies have highlighted the physiological significance of the power spectrum’s aperiodic (1/fχ) component (Donoghue et al., 2020; Lendner et al., 2020; Ostlund et al., 2021; Tröndle et al., 2022). In particular, the slope of the aperiodic component, given by the aperiodic exponent (χ), has been linked to the E/I balance. Higher exponents (steeper spectra) indicate a shift toward inhibition, and lower exponents toward excitation (Gao et al., 2017; Lendner et al., 2020). Moreover, E/I disruptions have been proposed as a common feature underlying different neuropsychiatric disorders (Ahmad et al., 2022; Sohal and Rubenstein, 2019). Accordingly, first EEG studies assessing the aperiodic component in several brain disorders have revealed changes in schizophrenia (Molina et al., 2020), ADHD (Ostlund et al., 2021; Robertson et al., 2019), and Alzheimer’s disease (Martínez-Cañada et al., 2023).
In the present study, we characterized the E/I balance in people with chronic pain non-invasively. To this end, we analyzed the aperiodic component of resting-state EEG recordings in a large cohort of people with chronic pain (N=149), compared it to healthy participants (N=115), and related it to pain intensity. We hypothesized that changes in the aperiodic exponent would particularly occur in the medial prefrontal cortex (mPFC), as this is a central hub for pain processing and E/I imbalances have been observed in this area (Bliss et al., 2016; Kummer et al., 2020; Shiers and Price, 2020; Tan and Kuner, 2021). We complemented this region-of-interest analysis with a whole-brain analysis of the aperiodic component. Furthermore, we investigated aperiodic exponents in subtypes of chronic pain.
The current study rigorously pursued open and reproducible science practices. We pre-registered the study (https://osf.io/xshmy), performed blind analyses with Bayesian hypothesis testing, and made all data and code openly available (Hardwicke and Wagenmakers, 2023; MacCoun and Perlmutter, 2015; Wagenmakers et al., 2018). Additionally, we conducted a multiverse analysis to assess the robustness of the results across different analytical choices (Simonsohn et al., 2020; Steegen et al., 2016). This is particularly relevant as aperiodic activity is an emerging metric for which the influence of several methodological decisions is not fully clear.
Together, the study is intended to represent a step forward in developing translational, non-invasive, and robust tools to better understand chronic pain mechanisms. Such insights might pave the way for developing clinically useful biomarkers and novel pharmacological and neuromodulatory treatments for chronic pain.
Results
To investigate the E/I balance in chronic pain, we analyzed the aperiodic exponent of the power spectrum in 149 people with different types of chronic pain and 115 healthy participants. The aperiodic component was characterized in eyes-closed resting-state EEG recordings following a pre-registered approach (https://osf.io/xshmy). Specifically, we addressed three research questions (Figure 1). We first studied how the aperiodic activity in the mPFC relates to a mixed sample of people with chronic pain (Q1). We further asked how aperiodic activity in the mPFC relates to different types of chronic pain (Q2). Finally, we assessed the relationship between aperiodic activity and chronic pain beyond the mPFC in a whole-brain analysis (Q3). For each question, we performed two pre-registered analyses. First, we analyzed whether aperiodic exponents differ between people with chronic pain and healthy participants. Second, we analyzed whether aperiodic exponents in people with chronic pain correlate with their average pain intensity.
We investigated the research questions in a two-step approach. We first estimated the aperiodic activity and conducted the two analyses according to the pre-registration. As a second step, we conducted a non-pre-registered multiverse analyses to test the robustness of the results, since the aperiodic component is an emerging metric, and it is unclear how different methodological decisions affect its estimation. Specifically, we performed specification curve analyses (Simonsohn et al., 2020), a type of multiverse analysis that includes the identification of theoretically and statistically valid ways of analyzing the data (specifications), the graphical display of evidence in a descriptive specification curve, and inference tests across specifications to reach an overall conclusion on the results. One specification thereby refers to one combination of analytical settings, i.e., one version of the analysis.
We identified five methodological decisions that could influence the estimation of aperiodic activity, each with two or three equally valid parameter settings (Table 1, Methods). The combination of all these parameter settings yielded 48 different specifications, including the specification of the pre-registered analyses. Thus, in each multiverse analysis, we estimated the aperiodic exponent in 48 different ways and repeated the corresponding pre-registered analysis. The result is a specification curve, depicting in the x-axis the specifications, sorted by increased effect size, in the y-axis the effect size, and color-coded the degree of evidence (BF10). We also quantified how likely the observed specification curve was, compared to a null distribution of curves indicating no effect with three inference tests (Simonsohn et al., 2020). We performed multiverse analyses for research questions Q1 and Q2. We did not perform multiverse analysis for Q3, as it would yield a higher number of results than what can be statistically addressed and presented (please refer to Methods).
How does aperiodic activity in the mPFC relate to chronic pain?
We first compared the aperiodic exponents in the mPFC between the entire sample of 149 people with chronic pain and 115 healthy participants. We controlled the potential confound of age in aperiodic exponents by regressing age out from all the studied variables. Average aperiodic exponents before age correction were 1.10±0.16 (mean ± std) for the healthy participants and 1.10±0.18 for the participants with chronic pain. The results provided moderate evidence against a difference between groups (Figure 2A, Cohen’s d=–0.011, BF10=0.136, Bayesian two-sided independent samples t-test). The specification curve analysis confirmed this finding (Figure 2B). All specifications indicated moderate evidence against a difference in aperiodic exponents between groups (all BF10 <0.22). Statistical inference on the specification curve indicated that the curve was not significantly different from a null distribution of specification curves (pmedian = 0.89, pshare = 1, paggregate = 0.95, please refer to the Methods for a description of the test statistics). Thus, both pre-defined and multiverse analysis provided robust evidence against a difference in aperiodic exponents in the mPFC between healthy participants and people with chronic pain.
Next, we correlated the aperiodic exponents in the mPFC of people with chronic pain with their average pain intensity. The results provided moderate evidence against a correlation (Figure 3A, Pearson’s R=0.006, BF10=0.103, Bayesian correlation test). The specification curve analysis confirmed this result (Figure 3B). All specifications indicated moderate evidence against a correlation between pain ratings and aperiodic exponents (all BF10 <0.26). Statistical inference on the specification curve indicated that the curve was not significantly different from a null distribution of specification curves (pmedian = 0.63, pshare = 1, paggregate = 0.74). Apart from the aperiodic exponent, the aperiodic offset is another parameter estimated when computing the aperiodic component. The aperiodic offset has a less clear physiological interpretation than the aperiodic exponent. Nevertheless, we also investigated whether it was altered in the mPFC in people with chronic pain or was related to pain intensity (Figure 3—figure supplement 1). We found moderate evidence against a difference in aperiodic offsets between people with chronic pain and healthy participants (Cohen’s d=–0.03, BF10=0.14). Furthermore, we found moderate evidence against a correlation between aperiodic offsets and average pain ratings (R=0.06, BF10=0.14). Moreover, we exploratively investigated the effects of age on aperiodic exponents in the mPFC in our dataset. The results are shown in Figure 3—figure supplement 2.
In summary, we found evidence against a difference in aperiodic activity in the mPFC between healthy participants and a diverse cohort of people with chronic pain. We also found evidence against a correlation between aperiodic exponents and pain intensity in people with chronic pain. Multiverse analyses confirmed the robustness of these findings across different analytical parameters.
How does aperiodic activity in the mPFC relate to different types of chronic pain?
The sample of people with chronic pain included participants with different types of chronic pain. Thus, effects specific to certain subtypes of pain might have been missed. Therefore, we repeated pre-registered and multiverse analyses for the two largest clinical subgroups, people with Chronic Back Pain (CBP, N=80) and people with Chronic Widespread Pain (CWP, N=33). We randomly drew a subsample of healthy participants of equal size for each subgroup, matched in age, gender, and dataset to which they belonged (see Methods). Average aperiodic exponents before age correction were 1.11±0.19 for the people with CBP and 1.02±0.13 for the people with CWP.
In people with CBP (Figure 4), we found moderate evidence against a difference between the mPFC aperiodic exponents of people with CBP and healthy participants (Cohen’s d=–0.12, BF10=0.23). The specification curve analysis confirmed this result across 46 specifications, with only two showing inconclusive evidence (0.17 < BF10<0.36). Statistical inference on the specification curve indicated that the curve was not significantly different from a null distribution of specification curves (pmedian = 0.59, pshare = 1, paggregate = 0.61). There was also moderate evidence against a correlation between the aperiodic exponent and pain intensity in people with CBP (Pearson’s R=–0.09, BF10=0.19). The specification curve analysis indicated that this result was robust, with 44 specifications showing evidence in favor of the null hypothesis and only four showing inconclusive evidence (0.14 < BF10<0.42). Statistical inference on the specification curve indicated that the curve was not significantly different from a null distribution of specification curves (pmedian = 0.23, pshare = 1, paggregate = 0.35).
In people with CWP (Figure 5), we found inconclusive evidence regarding a difference in the mPFC aperiodic exponent between people with CWP and healthy participants (Cohen’s d=0.26, BF10=0.41). The specification curve analysis revealed that approximately 1/3 of the specifications supported the null hypothesis (specifications 7–26, 0.25 < BF10<0.33), while 2/3 of the specifications showed inconclusive evidence (specifications 1–6 and 27–48, 0.34 < BF10<0.83). Statistical inference on the specification curve indicated that the curve was not significantly different from a null distribution of specification curves (pmedian = 0.31, pshare = 1, paggregate = 0.41). Furthermore, we found moderate evidence against a correlation between the aperiodic exponent and pain intensity in people with CWP (R=0.03, BF10=0.22). The specification curve analysis provided further evidence against a correlation in approximately half of the specifications (specifications 22–48, 0.22 < BF10<0.33). However, 20 specifications showed inconclusive evidence (specifications 2–21, 0.33 < BF10<1.1), and one specification indicated a negative correlation between aperiodic exponents and pain (specification 1, BF10=3.1). Statistical inference on the specification curve indicated that the curve was not significantly different from a null distribution of specification curves (pmedian = 0.25, pshare = 0.17, paggregate = 0.49).
In summary, we found moderate evidence against altered aperiodic exponents in the mPFC and against a correlation between aperiodic exponents and pain intensity in people with chronic back pain. However, we found inconclusive evidence for the role of aperiodic exponents in the mPFC in people with chronic widespread pain.
How does aperiodic activity beyond the mPFC relate to chronic pain?
So far, we focused on the medial prefrontal cortex due to previous evidence for its relevance to chronic pain (Bliss et al., 2016; Kummer et al., 2020; Shiers and Price, 2020; Tan and Kuner, 2021). However, pain is a network phenomenon involving multiple brain areas. We, therefore, extended our analyses beyond the mPFC and investigated aperiodic exponents at a whole-brain level (Figure 6).
We estimated aperiodic exponents using the pre-defined settings from the previous analyses in 100 points of the MNI template. Brain locations were the centroid regions of the 100-parcel version of the Schaefer atlas with 17 networks (Schaefer et al., 2018; Yeo et al., 2011). To assess differences in the aperiodic exponent between healthy participants and people with chronic pain, we performed independent sample two-sided t-tests between groups at every region. As there are no standards for correcting Bayesian statistics for multiple comparisons, we performed frequentist t-tests and adjusted the p-values with the resampling-based FDR correction (Yekutieli and Benjamini, 1999). We did not find a difference in aperiodic exponents between healthy participants and people with chronic pain in any of the 100 regions (all padj >0.05,–0.39<Cohen’s d<0.06). Raw values of the age-corrected aperiodic exponents of people with chronic pain and healthy participants are available in Figure 7.
To investigate the relationship between the aperiodic exponent and pain intensity in people with chronic pain beyond the mPFC, we performed a Pearson’s correlation between the two variables at each region of interest. Again, we did not observe any significant correlation after multiple comparisons correction with resampling-based FDR (all padj >0.05,–0.11<R < 0.22). Please note that no whole-brain multiverse analyses were conducted as the group comparisons and correlations would include 100 multiverse analyses each.
In summary, we did not observe any difference in aperiodic exponents between healthy participants and people with chronic pain at a whole brain level. We also did not observe a correlation between aperiodic exponents and pain intensity in any brain region beyond the mPFC.
Discussion
In the present study, we non-invasively characterized the balance between cerebral excitation and inhibition in people with chronic pain. To this end, we analyzed the aperiodic component of resting-state EEG data in a large cohort of people with chronic pain, compared it to healthy participants, and related it to pain intensity. Our results did not show differences in the aperiodic component between people with chronic pain and healthy participants nor relationships between the aperiodic component and pain intensity. These findings are strengthened by preregistration and Bayesian statistics, which allow the quantification and interpretation of both positive and negative results. Multiverse analyses confirmed the robustness of the results across a wide range of analytical decisions. In this way, the findings enhance the understanding of the brain mechanisms of chronic pain and can guide future research on developing non-invasive markers of cerebral excitability in chronic pain and other brain disorders.
The observed lack of a relationship between the aperiodic component and chronic pain can be due to several reasons. First, there is indeed no relationship between cerebral E/I and chronic pain. Considering the extensive evidence for pain-related changes in cerebral E/I in animals and humans, this explanation appears unlikely. In animal models of chronic pain, numerous studies have found alterations of E/I, mainly in the prefrontal cortex, but also in other cortical and subcortical brain areas (Bliss et al., 2016; Kummer et al., 2020; Shiers and Price, 2020; Tan and Kuner, 2021). In humans, magnetic resonance spectroscopy findings show abnormal glutamate and GABA signaling in cortical and subcortical regions in people with chronic pain (Peek et al., 2020; Zhao et al., 2017).
Second, there is a relationship between E/I and chronic pain, but it cannot be detected using EEG. This explanation appears reasonable since excitability changes in chronic pain might differ across cell types and brain regions (Bliss et al., 2016; Kummer et al., 2020; Shiers and Price, 2020; Tan and Kuner, 2021). Thus, possible E/I effects in chronic pain might be too subtle or too heterogeneous to be detectable at the scalp using EEG. However, EEG findings in brain disorders known to be associated with E/I alterations argue against this possibility (Martínez-Cañada et al., 2023; Molina et al., 2020; Ostlund et al., 2021; Robertson et al., 2019). Moreover, we have analyzed the aperiodic component exclusively in source space. However, source space analyses can yield imprecise results, particularly when default volume conduction models and electrode coordinates are used. A systematic exploration of different source space approaches and analyses in electrode space might address this issue.
Third, there is a relationship between E/I and chronic pain detectable with EEG, but aperiodic activity is an inappropriate measure of E/I. Although this explanation cannot be ruled out, a relationship between aperiodic activity and E/I has been supported by in-vivo electrophysiological recordings in rodents (Gao et al., 2017), computational models (Gao et al., 2017; Lombardi et al., 2017; Martínez-Cañada et al., 2023), the gabaergic drug propofol (Lendner et al., 2020), and neuromodulation studies (Martínez-Cañada et al., 2023). Cross-modal integration of findings from PET, MRS, and EEG studies would be desirable to strengthen the link between aperiodic EEG activity and E/I. Furthermore, future studies might also assess non-linear relationships between chronic pain and EEG measures of E/I. In addition, other EEG measures of E/I have been proposed (Ahmad et al., 2022) and could be explored in future studies.
Fourth, there is a relationship between E/I measured with the EEG aperiodic exponent and chronic pain, but our heterogeneous sample of people with chronic pain has hampered the detection of a small effect. It is plausible that E/I imbalances are more prominent in specific subtypes of chronic pain. For example, clinical observations suggest that changes in excitability are particularly relevant in chronic widespread pain (CWP) (Sluka and Clauw, 2016) and migraine (O’Hare et al., 2023), but possibly less in other chronic pain conditions. Our findings are compatible with this possibility, as they showed robust negative evidence in people with CBP, but partially inconclusive evidence in our small sample of people with CWP. Thus, future studies might include larger sample sizes and focus on specific subtypes of chronic pain, for example CWP or migraine. They might also include a standardized and more extensive assessment of clinical measures, as well as the analysis of EEG recordings with eyes open. Moreover, such studies should be designed and sufficiently powered to assess the effects of gender and medication on the aperiodic EEG component.
Considering the negative findings of the present study, does it make sense to further investigate the role of the aperiodic component of resting-state EEG in chronic pain? The possible gain and the necessary efforts might guide this decision. We consider the potential gain high. Finding a non-invasive measure of cortical excitability, a critical brain mechanism of chronic pain related to central sensitization (Woolf, 2011) and nociplastic pain (Kaplan et al., 2024), would significantly advance our understanding of the disease. Moreover, it could help to develop a broadly available, cost-efficient, and scalable biomarker. Such a biomarker could fulfill clinically valuable functions such as determining the risk/susceptibility to develop chronic pain or the prognosis of the disease and predicting treatment responses (Davis et al., 2017). A biomarker close to an essential neural mechanism of chronic pain, that is increased excitability, loss of inhibition, and central sensitization, is particularly promising. These arguments must be weighed against the efforts needed to obtain further insights into the role of the aperiodic EEG in chronic pain. EEG is broadly available, and the efforts to perform EEG recordings are decreasing with progress in EEG hardware with dry electrodes (Ng et al., 2022) and mobile use (Niso et al., 2023). Thus, even when aiming for large sample sizes, as opposed to other non-invasive measures of brain activity, the efforts are comparably low. Taken together the high possible gain and the comparatively low efforts, it might be worth further pursuing this line of research.
The present findings might also have implications for research beyond chronic pain. The coordination of excitatory and inhibitory brain activity is a fundamental property of brain function (Turrigiano and Nelson, 2004), and the E/I balance has been shown to play an important role in various brain functions and disorders. For instance, E/I imbalances have been frequently observed in autism (Sohal and Rubenstein, 2019) and schizophrenia (Liu et al., 2021). Research on the E/I balance in all these domains could benefit from methodological exchange, integration, and harmonization across fields and disorders. The present approach might represent an example of how a human EEG study can be designed to make positive as well as negative findings most valid and informative.
In conclusion, our study did not find a relationship between the aperiodic component of resting-state EEG and chronic pain, supported by high standards of open and reproducible science. These findings should guide, rather than discourage, further research on possible EEG biomarkers of chronic pain. Future studies might focus on specific chronic pain subtypes and aim at larger, multi-site samples. By pursuing this approach, we may identify a non-invasive and scalable marker of altered cerebral excitability in chronic pain. Eventually, this might help to develop biomarkers and even new therapeutic targets and endpoints aiming at restoring the E/I balance (Martínez-Cañada et al., 2023; van Bueren et al., 2023).
Methods
Study design
We investigated the relationship between E/I and pain using three public datasets acquired at the PainLabMunich, Technical University of Munich, Germany. The datasets included resting-state EEG data and clinical and demographic variables of 151 people with chronic pain and 115 healthy participants. Two datasets were already published and were re-analyzed. The third dataset was not analyzed so far and has been made publicly available for this project (https://osf.io/qgfma/). As each dataset was initially acquired for a different project, the EEG recording paradigm and the collected clinical and demographic variables slightly differed between datasets. We refer the reader to the original publications for further details regarding inclusion/exclusion criteria of participants, administered questionnaires, and technical characteristics of the EEG recordings (Heitmann et al., 2022; May et al., 2019; Ta Dinh et al., 2019; Tiemann et al., 2012).
This study was pre-registered in OSF Registries on 11 May 2023 (https://osf.io/xshmy) as a secondary data pre-registration. In the pre-registration, we hypothesized (1) differences in aperiodic activity in the mPFC between people with chronic pain and healthy participants and (2) a relationship between aperiodic activity in the mPFC and pain in people with chronic pain. We performed these two pre-registered analyses for our three research questions (Figure 1). We complemented each pre-registered analysis with a multiverse analysis to check the robustness of the results across different methodological decisions.
Participants
Dataset 1 contains cross-sectional data of 101 people with different chronic pain conditions and 88 healthy participants. The participants were recorded in three different cohorts. The first one included 20 people with chronic widespread pain and 22 healthy participants (Tiemann et al., 2012). The second one contained 34 people with chronic back pain (May et al., 2019). The third one contained 47 people with different types of chronic pain and 66 healthy participants (Ta Dinh et al., 2019). Dataset 2 is a longitudinal dataset comprising 50 participants with different types of chronic pain who underwent multimodal pain therapy (Heitmann et al., 2022). Only the first session (‘baseline’) is analyzed in the present study. Dataset 3 has not been previously analyzed or published and includes 27 healthy participants.
The gender ratio of both groups was similar (people with chronic pain: 98 females and 53 males, healthy: 74 females and 41 males). The total sample included participants of a wide age range (people with chronic pain: 86–18, healthy: 79–18) but mostly white men and women who were currently pursuing or previously acquired higher education. Thus, the study sample constitutes a specific sample not necessarily representative of Germany or any region.
Data acquisition
During data acquisition sessions, participants first completed clinical and demographic questionnaires. Then, brain activity was recorded using EEG during the resting state, i.e., participants were asked to stay relaxed and wakeful without doing any particular task.
EEG data was recorded with 64 pasive electrodes in the standard 10/20 positions (Easycap, Hersching, Germany) and a BrainAmp MR plus amplifier (Brain Products GmbH, Gilching, Germany). For most participants, two blocks of 5 min each were recorded, one with eyes closed and one with eyes open. This study only analyzed the eyes closed condition, available for all participants. During recording, electrodes were referenced to FCz and grounded at AFz. Data were obtained at a sampling frequency of 1000 Hz and were band-pass filtered online between 0.016 and 250 Hz.
The following demographic and clinical variables were used in this study: (1) Age of the participant at the recording date in years. (2) Average pain intensity, defined as the average pain intensity during the past four weeks obtained with a numerical rating scale from the PainDetect questionnaire (0=no pain, 10=maximum pain). In the Chronic Widespread Pain subset of Dataset 1 (N=20), this information was not available, and the variable ‘current pain at the time of the recording’ was used instead (0=no pain, 10=maximum pain). (3) Clinical diagnosis.
Data gathering, blinding, and preprocessing
The three datasets were merged into one, and participants’ IDs were randomly reassigned in the form sub-XXX, where XXX indicated a number between 001 and 266. In this way, researchers were blinded to the initial dataset to which the participant belonged. One participant in Dataset 2 was excluded because they had already been included in Dataset 1. Another participant of Dataset 2 was excluded due to technical problems during the recording. Therefore, the final sample size was 149 people with chronic pain and 115 healthy participants.
Next, two versions of the demographic data were prepared. In one version, participants’ group labels (chronic pain or healthy) were randomly re-assigned, keeping the ratio of the original labels. In the other version pain intensity scores were randomly re-assigned to people with chronic pain. All pre-registered analyses and statistical scripts were prepared using these modified versions of the demographic data. In this way, researchers were blinded during the computation and visualization of intermediate results, preventing p-hacking and confirmation bias (Hardwicke and Wagenmakers, 2023; MacCoun and Perlmutter, 2015). Only when the statistical and visualization scripts were fully developed were they executed with the original demographic data.
EEG data was automatically preprocessed in MATLAB (Mathworks, Natick, MA), with the default settings of the DISCOVER-EEG pipeline version 1.0.0. (Gil Ávila et al., 2023; Pernet et al., 2020) using the EEGLAB Toolbox (Delorme and Makeig, 2004). Briefly, raw data in the BIDS standard were loaded into EEGLAB and downsampled to 250 Hz. Line noise was removed with the CleanLine plug-in. Then, data were high-pass filtered at 0.5 Hz, and bad channels were rejected with the clean_rawdata function. Data was re-referenced to the average reference, and Independent Component Analysis with the ‘runica’ algorithm was performed. Independent components categorized as muscle or eye artifacts by the ICLabel classifier (Pion-Tonachini et al., 2019) were removed. Previously discarded channels were interpolated using spherical splines, and further time segments containing artifacts were removed using the Artifact Subspace Reconstruction method. Clean data was then saved to disk as a continuous segment and loaded into the Fieldtrip toolbox (Oostenveld et al., 2011) for further processing. Continuous data were segmented into 2 s or 5 s (see multiverse specifications in Table 1) epochs with 50% overlap. Then, segmented data were band-passed between 0.5 and 100.5 Hz for source localization. No recording was discarded based on preprocessing outcomes.
For subgroup analyses (Research Question 2), we defined two subgroups based on the clinical diagnosis: people with Chronic Back Pain (CBP, N=80) and people with Chronic Widespread Pain (CWP, N=33). To perform these analyses, we randomly drew a subsample of healthy participants of equal size for each subgroup, matched in age, gender, and dataset to which they belonged following an iterative process. An initial random sample of healthy participants containing the same number of participants per initial study as the people with chronic pain subgroup was chosen. Then differences in age and gender between the samples were tested with Bayesian independent samples t-tests. If moderate evidence in favor of the null hypothesis was observed (BF <0.33) in both age and gender, the healthy sample was kept. If not, a different random sample of healthy participants was chosen, and the procedure was repeated.
Aperiodic component computation
The aperiodic component is a recent measure of brain activity. Thus, how different analytical choices impact its calculation has not yet been characterized. For that reason, we performed the analyses using a two-step approach. First, we conducted the pre-registered analyses, for which we had pre-defined the steps and parameters to estimate the aperiodic activity. Second, we performed a multiverse analysis to investigate the robustness of the results across different methodological ways of computing the aperiodic activity.
Specifically, we performed a specification curve analysis (Simonsohn et al., 2020) in which we first identified a set of specifications, i.e., equally valid analytical options to estimate aperiodic activity. We then estimated aperiodic activity, performed the pre-defined statistical test for each specification, and visualized the results in a specification curve. We conducted three inference tests on the specification curve to assess its statistical significance (see Statistical Analyses). In Table 1, we list all the analytical decisions and their parameters, and describe their potential influence on the estimation of the aperiodic activity. The combination of all these analytical options yielded 48 different specifications, including the approach of the pre-registered analysis.
Source localization
To spatially localize brain activity, we performed a source reconstruction of the preprocessed, segmented, and band-passed filtered signals. We projected the data to source space using an array-gain Linear Constrained Minimum Variance (LCMV) beamformer. As a source model, we used the centroid regions of the 100-parcel version of the Schaefer atlas with 17 networks (Schaefer et al., 2018; Yeo et al., 2011). The lead field was built using a realistically shaped volume conduction model based on the Montreal Neurological Institute (MNI) template available in FieldTrip (standard_bem.mat) and the source model. Spatial filters were constructed with the covariance matrices of the band-pass filtered data and the described lead fields. A 5% regularization parameter was set to account for rank deficiencies in the covariance matrix. The dipole orientation was fixed to the direction of the maximum variance following the most recent recommendations (Westner et al., 2022).
We defined our primary region of interest, the medial prefrontal cortex, as all the parcels from the above-mentioned version of the Schaefer atlas labeled medial prefrontal cortex (‘PFCm’), medial posterior prefrontal cortex (‘PFCmp’), and anterior cingulate cortex (‘Cinga’). Four parcels comprised the mPFC by this definition (Figure 1).
Power spectrum computation
We reconstructed the virtual time series at each source location using the spatial filter and band-pass filtered sensor-level data. Power at each source location was then calculated between 1 and 100 Hz with a Fast Fourier Transform using the Fieldtrip function ft_freqanalysis and the ‘mtmfft’ method. We used Slepian multitapers with +/-1 Hz frequency smoothing for the pre-registered analyses. Average power spectra values are visualized in Supp. Fig. 3. For other specifications of the multiverse analysis, we used a single Hanning taper (Table 1).
To investigate aperiodic activity in the mPFC, we constructed a representative power spectrum of this region by averaging the power spectra of the four source locations that defined the mPFC. This representative power spectrum of the mPFC was used later to model aperiodic activity. We investigated an alternative way to summarize mPFC activity in the multiverse analysis by estimating aperiodic activity for each source location and then averaging aperiodic parameters in the mPFC, i.e., the aperiodic exponent (Table 1).
Power spectrum parametrization
Power spectra were parametrized into periodic and aperiodic components with a re-implementation in MATLAB of the ‘spectparam’ algorithm, formally known as FOOOF (Donoghue et al., 2020). The algorithm was reimplemented to integrate it with the multiverse analyses computation. The power spectrum was modeled in the log-log scale as a sum of an aperiodic component and N oscillatory peaks, each modeled individually by a Gaussian function. The aperiodic component, L, was modeled using a Lorentzian function.
Here, b denotes the aperiodic offset, k the ‘knee’ parameter, which controls for the kink in the aperiodic component, F frequency, and χ the aperiodic exponent. The power spectrum model with parameters b, k, and χ was fitted to the data in an iterative procedure as in Donoghue et al., 2020 using non-linear least squares solvers implemented with the Matlab ‘lsqcurvefit’ function.
For the pre-registered analyses, we modeled the power spectrum using the settings proposed in the ‘spectparam’ tutorial, that is modeling the power spectrum in the 2–40 Hz range and not including the knee parameter. In the multiverse analysis, we also explored two alternative frequency ranges: the 40–60 Hz range (Gao et al., 2017), and the 1–100 Hz range. We also investigated the influence of including the knee parameter in estimating the aperiodic component (Table 1). The rest of the algorithm parameters were not modified except for the minimum detectable peak width, which was increased from 0.5 Hz to 1 Hz to match twice the frequency resolution of our data. For a detailed description of these parameters, we refer the reader to the ‘spectparam’ documentation (https://fooof-tools.github.io/fooof/index.html) and the pre-registration.
Additionally, we extracted two measures of goodness of fit of the modeled power spectra: the mean average error (MAE) and the explained variance (R2). With the pre-registered settings, the power spectra of all participants could be satisfactorily modeled (MAE = 0.03 ± 0.01, R2=0.98 ± 0.01 [mean ± std]). Thus, no participant was excluded based on a poor model fit. In the multiverse analysis, in the eight specifications in which the power spectrum was modeled between 40 and 60 Hz and the knee parameter was estimated, the fitting algorithm did not converge in a minority of participants. These participants were excluded in the statistical analysis of the respective specification. The maximum number of discarded participants in a specification was 39.
Statistical analyses
Effect size
In the context of pre-registration, we conducted a sensitivity analysis in G*power (Faul et al., 2007) to determine the effect size detectable with a two-tailed independent samples t-test. With our pre-defined sample size, given that the data had already been acquired, an α error probability of 0.05, and a statistical power of 0.95, we could detect medium effect sizes (Cohen’s d=0.44).
Statistical models
We investigated two main topics in our three research questions. First, we investigated whether the aperiodic component differs between people with chronic pain and healthy participants. Second, we explored the relationship between the aperiodic exponent and pain intensity in people with chronic pain. To control for age effects, we used linear regression models to remove the influence of age from aperiodic exponents and pain intensity ratings. Thus, all statistical tests were performed on the residuals obtained from the regression models.
Research questions 1 and 2 focused on the medial prefrontal cortex. For these research questions, we performed a Bayesian two-tailed independent samples t-test to compare the age-corrected aperiodic exponents in the mPFC between people with chronic pain and healthy participants. We also performed a Bayesian linear correlation to explore the relationship between the age-corrected aperiodic exponent in the mPFC and age-corrected pain intensity in people with chronic pain. These tests were both part of our pre-defined analyses and applied to each specification of the multiverse analysis. The results of the multiverse analysis were gathered and visualized using a specification curve, depicting each specification’s effect size and Bayes Factor.
Research question 3 investigated aperiodic activity at the whole-brain level. To this end, we aimed to address the same questions as outlined above but focused on one hundred different brain locations. To our knowledge however, it is impossible to address the multiple comparison problem in a Bayesian framework. Therefore, we performed the frequentist version of the aforementioned statistical tests at each source location. Then, we adjusted the p-values with the resampling-based FDR correction (Yekutieli and Benjamini, 1999). Due to the nature of the question and the high number of statistical comparisons, we decided not to perform a multiverse analysis.
We implemented the statistical tests in R (R Development Core Team, 2021) using the package ‘BayesFactor’ for the Bayesian tests and ‘specr’ for the visualization of the specification curve.
Inference criteria
For research questions 1 and 2, statistical inference of the pre-defined analyses was based on Bayes Factors (BF10) obtained from each statistical test. Moderate evidence favoring the alternative hypothesis was inferred if BF10 >3, inconclusive evidence if 1/3 < BF10<3, and moderate evidence favoring the null hypothesis if BF10 <1/3. Bayes factors larger than 10 or smaller than 1/10 were considered strong evidence in favor or against the alternative hypothesis, respectively.
In the multiverse analysis, we integrated the statistical evidence from all the specifications by performing three inference tests based on a permutation approach, as proposed in the original publication (Simonsohn et al., 2020). Of note, these tests are based on frequentist statistics as, to the best of our knowledge, Bayesian approaches to specification curve inference do not exist yet. Thus, for each multiverse analysis, we generated a null distribution of 500 specification curves by randomly shuffling the group labels or the pain ratings and recomputing the analyses for all specifications. Then, we conducted three inference tests to assess the likelihood of obtaining the original specification curve under the null hypothesis of no effect. In each test a p-value was derived by comparing a test statistic of the original specification curve to the corresponding distribution of test statistics of the randomized specification curves. The considered test statistics were: (1) The median effect size, i.e., Cohen’s d or Pearson’s correlation coefficient depending on the analysis (pmedian). (2) The share of specifications indicating evidence for an effect, i.e., with BF10 >3, and the same direction of effect as the original curve (pshare). (3) Z-values associated with p-values summarized across all specifications according to Stouffer’s method (paggr). Note that for this test a frequentist version of t-test or correlation test was computed for each specification. The inference criterion for each of the three tests was a p-value <0.05.
For research question 3, we inferred statistical significance if the FDR-adjusted p-values were lower than 0.05.
Deviations from the pre-registration
During the pre-registration of research question 1, we proposed to perform a Bayesian ANCOVA with group (people with chronic pain/healthy) as a fixed factor and age as a covariate in JASP. However, it was impossible to perform this test programmatically for all specifications during the multiverse analysis. Therefore, we switched to independent samples t-tests, available in the R package ‘BayesFactor’ on age-corrected aperiodic exponents.
Additionally, we pre-registered another hypothesis that depended on a positive finding concerning pre-registered research question 1. As there were no differences in aperiodic exponents in the mPFC between people with chronic pain and healthy participants, that hypothesis was not tested, and whole-brain effects were investigated for research question 3.
Data availability
All the EEG raw data in BIDS format (Pernet et al., 2019), demographic and clinical variables, and code used in this work are available in the OSF project https://osf.io/qgfma/.
References
-
Synaptic plasticity in the anterior cingulate cortex in acute and chronic painNature Reviews. Neuroscience 17:485–496.https://doi.org/10.1038/nrn.2016.68
-
Neuronal oscillations in cortical networksScience 304:1926–1929.https://doi.org/10.1126/science.1099745
-
Brain imaging tests for chronic pain: medical, legal and ethical issues and recommendationsNature Reviews. Neurology 13:624–638.https://doi.org/10.1038/nrneurol.2017.122
-
EEGLAB: an open source toolbox for analysis of single-trial EEG dynamics including independent component analysisJournal of Neuroscience Methods 134:9–21.https://doi.org/10.1016/j.jneumeth.2003.10.009
-
Parameterizing neural power spectra into periodic and aperiodic componentsNature Neuroscience 23:1655–1665.https://doi.org/10.1038/s41593-020-00744-x
-
G*Power 3: a flexible statistical power analysis program for the social, behavioral, and biomedical sciencesBehavior Research Methods 39:175–191.https://doi.org/10.3758/bf03193146
-
BookLongitudinal Resting-State Electroencephalography in Chronic Pain Patients Undergoing Interdisciplinary Multimodal Pain TherapyPAIN.
-
Periodic and aperiodic neural activity displays age-dependent changes across early-to-middle childhoodDevelopmental Cognitive Neuroscience 54:101076.https://doi.org/10.1016/j.dcn.2022.101076
-
Deciphering nociplastic pain: clinical features, risk factors and potential mechanismsNature Reviews. Neurology 20:347–363.https://doi.org/10.1038/s41582-024-00966-8
-
The medial prefrontal cortex as a central hub for mental comorbidities associated with chronic painInternational Journal of Molecular Sciences 21:3440.https://doi.org/10.3390/ijms21103440
-
Structural plasticity and reorganisation in chronic painNature Reviews. Neuroscience 18:113.https://doi.org/10.1038/nrn.2017.5
-
A selective review of the excitatory-inhibitory imbalance in schizophrenia: underlying biology, genetics, microcircuits, and symptomsFrontiers in Cell and Developmental Biology 9:664535.https://doi.org/10.3389/fcell.2021.664535
-
BookPrefrontal Gamma Oscillations Reflect Ongoing Pain Intensity in Chronic Back Pain PatientsHum Brain Mapp.
-
Power-law scaling in the brain surface electric potentialPLOS Computational Biology 5:e1000609.https://doi.org/10.1371/journal.pcbi.1000609
-
Memantine effects on electroencephalographic measures of putative excitatory/inhibitory balance in schizophreniaBiological Psychiatry. Cognitive Neuroscience and Neuroimaging 5:562–568.https://doi.org/10.1016/j.bpsc.2020.02.004
-
Excitation-inhibition imbalance in migraine: from neurotransmitters to brain oscillationsInternational Journal of Molecular Sciences 24:10093.https://doi.org/10.3390/ijms241210093
-
FieldTrip: Open source software for advanced analysis of MEG, EEG, and invasive electrophysiological dataComputational Intelligence and Neuroscience 2011:156869.https://doi.org/10.1155/2011/156869
-
Behavioral and cognitive correlates of the aperiodic (1/f-like) exponent of the EEG power spectrum in adolescents with and without ADHDDevelopmental Cognitive Neuroscience 48:100931.https://doi.org/10.1016/j.dcn.2021.100931
-
Functional MRS studies of GABA and glutamate/Glx - A systematic review and meta-analysisNeuroscience and Biobehavioral Reviews 144:104940.https://doi.org/10.1016/j.neubiorev.2022.104940
-
SoftwareR: A language andenvironment for statistical computingR Foundation for Statistical Computing, Vienna, Austria.
-
EEG power spectral slope differs by ADHD status and stimulant medication exposure in early childhoodJournal of Neurophysiology 122:2427–2437.https://doi.org/10.1152/jn.00388.2019
-
Specification curve analysisNature Human Behaviour 4:1208–1214.https://doi.org/10.1038/s41562-020-0912-z
-
Increasing transparency through a multiverse analysisPerspectives on Psychological Science 11:702–712.https://doi.org/10.1177/1745691616658637
-
Neocortical circuits in pain and pain reliefNature Reviews. Neuroscience 22:458–471.https://doi.org/10.1038/s41583-021-00468-2
-
Homeostatic plasticity in the developing nervous systemNature Reviews. Neuroscience 5:97–107.https://doi.org/10.1038/nrn1327
-
Bayesian inference for psychologyPart I: Theoretical Advantages and Practical Ramifications. Psychon Bull Rev 25:35–57.https://doi.org/10.3758/s13423-017-1343-3
-
Resampling-based false discovery rate controlling multiple test procedures for correlated test statisticsJournal of Statistical Planning and Inference 82:171–196.https://doi.org/10.1016/S0378-3758(99)00041-5
-
The organization of the human cerebral cortex estimated by intrinsic functional connectivityJournal of Neurophysiology 106:1125–1165.https://doi.org/10.1152/jn.00338.2011
Article and author information
Author details
Funding
Deutsche Forschungsgemeinschaft (SFB1158)
- Markus Ploner
Deutsche Forschungsgemeinschaft (PL321/14-1)
- Markus Ploner
Technical University of Munich (Neurotech)
- Markus Ploner
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The study was supported by the TUM Innovation Network Neurotechnology in Mental Health (NEUROTECH) and the Deutsche Forschungsgemeinschaft (PL 321/14–1, SFB1158).
Ethics
The study was a pre-registered analysis of public datasets.
Version history
- Preprint posted:
- Sent for peer review:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Version of Record published:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.101727. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2024, Gil Avila et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 283
- views
-
- 4
- downloads
-
- 0
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Neuroscience
Data-driven models of neurons and circuits are important for understanding how the properties of membrane conductances, synapses, dendrites, and the anatomical connectivity between neurons generate the complex dynamical behaviors of brain circuits in health and disease. However, the inherent complexity of these biological processes makes the construction and reuse of biologically detailed models challenging. A wide range of tools have been developed to aid their construction and simulation, but differences in design and internal representation act as technical barriers to those who wish to use data-driven models in their research workflows. NeuroML, a model description language for computational neuroscience, was developed to address this fragmentation in modeling tools. Since its inception, NeuroML has evolved into a mature community standard that encompasses a wide range of model types and approaches in computational neuroscience. It has enabled the development of a large ecosystem of interoperable open-source software tools for the creation, visualization, validation, and simulation of data-driven models. Here, we describe how the NeuroML ecosystem can be incorporated into research workflows to simplify the construction, testing, and analysis of standardized models of neural systems, and supports the FAIR (Findability, Accessibility, Interoperability, and Reusability) principles, thus promoting open, transparent and reproducible science.
-
- Neuroscience
Recent studies suggest that calcitonin gene-related peptide (CGRP) neurons in the parabrachial nucleus (PBN) represent aversive information and signal a general alarm to the forebrain. If CGRP neurons serve as a true general alarm, their activation would modulate both passive nad active defensive behaviors depending on the magnitude and context of the threat. However, most prior research has focused on the role of CGRP neurons in passive freezing responses, with limited exploration of their involvement in active defensive behaviors. To address this, we examined the role of CGRP neurons in active defensive behavior using a predator-like robot programmed to chase mice. Our electrophysiological results revealed that CGRP neurons encode the intensity of aversive stimuli through variations in firing durations and amplitudes. Optogenetic activation of CGRP neuron during robot chasing elevated flight responses in both conditioning and retention tests, presumably by amyplifying the perception of the threat as more imminent and dangerous. In contrast, animals with inactivated CGRP neurons exhibited reduced flight responses, even when the robot was programmed to appear highly threatening during conditioning. These findings expand the understanding of CGRP neurons in the PBN as a critical alarm system, capable of dynamically regulating active defensive behaviors by amplifying threat perception, ensuring adaptive responses to varying levels of danger.