Abstract
Background
Non-invasive painless signaling therapy (NPST) is an electro-cutaneous treatment that converts endogenous pain information into synthetic non-pain information. This study explored whether pain improvement by NPST in failed back surgery syndrome (FBSS) patients is related to cerebral modulation.
Methods
Electroencephalography (EEG) analysis was performed in 11 patients with FBSS. Subjects received daily NPST for 5 days. Before the first treatment, patients completed the Brief Pain Inventory (BPI) and Beck Depression Inventory and underwent baseline EEG. After the final treatment, they responded again to the BPI, reported the percent pain improvement (PPI), and then underwent post-treatment EEG. If the PPI grade was zero, they were assigned to the ineffective group, while all others were assigned to the effective group. We used standardized low-resolution brain electromagnetic tomography (sLORETA) to explore the EEG current-source distribution (CSD) associated with pain improvement by NPST.
Results
The 11 participants had a median age of 67.0 years, and 63.6% were female. The sLORETA images revealed a beta-2 CSD increment in 12 voxels of the right anterior cingulate gyrus (ACG) and the right medial frontal area. The point of maximal CSD changes was in the right ACG. The alpha band CSD increased in 2 voxels of the left transverse gyrus.
According to one study, the global point prevalence of low back pain was 9.4% [1]. In adults, the prevalence of chronic back pain is on the rise [2], and as a result the number of surgeries is also increasing [3]. However, surgery may not relieve the pain. Failed back surgery syndrome (FBSS) is one such outcome [4]. For the treatment of FBSS, a type of chronic pain, the first recommendation is conservative management. Surgical consideration is the final stage of the treatment algorithm, but it does not always resolve the pain satisfactorily and involves risks [5,6]. Among the modalities of conservative management, drug therapy is the most accessible. However, more than half of patients with chronic pain combine multiple drugs [7]. Such medication polytherapy increases the risk of side effects. Non-invasive painless signaling therapy (NPST) is a treatment modality that provides non-pain information to the cutaneous nerves using patient-specific electrostimulation to control the cerebral perception of pain information. NPST has the advantage of being a minimally invasive non-pharmacological treatment [8-10].
The brain interprets the pain stimuli transmitted from the peripheral nervous system and integrates this interpretation with other stored information. When perceiving pain, the brain employs a process that includes selective attention and intensity coding [11-16]. The pain perceived by the brain may have a different intensity than pain at the peripheral level; pain perception may even occur in the brain in response to non-painful peripheral stimuli [17-20]. The pain matrix is the cerebral network involved in pain perception [20,21]. The somatosensory cortex, insula, and anterior cingulate gyrus (ACG) participate in the pain matrix. When the pain matrix is activated, the brain senses and focuses on sensory input from the peripheral system. The brain then integrates the stimuli with the multimodal somatic representation in the cerebral cortex and interprets the integrated information to convert it into a conscious pain perception [20].
Electrical source imaging (ESI) shows the cerebral activity through electroencephalography (EEG) current-source analysis on magnetic resonance imaging (MRI) matched with a human head model [22-24]. To assess the neurophysiological state and treatment response in cases of chronic pain, the generally accepted test modality has been functional MRI (fMRI). fMRI measures neural brain activity indirectly using the hemodynamic response [25-28]. ESI and fMRI complement each other [29]. Although its spatial resolution is inferior to that of fMRI, ESI has the advantage of mapping cerebral neuronal activity directly and at high temporal resolution [22-24,30,31]. Research using ESI to investigate cerebral electrophysiology associated with chronic pain remains rare.
The effects of NPST might be associated with functional modulation of cerebral structures involving pain perception, such as the pain matrix. Our goal was to identify the functional changes in the brain associated with pain improvement by NPST in FBSS patients. For that purpose, using EEG current-source analysis as a frequency domain analysis, we observed the current-source distribution (CSD) difference between FBSS patients with pain improvement by NPST and those with no therapeutic effect.
The Institutional Review Board (IRB) of Gyeongsang National University Hospital approved the present study; the IRB number was 2018-08-016-001. We also registered the study with the World Health Organization International Clinical Trials Registry Platform (WHO ICTRP) through the Clinical Research Information Service, which enables the use of the WHO ICTRP in Korea; the registration number was KCT0003347.
The Department of Rehabilitation Medicine at Gyeongsang University Hospital recruited patients diagnosed with FBSS from November 2018 to June 2019. Patient age was limited to 30-70 years. To avoid difficulties related to cooperation with treatment and evaluation, patients with intellectual disabilities or severe medical, neurological, or psychiatric disorders were excluded. Patients with various conditions that may affect pain sensitivity were also excluded, as were those with a history of ischemic heart disease within 6 months, mesenteric nerve block within 4 weeks, pregnancy or lactation, scars on the back, or chronic use of antiepileptic drugs. Finally, because depressive symptoms can affect pain and brain activity, patients whose degree of depressive symptoms as evaluated by the Beck Depressive Inventory (BDI) was moderate or severe were excluded.
After agreeing to participate, enrolled patients visited the hospital every weekday, Monday through Friday, for one week. On Monday, they completed the Brief Pain Inventory (BPI) to determine the baseline pain level and the BDI to assess depressive symptoms. Then, we recorded the patients’ first baseline EEGs (time point 1, T1). Next, on the same day, patients received the first session of NPST. They were treated with NPST for 5 days in a row, Monday through Friday. Immediately after the fifth and last therapy session on Friday, they completed the follow-up BPI and reported the percent pain improvement (PPI). After the questionnaires were completed, we recorded the second, post-treatment EEGs (time point 2, T2) for comparison with the baseline T1 EEG.
We then divided the patients into two groups: effective and ineffective groups. The effective group included patients for whom NPST resulted in reduced pain, and the ineffective group included those for whom there was no therapeutic effect. The therapeutic effect was evaluated based on the PPI value reported just after the last (Friday) session of NPST. When evaluating the PPI in terms of grades, a grade of 1 or higher was regarded effective, and a grade of 0 as ineffective.
The BDI is a self-reported inventory and a useful tool for detecting depressive symptoms. It requires that subjects rate their condition during the last week using 21 items describing depressive symptoms. It is a Likert scale-based inventory where each item is given a score from 0 to 3. The total BDI score combines the 21 items, yielding an overall score of 0-63. Based on the total score, the degrees of depressive symptoms are described as follows: 0-16, normal; 17-20, mild; 21-23, moderate; and 24-63, severe [32-34].
The BPI was used to evaluate pain severity [35-37]. It is a Likert scale-based inventory with 11 items that evaluate pain during the last 24 hours. Each of the 11 items is given a score ranging from 0 to 10. Four items address pain severity, and the other 7 evaluate interference to living functions. In this study, the 4 BPI items addressing pain severity (PSIs-BPI) were used to assess pain severity at two time points: T1 and T2. The BPI guides patients to report their pain intensity on a scale from 0 to 10. A score of 0 indicates “no pain,” and a score of 10 means “pain as bad as you can imagine.” Total scores on the four-item PSIs-BPI could range from 0-40.
The NPST was administered using a Pain Jammer ENS-1140 device (ENS System, Bucheon, Korea). The device uses alternating current stimulation, and the maximum current available for treatment was 5.5 mA. The treatment was administered to the area in the low back and buttock region where the patient complained of the most severe pain. The therapist placed four stimulation electrodes along the lines of pain and numbness at the designated treatment location.
After determining the positions of the stimulation electrodes for NPTS, at those spots, we discovered the maximum tolerance intensity of the electric stimulation. Stimulation was started at zero intensity and was then gradually raised until a patient felt a sting or burning sensation. At that point, we added more stimuli with that intensity. If the patient perceived the stimuli as a tolerable tingling sensation by adaptation, we slightly increased the stimulation intensity again. When the patient could not tolerate any more for a certain level of intensity, we slightly reduced the intensity and found the maximum tolerance intensity. NPST was performed using the maximum tolerance intensity, and each treatment session lasted 50 minutes.
During EEG recording, patients closed their eyes and remained awake in a calm state. Technicians recorded the EEGs using two digital EEG machines; the recording time was 50 minutes. Two machines were used: a Comet (Grass Telefactor, West Warwick, RI) or a Beehive Horizon (Grass Telefactor). For EEG recording, 25 electrodes were placed on the scalp referencing the International 10-20 system. The positions of the electrodes were as follows: Fp1, F3, C3, P3, O1, Fp2, F4, C4, P4, O2, F7, T7, P7, F8, T8, P8, F9, P9, T9, F10, P10, T10, Fz, Cz, and Pz. For the EEG recording, the filter setting was 1-70 Hz with a notch filter of 60 Hz, and the sampling rate was 400 Hz.
We extracted the appropriate EEG segments from the recorded EEGs. Extraction of the EEG segments involved use of the brain electrical source analysis (BESA, v.7.0; MEGIS, Grafelfing, Germany) data review and processing module. The EEG segment duration was 5 seconds, and it was extracted from EEG records with occipital alpha activity and without artifacts. Each patient underwent EEG twice (T1 and T2), and each EEG provided five segments, yielding 10 EEG segments per patient. Among these, 5 were EEG segments at T1, and the other 5 at T2.
The EEG segments selected through preprocessing were used as data to identify CSD associated with pain improvement by NPST. The current-source analysis proceeded in two steps. The first step involved subtracting the EEG at T1 from the EEG at T2 to obtain the CSD changes induced by NPST in the effective and ineffective groups. The second step entailed subtracting the CSD changes induced by NPST in the ineffective group from those in the effective group to obtain the CSD associated with pain improvement (CSD of effective group [T2 – T1] – CSD of ineffective group [T2 – T1]).
We used standardized low-resolution brain electromagnetic tomography (sLORETA; Key Institute for Brain-Mind Research, Zurich, Switzerland) software to calculate and describe the cerebral CSD. The sLORETA expresses EEG CSD in the brain using a distributed source model. It finds CSD in a three-shell spherical head model registered to the Talairach human brain atlas, converts the CSD found into Montreal Neurological Institute (MNI) coordinates, and represents the transformed coordinates in a digitized MRI from the Brain Imaging Center, MNI. Thus, the sLORETA images obtained through EEG analysis using sLORETA are three-dimensional functional images. The images consisted of 6,430 cortical gray voxels, and their spatial resolution was 5 mm [38].
By performing the current-source analysis, we obtained ESIs in the frequency domain. The seven frequency domains obtained were as follows; 1-3 Hz, delta; 4-7 Hz, theta; 8-12 Hz, alpha; 3-18 Hz, beta-1; 19-21 Hz, beta-2; and 22-30 Hz, beta-3. Statistical significance set at P < 0.05 (two tailed) expressed the suprathreshold voxels in the sLORETA images.
We performed univariate analyses using Mann–Whitney tests for continuous variables and Fisher’s exact tests for categorical variables to compare the clinical features between the effective and ineffective groups. R software (version 3.5.3; R Foundation for Statistical Computing, Vienna, Austria) was used for these statistical analyses [39].
During the study period, 11 of 45 FBSS patients were enrolled in this study (Fig. 1, Table 1). The 11 enrolled participants had a median age of 67.0 years, and 7 (63.6%) were female. The enrolled FBSS patients had undergone spinal surgery between one and three times. Nine patients had undergone operation twice, and the remaining two patients had undergone surgery three times and once, respectively. Three spinal diseases were the targets of the operations. These include herniated lumbar disc in nine patients, spinal stenosis in seven, and spondylolisthesis in one (Table 1).
Seven of the 11 were classified as belonging to the effective group, and the remaining 4 made up the ineffective group. Table 2 presents a comparison of patients in the effective group with those in the ineffective group. Significant differences were observed between the T1–T2 PSIs-BPI scores and in PPI, but no other differences were found (Table 2). The median PPI was 40.0% and 7.5% in the effective and ineffective groups, respectively (P < 0.001). In terms of the grade of pain improvement, three patients in the effective group had grade 1 improvement, two had grade 2, and the remaining two had grades 3 and 4. All four patients in the ineffective group had grade 0 (P = 0.030). The median improvement in PSIs-BPI score (T1 – T2) was 6.0 points in the effective group and 0.4 points in the ineffective group (P < 0.001).
In the sLORETA images using a two-tailed test with a significance level of 0.05, the CSD of the beta-2 and alpha band significantly increased. The point of maximal CSD change was in the right ACG (Brodmann’s area [BA] 32, MNI coordinates: x = 10, y = 25, z = 25) (Table 3, Fig. 2). No significant CSD changes occurred in the other frequency domains.
The beta-2 CSD increment was observed in 12 suprathreshold voxels of the right hemisphere located in BA 32, the medial part of BA 24, and BA 9 (Table 3, Fig. 2). The three BA locations were the following: BA 32, the dorsolateral ACG; BA 24, the ventral ACG; and BA9, the medial- and dorsolateral frontal area [40]. Of the 12 suprathreshold voxels, 9 were in the ACG, and 3 in the medial frontal gyrus. The CSD of the alpha band increased in BA 41 in the left hemisphere, with two suprathreshold voxels (Table 3, Fig. 3). BA 41 is a transverse gyrus corresponding to the posterosuperior part of the superior temporal gyrus [40].
This study used the ineffective group as controls. The EEG of each patient was recorded twice, at T1 and T2. The EEG analysis design for sLORETA image acquisition was the following: CSD of the effective group (T2 – T1) – CSD of the ineffective group (T2 – T1). The CSD difference between the two groups was the primary finding of this study, so the changes in CSD were considered to indicate cerebral activity linked to pain reduction by NPST. As elements of that observation, the CSD of beta-2 and alpha significantly increased in certain cerebral regions. Compared to the alpha band, the CSD change was substantial and broad in the beta-2 band. The CSD increment locations of the beta-2 band were the right ACG gyrus and the right medial frontal area. The CSD increment of the alpha band occurred in the left transverse gyrus. Among the three areas with increased activity, the broadest and most marked increase was in the right ACG. Because the ACG is a component of the pain matrix, the pain matrix could be functionally associated with the therapeutic NPST effect.
The ACG participates in two functional pain perception networks: one of these detects pain information, and the other perceives pain intensity. When these two networks work together as an integrative functional map, people experience subjective pain [11-16]. Because of this pain-perceiving network involving the ACG, people sense subjective pain rather than the actual intensity of a noxious stimulus [11].
In the present study, increased activity in the right ACG was associated with pain reduction by NPST. Therefore, it seems that increasing ACG activity may cure chronic pain. It might be possible to correctly code the pain intensity in patients with chronic pain if NPST increases ACG activity. A previous study reported that daily life disruption due to pain declined when the dorsal ACG was stimulated with deep brain stimulation in patients with thalamic pain syndrome [13].
Functional representations associated with pain in the medial frontal cortex are well known. For the representations, pain is localized in the anterior midcingulate cortex, negative emotion in the ventromedial prefrontal cortex, and cognitive control in the dorsal midcingulate cortex. Furthermore, these three representations have a meaningful correlation in the dorsomedial frontal cortex [41]. In the present study, we observed an increase in the beta-2 CSD in three suprathreshold voxels of the right medial frontal gyrus. The increase in activity may be significant, considering the functional representations in the medial frontal cortex. Our findings indicate that increased activity was related to pain reduction via NPST.
A pathophysiological alteration in the superior temporal gyrus may occur in patients with postherpetic neuralgia (PHN) [42]. Previous research has revealed that regional homogeneity of the brain increased in PHN patients with effective pain control in some areas, including the right superior temporal gyrus, indicating increased local functional connectivity [43,44]. We also observed a CSD change in the superior temporal gyrus, although it was in the opposite hemisphere, and this change was associated with pain reduction in FBSS. Chronic pain may associate with the abnormality of the superior temporal gyrus in both the right and left hemispheres. It seems that pain reduction may result if the areas obtain functional improvement.
The sLORETA model is a distributed ESI model and is free from localization errors compared to the previous version, the LORETA, but it still has spatial resolution limitations [38]. Various factors, such as the number and distribution of electrodes, head models, and EEG-MRI co-registration, affect ESI accuracy [22,45]. The localization accuracy of ESI increases as the number of scalp electrodes increases, but accuracy does not continue to improve at more than 89 electrodes. Uniform sampling of the full head surface is required for source localization to work correctly. The accuracy of source localization is comparable between 19 and 46 electrodes when obtaining a LORETA image based on the three-shell spherical head model if the electrodes are distributed across the entire head [22].
The locations of scalp electrodes must be correct to increase the accuracy of source localization in ESI. The individual EEG-MRI co-registration makes it possible to overcome interindividual variation in electrode positions [22]. Selection of head models is another crucial aspect of EEG-MRI co-registration. The spherical model with uniform conductivity is the simplest head model. It is easy to apply and allows rapid calculation, but the human head is not spherical and is not spatially uniform. Realistic head models increase ESI accuracy but take longer to compute [22,45].
ESI serves as an adjunctive modality that finds the epileptogenic zone [46]. Precise localization is required when performing ESI for presurgical assessment, and the ESI used in the present study (25 electrodes and based on a 3-shell spherical head model) would be insufficient for this task. However, our goal was to investigate cerebral structures associated with NPST-induced pain relief in FBSS patients. The resolution of the sLORETA imaging obtained in this study may be sufficient to discriminate among the brain structures. In any case, the findings of this study help clarify the neurophysiological mechanisms of FBSS pain.
This study involved some limitations. The first weakness was the small number of enrolled patients. However, the findings are of value because we observed direct neural activity using ESI. Few researchers have observed ESI associated with the therapeutic effect of NPST in FBSS. The second limitation was the short-term follow-up, which made it difficult to know how long the observed CSD changes might persist. More research is needed with larger populations of patients with FBSS to explore whether our results can be reproduced, and long-term follow-up research is needed to further evaluate the long-term effects of NPST. Another drawback is that this is an observational study rather than an experimental one, so that it is unable to examine the causal relationship between pain reduction and CSD changes. The interpretation that there is an association between them, on the other hand, maybe appropriate. The present authors excluded patients taking antiepileptic drugs chronically from this study. This exclusion could be a drawback if this study aimed to observe the effectiveness of NPST in FBSS patients. Rather than that, the purpose of this study was to research cerebral functional changes associated with NPST-induced pain improvement in FBSS patients.
Certain brain structures, such as the pain matrix, serve to detect and perceive pain. The brain also has the plasticity to adapt to alterations of afferent peripheral stimulus information. NPST modifies such afferent information. In this study, the pain reduction induced by NPST in patients with FBSS was associated with increased cerebral activity, mainly in the right ACG. The ACG is a cerebral area that is included in the pain matrix. Thus, NPST may generate therapeutic effects by affecting the pain matrix.
Notes
REFERENCES
1. Hoy D, March L, Brooks P, Blyth F, Woolf A, Bain C, et al. 2014; The global burden of low back pain: estimates from the Global Burden of Disease 2010 study. Ann Rheum Dis. 73:968–74. DOI: 10.1136/annrheumdis-2013-204428. PMID: 24665116.
2. Smith M, Davis MA, Stano M, Whedon JM. 2013; Aging baby boomers and the rising cost of chronic back pain: secular trend analysis of longitudinal Medical Expenditures Panel Survey data for years 2000 to 2007. J Manipulative Physiol Ther. 36:2–11. DOI: 10.1016/j.jmpt.2012.12.001. PMID: 23380209. PMCID: PMC3596008.
3. Rajaee SS, Bae HW, Kanim LE, Delamarter RB. 2012; Spinal fusion in the United States: analysis of trends from 1998 to 2008. Spine (Phila Pa 1976). 37:67–76. DOI: 10.1097/BRS.0b013e31820cccfb. PMID: 21311399.
4. Harvey AM. 1995; Classification of chronic pain-descriptions of chronic pain syndromes and definitions of pain terms. Clin J Pain. 11:163. https://journals.lww.com/clinicalpain/Citation/1995/06000/Classification_of_Chronic_Pain_Descriptions_of.24.aspx. DOI: 10.1097/00002508-199506000-00024. PMID: 3461421.
5. Diebo BG, Passias PG, Marascalchi BJ, Jalai CM, Worley NJ, Errico TJ, et al. 2015; Primary versus revision surgery in the setting of adult spinal deformity: a nationwide study on 10,912 patients. Spine (Phila Pa 1976). 40:1674–80. DOI: 10.1097/BRS.0000000000001114. PMID: 26267823.
6. Baber Z, Erdek MA. 2016; Failed back surgery syndrome: current perspectives. J Pain Res. 9:979–87. DOI: 10.2147/JPR.S92776. PMID: 27853391. PMCID: PMC5106227.
7. Freynhagen R, Bennett MI. 2009; Diagnosis and management of neuropathic pain. BMJ. 339:b3002. DOI: 10.1136/bmj.b3002. PMID: 19675082.
8. Ricci M, Pirotti S, Scarpi E, Burgio M, Maltoni M, Sansoni E, et al. 2012; Managing chronic pain: results from an open-label study using MC5-A Calmare®device. Support Care Cancer. 20:405–12. DOI: 10.1007/s00520-011-1128-6. PMID: 21394458.
9. Smith TJ, Coyne PJ, Parker GL, Dodson P, Ramakrishnan V. 2010; Pilot trial of a patient-specific cutaneous electrostimulation device (MC5-A Calmare®) for chemotherapy-induced peripheral neuropathy. J Pain Symptom Manage. 40:883–91. DOI: 10.1016/j.jpainsymman.2010.03.022. PMID: 20813492. PMCID: PMC4383258.
10. Shannon CE. 1948; A mathematical theory of communication. Bell Syst Tech J. 27:379–423. https://ieeexplore.ieee.org/abstract/document/6773024. DOI: 10.1002/j.1538-7305.1948.tb01338.x.
11. Craig AD, Reiman EM, Evans A, Bushnell MC. 1996; Functional imaging of an illusion of pain. Nature. 384:258–60. DOI: 10.1038/384258a0. PMID: 8918874.
12. Foltz EL, White LE Jr. 1962; Pain "relief" by frontal cingulumotomy. J Neurosurg. 19:89–100. DOI: 10.3171/jns.1962.19.2.0089. PMID: 13893868.
13. Levi V, Cordella R, D'Ammando A, Tringali G, Dones I, Messina G, et al. 2019; Dorsal anterior cingulate cortex (ACC) deep brain stimulation (DBS): a promising surgical option for the treatment of refractory thalamic pain syndrome (TPS). Acta Neurochir (Wien). 161:1579–88. DOI: 10.1007/s00701-019-03975-5. PMID: 31209628.
14. Peyron R, García-Larrea L, Grégoire MC, Costes N, Convers P, Lavenne F, et al. 1999; Haemodynamic brain responses to acute pain in humans: sensory and attentional networks. Brain. 122(Pt 9):1765–80. DOI: 10.1093/brain/122.9.1765. PMID: 10468515.
15. Sikes RW, Vogt BA. 1992; Nociceptive neurons in area 24 of rabbit cingulate cortex. J Neurophysiol. 68:1720–32. DOI: 10.1152/jn.1992.68.5.1720. PMID: 1479441.
16. Vaccarino AL, Melzack R. 1989; Analgesia produced by injection of lidocaine into the anterior cingulum bundle of the rat. Pain. 39:213–9. DOI: 10.1016/0304-3959(89)90008-0. PMID: 2594399.
17. Baumgärtner U, Tiede W, Treede RD, Craig AD. 2006; Laser-evoked potentials are graded and somatotopically organized anteroposteriorly in the operculoinsular cortex of anesthetized monkeys. J Neurophysiol. 96:2802–8. DOI: 10.1152/jn.00512.2006. PMID: 16899640.
18. Hofbauer RK, Fiset P, Plourde G, Backman SB, Bushnell MC. 2004; Dose-dependent effects of propofol on the central processing of thermal pain. Anesthesiology. 100:386–94. DOI: 10.1097/00000542-200402000-00031. PMID: 14739816.
19. Lee MC, Mouraux A, Iannetti GD. 2009; Characterizing the cortical activity through which pain emerges from nociception. J Neurosci. 29:7909–16. DOI: 10.1523/JNEUROSCI.0014-09.2009. PMID: 19535602. PMCID: PMC6665631.
20. Legrain V, Iannetti GD, Plaghki L, Mouraux A. 2011; The pain matrix reloaded: a salience detection system for the body. Prog Neurobiol. 93:111–24. DOI: 10.1016/j.pneurobio.2010.10.005. PMID: 21040755.
21. Peyron R, Laurent B, García-Larrea L. 2000; Functional imaging of brain responses to pain. A review and meta-analysis (2000). Neurophysiol Clin. 30:263–88. DOI: 10.1016/S0987-7053(00)00227-6. PMID: 11126640.
22. Michel CM, Murray MM, Lantz G, Gonzalez S, Spinelli L, Grave de Peralta R. 2004; EEG source imaging. Clin Neurophysiol. 115:2195–222. DOI: 10.1016/j.clinph.2004.06.001. PMID: 15351361.
23. He B, Sohrabpour A, Brown E, Liu Z. 2018; Electrophysiological source imaging: a noninvasive window to brain dynamics. Annu Rev Biomed Eng. 20:171–96. DOI: 10.1146/annurev-bioeng-062117-120853. PMID: 29494213. PMCID: PMC7941524.
24. Gross J, Kujala J, Hämäläinen M, Timmermann L, Schnitzler A, Salmelin R. 2001; Dynamic imaging of coherent sources: studying neural interactions in the human brain. Proc Natl Acad Sci U S A. 98:694–9. DOI: 10.1073/pnas.98.2.694. PMID: 11209067. PMCID: PMC14650.
25. Davis KD, Flor H, Greely HT, Iannetti GD, Mackey S, Ploner M, et al. 2017; Brain imaging tests for chronic pain: medical, legal and ethical issues and recommendations. Nat Rev Neurol. 13:624–38. DOI: 10.1038/nrneurol.2017.122. PMID: 28884750.
26. Tracey I, Johns E. 2010; The pain matrix: reloaded or reborn as we image tonic pain using arterial spin labelling. Pain. 148:359–60. DOI: 10.1016/j.pain.2009.11.009. PMID: 20080346.
27. Bélanger M, Allaman I, Magistretti PJ. 2011; Brain energy metabolism: focus on astrocyte-neuron metabolic cooperation. Cell Metab. 14:724–38. DOI: 10.1016/j.cmet.2011.08.016. PMID: 22152301.
28. Logothetis NK. 2003; The underpinnings of the BOLD functional magnetic resonance imaging signal. J Neurosci. 23:3963–71. DOI: 10.1523/JNEUROSCI.23-10-03963.2003. PMID: 12764080. PMCID: PMC6741096.
29. Nguyen T, Potter T, Karmonik C, Grossman R, Zhang Y. 2018; Concurrent EEG and functional MRI recording and integration analysis for dynamic cortical activity imaging. J Vis Exp. 136:56417. DOI: 10.3791/56417. PMID: 30010646. PMCID: PMC6102018.
30. Seeber M, Cantonas LM, Hoevels M, Sesia T, Visser-Vandewalle V, Michel CM. 2019; Subcortical electrophysiological activity is detectable with high-density EEG source imaging. Nat Commun. 10:753. DOI: 10.1038/s41467-019-08725-w. PMID: 30765707. PMCID: PMC6376013.
31. Ramon C, Schimpf PH, Haueisen J. 2006; Influence of head models on EEG simulations and inverse source localizations. Biomed Eng Online. 5:10. DOI: 10.1186/1475-925X-5-10. PMID: 16466570. PMCID: PMC1389789.
32. Beck AT, Ward CH, Mendelson M, Mock J, Erbaugh J. 1961; An inventory for measuring depression. Arch Gen Psychiatry. 4:561–71. DOI: 10.1001/archpsyc.1961.01710120031004. PMID: 13688369.
33. Kwon OY, Park SP. 2013; Frequency of affective symptoms and their psychosocial impact in Korean people with epilepsy: a survey at two tertiary care hospitals. Epilepsy Behav. 26:51–6. DOI: 10.1016/j.yebeh.2012.10.020. PMID: 23207517.
34. Rhee MK, Lee YH, Park SH, Sohn CH, Chung YC, Hong SK, et al. 1995; A standardization study of Beck Depression Inventory 1 - Korean version (K-BDI): reliability and factor analysis. Korean J Psychopathol. 4:77–95. http://www.riss.kr/search/detail/DetailView.do?p_mat_type=1a0202e37d52c72d&control_no=6860a27c5da18618ffe0bdc3ef48d419.
35. Cleeland CS, Ladinsky JL, Serlin RC, Nugyen CT. 1988; Multidimensional measurement of cancer pain: comparisons of US and Vietnamese patients. J Pain Symptom Manage. 3:23–7. DOI: 10.1016/0885-3924(88)90134-0. PMID: 3351345.
36. Cleeland C. 1991; Research in cancer pain. What we know and what we need to know. Cancer. 67(3 Suppl):823–7. DOI: 10.1002/1097-0142(19910201)67:3+<823::AID-CNCR2820671412>3.0.CO;2-S. PMID: 1986852.
37. Yun YH, Mendoza TR, Heo DS, Yoo T, Heo BY, Park HA, et al. 2004; Development of a cancer pain assessment tool in Korea: a validation study of a Korean version of the brief pain inventory. Oncology. 66:439–44. DOI: 10.1159/000079497. PMID: 15452372.
38. Pascual-Marqui RD. 2002; Standardized low-resolution brain electromagnetic tomography (sLORETA): technical details. Methods Find Exp Clin Pharmacol. 24 Suppl D:5–12. PMID: 12575463.
39. R Core Team. 2020. R: a language and environment for statistical computing [Internet]. R Foundation for Statistical Computing;Vienna: Available at: https://www.R-project.org/.
40. Simić G, Hof PR. 2015; In search of the definitive Brodmann's map of cortical areas in human. J Comp Neurol. 523:5–14. DOI: 10.1002/cne.23636. PMID: 24889330.
41. Kragel PA, Kano M, Van Oudenhove L, Ly HG, Dupont P, Rubio A, et al. 2018; Generalizable representations of pain, cognitive control, and negative emotion in medial frontal cortex. Nat Neurosci. 21:283–9. DOI: 10.1038/s41593-017-0051-7. PMID: 29292378. PMCID: PMC5801068.
42. Zhang Y, Yu T, Qin B, Li Y, Song G, Yu B. 2016; Microstructural abnormalities in gray matter of patients with postherpetic neuralgia: a diffusional kurtosis imaging study. Pain Physician. 19:E601–11. DOI: 10.36076/ppj/2019.19.E601. PMID: 27228526.
43. Jiang L, Zuo XN. 2016; Regional homogeneity: a multimodal, multiscale neuroimaging marker of the human connectome. Neuroscientist. 22:486–505. DOI: 10.1177/1073858415595004. PMID: 26170004. PMCID: PMC5021216.
44. Zhang Y, Cao S, Yuan J, Song G, Yu T, Liang X. 2020; Functional and structural changes in postherpetic neuralgia brain before and six months after pain relieving. J Pain Res. 13:909–18. DOI: 10.2147/JPR.S246745. PMID: 32440196. PMCID: PMC7210030.
45. Koles ZJ. 1998; Trends in EEG source localization. Electroencephalogr Clin Neurophysiol. 106:127–37. DOI: 10.1016/S0013-4694(97)00115-6. PMID: 9741773.
46. Sharma P, Seeck M, Beniczky S. 2019; Accuracy of interictal and ictal electric and magnetic source imaging: a systematic review and meta-analysis. Front Neurol. 10:1250. DOI: 10.3389/fneur.2019.01250. PMID: 31849817. PMCID: PMC6901665.