ORIGINAL RESEARCH article
Increased Low- and High-Frequency Oscillatory Activity in the Prefrontal Cortex of Fibromyalgia Patients
- 1Neuroscience Research Institute, Seoul National University College of Medicine, Seoul, South Korea
- 2Department of Brain and Cognitive Sciences, Seoul National University College of Natural Sciences, Seoul, South Korea
- 3Department of Neurosurgery, Seoul National University Hospital, Seoul, South Korea
Recent human neuroimaging studies have suggested that fibromyalgia (FM), a chronic widespread pain disorder, exhibits altered thalamic structure and function. Since the thalamus has extensive reciprocal connection with the cortex, structural and functional thalamic alterations in FM might be linked to aberrant thalamocortical oscillation. This study investigated the presence of abnormal brain rhythmicity in low- and high-frequency bands during resting state in patients with FM and their relationship to clinical pain symptom. Spontaneous magnetoencephalography (MEG) activity was recorded in 18 females with FM and 18 age- and sex-matched healthy control (HC) subjects. The most remarkable finding was that FM patients had general increases in theta, beta and gamma power along with a slowing of the dominant alpha peak. Increased spectral powers in the theta-band were primarily localized to the left dorsolateral prefrontal (DLPFC) and orbitofrontal cortex (OFC). Beta and gamma over-activation were localized to insular, primary motor and primary and secondary somatosensory (S2) cortices, as well as the DLPFC and OFC. Furthermore, enhanced high-frequency oscillatory activities in the DLPFC and OFC were associated with higher affective pain scores in patients with FM. Our results demonstrate that FM patients feature enhanced low- and high-frequency oscillatory activity in the brain areas related to cognitive and emotional modulation of pain. Increased low- and high-frequency activity of the prefrontal cortex may contribute to persistent perception of pain in FM. Therapeutic intervention based on manipulating neural oscillation to restore normal thalamocortical rhythmicity may be beneficial to pain relief in FM.
Fibromyalgia (FM) is a chronic pain disorder characterized by the widespread pain and tenderness, and which is often accompanied by affective and cognitive symptoms (Wolfe et al., 1990, 2010; Bartley et al., 2009; Schmidt-Wilcke and Clauw, 2011). Although the underlying cause of the symptoms of FM remains elusive, recent neuroimaging findings suggest that abnormally increased cortical excitability and dysfunctional endogenous pain modulation are important in maintenance of FM pain (Jensen et al., 2009; Mhalla et al., 2010; Woolf, 2011; de Tommaso et al., 2011; Dailey et al., 2013; Villamar et al., 2013; Kim D. J. et al., 2014; Loggia et al., 2014; López-Solà et al., 2014; Pujol et al., 2014; Choi et al., 2015; Foerster et al., 2015; Lim et al., 2015).
Of particular note, the thalamus is the region of the brain that mediates transmission of sensory and pain signals to the cortex. Single-photon-emission computed tomography studies reported reduced resting regional cerebral blood flow in the thalamus in FM patients compared with healthy control (HC) subjects (Mountz et al., 1995; Kwiatek et al., 2000). Abnormal thalamic activity in FM was also reported in functional magnetic resonance imaging (fMRI) studies during a painful pressure stimulus. Patients with FM showed no thalamic increases (Gracely et al., 2002) or deactivation of the thalamus as opposed to increased activity in HC subjects (Jensen et al., 2009). A recent study reported the reduced functional connectivity between the thalamus and the lateral orbitofrontal cortex (OFC) in FM patients (Jensen et al., 2012b). FM has been also associated with altered brain structure in the thalamus in the form of decreased gray matter in the left posterior thalamus (Schmidt-Wilcke et al., 2007) and lower fractional anisotropy in the bilateral posterior thalami (Lutz et al., 2008). Since the thalamus has extensive reciprocal connection with the cortex (Jones, 2001), it seems likely that structural and functional thalamic alterations in FM might be linked to aberrant thalamocortical oscillation.
Thalamocortical dysrhythmia (TCD) has been proposed as the underlying mechanism of chronic neuropathic pain and other disorders (Llinás et al., 1999, 2005; Schulman et al., 2005, 2011). TCD in chronic pain is characterized by abnormal oscillatory activity in the form of increased spectral power and a dominant peak shifted towards a lower frequency (Llinás et al., 2005; Sarnthein et al., 2006; Stern et al., 2006; Jensen et al., 2013; Vuckovic et al., 2014). Spontaneous rhythmic activity measured by magnetoencephalography (MEG) or electroencephalography (EEG) shows enhanced theta and beta activity in patients with neuropathic pain compared to HCs (Sarnthein et al., 2006; Stern et al., 2006; Walton et al., 2010). Theta and beta bands over-activations in patients with neuropathic pain are localized to multiple area of the cortical pain matrix, including the prefrontal, anterior cingulate, and insular cortices as well as primary (S1) and secondary somatosensory (S2) cortices (Stern et al., 2006). These over-activations in the cortical pain matrix can be reduced after a therapeutic lesion in the central lateral nucleus of the thalamus (Stern et al., 2006). Therefore, TCD may perpetuate and exacerbate the chronic pain condition when abnormal rhythmic activity is present in the cortical pain network. Especially, high-frequency over-activation, referred to as an edge effect, can lead to pain symptoms (Llinás et al., 1999, 2005; Walton et al., 2010). However, to our knowledge, an abnormality in spontaneous brain activity in patients with FM has not been investigated.
The present study aimed to investigate the presence of abnormal brain rhythmicity in low- and high-frequency bands during resting state in patients with FM. By applying our MEG source localization approach to the cortical surface model (Kim and Chung, 2008; Kim J. S. et al., 2014), we further aimed to identify the cortical sources that generate different oscillatory activity between FM patients and HC subjects. We hypothesized that FM patients have increased spontaneous brain activity in both low- and high-frequency bands, and that these over-activations are related to clinical pain severity in FM patients.
Patients and Methods
The study participants were recruited from the outpatient clinics of the Rheumatology Departments of the Seoul National University Hospital and Hallym University Sacred Heart Hospital. Eligibility criteria for patients with FM were: (1) meeting the American College of Rheumatology 1990 criteria for primary FM (Wolfe et al., 1990); (2) a duration of widespread pain of at least 3 months but less than 10 years; (3) experiencing pain intensity of at least 40 on a 0–100 mm pain visual analog scale (VAS) over the past week; (4) age of 30–60 years; (5) female; (6) right-handed (Oldfield, 1971); and (7) willing to stop taking medications that might affect brain electrophysiology, such as analgesics, antidepressants, and anticonvulsants, at least 3 days before the assessments. Patients were excluded if they had: (1) secondary FM associated with inflammatory arthritis; (2) history of substance abuse; (3) signs of peripheral neuropathy or concomitant acute pain in the upper extremities; (4) hearing loss or use of hearing aids; (5) pregnancy or breastfeeding; or (6) contradictions for MEG or MRI assessments. Age-, gender-, and education-matched HC subjects were recruited by local advertisement. Exclusion criteria were the same as for patients. The study protocol was approved by the Institutional Review Boards at Seoul National University Hospital and Hallym University Sacred Heart Hospital and was conducted in compliance with the Declaration of Helsinki. All participants (19 FM patients and 21 HC subjects) provided written informed consent.
The Beck Depression Inventory (BDI; Beck et al., 1961), and Beck Anxiety Inventory (BAI; Beck et al., 1988) were assessed. The Fibromyalgia Impact Questionnaire (FIQ; Burckhardt et al., 1991) was used to assess the composite impact of FM symptoms. The sensory and affective components of pain were assessed using the short-form McGill Pain Questionnaire (SF-MPQ; Melzack, 1987). The questionnaires were performed on the same day as the MEG recordings, except for nine subjects. For nine subjects, the questionnaires were completed at the first visit and resting MEG data was acquired at the second visit. The mean ± standard deviation (SD) interval between sessions was 12.3 ± 5.7 days for FM (n = 6) and 18.3 ± 4.2 days for HC subjects (n = 3). The demographic and clinical characteristics of the participants are presented in Table 1.
Subjects were asked to refrain from consuming caffeinated beverages before the recording to avoid the caffeine-induced theta power decrease (Landolt et al., 2004). Spontaneous brain activity was consecutively recorded with subjects in a resting state with their eyes-open and resting state with their eyes-closed. In the eyes-open condition, the subjects were instructed to keep their gaze fixed on a cross in the center of the front wall. Subjects were instructed to relax, but to stay alert during the recording. Subject’s alertness during the experiment was checked by self-report after the MEG recording. The recording time of each condition was 270 s. The experimenter (M.L) monitored the subject’s head online using a video camera throughout the recordings. All subjects were carefully immobilized and did not markedly move their heads during the recordings. We focused our analysis on the eyes-closed condition (Sarnthein et al., 2006), as this condition is a relatively straightforward means to standardize the experiment between healthy and patient groups, and can minimize eye movements, such as eye blinks, which can affect frontal brain activity (van Diessen et al., 2015).
MEG and MRI data were acquired at Seoul National University Hospital. Spontaneous MEG signals were recorded with a VectorViewTM 306-channel whole-head neuromagnetometer (Elekta Neuromag, Helsinki, Finland) in a magnetically shielded room. The device contains 102 identical triple sensors, comprising two orthogonal planar gradiometers and one magnetometer. During the MEG recordings, the subjects sat comfortably under the helmet-shaped sensor array and were asked to keep their heads as still as possible. The exact head position with respect to the sensors was determined by measuring the magnetic signals produced by currents delivered to four head position indicator coils placed at known sites on the scalp. The two groups did not differ in their head position relative to the MEG sensors. The location of head position indicator coils and three anatomical landmarks, the nasion and two preauricular points, were measured by using a FASTRAKTM three-dimensional digitizer (Polhemus, Colchester, VT, USA), to allow alignment of the MEG and MRI coordinate systems. The x-axis passed through the two preauricular points with the positive direction to the right. The positive y-axis passed through the nasion, and the z-axis pointed upward. The MEG signals were band-pass filtered at 0.1–300 Hz and sampled at 1 kHz. To reduce environmental and biological noise, we applied the spatiotemporal signal space separation method using MaxFilter software version 2.2.10 (Elekta Neuromag; Taulu and Simola, 2006; Lim et al., 2014, 2015; Choi et al., 2015).
T1-weighted brain MRI was acquired for each subject using a Magnetom TrioTim 3T scanner (Siemens, Erlangen, Germany). Parameters were: sagittal acquisition with a 256 × 256 matrix; field of view = 250 mm; voxel size = 1 mm × 1 mm × 1 mm; slice thickness = 1.0 mm with no gap; repetition time/echo time = 1670/1.89 ms; flip angle = 9°; 1 excitation. One HC subject whose MRI examination could not be completed due to claustrophobia was excluded from further analysis. We found no group differences in whole-brain volume measures in our previous study (Kim D. J. et al., 2014).
We applied the constrained Laplacian anatomic segmentation using a proximity algorithm to reconstruct the inner and outer interfaces of the gray matter with 40,962 vertices for each hemisphere (Kim et al., 2005). A mid-surface was generated by averaging the corresponding points of the inner and outer interfaces. The vertices consisting of the mid surface were down-sampled to 5 mm intervals between points. The number of source points was about 8200 in the whole-brain model after the down-sampling. The source images were reconstructed on the down-sampled points. A linear interpolation was applied to display the source imaging on the cortical surface.
The solution of the forward problem has a linear relationship between the current sources and the measured signals:
where B is the vector of the measured signals, S is the vector of current sources, A is the lead field matrix, which describes the sensitivity pattern of each MEG sensor (magnetometers and gradiometers), and n is a noise vector. A spherical head model was used for source modeling (Hämäläinen et al., 1993). The center of the sphere was fit from manually selected inner skull points.
We used a standardized low-resolution brain electromagnetic tomography (sLORETA) algorithm (Pascual-Marqui, 2002) to localize the sources of oscillatory activity in various frequency bands: theta (5–8 Hz), alpha (8–13 Hz), beta (13–30 Hz), and gamma (30–48 Hz). The modified pseudo-statistics of sLORETA were used to analyze the absolute activation at each point (Kim and Chung, 2008; Kim J. S. et al., 2014). The weight matrix at the j-th point, Wj, is given by:
where Wj is a row vector, I is the identity matrix, and α is the regularization parameter. For any matrix M, M+ denotes its Moore-Penrose pseudo-inverse (Rao and Mitra, 1973). The denominator indicates the SD of the estimated current density.
The current density at the j-th point, Sj, is given by:
The signal at each polygonal point consisted of two orthogonal source vectors tangential to the surface of the spherical head model. Since two source activities are mutually orthogonal, source power at each location can be determined by vector addition.
The MEG signals were segmented into 4 s epochs without overlap. Epochs with amplitudes exceeding 3000 fT/cm for MEG channel or 150 μV for electrooculogram were excluded from the average. The data of one FM and two HC subjects were excluded from further analyses because MEG signals were heavily contaminated by MEG sensor noise or eye movements. The mean ± SD artifact free epochs were 61.4 ± 1.5 and 60.8 ± 0.7 (FM and HC, respectively). All artifact free epochs were extracted through a Hanning window. Power spectral analysis was performed with a fast Fourier transform of 4096 points. The real and imaginary parts of fast Fourier transform at each source point were estimated separately using the spatial filter matrix described in the previous section of this article. The details of the methods were described in previous studies (Jensen and Vanni, 2002; Kim and Chung, 2008; Kim J. S. et al., 2014). The Fourier transformed signals Sk(f) were calculated for each segment, k, with respect to the frequency of f. The real Sk(f)Re and the imaginary part Sk(f)Im of the Fourier-transformed signal were then applied to the sLORETA filter resulting in the current distributions, Qk(f)Re and Qk(f)Im, respectively. The absolute current estimates at the j-th source point for all the segments were averaged in the source space as follows:
Finally, we divided the frequency spectrum into four frequency bands.
Individual surface models were nonlinearly transformed to the template using a surface registration algorithm (Robbins et al., 2004). The spatially normalized sLORETA images of the two groups were compared at each vertex using general linear models adjusting for age. We performed a permutation tests with a P value of 0.01 to determine whether significant differences were not due to chance in multiple comparisons (Bullmore et al., 1999). Subjects were randomly assigned to groups across 10,000 new randomized analyses at each vertex, and the number of significant results (i.e., the power spectral density at any vertex that differed significantly between groups at a threshold of P < 0.01) that occurred in the real test for group differences were compared with the null distribution of significant results that occurred by chance (Kim J. S. et al., 2014).
The regional power spectra were extracted from the locations showing the maximal difference of theta power in the dorsolateral prefrontal (DLPFC) and OFC areas. Relationships between clinical pain symptoms of FM patients and the spectral power in the high-frequency (beta and gamma) oscillations were tested using Pearson’s correlation coefficient. P values < 0.0125 were considered to indicate significance after Bonferroni correction for multiple comparisons. All group data are presented as mean ± SD.
Power Spectral Analysis
The individual global power spectra for the FM patients and HC subjects are superimposed in Figures 1A,B, respectively. Figure 1C shows the grand average power spectra in the resting state of the FM patients and HC subjects. Notably, FM patients had general increases in theta, beta and gamma power and a dominant alpha peak shifted towards a lower frequency. The mean dominant alpha peak of the FM group (9.6 ± 0.6 Hz) was lower than that of the HC group (10.3 ± 0.8 Hz; P = 0.005; Figure 1D).
Figure 1. Resting state oscillatory brain activity. Superposition of individual mean power spectra of spontaneous brain activity recorded from all MEG gradiometer sensor pairs in FM patients (A) and HC subjects (B). Grand averaged power spectra of FM patients (red) and HC subjects (blue) (C). Box plots of the dominant peak frequency in the alpha band (D). On each box, the horizontal line is the median and the edges of the box are the 25th and 75th percentiles. Whiskers extend to the highest and lowest values within 1.5 times the interquartile range. The filled circles represent outliers. P value was determined by the two-tailed t-test for independent samples. FM, fibromyalgia; HC, healthy control.
Spontaneous peak over-activations within the theta (5–8 Hz), alpha (8–13 Hz), beta (13–30 Hz), and gamma (30–48 Hz) frequency bands were commonly localized to the left DLPFC and OFC (Figure 2). In the beta frequency band, peak over-activations were found in the anterior insular cortex (AIC) and primary motor cortex (M1), as well as S1 and S2 in the left hemisphere. Peak over-activations in the left AIC, M1, and S1 were also persistent in the gamma frequency band.
Figure 2. Source imaging of the significant power differences in various frequency bands between FM patients and HC subjects. DLPFC, dorsolateral prefrontal cortex; OFC, orbitofrontal cortex; AIC, anterior insular cortex; M1, primary motor cortex; S1, primary somatosensory cortex; S2, secondary somatosensory cortex.
Clinical Symptom Correlation
Figures 3A,D depict the power spectrum of the left DLPFC and OFC source, respectively. These areas are involved in cognitive-affective modulation of pain. Interestingly, abnormally increased rhythmicity was seen at not only low frequency, but also high frequency. Thus, we tested the relationship between over-activation in the high-frequency band in these areas and affective pain intensity in patients with FM. Both mean beta and gamma power in the DLPFC was positively associated with clinical ratings in the affective dimension of pain (r = 0.684, P = 0.002; r = 0.628, P = 0.005; Figures 3B,C). In addition, higher mean beta power in the OFC was also associated with higher score in the affective dimension of pain (r = 0.625, P = 0.006; Figure 3E). After correcting for multiple comparisons, higher mean gamma power in the OFC was marginally associated with increased affective pain intensity (r = 0.475, P = 0.047; Figure 3F). In our exploratory analysis, we did not observe any significant correlation between theta and alpha powers in the DLPFC or OFC and affective pain score in patients with FM (all Ps > 0.05).
Figure 3. Power spectrums of the prefrontal source activity (A,D) and relationship between affective pain intensity and high-frequency power of prefrontal regions in patients with FM (B,C,E,F). The yellow dots in the cortical surface indicate the locations of the maximal difference in theta power in the DLPFC (A) and OFC (D) areas. Red and blue lines indicate grand averaged source power spectra for FM patients and HC subjects, respectively. Mean beta (B,E) and gamma (C,F) powers were calculated from the individual DLPFC (A) and OFC (D) source power spectra in FM patients. DLPFC, dorsolateral prefrontal cortex; OFC, orbitofrontal cortex; SF-MPQ, short-form McGill Pain Questionnaire; FM, fibromyalgia; HC, healthy control.
To our knowledge, this is the first study to investigate changes of spontaneous neuromagnetic activity during the resting state and their relationship with clinical pain symptoms in patients with FM. FM patients displayed general increases in theta, beta and gamma power, along with a slowing of the dominant alpha peak. Spontaneous oscillatory over-activations were generated primarily in the prefrontal, insular, S1, S2, and M1 areas. Increased high-frequency oscillatory activity in the DLPFC and OFC was associated with higher affective pain scores in patients with FM.
Our findings of abnormally increased low-frequency activity in FM patients are partly in agreement with TCD as a pathophysiological model for chronic neuropathic pain (Llinás et al., 1999, 2005; Schulman et al., 2005; Saab, 2012; Leblanc et al., 2014). A prior resting EEG study showed that neuropathic pain patients exhibited higher spectral power over the frequency range of 2–25 Hz with a dominant peak shifted towards a lower frequency (Sarnthein et al., 2006). The authors also reported reduced pain symptom accompanied by decreases in theta EEG power after therapeutic lesion in the thalamus (central lateral thalamotomy) in a subgroup of patients. This result suggests that both neuropathic pain and low-frequency oscillations are presumably determined by thalamocortical loops (Sarnthein et al., 2006). This notion is supported by the evidence of high theta coherence between EEG and thalamic local field potentials from central lateral nucleus in neuropathic pain patients (Sarnthein and Jeanmonod, 2008). In addition, a recent resting MEG study also revealed a leftward shift of the dominant peak in the power spectrum and higher mean spectral power ratio (4–8 Hz/8–12 Hz) in patients with complex regional pain syndrome (Walton et al., 2010). Therefore, although the resting state MEG activity does not directly reflect thalamus activity, our results indicate that the pathophysiology of FM may share some characteristics with TCD.
The most obvious characteristic is the enhanced spontaneous oscillatory activity in not only the theta but also the beta and gamma bands of FM patients compared to HC subjects. In the framework of TCD (Llinás et al., 1999, 2005), the deafferentation of excitatory inputs on thalamic cells results in cell membrane hyperpolarization. This in turn generates low-threshold calcium spike bursts at theta frequency by deinactivation of calcium T-channels (Llinás and Jahnsen, 1982). Low-threshold calcium spike bursts of thalamic relay neurons exert an oscillatory activity in the thalamocortical loops (Sarnthein and Jeanmonod, 2008) and lead to over-production of theta oscillatory activity in the cortex (Sarnthein et al., 2006; Stern et al., 2006; Walton et al., 2010). Since the DLPFC and OFC receive projections from the mediodorsal nucleus of the thalamus (Klein et al., 2010), theta frequency activity in the mediodorsal thalamus could lead to theta over-production in the DLFPC and OFC. In the final step, thalamocortical modules in theta mode disrupt lateral inhibition on other cortical modules and these modules are over-activated at high frequencies (beta and gamma), which constitutes the edge effect (Llinás et al., 1999, 2005). This edge effect has been suggested to be responsible for the positive symptoms of chronic pain patients by acting in the cortical pain matrix (Schulman et al., 2005; Stern et al., 2006; Walton et al., 2010).
In this study, both low- and high-frequency over-activations in FM patients were commonly generated in the DLPFC and OFC regions. Our results are partly consistent with a recent resting state functional MRI study showing that FM patients had increased spectral power at 0.01–0.08 Hz, which was thought to be related to brain activity fluctuation, in several brain regions associated with pain processing including the DLPFC (Kim et al., 2013). DLPFC activity has been implicated in the cognitive and attentional processing of pain (Lorenz et al., 2003; Wiech et al., 2008). In addition, neuronal gamma oscillation has been suggested to be involved in attentional effects of pain (Hauck et al., 2007; Tiemann et al., 2010). Thus, we reasoned that over-activation of high frequency power, including the beta and gamma bands, in the DLPFC during the resting state may at least reflect sustained attention to spontaneous pain in patients with FM. The present results of a positive relationship between affective pain symptom and high-frequency power in the DLPFC support this interpretation.
Since the OFC is directly connected to the limbic system including amygdala, it is critically involved in regulation of emotion (Kringelbach and Rolls, 2004; Kringelbach, 2005). There is interesting evidence that activation of the OFC is related to the affective aspects of touch; e.g., pleasant and painful touch, while the S1 is activated more by neutral touch (Rolls et al., 2003). A prior study in patients with peripheral nerve injury revealed that brush-evoked allodynia activated the OFC and AIC, possibly due to the stronger emotional load of neuropathic pain (Witting et al., 2006). Based on these findings, we speculate that abnormally increased high-frequency power in the OFC may indicate emotional and affective dysregulation during resting state in FM. This hyperexcitation state would result in an overall decrease in brain resources for processing in emotional modulation. In this context, a recent study reported that a deficit in modulation of painful stimuli in a positive emotional context is related to less brain activation, including in the OFC, for FM patients (Kamping et al., 2013). FM has been also associated with structural change; i.e., increase in gray matter in the left OFC (Schmidt-Wilcke et al., 2007). Thus, the observed increases in spectral power in the OFC could be a consequence of its increased gray matter volume. Given a finite degree of synchronous firing, a higher cell population would lead to a stronger signal. Taken together, the functional and structural abnormalities of the OFC may contribute to altered emotional regulation of pain in FM (González-Roldán et al., 2013; Fallon et al., 2015).
The results of beta and gamma over-activation in the S1, S2 and M1 further support the view of altered excitability of the sensorimotor system in patients with FM (Mhalla et al., 2010, 2011; Lim et al., 2015). Recent studies suggested that beta and gamma oscillatory activity in the sensorimotor cortex have an important role in encoding of pain. Previous studies showed that painful stimuli suppress spontaneous beta-band oscillation in the sensorimotor cortex (Ploner et al., 2006; Kirveskari et al., 2010). In addition, pain induced gamma-band activity in the S1 (Gross et al., 2007; Hauck et al., 2013) and sensorimotor (Hauck et al., 2007) cortex was positively associated with the subjective perception of pain. Therefore, further studies are needed to investigate how the neuronal oscillatory response to painful stimuli in the sensorimotor cortex influence changes in the perception and processing of pain in FM. For example, a recent EEG study that reported beta-band suppression in the ipsilateral S1 and S2 regions during mechanical brush-evoked allodynia in FM patients indicated altered oscillatory processing of innocuous somatosensory input, which may then affect clinical symptom severity (Fallon et al., 2013).
Recent studies using multiple brain imaging modalities have demonstrated that chronic neuropathic pain is associated with structural and functional thalamic alterations, which may result in development of TCD (Henderson et al., 2013; Gustin et al., 2014; Youssef et al., 2014). Specifically, patients with painful trigeminal neuropathy display a reduction in gray matter volume in the somatosensory thalamus (ventroposterior nucleus) and reduced blood flow in the thalamic reticular nucleus (Henderson et al., 2013). The findings suggested that reduction of thalamic GABA level, likely as a consequence of decreased thalamic reticular nucleus activity, leads to altered thalamocortical rhythmic oscillation, which in turn results in the perception of persistent pain (Henderson et al., 2013). Therefore, we reasoned that abnormal spontaneous oscillatory activity in FM patients may be due to structural (Schmidt-Wilcke et al., 2007; Lutz et al., 2008) and functional (Mountz et al., 1995; Kwiatek et al., 2000; Gracely et al., 2002; Jensen et al., 2009, 2012a,b; Garcia-Larrea and Peyron, 2013) changes in the thalamus.
Several limitations should be considered in interpreting our results. First, since we did not measure thalamic activity, the present study could not provide direct evidence of thalamic dysregulation of cortical activity in FM. Second, although patients stopped all medications 3 days before the MEG recording, this may not have been sufficient time to exclude the long-term effect of medication on cortical excitability. On the other hand, this medication-discontinuation could have led to withdrawal symptoms. We confirmed that every patient was in good condition before and after the experiments and related problem was not reported by every patient. Even if the patients had symptoms from the discontinuation, these may have been minimal or not to the extent that could be clearly differentiated. Third, our source imaging results were lateralized to the left hemisphere. We have no clear explanation for this finding. One possibility is that abnormal thalamocortical low frequency activity triggered in the left prefrontal regions may cause less collateral inhibition on neighboring cortical modules, thereby promoting over-production of high frequency activity (Llinás et al., 1999). Further analysis of intra- and inter-hemispheric cross-frequency coupling (Adamchic et al., 2014; De Ridder et al., 2015) would facilitate the characterization of pathological oscillatory dynamics across brain networks in FM. Lastly, the source level findings would benefit from a more advanced forward and inverse modeling approach. In this study, our use of the spherical head model is a kind of compromise between precision and practicability. However, the oversimplified single-sphere model might have resulted in localization errors and better accuracy could have been obtained with a more realistic head model, especially for frontal areas (Stenroos et al., 2014). This means that there is some uncertainty whether the activity really originates from DLPFC and OFC. There is a pressing need for future studies which use a more realistic head model in order to confirm the validity of our results. In addition, although sLORETA could result in robust source regions, the source image can be blurred. The use of beamforming techniques could help identify frequency band-specific neural sources more precisely (Hillebrand et al., 2005).
In summary, this study demonstrates that FM patients present enhanced low- and high-frequency oscillatory activity in the brain areas related to cognitive and emotional modulation of pain. Increased low- and high-frequency activity of the prefrontal cortex may contribute to persistent perception of pain in FM. Therapeutic intervention based on manipulating neural oscillation to restore normal thalamocortical rhythmicity may be beneficial to pain relief in FM.
All authors were involved in revising the article critically for important intellectual content. CKC had full access to all of the data in the study and takes responsibility for the integrity of the data and accuracy of the data analysis. ML, JSK, DJK, and CKC: study conception and design, acquisition of data, and analysis and interpretation of data. ML: drafting the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This study was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIP; NRF-2014R1A2A1A11049662) and Basic Science Research Program through the NRF funded by the Ministry of Education (NRF-2014R1A6A3A01009347).
Beck, A. T., Epstein, N., Brown, G., and Steer, R. A. (1988). An inventory for measuring clinical anxiety: psychometric properties. J. Consult. Clin. Psychol. 56, 893–897. doi: 10.1037/0022-006x.56.6.893
Bullmore, E. T., Suckling, J., Overmeyer, S., Rabe-Hesketh, S., Taylor, E., and Brammer, M. J. (1999). Global, voxel and cluster tests, by theory and permutation, for a difference between two groups of structural MR images of the brain. IEEE Trans. Med. Imaging 18, 32–42. doi: 10.1109/42.750253
Choi, W., Lim, M., Kim, J. S., Kim, D. J., and Chung, C. K. (2015). Impaired pre-attentive auditory processing in fibromyalgia: a mismatch negativity (MMN) study. Clin. Neurophysiol. 126, 1310–1318. doi: 10.1016/j.clinph.2014.10.012
Dailey, D. L., Rakel, B. A., Vance, C. G., Liebano, R. E., Amrit, A. S., Bush, H. M., et al. (2013). Transcutaneous electrical nerve stimulation reduces pain, fatigue and hyperalgesia while restoring central inhibition in primary fibromyalgia. Pain 154, 2554–2562. doi: 10.1016/j.pain.2013.07.043
de Tommaso, M., Federici, A., Santostasi, R., Calabrese, R., Vecchio, E., Lapadula, G., et al. (2011). Laser-evoked potentials habituation in fibromyalgia. J. Pain 12, 116–124. doi: 10.1016/j.jpain.2010.06.004
Fallon, N., Chiu, Y. H., Li, X., Nurmikko, T. J., and Stancak, A. (2013). Ipsilateral cortical activation in fibromyalgia patients during brushing correlates with symptom severity. Clin. Neurophysiol. 124, 154–163. doi: 10.1016/j.clinph.2012.06.014
Fallon, N., Li, X., Chiu, Y., Nurmikko, T., and Stancak, A. (2015). Altered cortical processing of observed pain in patients with fibromyalgia syndrome. J. Pain 16, 717–726. doi: 10.1016/j.jpain.2015.04.008
Foerster, B. R., Nascimento, T. D., Deboer, M., Bender, M. A., Rice, I. C., Truong, D. Q., et al. (2015). Excitatory and inhibitory brain metabolites as targets of motor cortex transcranial direct current stimulation therapy and predictors of its efficacy in fibromyalgia. Arthritis Rheumatol. 67, 576–581. doi: 10.1002/art.38945
González-Roldán, A. M., Muñoz, M. A., Cifre, I., Sitges, C., and Montoya, P. (2013). Altered psychophysiological responses to the view of others’ pain and anger faces in fibromyalgia patients. J. Pain 14, 709–719. doi: 10.1016/j.jpain.2013.01.775
Gracely, R. H., Petzke, F., Wolf, J. M., and Clauw, D. J. (2002). Functional magnetic resonance imaging evidence of augmented pain processing in fibromyalgia. Arthritis Rheum. 46, 1333–1343. doi: 10.1002/art.10225
Gross, J., Schnitzler, A., Timmermann, L., and Ploner, M. (2007). Gamma oscillations in human primary somatosensory cortex reflect pain perception. PLoS Biol. 5:e133. doi: 10.1371/journal.pbio.0050133
Gustin, S. M., Wrigley, P. J., Youssef, A. M., Mcindoe, L., Wilcox, S. L., Rae, C. D., et al. (2014). Thalamic activity and biochemical changes in individuals with neuropathic pain after spinal cord injury. Pain 155, 1027–1036. doi: 10.1016/j.pain.2014.02.008
Hämäläinen, M., Hari, R., Ilmoniemi, R. J., Knuutila, J., and Lounasmaa, O. V. (1993). Magnetoencephalography–theory, instrumentation and applications to noninvasive studies of the working human brain. Rev. Mod. Phys. 65, 413–497. doi: 10.1103/revmodphys.65.413
Hauck, M., Lorenz, J., and Engel, A. K. (2007). Attention to painful stimulation enhances gamma-band activity and synchronization in human sensorimotor cortex. J. Neurosci. 27, 9270–9277. doi: 10.1523/jneurosci.2283-07.2007
Hauck, M., Metzner, S., Rohlffs, F., Lorenz, J., and Engel, A. K. (2013). The influence of music and music therapy on pain-induced neuronal oscillations measured by magnetencephalography. Pain 154, 539–547. doi: 10.1016/j.pain.2012.12.016
Henderson, L. A., Peck, C. C., Petersen, E. T., Rae, C. D., Youssef, A. M., Reeves, J. M., et al. (2013). Chronic pain: lost inhibition? J. Neurosci. 33, 7574–7582. doi: 10.1523/JNEUROSCI.0174-13.2013
Hillebrand, A., Singh, K. D., Holliday, I. E., Furlong, P. L., and Barnes, G. R. (2005). A new approach to neuroimaging with magnetoencephalography. Hum. Brain Mapp. 25, 199–211. doi: 10.1002/hbm.20102
Jensen, K. B., Kosek, E., Petzke, F., Carville, S., Fransson, P., Marcus, H., et al. (2009). Evidence of dysfunctional pain inhibition in Fibromyalgia reflected in rACC during provoked pain. Pain 144, 95–100. doi: 10.1016/j.pain.2009.03.018
Jensen, K. B., Kosek, E., Wicksell, R., Kemani, M., Olsson, G., Merle, J. V., et al. (2012a). Cognitive behavioral therapy increases pain-evoked activation of the prefrontal cortex in patients with fibromyalgia. Pain 153, 1495–1503. doi: 10.1016/j.pain.2012.04.010
Jensen, K. B., Loitoile, R., Kosek, E., Petzke, F., Carville, S., Fransson, P., et al. (2012b). Patients with fibromyalgia display less functional connectivity in the brain’s pain inhibitory network. Mol. Pain 8:32. doi: 10.1186/1744-8069-8-32
Jensen, M. P., Sherlin, L. H., Gertz, K. J., Braden, A. L., Kupper, A. E., Gianas, A., et al. (2013). Brain EEG activity correlates of chronic pain in persons with spinal cord injury: clinical implications. Spinal Cord 51, 55–58. doi: 10.1038/sc.2012.84
Kamping, S., Bomba, I. C., Kanske, P., Diesch, E., and Flor, H. (2013). Deficient modulation of pain by a positive emotional context in fibromyalgia patients. Pain 154, 1846–1855. doi: 10.1016/j.pain.2013.06.003
Kim, J. S., and Chung, C. K. (2008). Language lateralization using MEG beta frequency desynchronization during auditory oddball stimulation with one-syllable words. Neuroimage 42, 1499–1507. doi: 10.1016/j.neuroimage.2008.06.001
Kim, J. Y., Kim, S. H., Seo, J., Kim, S. H., Han, S. W., Nam, E. J., et al. (2013). Increased power spectral density in resting-state pain-related brain networks in fibromyalgia. Pain 154, 1792–1797. doi: 10.1016/j.pain.2013.05.040
Kim, D. J., Lim, M., Kim, J. S., Son, K. M., Kim, H. A., and Chung, C. K. (2014). Altered white matter integrity in the corpus callosum in fibromyalgia patients identified by tract-based spatial statistical analysis. Arthritis Rheumatol. 66, 3190–3199. doi: 10.1002/art.38771
Kim, J. S., Shin, K. S., Jung, W. H., Kim, S. N., Kwon, J. S., and Chung, C. K. (2014). Power spectral aspects of the default mode network in schizophrenia: an MEG study. BMC Neurosci. 15:104. doi: 10.1186/1471-2202-15-104
Kim, J. S., Singh, V., Lee, J. K., Lerch, J., Ad-Dab’bagh, Y., Macdonald, D., et al. (2005). Automated 3-D extraction and evaluation of the inner and outer cortical surfaces using a Laplacian map and partial volume effect classification. Neuroimage 27, 210–221. doi: 10.1016/j.neuroimage.2005.03.036
Klein, J. C., Rushworth, M. F., Behrens, T. E., Mackay, C. E., De Crespigny, A. J., D’arceuil, H., et al. (2010). Topography of connections between human prefrontal cortex and mediodorsal thalamus studied with diffusion tractography. Neuroimage 51, 555–564. doi: 10.1016/j.neuroimage.2010.02.062
Kringelbach, M. L., and Rolls, E. T. (2004). The functional neuroanatomy of the human orbitofrontal cortex: evidence from neuroimaging and neuropsychology. Prog. Neurobiol. 72, 341–372. doi: 10.1016/j.pneurobio.2004.03.006
Kwiatek, R., Barnden, L., Tedman, R., Jarrett, R., Chew, J., Rowe, C., et al. (2000). Regional cerebral blood flow in fibromyalgia: single-photon-emission computed tomography evidence of reduction in the pontine tegmentum and thalami. Arthritis Rheum. 43, 2823–2833. doi: 10.1002/1529-0131(200012)43:12<2823::aid-anr24>3.0.co;2-e
Landolt, H. P., Rétey, J. V., Tönz, K., Gottselig, J. M., Khatami, R., Buckelmüller, I., et al. (2004). Caffeine attenuates waking and sleep electroencephalographic markers of sleep homeostasis in humans. Neuropsychopharmacology 29, 1933–1939. doi: 10.1038/sj.npp.1300526
Leblanc, B. W., Lii, T. R., Silverman, A. E., Alleyne, R. T., and Saab, C. Y. (2014). Cortical theta is increased while thalamocortical coherence is decreased in rat models of acute and chronic pain. Pain 155, 773–782. doi: 10.1016/j.pain.2014.01.013
Lim, M., Kim, J. S., Kim, M., and Chung, C. K. (2014). Ascending beta oscillation from finger muscle to sensorimotor cortex contributes to enhanced steady-state isometric contraction in humans. Clin. Neurophysiol. 125, 2036–2045. doi: 10.1016/j.clinph.2014.02.006
Lim, M., Roosink, M., Kim, J. S., Kim, D. J., Kim, H. W., Lee, E. B., et al. (2015). Disinhibition of the primary somatosensory cortex in patients with fibromyalgia. Pain 156, 666–674. doi: 10.1097/j.pain.0000000000000096
Llinás, R. R., Ribary, U., Jeanmonod, D., Kronberg, E., and Mitra, P. P. (1999). Thalamocortical dysrhythmia: a neurological and neuropsychiatric syndrome characterized by magnetoencephalography. Proc. Natl. Acad. Sci. U S A 96, 15222–15227. doi: 10.1073/pnas.96.26.15222
Llinás, R. R., Urbano, F. J., Leznik, E., Ramirez, R. R., and van Marle, H. J. (2005). Rhythmic and dysrhythmic thalamocortical dynamics: GABA systems and the edge effect. Trends Neurosci. 28, 325–333. doi: 10.1016/j.tins.2005.04.006
Loggia, M. L., Berna, C., Kim, J., Cahalan, C. M., Gollub, R. L., Wasan, A. D., et al. (2014). Disrupted brain circuitry for pain-related reward/punishment in fibromyalgia. Arthritis Rheumatol. 66, 203–212. doi: 10.1002/art.38191
López-Solà, M., Pujol, J., Wager, T. D., Garcia-Fontanals, A., Blanco-Hinojo, L., Garcia-Blanco, S., et al. (2014). Altered functional magnetic resonance imaging responses to nonpainful sensory stimulation in fibromyalgia patients. Arthritis Rheumatol. 66, 3200–3209. doi: 10.1002/art.38781
Lutz, J., Jager, L., De Quervain, D., Krauseneck, T., Padberg, F., Wichnalek, M., et al. (2008). White and gray matter abnormalities in the brain of patients with fibromyalgia: a diffusion-tensor and volumetric imaging study. Arthritis Rheum. 58, 3960–3969. doi: 10.1002/art.24070
Mhalla, A., Baudic, S., Ciampi de Andrade, D., Gautron, M., Perrot, S., Teixeira, M. J., et al. (2011). Long-term maintenance of the analgesic effects of transcranial magnetic stimulation in fibromyalgia. Pain 152, 1478–1485. doi: 10.1016/j.pain.2011.01.034
Mountz, J. M., Bradley, L. A., Modell, J. G., Alexander, R. W., Triana-Alexander, M., Aaron, L. A., et al. (1995). Fibromyalgia in women. Abnormalities of regional cerebral blood flow in the thalamus and the caudate nucleus are associated with low pain threshold levels. Arthritis Rheum. 38, 926–938. doi: 10.1002/art.1780380708
Pujol, J., Macià, D., Garcia-Fontanals, A., Blanco-Hinojo, L., López-Solà, M., Garcia-Blanco, S., et al. (2014). The contribution of sensory system functional connectivity reduction to clinical pain in fibromyalgia. Pain 155, 1492–1503. doi: 10.1016/j.pain.2014.04.028
Rolls, E. T., O’Doherty, J., Kringelbach, M. L., Francis, S., Bowtell, R., and McGlone, F. (2003). Representations of pleasant and painful touch in the human orbitofrontal and cingulate cortices. Cereb. Cortex 13, 308–317. doi: 10.1093/cercor/13.3.308
Sarnthein, J., Stern, J., Aufenberg, C., Rousson, V., and Jeanmonod, D. (2006). Increased EEG power and slowed dominant frequency in patients with neurogenic pain. Brain 129, 55–64. doi: 10.1093/brain/awh631
Schmidt-Wilcke, T., Luerding, R., Weigand, T., Jurgens, T., Schuierer, G., Leinisch, E., et al. (2007). Striatal grey matter increase in patients suffering from fibromyalgia—a voxel-based morphometry study. Pain 132, S109–S116. doi: 10.1016/j.pain.2007.05.010
Schulman, J. J., Cancro, R., Lowe, S., Lu, F., Walton, K. D., and Llinás, R. R. (2011). Imaging of thalamocortical dysrhythmia in neuropsychiatry. Front. Hum. Neurosci. 5:69. doi: 10.3389/fnhum.2011.00069
Schulman, J. J., Ramirez, R. R., Zonenshayn, M., Ribary, U. R. S., and Llinás, R. R. (2005). Thalamocortical dysrhythmia syndrome: MEG imaging of neuropathic pain. Thalamus Relat. Syst. 3:33. doi: 10.1017/s1472928805000063
Stenroos, M., Hunold, A., and Haueisen, J. (2014). Comparison of three-shell and simplified volume conductor models in magnetoencephalography. Neuroimage 94, 337–348. doi: 10.1016/j.neuroimage.2014.01.006
Stern, J., Jeanmonod, D., and Sarnthein, J. (2006). Persistent EEG overactivation in the cortical pain matrix of neurogenic pain patients. Neuroimage 31, 721–731. doi: 10.1016/j.neuroimage.2005.12.042
van Diessen, E., Numan, T., Van Dellen, E., Van Der Kooi, A. W., Boersma, M., Hofman, D., et al. (2015). Opportunities and methodological challenges in EEG and MEG resting state functional brain network research. Clin. Neurophysiol. 126, 1468–1481. doi: 10.1016/j.clinph.2014.11.018
Villamar, M. F., Wivatvongvana, P., Patumanond, J., Bikson, M., Truong, D. Q., Datta, A., et al. (2013). Focal modulation of the primary motor cortex in fibromyalgia using 4×1-ring high-definition transcranial direct current stimulation (HD-tDCS): immediate and delayed analgesic effects of cathodal and anodal stimulation. J. Pain 14, 371–383. doi: 10.1016/j.jpain.2012.12.007
Vuckovic, A., Hasan, M. A., Fraser, M., Conway, B. A., Nasseroleslami, B., and Allan, D. B. (2014). Dynamic oscillatory signatures of central neuropathic pain in spinal cord injury. J. Pain 15, 645–655. doi: 10.1016/j.jpain.2014.02.005
Witting, N., Kupers, R. C., Svensson, P., and Jensen, T. S. (2006). A PET activation study of brush-evoked allodynia in patients with nerve injury pain. Pain 120, 145–154. doi: 10.1016/j.pain.2005.10.034
Wolfe, F., Clauw, D. J., Fitzcharles, M. A., Goldenberg, D. L., Katz, R. S., Mease, P., et al. (2010). The American college of rheumatology preliminary diagnostic criteria for fibromyalgia and measurement of symptom severity. Arthritis Care Res. (Hoboken) 62, 600–610. doi: 10.1002/acr.20140
Wolfe, F., Smythe, H. A., Yunus, M. B., Bennett, R. M., Bombardier, C., Goldenberg, D. L., et al. (1990). The American college of rheumatology 1990 criteria for the classification of fibromyalgia. Report of the multicenter criteria committee. Arthritis Rheum. 33, 160–172. doi: 10.1002/art.1780330203
Youssef, A. M., Gustin, S. M., Nash, P. G., Reeves, J. M., Petersen, E. T., Peck, C. C., et al. (2014). Differential brain activity in subjects with painful trigeminal neuropathy and painful temporomandibular disorder. Pain 155, 467–475. doi: 10.1016/j.pain.2013.11.008
Keywords: fibromyalgia, thalamocortical dysrhythmia, prefrontal cortex, gamma oscillation, pain, magnetoencephalography
Citation: Lim M, Kim JS, Kim DJ and Chung CK (2016) Increased Low- and High-Frequency Oscillatory Activity in the Prefrontal Cortex of Fibromyalgia Patients. Front. Hum. Neurosci. 10:111. doi: 10.3389/fnhum.2016.00111
Received: 20 November 2015; Accepted: 29 February 2016;
Published: 14 March 2016.
Edited by:Tetsuo Kida, National Institute for Physiological Sciences, Japan
Reviewed by:Bernadette Van Wijk, University College London, UK
Angela Cristina Valle, University of São Paulo, Brazil
Copyright © 2016 Lim, Kim, Kim and Chung. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution and reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chun Kee Chung, email@example.com
† These authors have contributed equally to this work.