Korean J Pain 2023; 36(1): 113-127
Published online January 1, 2023 https://doi.org/10.3344/kjp.22225
Copyright © The Korean Pain Society.
Matheus Soldatelli1,2 , Álvaro de Oliveira Franco2 , Felipe Picon2,3,4 , Juliana Ávila Duarte1,5 , Ricardo Scherer2 , Janete Bandeira2 , Maxciel Zortea2 , Iraci Lucena da Silva Torres6 , Felipe Fregni7 , Wolnei Caumo1,2,7,8
1Post-Graduate Program in Medical Sciences, School of Medicine, Universidade Federal do Rio Grande do Sul (UFRGS), Porto Alegre, Brazil, 2Laboratory of Pain and Neuromodulation at Hospital de Clínicas de Porto Alegre (HCPA), Porto Alegre, Brazil, 3Department of Psychiatry, Faculdade de Medicina, UFRGS, Porto Alegre, Brazil, 4ADHD Outpatient Program, HCPA, Porto Alegre, Brazil, 5Department of Internal Medicine, UFRGS, Porto Alegre, Brazil, 6Laboratory of Pharmacology in Pain and Neuromodulation: Pre-clinical Investigations, Experimental Research Center, HCPA, Porto Alegre, Brazil, 7Pain and Palliative Care Service, HCPA, Porto Alegre, Brazil, 8Laboratory of Neuromodulation and Center for Clinical Research Learning, Physics and Rehabilitation Department, Spaulding Rehabilitation Hospital, Boston, MA, USA
Correspondence to:Wolnei Caumo
Laboratory of Pain and Neuromodulation at Hospital de Clínicas de Porto Alegre (HCPA), Rua Ramiro Barcelos, 2350 - Rio Branco, Porto Alegre 90035-003, Brazil
Tel: +55-51-3359-8083, Fax: +55-51-3359-8083, E-mail: wcaumo@hcpa.edu.br
*These authors contributed equally to this work.
Handling Editor: Young Hoon Kim
Received: July 4, 2022; Revised: November 9, 2022; Accepted: November 9, 2022
This is an open-access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/), which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
Background: Resting-state functional connectivity (rs-FC) may aid in understanding the link between pain-modulating brain regions and the descending pain modulatory system (DPMS) in fibromyalgia (FM). This study investigated whether the differences in rs-FC of the primary somatosensory cortex in responders and non-responders to the conditioned pain modulation test (CPM-test) are related to pain, sleep quality, central sensitization, and the impact of FM on quality of life.
Methods: This cross-sectional study included 33 females with FM. rs-FC was assessed by functional magnetic resonance imaging. Change in the numerical pain scale during the CPM-test assessed the DPMS function. Subjects were classified either as non-responders (i.e., DPMS dysfunction, n = 13) or responders (n = 20) to CPM-test. A generalized linear model (GLM) and a receiver operating characteristic (ROC) curve analysis were performed to check the accuracy of the rs-FC to differentiate each group.
Results: Non-responders showed a decreased rs-FC between the left somatosensory cortex (S1) and the periaqueductal gray (PAG) (P < 0.001). The GLM analysis revealed that the S1-PAG rs-FC in the left-brain hemisphere was positively correlated with a central sensitization symptom and negatively correlated with sleep quality and pain scores. ROC curve analysis showed that left S1-PAG rs-FC offers a sensitivity and specificity of 85% or higher (area under the curve, 0.78, 95% confidence interval, 0.63–0.94) to discriminate who does/does not respond to the CPM-test.
Conclusions: These results support using the rs-FC patterns in the left S1-PAG as a marker for predicting CPM-test response, which may aid in treatment individualization in FM patients.
Keywords: Central Nervous System Sensitization, Chronic Pain, Fibromyalgia, Functional Neuroimaging, Magnetic Resonance Imaging, Neural Pathways, Pain Perception, Periaqueductal Gray, Psychophysics, Somatosensory Cortex.
Fibromyalgia (FM) is a chronic pain condition that comprises a widespread pain index (WPI) and emotional distress, fatigue, sleep disorders, depression, catastrophizing thinking related to pain, and cognitive compromise. Although its pathophysiology remains unclear, several studies have found altered activations of several brain regions, such as the thalamic nuclei, somatosensory cortices, anterior cingulate, and the insular and prefrontal cortices [1–3]. A review of structural and functional rearrangements in chronic pain indicates more active neural networks in pain, including the thalamus, anterior cingulate cortex (ACC), primary and secondary somatosensory cortices (S1 and S2, respectively), and the insular cortex [4].
According to functional magnetic resonance imaging (fMRI) data, the ventrolateral-periaqueductal gray matter (PAG) has indirect connections with central lateral and medial pathways, the ACC, and upper pons/medulla [5]. The PAG relates to somatosensory and affective-attentional pain components [2,6], with downward projections to spinal dorsal horn neurons to modulate pain transmission information [7]. These pathways constitute the descending pain modulatory system (DPMS), whose dysfunction is frequent in chronic pain syndromes [8].
The conditioned pain modulation test (CPM-test) assesses the function of the DPMS based on the "pain suppresses pain" phenomenon [9]. The CPM-test activates a cortical spinal-bulb-spinal loop responsible for diffuse noxious inhibitory control. The stimulation parameters, test sites, and study population influence the reliability of the CPM-test and conditioning stimuli [10]. A systematic review found fifty percent of intersession reliability measured by intraclass correlation coefficient (ICC) rated as good (ICC 0.60 to 0.75) or excellent (ICC > 0.75) [10]. The CPM-test permits the evaluation of the DPMS function. Its dysfunction in FM has been correlated with the serum biomarkers of neuroplasticity, such as the brain-derived neurotrophic factor and S100-B protein [11]. Furthermore, the dysfunction of DPMS is related to hyperinhibition at the cortical level and increased scores in the Central Sensitization Inventory (CSI) [12].
Resting-state functional magnetic resonance imaging (rs-fMRI) allows the evaluation of the resting-state functional connectivity (rs-FC) between different brain regions. This is a practical approach when investigating the existence of a specific brain’s functional signature and coherence patterns associated with specific patient groups across chronic pain conditions [13]. Regions with synchronous activity tend to correlate when blood oxygen level-dependent (BOLD) activity is measured. However, studies have found no information on the intrinsic FC of pain processing areas with the severity of dysfunction of the DPMS.
FM is a disease associated with central sensitization syndrome (CSS), likely related to a maladaptive function of the neural networks involved in pain processing. However, a gap in the literature persists regarding the understanding about the connection between pain processing brain networks and the modulation of the pain by the DPMS. Thus, the authors aimed to answer the following questions: (i) to examine whether there are differences in the rs-FC pattern between the primary somatosensory cortex and pain processing areas in patients who do/do not respond to the CPM-test, and (ii) to evaluate whether the differences in primary somatosensory cortex rs-FC patterns in patients who do/do not respond to the CPM-test are related to pain, sleep quality, central sensitization, and the impact of FM on quality of life. It was hypothesized that the connection between the networks involved in the pain neuromatrix and the DPMS might be a suitable marker for identifying patients with more severe clinical symptoms.
The Institutional Review Board (IRB, CAAE 2018-0353) at the Hospital de Clínicas de Porto Alegre (HCPA) approved the protocol of this cross-sectional study, which was conducted according to the Declaration of Helsinki. Participants provided verbal and written informed consent before participating and did not receive payment in exchange for their participation. Recruitment was undertaken from January 2018 to December 2019.
We included 33 adult females, all right-handed, ages 18 to 65, who could read and write. Patients from the Pain Outpatient Clinic of the HCPA were invited to the study. Other participants became aware of the study through local newspaper publicity. All included subjects met the diagnostic criteria of FM by the standard assessment protocol criteria for the diagnosis of FM—according to the American College of Rheumatology (ACR) 2016 criteria [14]—which was applied by senior physicians with more than fifteen years of experience in chronic pain care. Furthermore, they had to demonstrate daily disability for routine activities due to FM for the three months preceding enrollment. In addition, they needed to report a score of at least six on the numerical pain scale (NPS, 0–10).
We did not include people who had a history of rheumatoid arthritis, lupus, or any other autoimmune, neurologic, or oncological disease, any uncompensated clinical disease (
A prior sample size was estimated, expecting a large effect size (ES, f squared, f2 0.35) for multiple regression analysis, allowing for five predictors. For error types I and II, which had values respectively of 0.05 and 0.20, the estimation showed that the sample size should be 31 patients (Post-hoc Statistical Power Calculator for Hierarchical Multiple Regression: https://www.danielsoper.com/statcalc/calculator.aspx?id=17) [15]. To guarantee the study’s power, it was decided to include 33 patients. However, the outcome variable (
All psychological tests used in this study have been validated for the target population. The following instruments were used to assess psychological symptoms and sleep quality: the Central Sensitization Inventory [12], the Beck Depression Inventory-II [16], the Pain Catastrophizing Scale [17], the State-Trait Anxiety Inventory [18], and the Pittsburgh Sleep Quality Index [19].
The dependent variables were the rs-FC in the left S1 and the PAG. The main interest factor was the DPMS function, which was assessed through a CPM-test—expressed by the difference in score on an NPS (0–10) produced by the Quantitative Sensory Testing (QST) before and concurrently with the CPM-test induced by immersion of the dominant hand in cold water (zero to 1°C). The CPM-test permits the classification of non-responders and responders. Independent variables included age, years of education, pain catastrophizing, pain-related disability, quality of sleep, pain scores, anxiety, number of psychiatric diagnoses, use of analgesics, and heat pain threshold. The timeline of assessments is shown in Fig. 1.
Structural and functional images were acquired with a 3T scanner (Ingenia 3.0T; Phillips, Best, Netherlands) using a 15-channel head coil. The rs-fMRI was a single-shot T2*-weighted fast-field echo, echo-planar imaging sequence (repetition time [TR] = 2,000 ms, repetition time [TE] = 30 ms, Matrix = 80 × 80, field of view [FOV] = 240 mm, flip angle = 90, 3 × 3 × 3.5 mm voxel size, 36 slices in ascending order with 0.35 mm gap) with 300 volumes, totaling 10 minutes. Anatomical references were acquired using a T1-weighted gradient-echo sequence (TR/TE = 8.5/3.9 sec, Matrix/FOV = 240/240 mm, flip angle = 8°, 0.94 mm isotropic voxel size, 200 sagittal slices).
All rs-fMRI processing was done using CONN18 (www.nitrc.org/projects/conn) [20], which is a toolbox that uses SPM12 (www.fil.ion.ucl.ac.uk/spm) commands over MATLAB (MathWorks, Natick, MA) [21]. The pipeline starts with preprocessing steps (segmenting anatomical volumes in gray matter, white matter, and cerebrospinal fluid, realigning and unwarping, normalizing, and smoothing the functional volumes). CONN’s default denoising pipeline combines two general steps: linear regression of potential confounding effects in the BOLD signal and temporal band-pass filtering. The BOLD signal was extracted from each region of interest (ROI) based on the Harvard-Oxford Atlas predefined within CONN standard settings. Other ROIs were added based on the relevant areas related to the pain connectivity literature (described in more detail below). Factors identified as potential confounding effects on the estimated BOLD signal were estimated and removed separately for each voxel, subject, and functional run/session using Ordinary Least Squares regression to project each BOLD signal time series to the subspace orthogonal to all potential confounding effects. The temporal preprocessing started with a regression out of the BOLD signal to control motion artifacts and residual physiological noise, then controlling the confounding variables (first-level covariates, the BOLD signal from white matter, cerebrospinal fluid), realignment parameters, scrubbing, and the CompCor approach [22]. Twelve potential noise components were defined from the estimated subject-motion parameters to minimize motion-related BOLD variability, including three translation parameters and three rotation parameters plus their associated first-order derivatives. A variable number of noise components (one for each identified outlier scan during the outlier identification preprocessing step) were used as potential confounding effects to remove any influence of these outlier scans on the BOLD signal. Temporal frequencies below 0.008 Hz or above 0.09 Hz were removed from the BOLD signal to focus on slow-frequency fluctuations while minimizing the influence of physiological noise, noise from head motion, and other noise sources [20].
Then, the BOLD signal from each individual was filtered (0.008–0.09 Hz). Filtering was implemented using a discrete cosine transform windowing operation to minimize border effects, having been performed after regression to avoid any frequency mismatch in the nuisance regression procedure. Each rs-FC correlation measure was calculated from all ROIs concerning each other in an FC matrix. Finally, after extracting the average time course from each ROI, rs-FC was estimated from Fisher Z scores.
We applied ART (Artifact detection/identification toolbox) within the CONN [20,22] environment to control for possible spurious correlations brought by head movement artifacts. A conservative threshold was chosen within CONN that sets the subject-motion threshold at 0.5 mm.
Spherical ROIs with MNI (Montreal Neurological Institute) coordinates [23] were chosen as predefined ROIs based on the Harvard-Oxford Atlas using the CONN-fMRI Functional connectivity toolbox [24]. The radius used was 6 mm for the PAG and 10 mm for the remaining areas. ROIs were acquired from prior studies as seeds. Together, they were used to create rs-FC maps of specific networks related to chronic pain and FM for each subject, with rs-FC being measured using the following source seeds (see Table 1 for coordinates): PAG [6], left (L) and right (R) primary and secondary somatosensory cortex (L_S1, R_S1, L_S2, R_S2) [25], left and right primary motor cortex (L_M1, R_M1) [26], left and right dorsolateral and ventromedial prefrontal cortexes (L_DLPFC, R_DPLFC, L_MPFC, R_MPFC) [27,28], left and right anterior and posterior insular cortexes (L_aIC, L_pIC, R_aIC, R_pIC), left and right ACC (L_ACC, R_ACC) [25], left and right ventrolateral and mediodorsal thalamus (L_VLTh, R_VLTh, L_MDTh, R_MDTh) [29,30], left and right hippocampus (L_Hippo, R_Hippo) [31], amygdalae (L_Amyg, R_Amyg) [32], and nucleus accumbens (L_NA, R_NA) [33]. After calculating the subjects’ FC map, the second-level comparison was performed through an ROI-to-ROI analysis. The threshold for bidirectional explorations of FC was
Table 1 Results from functional connectivity analysis (Z score threshold > 3.96 and seed-level-corrected-FDR
Region | Zmax | p-unc | p-FDR | |||
---|---|---|---|---|---|---|
PAG | 0.0049 | 0.1417 | 0 | –32 | –10 | |
L_S1 | –4.46 | 0.0001 | 0.0401 | –50 | –28 | 44 |
R_S1 | 0.94 | 0.4968 | 0.8475 | 48 | –28 | 46 |
L_S2 | 1.74 | 0.1451 | 0.4002 | –64 | –26 | 18 |
R_S2 | 1.89 | 0.1139 | 0.4002 | 50 | –26 | 16 |
L_M1 | 2.00 | 0.0948 | 0.4002 | –32 | –16 | 58 |
R_M1 | 1.71 | 0.1518 | 0.4002 | 34 | –12 | 74 |
L_DLPFC | 2.83 | 0.0256 | 0.2471 | –39 | 29 | 26 |
R_DLPFC | 0.21 | 0.9789 | 0.9789 | 42 | 35 | 26 |
L_MPFC | 0.61 | 0.7406 | 0.9338 | –12 | 60 | –8 |
R_MPFC | 0.71 | 0.6603 | 0.9244 | 22 | 58 | –10 |
L_aIC | 1.62 | 0.1750 | 0.4230 | –32 | 16 | 6 |
L_pIC | 1.02 | 0.4445 | 0.8475 | –39 | –15 | 1 |
R_aIC | 3.26 | 0.0133 | 0.1933 | 32 | 16 | 6 |
R_pIC | 0.54 | 0.7989 | 0.9510 | 39 | –15 | 8 |
L_aACC | 1.10 | 0.3910 | 0.8099 | –2 | 14 | 30 |
R_aACC | 2.09 | 0.0823 | 0.4002 | 2 | 12 | 32 |
L_VLTh | 0.98 | 0.4701 | 0.8475 | –13 | –17 | 1 |
R_VLTh | 1.83 | 0.1247 | 0.4002 | 14 | –17 | 1 |
L_MDTh | 0.66 | 0.7013 | 0.9244 | –11 | –22 | 3 |
R_MDTh | 0.83 | 0.5726 | 0.9225 | 9 | –19 | 3 |
L_Hippo | 0.74 | 0.6427 | 0.9244 | –22 | –16 | –20 |
R_Hippo | 1.96 | 0.1020 | 0.4002 | 22 | –20 | –16 |
L_Amyg | 0.51 | 0.8198 | 0.9510 | –32 | 2 | –16 |
R_Amyg | 1.38 | 0.2570 | 0.5733 | 32 | 2 | 16 |
L_NA | 0.68 | 0.6858 | 0.9244 | –8 | 12 | 1 |
R_NA | 0.22 | 0.9761 | 0.9789 | 11 | 11 | 1 |
ROI: region of interest, FDR: false discovery rate, PAG: periaqueductal gray, p-unc: uncorrected
1) For the FM diagnosis, the ACR [14] criteria were applied by a physician. The scale for FM symptoms ranges from 0 to 3 and is composed of the WPI and the modified severity scale (SS scale) [14].
2) The
3) The
4) The
5) The
6) The
7) The
1) The
2) The
1) Demographic data and medical comorbidities were evaluated using a standardized questionnaire. Psychiatric diagnoses were assessed by the Mini-International Neuropsychiatric Interview (MINI) [37]. Weekly analgesic intake in the last three months was recorded. They were dichotomized as < 4 times per week or > 4 times per week.
2) The visual analog scale (VAS) is a visual scale for pain score assessment based on millimeters, which ranges from "no pain" (zero) to “the worst possible pain" (100 mm). Patients had to score their worst pain during the last 24 hours.
Subjects were chosen based on pre-set criteria evaluated by two researchers with clinical experience in outpatient treatment for chronic pain. Two evaluators specifically trained for performing all assessments applied standardized protocols (
Categorical and continuous variables were summarized using conventional descriptive statistics. Continuous variables were compared with
Table 2 Demographic, clinical characteristics, psychological, and pain measures related to the pain of the study sample
Variable | Non-responders (n = 13) | Responders (n = 20) | |
---|---|---|---|
Age (yr) | 52.11 ± 6.67 | 48.00 ± 7.98 | 0.980 |
Education (yr) | 11.44 ± 4.09 | 12.55 ± 5.29 | 0.643 |
Time of diagnosis (yr) | 6.78 ± 0.44 | 6.95 ± 0.22 | 0.090 |
Smoking (Yes/No) | 2/11 | 7/13 | 0.200 |
Clinical comorbidity (Yes/No) | 10 (76.9)/3 (23.1) | 14 (70.0)/6 (30.0) | 0.493 |
Hypertension (Yes) | 9 | 4 | |
Cardiac ischemia (Yes) | 2 | 0 | |
Diabetes (Yes) | 0 | 2 | |
Hypothyroidism (Yes) | 5 | 2 | |
Asthma (Yes) | 0 | 1 | |
Other (Yes) | 2 | 0 | |
Number of days analgesic used per week in the last three months (< 4 times/> 4 times)b | 5 (38.5)/8 (61.5) | 10 (50.0)/10 (50.0) | 0.511 |
Acetaminophen/Dipirone (Yes) | 9 | 8 | |
Dorflex (Yes) orphenadrine citrate + dipirone | 7 | 6 | |
NSAID (Yes) | 4 | 3 | |
Opioid (Yes) | 4 | 3 | |
Psychiatric disorder according to MINI (Yes/No)a | 11 (84.6)/2 (15.4) | 10 (50.0)/10 (50.0) | 0.041 |
Anxiety | 4 | 5 | |
Major depression | 10 | 12 | |
Panic Syndrome | 1 | 4 | |
Obsessive-compulsive disorder | 0 | 1 | |
Use of central nervous system medication (Yes/No)b | 11 (84.6)/2 (15.4) | 10 (50.0)/10 (50.0) | 0.040 |
Tricyclic or dual inhibition of monoamine reuptake antidepressants | 6 | 7 | |
Selective serotonin reuptake inhibitor antidepressant | 5 | 4 | |
Anticonvulsant drugs | 3 | 4 | |
Benzodiazepine | 2 | 2 | |
Pain measures and severity of symptoms | |||
Visual Analogue Scale (0–10 cm) | 8.63 ± 0.97 | 8.25 ± 1.25 | 0.206 |
Widespread pain index | 13.83 ± 1.99 | 12.63 ± 1.92 | 0.655 |
Symptom severity scale | 9.58 ± 1.83 | 10.05 ± 1.62 | 0.790 |
American College Rheumatology 2016 | 22.22 ± 3.86 | 22.60 ± 2.42 | 0.812 |
Heat pain threshold – Celsius degree | 44.60 ± 3.33 | 45.78 ± 3.25 | 0.381 |
Central sensitization symptoms, sleep quality, and psychological symptoms | |||
Central Sensitization Inventory | 61.23 ± 15.54 | 64.30 ± 14.39 | 0.566 |
Pittsburgh Sleep Quality Index | 12.22 ± 3.88 | 13.05 ± 3.69 | 0.145 |
Trait-Anxiety Inventory (reduced version) | 33.00 ± 5.70 | 27.15 ± 9.93 | 0.110 |
State-Anxiety Inventory (reduced version) | 28.67 ± 5.93 | 25.75 ± 6.43 | 0.265 |
Beck Depression Inventory | 21.33 ± 9.53 | 23.55 ± 12.04 | 0.808 |
Pain Catastrophizing Scale | 38.61 ± 7.17 | 30.45 ± 11.35 | 0.028 |
Values are presented as mean ± standard deviation, number only, or number (%).
For data analysis, analgesic use was included as a dichotomous variable (analgesics less than four days per week or more than four days per week). For patients with chronic pain, the rescue analgesic use changes each week depending on their pain level.
NSAID: nonsteroidal anti-inflammatory drug, MINI: mini-international Neuropsychiatric Interview.
aPatients could have none or more than one psychiatric disorder.
bSome patients were using more than one type of drug.
Table 3 Correlation among pain catastrophizing, depressive symptoms, central sensitization score, sleep quality, and disability due to pain (n = 33)
Variables | Responders to CPM-test (n = 20) | Non-responders to CPM-test (n = 13) | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
(1) | (2) | (3) | (4) | (5) | (1) | (2) | (3) | (4) | (5) | ||
Fibromyalgia Impact Questionnaire (1) | 1 | 1 | |||||||||
Pain Catastrophizing Scale (2) | –0.07 | 1 | 0.27 | 1 | |||||||
Beck Depression Inventory (3) | 0.03 | 0.64** | 1 | 0.25 | 0.49 | 1 | |||||
Pittsburgh Sleep Quality Index (4) | 0.21 | 0.58** | 0.59** | 1 | 0.09 | –0.13 | 0.27 | 1 | |||
Central Sensitization Inventory (5) | –0.12 | 0.49* | 0.43 | 0.54* | 1 | 0.74** | 0.46 | 0.49 | 0.20 | 1 | |
Left S1-PAG rs-FC (6) | –0.13 | –0.11 | –0.12 | –0.07 | 0.37 | –0.40 | 0.06 | –0.03 | –0.32 | –0.08 |
CPM-test: conditioned pain modulation test, rs-FC: resting-state functional connectivity, S1: primary somatosensory cortex, PAG: periaqueductal gray.
*
The cutoff values with the highest Youden index, sensitivity and specificity values were plotted and visually inspected for a plausible cutoff point with a higher value. Multiple comparisons were adjusted using the Bonferroni test, and two-tailed tests were used. A type I error of 5% was accepted. The authors used the software IBM SPSS Statistics for Windows, version 22.0 (IBM Corp., Armonk, NY) for the statistical analysis.
We included 33 females with FM out of 50 screened candidates. The clinical and sociodemographic descriptions of the final sample are given in detail in Table 2. Six patients did not fulfill the diagnostic criteria for FM or presented pain levels on most days for the last three months. Five had another clinical diagnosis defined as an exclusion criterion. In addition, six had contraindications for fMRI scanning. Compared to responder subjects (n = 20), non-responders (n = 13) to the CPM-test had a higher prevalence of psychiatric disorders according to the MINI, and higher use of central nervous system medication. Besides, they showed higher levels of pain catastrophizing.
The results from rs-FC analysis (Z score threshold > 3.96 and seed-level-corrected-FDR
The ROI-to-ROI rs-FC analysis comparing non-responders to responders is shown in Fig. 2. Non-responders showed a decreased rs-FC between the left S1 and the PAG, as shown by the blue line, compared to the responders. This finding was significant, using a threshold of
Spearman correlation analysis was used to look at the relationship between the left S1-PAG rs-FC (the outcome) between groups of responders and non-responders to the CPM-test with the following covariates: impact of FM symptoms on quality of life, pain catastrophizing, depressive symptoms, central sensitization score, and sleep quality. Non-responders to the CPM-test showed a moderately positive correlation between the left S1-PAG rs-FC and pain catastrophizing, central sensitization score, and quality of life scores. Data are presented in Table 3.
The left S1-PAG rs-FC marginal means adjusted by GLM and standard error in responders
The S1-PAG rs-FC in the left-brain hemisphere was negatively correlated with clinical symptoms, such as sleep quality, disability due to pain, and the severity of pain. In contrast, the severity of CS was positively correlated with the left S1-PAG rs-FC. This result indicates that increased FC was associated with more severe symptoms of CSS. The adjusted ES according to the CMP-test group on the left S1-PAG rs-FC, considering the final GLM, with all variables, was large (the χ2 = 10.14; ES = 0.56). In contrast, the ES of clinical symptoms on the S1-PAG rs-FC was moderate (Table 4).
Table 4 Generalized linear model analysis to assess the relationship of the descending pain modulatory system function according to responders (n = 20) and non-responders (n = 13) to the conditioned pain modulation test (CPM-test) on the clinical symptoms and left S1-PAG rs-FC (n = 33)
Outcome: Left S1-PAG rs-FC | B | SE | 95% CI | Wald χ2 | df | Effect size | |
---|---|---|---|---|---|---|---|
(Intercept) | 0.335 | 0.1816 | (–0.021 to 0.690) | 3.394 | 1 | 0.065 | |
Responders to CPM-test | 0.145 | 0.0457 | (0.056 to 0.240) | 10.109 | 1 | 0.001 | 0.58 |
Non-responders to CPM-test | reference | ||||||
Pittsburgh Sleep Quality Index | –0.014 | 0.0058 | (–0.026 to –0.003) | 6.162 | 1 | 0.013 | 0.45 |
Fibromyalgia Impact Questionnaire | –0.002 | 0.0010 | (–0.004 to 0.001) | 4.417 | 1 | 0.036 | 0.38 |
Pain severity on Numerical Pain Scale (0–10) | –0.011 | 0.0052 | (–0.021 to 0.001) | 4.395 | 1 | 0.036 | 0.38 |
Central Sensitization Inventory | 0.003 | 0.0012 | (0.001 to 0.005) | 5.737 | 1 | 0.017 | 0.44 |
Primary outcome – generalized linear model analysis to compare responders and non-responders.
The Cramerʼs V was used as a measure of effect size for chi-square tests. The size effect was interpreted as follows: Standards for interpreting Cramerʼs V proposed by Cohen (1988) are the following. DF (degrees of freedom) = 1 (0.10 = small effect) (0.30 = medium effect) (0.50 = large effect). https://www.campbellcollaboration.org/escalc/html/EffectSizeCalculator-R5.php.
CPM-test: conditioned pain modulation test, rs-FC: resting-state functional connectivity, S1: primary somatosensory cortex, PAG: periaqueductal gray, B: regression coefficient, SE: standard error, CI: confidence interval, χ2: Wald chi-square, df: degrees of freedom.
ROC analysis applying the Youden index extension showed that non-responders could be distinguished from responders using the cutoff point on the left S1-PAG rs-FC set to –0.24, offering a sensitivity of 95% and a specificity of 85% or higher (AUC 0.78, 95% CI, 0.63 to 0.94). The non-responders had a negative mean value in the left S1-PAG FC. The ROC curves are displayed in Fig. 3.
This study revealed that the rs-FC between the left S1 and the PAG is significantly related to the dysfunction of the DPMS. These results show that the left S1-PAG rs-FC was negatively linked to a lower quality of life, worse sleep quality, and more severe pain in women with FM. Conversely, it was positively associated with central sensitization. Furthermore, they revealed that left S1-PAG rs-FC might be able to distinguish patients with a failure of the DPMS according to groups of responders and non-responders to the CPM-test. To the best of the authors’ knowledge, this is the first study that suggests a dysfunction in the cortical area involved in the sensorial-discriminative component of pain (
Although the study design hinders us from causally explaining these associations, it is plausible that they could be a compensatory response to the persistent hyperexcitability related to chronic pain adaptation. Hence, this could result in the dysregulation of cortical function and its connection with PAG. This hypothesis supports the PAG neurons’ physiological processes, which play a critical role in autonomic, motivated behavior, cortical motor, and perception networks. Besides, it is the primary control center for descending pain modulation [38]. Therefore, the increase in connectivity might indicate "pain sensitization" rather than "pain intensity". Given the earlier studies’ results, the increased rs-FC in the left S1-PAG is possibly related to CSS [11].
Our findings are aligned with previous studies in FM that found higher rs-FC between the insula and interconnected networks of the brain’s default mode network (DMN), including MPFC, posterior cingulate cortex (PCC), precuneus, inferior parietal lobule, hippocampal formation, and lateral temporal cortex [13,39]. Likewise, an increase in connectivity between the insula and the DMN has been documented in several pain conditions, including chronic low back pain, osteoarthritis, and complex regional pain syndrome [13]. In contrast, a previous study found a decreased rs-FC of the PAG in areas associated with motor/executive, DMN, and premotor cortex in FM patients compared to healthy controls [6].
Our results, therefore, add to the evidence about increased rs-FC between brain areas involved in pain processing [40,41]. Its significance highlights that the left S1-PAG rs-FC is positively correlated with the symptoms of central sensitization. This implies that increased FC between these brain regions may point to decreased function in these neural networks. Based on the weakening of cortical inhibition, the authors think the shift in the left S1-PAG may signify an up-regulation phenomenon of the intracortical inhibitory networks, as evidenced by an increase in short intracortical inhibition [42]. Even though the exact mechanism of central sensitization is not fully understood, changes in transmission and weak synapses have been suggested as possible causes of the symptoms of CSS. So, the rs-FC might help to comprehend how unfavorable interactions among cortical areas associated with sensory discrimination, motivation, emotion, motor, attention, arousal, and response selection might arise with symptoms of CSS [43]. One important thing to remember about this approach is that the correlation analysis used to measure how different brain areas work together is limited because it cannot disclose the type of relationship, causality, or direction of causation. According to earlier studies, the dynamic processes of rs-FC can change in rest, compared to stimulus, in the way FM displays a substantial imbalance in the connectivity within the pain network during rest [39], as well as the descending pain inhibition might have a ceiling effect that lessens thalamic activity [41]. Considering this, it is possible that the left S1-PAG rs-FC results from an increased persistent excitatory input of pain signals between important antinociceptive locations, such as the brainstem. While several cortical regions have been proposed as possible sources of faulty pain inhibition in FM, the CSS scores reveal abnormal processes in the PFC and motor cortex [44].
Also, it is known that the connectivity measurements of the forebrain regions may only be indirectly related to these brain areas. The relationship between rs-FC and the various target areas involved in pain processing and in the severity of clinical symptoms is an active area of research. More longitudinal studies are needed to understand these complex relationships fully. In this situation, connectivity measures must be seen as proxies; yet this does not limit the use of connectivity metrics to identify FM patients with more severe clinical symptoms. Overall, these results show that the efficiency of DPMS and cortical dysfunction are related, and they also show the importance of looking at the rs-FC as a possible indicator of disease severity and at its relationship with clinical symptoms (
This study found, in the discriminative analysis using ROC analysis to screen non-responders to the CPM-test, that the cutoff value of 0.24 gave a good profile in terms of sensitivity, specificity, and an AUC of 0.78 for finding people with more DPMS dysfunction. These results suggest that the abnormal rs-FC in these neural networks may be a marker of the dysfunction of the neurobiological systems underpinning FM symptoms. From a theoretical point of view, these results suggest that a maladaptive function of these areas is part of the cascade of events connected to the inefficiency of the descending endogenous pain-modulation system. At the same time, the leftward asymmetry in FC with the left S1-PAG neurons could reflect differences in integration and specialization in pain processing, or be a severity marker of dysfunction in the cortical processing of pain. According to a previous study, the left amygdala lateralization in pain processing is responsible for the main component of the affective-emotional pain pathway related to the DPMS [45]. Thus, the asymmetry in the rs-FC left S1-PAG found in the current study may explain this result. Many networks are strongly lateralized, such as language and visuals in the left and right hemispheres; the sensorimotor networks’ laterality is as well present and progressively being described [46]. The authors are conscious that more results are needed to allow some conclusion in this regard, but the current results open a new avenue to investigate lateralization of pain processing areas.
We found a moderate correlation between pain catastrophizing and the left S1-PAG rs-FC in univariate analysis. However, this correlation lost statistical significance in the multivariate analysis (Table 4). According to the literature, multiple cortical areas’ FC are involved with pain catastrophizing thinking, such as the S1, S2, DLPFC, mPFC, and ACC, the anterior insula, thalamus, and DMN (mPFC- PCC/precuneus) [47]. In light of these mixed results, further studies using a similar paradigm with a large sample size are needed to clarify these differences. Data suggests that pain catastrophizing may contribute to individual differences and susceptibility to subsequent chronic pain disorders. According to an earlier study, cognitive-behavioral therapy to reduce pain catastrophizing reduced FC between S1 and the anterior/medial insular cortex [48].
The negative correlation of quality-of-life scores, sleep quality, and severity of pain with the left S1-PAG rs-FC provides input to comprehend the relationship between the dysfunction of cortical processing in target areas involved in pain processing and its connection with a site with a significant role in the DPMS. It is conceivable that the increased left S1-PAG rs-FC reflects an adaptation to counter-regulate cortical hyperactivation due to DMPS inefficiency. This idea finds some ground in a previous study on FM with fMRI that showed more significant activity than controls in response to nonpainful warm stimuli in the prefrontal, supplemental motor, insular, and anterior cingulate cortices [49]. Also, it agreed with previous studies that the severity of clinical symptoms, including pain, is linked to the dysfunction of the DPMS [11]. So, it aligned with emerging evidence suggesting that FM pain is mediated by central nervous system hyperexcitability rather than only peripheral pathology.
Some points were addressed concerning the study design. First, although the literature supports selecting ROIs to study FC, they might not correspond precisely to the broad anatomical spectrum of pain-related areas. Hence, this limitation should be considered in interpreting current findings, especially regarding PAG, due to its small representation area. In addition, to improve the precision of the measure, a ten-minute-long resting-state scan was performed to overcome technical issues related to data acquisition. In contrast, most previous studies used a 5- or 7-min scan. Second, there was no control group because the authors’ main interest was comparing FM groups according to their DPMS dysfunction. Third, psychiatric disorders are potentially confounding in cortical processing, as well as medications such as antidepressants, analgesics, mood stabilizers, and antipsychotics. Fourth, only females were included due to the higher prevalence of FM in women, which is related to sex differences in pain processing and in the functioning of the DPMS [50]. Fifth, it was observed that the cutoff point of –0.24 showed the best equilibrium between sensitivity and specificity to identify non-responders. It indicates that the values are quite spread out among non-responders to the CPM-test. In the context of diagnosis, it appears that it identifies subjects with a more severe reduction in the FC indexed to the left S1-PAG measures. Finally, FM exhibits substantial variability in the somatic and cognitive symptoms. Thus, this heterogeneity translates into fMRI brain phenomenology, limiting the precision and specificity of making generalized interpretations for other pain conditions.
The rs-FC patterns of the left S1-PAG may help to identify those patients with more severe inhibitory dysfunction of the DPMS. In addition, they support the hypothesis that the DPMS and the somatosensory cortex are involved in maladaptive neuroplasticity processes. Overall, they shed light on how brain function relates to the pathophysiology of primary chronic pain. So, they open up new ways to customize top-down therapy approaches to improve the effectiveness of the DPMS and cortical areas that play a key role in pain processing.
The datasets supporting the findings of this study are available from the corresponding author upon reasonable request.
No potential conflict of interest relevant to this article was reported.
The following Brazilian funding agencies supported the present research:
1. Committee for the Development of Higher Education Personnel—CAPES (PROEX). To Wolnei Caumo.
2. National Council for Scientific and Technological Development - CNPq (grant to WC number: 420826/2018-1; 302688/2017-0).
3. Postgraduate Research Group at the Hospital de Clínicas de Porto Alegre, to Wolnei Caumo (FIPE, project number 2018-0353).
4. Foundation for the Support of Research at the Rio Grande do Sul (FAPERGS)Array of state of Rio Grande do Sul, Brazil (SEARS) n. 03/2017 (PPSUS), to Wolnei Caumo. number: 17/2551-0001.
5. Brazilian Innovation Agency (FINEP [Financiadora de Estudos e Projetos]); process number 1245/13 to Iraci LS Torres and Wolnei Caumo).
Wolnei Caumo agrees to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.
Matheus Soldatelli: Writing/manuscript preparation; Álvaro de Oliveira Franco: Writing/manuscript preparation; Felipe Picon: Computation; Juliana Ávila Duarte: Investigation; Ricardo Scherer: Investigation; Janete Bandeira: Investigation; Maxciel Zortea: Investigation; Iraci Lucena da Silva Torres: Writing/manuscript preparation; Felipe Fregni: Writing/manuscript preparation; Wolnei Caumo: Supervision.