Document Type : Original Research

Authors

1 Department of Radiology, School of Para-Medicine, Hormozgan University of Medical Sciences (HUMS), Bandar Abbas, Iran

2 Department of Medical Physics and Biomedical Engineering, School of Medicine, Tehran University of Medical Sciences (TUMS), Tehran, Iran

3 PERFORM Preventive Medicine and Health Care Centre, Concordia University, Montreal, Quebec, Canada

4 Department of Sports Medicine, School of Medicine, Tehran University of Medical Sciences (TUMS), Tehran, Iran

5 Spine Center of Excellence, Yas Hospital, Tehran University of Medical Sciences, Tehran, Iran

6 Department of Neurosurgery, Shariati hospital, Tehran University of Medical Sciences, Tehran, Iran

7 Department of Epidemiology and Biostatistics, School of Public Health, Tehran University of Medical Sciences (TUMS), Tehran, Iran

8 Department of Radiology, School of Medicine, Tehran University of Medical Sciences (TUMS), Tehran, Iran

10.31661/jbpe.v0i0.2204-1481

Abstract

Background: Non-specific chronic low back pain (CLBP) is a common painful condition and is responsible for different physical disorders. Despite alternative therapies, patients still suffer from persistent pain. Repetitive transcranial magnetic stimulation (rTMS) has provided much evidence of pain reduction, but results have not been examined deeply in CLBP symptoms.
Objective: The analgesic effect of rTMS in non-specific CLBP patients was evaluated by the amplitude of low-frequency fluctuation (ALFF) analysis in resting-state fMRI.
Material and Methods: In this experimental study, fifteen non-specific CLBP participants (46.87±10.89 years) received 20 Hz rTMS over the motor cortex. The pain intensity and brain functional scan were obtained during pre and post-stimulation for all participants. The ALFF maps of the brain in two scan sessions were identified and the percentage of pain reduction (PPR%) was determined using paired t-test. Also, correlation analysis was used to find a relationship between ALFFs and pain intensity.
Results: Pain intensity was significantly reduced after induced-rTMS in non-specific CLBP (36.22%±13.28, P<0.05). Positive correlation was found between ALFF in the insula (INS) and pain intensity (rpre-rTMS=0.59, rpost-rTMS=0.58) while ALFF in medial prefrontal cortex (mPFC) and pain intensity had negatively correlated (rpre-rTMS=-0.54, rpost-rTMS=-0.56) (P<0.05). ALFF increased in mPFC while INS, thalamus (THA), and supplementary motor area (SMA) showed decremental ALFF followed by rTMS. 
Conclusion: This study demonstrated that ALFF in INS, THA, mPFC, and SMA is associated with CLBP symptoms and analgesic effects of rTMS. ALFF potentially seems to be a proper objective neuroimaging parameter to link spontaneous brain activity with pain intensity in non-specific CLBP patients.

Highlights

Mahboubeh Masoumbeigi (Google Scholar)

Nader Riyahi Alam (Google Scholar)

Keywords

Introduction

Chronic low back pain (CLBP) with no clear clinical evidence which is known as non-specific is a highly prevalent pain condition and is a leading cause of temporary or permanent physical and even psychological disorders worldwide [ 1 ]. Despite alternative therapies such as pharmacological prescriptions, exercise programs, bed rest, acupuncture, spinal manipulation, or even surgery and injection, patients with CLBP still suffer from persistence or recurrence of pain and movement problems as serious difficulties in daily life [ 2 - 5 ]. In this context, high frequency repetitive transcranial magnetic stimulation (rTMS) as an innovative, safe, and noninvasive alternative technique has been recently proposed for possible use in chronic pain treatment [ 6 ]. In the TMS method, an electric current is rapidly changed and flows through an electromagnetic coil that has been used to produce short magnetic pulses that pass easily and painlessly through the brain. These pulses induce an electric field that results in neuron depolarization and cortical stimulation that has been applied to modulate altered brain activities to alleviate pain [ 7 , 8 ].

Different studies showed that a single session of high-frequency rTMS over the primary motor cortex (M1) results in pain relief for several days (up to a week) after stimulation that the optimal effect is delayed by 2-4 days [ 9 - 12 ]. Nevertheless, the neural mechanisms that contribute to the pain reduction effects of rTMS for patients with pain are still not completely understood. However, the modulation of neuronal plasticity and cortical excitability after induced high-frequency rTMS has been recently observed in functional neuroimaging studies [ 13 , 14 ]. In this regard, the development of imaging techniques such as resting-state functional magnetic resonance imaging (rs-fMRI) can be used as a noninvasive method to investigate neural mechanisms involved with rTMS in humans [ 11 ].

Recently, accumulating evidence suggests that low-frequency (0.01-0.08 Hz) spontaneous blood oxygen level-dependent (BOLD) fluctuations in resting-state fMRI (rs-fMRI) provide physiological insights into neural spontaneous activity [ 15 , 16 ]. The amplitude of low-frequency fluctuations (ALFFs) as a trustworthy method identifies the voxel-level spontaneous activity in the brain region. ALFF can measure the amplitude values of low-frequency oscillations (LFOs) of rs-fMRI time courses in the cortical regions and focus on the neural processes and activity in different brain regions rather than time-series similarities among regions [ 17 , 18 ]. Importantly, disrupted LFOs amplitude has been reported in various clinical disorders including attention-deficit hyperactivity disorder (ADHD), obsessive-compulsive disorder (OCD), Alzheimer’s disease, schizophrenia, depressive disorders, and different pain symptoms [ 15 , 16 , 19 - 25 ]. Specifically, abnormalities in spontaneous brain activity in CLBP patients have been revealed in widespread brain regions such as the medial prefrontal cortex (mPFC), supplementary motor area (SMA), anterior cingulate cortex (ACC), precuneus (PC), amygdala (AMY), thalamus (THA), and insula (INS) [ 18 ].

Although there is a study to evaluate the effect of the acupuncture technique on the low-frequency BOLD signal oscillation response in the insular in patients with CLBP [ 26 ], no study has not been existed to investigate ALFF changes followed by rTMS application in different brain regions of CLBP patients. Therefore, the present study for the first time aimed to determine the regional spontaneous BOLD signal fluctuations and pain alterations followed by rTMS in non-specific CLBP patients.

Material and Methods

Participants

Fifteen (7 males, 8 females) right-handed non-specific CLBP patients, originally recruited from Tehran University of Medical Sciences (TUMS) hospitals, participated in this experimental study. The inclusion criteria were that all patients would be aged between 20 and 60 years with a disease duration of more than three months, and also experienced pain intensity of at least 4 (Visual Analog Scale/Score (VAS)=0-10) in the week before the first rTMS session. Patients were excluded if they had potential contraindications for the use of TMS or MRI such as a history of seizures, claustrophobia, brain trauma, brain surgery, a pacemaker, or other metallic implants.

Ten right-handed healthy control (HC) voluntaries (4 males, 6 females) also participated in this study. The inclusion criteria were that all subjects would be aged between 20 and 60 years without any history of chronic pain, and also be able to undergo MRI. Written informed consent was obtained from all participants including patients and HCs.

Study Paradigm

As shown in Figure 1, this was a pre-post study design. Resting-state fMRI was designed to evaluate the effect of brain stimulation of rTMS as a state-of-the-art pain relief method on spontaneous neural activity in non-specific CLBP patients. Two fMRI scanning sessions, before and after induce of the brain stimulation, were obtained. This study focused on the low frequency of BOLD signals. At each study session (i.e., before and after the brain stimulation), subjects were given questionnaires by Roland–Morris (R-M) to assess pain-related disability, also the conventional method of VAS to determine the pain intensity.

Figure 1. Experimental protocol: The experiment consisted of two sessions of fMRI before and 2 days after rTMS for every subject. Pain intensity and disability score were assessed every session. (CLBP: Chronic low back pain, fMRI: Functional magnetic resonance imaging, rTMS: Repetition transcranial magnetic stimulation, ALFF: Amplitude of low-frequency fluctuation)

Brain Stimulation

All participants underwent a brain rTMS stimulation in the first session of the study. The rTMS was performed via an eight-shaped coil (MC-B65-HO butterfly shape coil MagVenture) connected to a magnetic stimulator (MagPro X100, MagVenture) placed over the motor cortex corresponding to the M1 brain region. The stimulation parameters provided a frequency of 20 Hz with a field intensity of 95% of motor threshold (MT) and trains of 40 pulses (in 2 sec), followed by a 28 sec rest period. It should be noted that MT was defined as the lowest intensity that yielded motor evoked potentials (MEPs) in 50% of consecutive trials for each participant obtained with the coil located on the hand area of the left hemisphere under TMS conditions until observed the muscle contraction of the right hand.

Pain Evaluation

To measure the pain intensity experienced by the CLBP patients, a VAS was used. The VAS comprised a horizontal line of 10-cm which no pain (0) indicates on the left side (start point of the line) and the strongest pain imaginable (10) at the right side (the end point of the line). A ruler was used to measure the pain intensity by the measured length (in centimeters) from the start point to the point marking a cross on the VAS by each patient. The pain intensity was again determined within 2 days of completing the rTMS session using the same procedures described above. The percentage of pain reduction (PPR) was obtained from the measurement of VAS scores at pre and post rTMS sessions according to the following equation [ 10 ]:

PPR=VASpre-VASpostVASpre×100 (1)

Disability Evaluation

To measure the disability scale, a self-administered questionnaire of Roland-Morris (RM) with a 24-point scale was used in which greater levels of disability were linked directly with the higher numbers on it. The RM is sensitive to change over time for low back pain patients. Participants were asked to mark next to every appropriate statement on a 24-point scale. To get a disability score for non-specific CLBP patients, add up all the marked statements. The disability was again determined within 2 days of completing the rTMS session using the same procedures described above. The percentage of disability change (PDC) was calculated from the Roland–Morris (RM) questionnaire scores measured at pre and post rTMS sessions according to the following equation:

PDC=RMpre-RMpostRMpre×100 (2)

Imaging Data Acquisition

MRI data were acquired using a 3.0 Tesla Siemens, MAGNETRON Prisma (Munich, Germany) scanner with a standard 64-channel radio-frequency head coil. For each subject, a high-resolution T1-weighted structural image for registration purposes was collected using a magnetization-prepared rapid gradient-echo imaging (MPRAGE) with repetition time/echo time=2300/2.32 ms; FOV=240 mm×240 mm; matrix size=256×256; voxel size =0.9×0.9×0.9 mm3 and flip angle=8°. During the fMRI scan, all participants were asked to lie in a supine position in the scanner with their eyes closed, without falling asleep, and not to think about anything in particular. Their heads were fixed to minimize head motion using pillows and foam pads. T2*-weighted resting fMRI scans were acquired using a gradient recalled echo echo-planar imaging (GRE-EPI) sequence with repetition time=3000 ms; echo time=40 ms; flip angle=90°; number of slices=42; field of view=220×220 mm2 without gap, matrix size=64×64; voxel size=3.4×3.4×3 mm3, slice thickness=3 mm; and 200 volumes.

fMRI Signal Analysis

Resting-state fMRI data were preprocessed using the Data Processing & Analysis for Brain Imaging (DPABI v6.0, http://www.rfmri.org/dpabi) toolbox for SPM12 (Wellcome Department of Cognitive Neurology. London, UK) running in MATLAB 2017b (Mathworks, Inc., Natick, MA, USA) with these steps: for stabilization, the first 10 time points were discarded. Then slice timing correction, motion correction, registration, and normalization in Montreal Neurological Institute [MNI] 152 space, and spatial smoothing with 6 mm FWHM were respectively performed. It should be mentioned that the exclusion criteria for the max head motion in the motion correction step were 3 mm for translation and 3° for rotation motion. Also, in the normalization step, T1-weighted high-resolution images were registered to the mean fMRI images, and the resulting images were segmented and transformed into standard MNI spaces using the Diffeomorphic Anatomical Registration Through Exponentiated Lie Algebra toolbox (DARTEL) [ 17 ].

ALFF Computing

For ALFF analysis, a preprocessed BOLD time series x(i)=x(ti), i=1, 2, …, N with N time points can be written in the frequency domain as:

xi=1Nk=1NXke2πjNki, i=1,2,…,N (3)

Where the spectral amplitude terms, Xk=x(fk) with k=1, 2, …, N are the discrete Fourier transformed coefficients. Then, the power spectrum (PS) can be calculated at the frequency band of 0.01<f<0.08 Hz to reduce the low-frequency drift effects and high-frequency physiological noise in each voxel.

PS(0.01,0.08)=1Nk=0.010.08|Xk|2 (4)

Since each term in Eq. 4 is proportional to the square of the amplitude at the frequency k, ALFF can be computed using the following equation (i.e., the square root of the PS density at each frequency component) [ 17 , 21 , 22 ]:

ALFF=1Nk=0.010.08|Xk| (5)

Statistical Analyses

Individual-level values were analyzed using the Statistical Package for Social Science (SPSS, v. 22.0, IBM Corporation, Armonk, NY). The characteristics variables were performed using descriptive analysis (mean±SD) for quantitative variables. Percentage of pain reduction and percentage of disability changes were compared within the/a patient group using paired t-test analysis, with the statistical significance levels set at P≤0.05.

DPABI V6.0 was used for statistical analysis on ALFF. Different types of t-test statistical analysis were used to compare the values in groups. Multiple comparison corrections were performed based on Gaussian random field (GRF) theory (voxel wise P<0.01, cluster-wise P<0.05, and two-tailed test).

Finally, the correlation analysis was used to explore the relationship between ALFF values in brain areas and VAS scores in pre and post-stimulation sessions.

Results

Table 1, shows the descriptive statistics for general demographics data and pain information of CLBP participants and healthy control.

Characteristics Age (year) Female/Men (%) Pain Duration (year) *BMI (Kg/m2)
Non-specific CLBP (N=15) 44.87 (10.89) 54/46 1.88 (1.10) 27.39 (4.59)
HC (N=10) 40.88 (8.73) 60/40 --- 23.70 (5.46)
N: Number of participants, CLBP: Chronic low back pain, HC: Health control, BMI: Body mass index, *Normal BMI: 18.5-25, more than 25 is overweight, less than 18.5 is underweight.
Table 1.Demographics for non-specific CLBP (Chronic low back pain) patients and HCs (Health control). The values presented are ‘mean (standard deviation)’ for characteristics.

The measured pain intensity and disability score for non-specific CLBP patients in every session are presented in Table 2. The mean pain intensity changed significantly from 7.25±0.99 in the first session to 4.60±1.03 in the second session (i.e, 2 days after brain stimulation using rTMS) in the patient group (t=6.80, P<0.001). The mean percentage of pain reduction (PPR) was 36.23%±13.28. The mean disability score was 15.58±3.40 before brain stimulation while it was 13.25±3.46 after stimulation (t=12.72, P<0.001). The mean percentage of disability change (PDC) was 15.89%±7.82.

Non-specific CLBP (N=15) VAS Score (#) Mean±SD R-M Score (#) Mean±SD
Pre Stimulation 7.25±0.99 15.58±3.40
Post Stimulation 4.60±1.03 13.25±3.46
PPR% or PDC% 36.22±13.28 15.89±7.82
CLBP: Chronic low back pain, VAS: Visual analysis scale, R-M: Roland-Morris, SD: Standard deviation, PPR: Percentage of pain reduction, PDC: Percentage of disability change
Table 2.Measured pain and disability results before and after rTMS (Repetition transcranial magnetic stimulation) stimulation in CLBP (Chronic low back pain) patients.

Figure 2 shows the ALFF analysis process from non-specific CLBP patients in rs-fMRI data of the pre and post-rTMS stimulation sessions. Figure 2 (A) demonstrates the BOLD signal was extracted from a 6 mm diameter ROI (e.g., INS (L)) centered on the peak coordinates (-47, 11, -5) during the scan time. The time series were transformed into the frequency domain with a Fast Fourier Transform (FFT) (Figure 2 (B)) and the power spectrum was obtained (Figure 2 (C)). A filtered power spectrum in the frequency band of 0.01 Hz to 0.08 Hz was used to measure ALFF (Figure 2 (D)). The power of a frequency spectrum is proportional to the square of the amplitude of this frequency component, so, the square root was calculated at the frequency band (0.01–0.08 Hz) to measure the average of the square root as ALFF in this region (Figure 2 (E)).

Figure 2. ALFF analysis for pre (blue) and post (red) rTMS sessions: (A) Time series of an individual ROI (INS (L)), (B) transformed signal by FFT, (C) power spectrum of the signal, (D) filtered power spectrum in 0.01 to 0.08, (E) power spectrum between 0.01 and 0.08 Hz, i.e., ALFF. rTMS: Repetition transcranial magnetic stimulation, ROI: region of interest, INS: Insula, FFT: Fast Fourier Transform, ALFF: Amplitude of low-frequency fluctuation

Accordingly, Figure 3 and Table 3 show the ALFF analysis results in the HC group and non-specific CLBP patients before and after brain stimulation of rTMS.

Figure 3. ALFF maps: Axial view of the results from the ALFF analysis (averaged across subjects) in (A) healthy control group (B) patient group at the first session before rTMS and (C) patient group at second session in 2-4 days after rTMS session (P<0.05). Warm colors indicate increased ALFF in different brain regions. ALFF: Amplitude of low-frequency fluctuation, rTMS: Repetition transcranial magnetic stimulation, INS: Insula, mPFC: Medial prefrontal cortex, THA: Thalamus, SMA: Supplementary motor area, PC: Precuneus

Region Peak MNI T-value T-value T-value
X Y Z HC Pre-stimulation Post-stimulation
INS (L) -47 11 -5 3.89 5.44 3.80
INS (R) 46 15 -5 3.81 5.02 3.83
SMA (L) 0 23 65 3.91 5.06 3.71
SMA (R) 3 23 65 3.98 4.38 3.64
THA (L) -12 -14 20 4.63 4.52 3.79
THA (R) 9 -6 15 5.41 6.37 3.90
PC (L) -6 -55 70 4.33 4.02 2.81
PC (R) 5 -60 70 5.55 3.05 3.73
mPFC 1 55 -3 4.81 4.69 5.75
MNI: Montreal neurological institute, HC: Health control, INS: Insula, SMA: Supplementary motor area, THA: Thalamus, PC: Precuneus, Mpfc: Medial prefrontal cortex
Table 3.ALFF (Amplitude of low-frequency fluctuation) analysis results derived in HC (Health control) group and non-specific CLBP (Chronic low back pain) patients.

Subsequently, paired t-test analysis was performed to find ALFF differences before and after rTMS in non-specific CLBP patients, and a two-sample t-test was used to find differences between patients and HCs (Figure 4).

Figure 4. ALFF group analysis. (A) Within-group analyses of the ALFF in non-specific CLBP patients before and after rTMS. (B) Between-group analyses of the ALFF between patients and HC groups (two-tailed, voxel-level P<0.01; GRF correction, cluster-level P<0.05). CLBP: Chronic low back pain, HC: Health control, ALFF: Amplitude of low-frequency fluctuation, rTMS: Repetition transcranial magnetic stimulation, GRF: Gaussian random field

The correlation analysis was used to explore the logical potential relationship between the pain intensity in pre/post-stimulation sessions and the ALFF values in brain regions. As shown in Figure 5, the Pearson correlation coefficient between VAS and ALFF in mPFC in the pre-rTMS session and the post-rTMS session was -0.54 and -0.56, respectively. In addition, the Pearson correlation coefficient between VAS and ALFF in the INS (L) in the pre-rTMS session is 0.59 whereas this is 0.58 in the post-rTMS session (Figure 5).

Figure 5. Correlation analysis of ALFF and pain intensity. (A) The negative correlation between pain intensity (i.e., VAS) and ALFF in mPFC before rTMS stimulation (Pearson correlation, r=-0.54), (B) the negative correlation between VAS and ALFF in mPFC after rTMS stimulation (Pearson correlation, r=-0.56), (C) the positive correlation between pain intensity (i.e., VAS) and ALFF in INS (L) before rTMS stimulation (Pearson correlation, r=0.59), (D) the negative correlation between VAS and ALFF in INS (L) after rTMS stimulation (Pearson correlation, r=0.58). ALFF: Amplitude of low-frequency fluctuation, mPFC: Medial prefrontal cortex, INS (L): Left insula, rTMS: Repetition transcranial magnetic stimulation, VAS: Visual analysis scale.

Discussion

The major aim of this study was to investigate the differences in the regional spontaneous activity of the brain in non-specific CLBP patients before and after rTMS brain stimulation by using ALFF analysis. Although disrupted low-frequency oscillation amplitude has been reported in various clinical disorders relative to healthy control [ 19 - 22 , 25 ], neural modulation effects investigation using ALFF analysis is a new approach in neuroimaging chronic pain symptoms studies. Therefore, in the present study, the ALFF analysis was applied to evaluate the alterations of spontaneous neural activity in non-specific CLBP patients after induced rTMS.

The results of rs-fMRI imaging data analysis in the present study after induced-rTMS compared with pre-stimulation in non-specific CLBP patients, initially demonstrated changes in spontaneous brain activity of the bilateral INS, bilateral THA, bilateral SMA, bilateral PC, and mPFC. Although the mechanism of brain changes after rTMS has not been clearly known yet, some previous imaging studies demonstrated that M1 stimulation by rTMS resulted in the brain activity alteration not only in the motor system but also in the regions which are related to pain processing and modulation [ 13 , 27 ]. Importantly, in the present study, decreased ALFF was observed in bilateral INS, bilateral THA, and bilateral SMA after using high-frequency rTMS for patients suffering from CLBP. In contrast, ALFF of mPFC increased after brain stimulation in the patient group (Table 3). Besides, THA showed a notable difference in pre and post-rTMS brain stimulation sessions in CLBP patients (Figure 4). Specifically, INS, THA, and prefrontal cortices have been introduced as essential brain areas in the pain matrix that disrupted activation in these regions is associated with the experience of pain [ 28 ]. Therefore, in line with recent studies, these brain regions that are sensitive to pain intensity changes may have the potential for measurement of treatment response in CLBP patients [ 18 ]. In this regard, consistent with our findings, decreased spontaneous neural activity in INS simultaneous with decreased pain intensity in CLBP patients after ankle acupuncture was observed [ 26 ]. Despite applying different pain reduction methods in that study (using acupuncture) with our study (using rTMS), it seems that INS is a major brain region in response to the sense of pain in patients suffering from CLBP. However, no study exists to investigate spontaneous neural activity changes followed by rTMS in non-specific CLBP patients.

Another major concern in our study is that the spontaneous activation of brain regions demonstrated differences between CLBP patients and healthy control. Actually, CLBP patients showed the higher amplitude of low-frequency oscillation values in bilateral SMA and bilateral INS relative to HC while the lower ALFF was observed in bilateral PC and mPFC (Table 3). Some previous studies also found the augmented activation in the SMA in chronic low back pain patients compared with HC, thereby, SMA has been suggested to discriminate CLBP subjects from HCs [ 18 ]. Similarly, in our study, SMA showed a significant difference in non-specific CLBP compared with HC (Figure 4). Moreover, the mPFC and PC are important regions in the default mode network (DMN) and the INS is an essential region in the salience network (SN). DMN is disrupted in CLBP and involved in the spontaneous disengagement of pain attention and is suppressed in pain presence [ 18 , 29 ]. In addition, decreased ALFF was observed in the PC in CLBP patients compared with HC which is under another low back pain study. However, they induced acute LBP by injection of hypertonic saline into the back muscles of the healthy subjects [ 24 ]. Despite the difference in the kind of LBP (being acute or chronic), the decreased spontaneous BOLD signal in PC was observed in both conditions. On the other side, the salience network (SN) as one of the important brain responsive regions in CLBP which monitors sensory input alterations, is abnormal in CLBP [ 26 , 29 ].

On the clinical side, the result of the current study in agreement with plenty of research has shown to decrease significantly in VAS pain scores after a single session of high-frequency rTMS in chronic pain patients [ 9 , 10 , 12 , 30 ]. Indeed, the mean analgesic effect, measured with the VAS, in CLBP participants was 36.22% after high-frequency rTMS of M1 stimulation. This result is completely consistent with other reports that are expected to decrease pain intensity by 20% to 45% after active brain stimulation using high-frequency rTMS of M1 [ 7 , 12 , 13 , 31 , 32 ]. So, this result also supports the idea of the beneficial role of rTMS in pain management.

Furthermore, the results of the present study demonstrated a significant correlation between ALFF of mPFC and pain intensity in CLBP patients (P<0.05) (Figure 5). Pearson correlation analysis revealed an indirect correlation between mPFC ALFF values and VAS scores in CLBP patients in the pre and post-stimulation of -0.54 and -0.56, respectively. Recent studies also found the ALFF values in default mode network (DMN) regions are negatively correlated to the pain intensity in CLBP patients and also in patients with discogenic low back and leg pain [ 17 , 18 ]. Therefore, under these studies, the correlation analysis showed a possible relationship between the spontaneous brain activity of mPFC as the main part of DMN and the pain intensity of CLBP in two sessions of imaging before and after induced-rTMS. This result highlighted the role of frontal regions of the brain in the cognitive dimension of pain that is reported in some previous research [ 33 ]. Moreover, a direct correlation between the ALFF in the INS (L) and pain intensity in pre and post-stimulation sessions was obtained in the present study (rpre-rTMS=0.59, rpost-rTMS=0.58) (Figure 5). Inconsistent with our finding, a positive correlation between the VAS score of CLBP patients and ALFF values in INS (L) was explored in some previous studies [ 18 , 26 ]. The results support the idea that INS is an important part of the pain matrix and also have a multidimensional role in different aspect of pain including sensory, cognitive, and emotional [ 24 , 26 , 33 - 35 ].

Conclusion

In conclusion, altered spontaneous neural activity in different brain regions including the insula, thalamus, medial prefrontal cortex, precuneus, and the supplementary motor area may provide evidence to introduce ALFF as a non-invasive neuro biomarker to monitor CLBP symptoms. The results of correlation analyses potentially open up a new approach to link ALFF as a proper objective neuroimaging parameter with the VAS score as a classic subjective pain intensity evaluation.

Acknowledgment

The authors thank with grateful appreciation for the technical support from the staff at National Brain Mapping Lab (NBML), Tehran, Iran. This study was part of the Ph.D. thesis of Mahboubeh Masoumbeigi.

Authors’ Contribution

M. Masoumbeigi performed the data acquisition, fMRI analysis, interpretation of data, and drafting of the work as the principal author. N. Riyahi Alam designed the main conception of this work and approved the final version to be published. R. Kordi and M. Rostami contributed to the recruitment and clinical assessment of the patients. A. Rahimi foroushani analyzed the results statistically. AH. Jafari and H. Hashemi revised the manuscript critically for important intellectual content. A. Ebrahimpour contributed to the data acquisition. All the authors read, modified, and approved the final version of the manuscript.

Ethical Approval

The protocol of the human study was approved by the local ethical committee, Tehran University of Medical Sciences (TUMS), Tehran, Iran (Approval number: IR.TUMS.MEDICINE.REC.1397.957).

Informed Consent

Written informed consent was obtained from all participants including patients and HCs.

Funding

The research leading to these results received funding from the research chancellor of Tehran University of Medical Sciences (TUMS) Tehran, Iran (Grant number: 97-03-30-40194).

Conflict of Interest

None

References

  1. Krismer M, Van Tulder M. Strategies for prevention and management of musculoskeletal conditions. Best Pract Res Clin Rheumatol. 2007; 21(1):77-91. DOI | PubMed
  2. Oliveira CB, Maher CG, Pinto RZ, Traeger AC, Lin CC, et al. Clinical practice guidelines for the management of non-specific low back pain in primary care: an updated overview. Eur Spine J. 2018; 27(11):2791-803. DOI | PubMed
  3. Machado LA, Kamper SJ, Herbert RD, Maher CG, McAuley JH. Analgesic effects of treatments for non-specific low back pain: a meta-analysis of placebo-controlled randomized trials. Rheumatology (Oxford).. 2009; 48(5):520-7. DOI | PubMed
  4. Ghafouri M, Teymourzadeh A, Nakhostin-Ansari A, Sepanlou SG, Dalvand S, et al. Prevalence and predictors of low back pain among the Iranian population: Results from the Persian cohort study. Ann Med Surg (Lond).. 2022; 74:103243. Publisher Full Text | DOI | PubMed
  5. Noormohammadpour P, Mansournia MA, Koohpayehzadeh J, Asgari F, Rostami M, Rafei A, Kordi R. Prevalence of Chronic Neck Pain, Low Back Pain, and Knee Pain and Their Related Factors in Community-Dwelling Adults in Iran: A Population-based National Study. Clin J Pain. 2017; 33(2):181-7. DOI | PubMed
  6. Lefaucheur JP, Nguyen JP. A practical algorithm for using rTMS to treat patients with chronic pain. Neurophysiol Clin. 2019; 49(4):301-7. DOI | PubMed
  7. Galhardoni R, Correia GS, Araujo H, Yeng LT, Fernandes DT, et al. Repetitive transcranial magnetic stimulation in chronic pain: a review of the literature. Arch Phys Med Rehabil. 2015; 96(4 Suppl):S156-72. DOI | PubMed
  8. Fregni F, Freedman S, Pascual-Leone A. Recent advances in the treatment of chronic pain with non-invasive brain stimulation techniques. The Lancet Neurology. 2007; 6(2):188-91. DOI
  9. Lefaucheur JP, Drouot X, Menard-Lefaucheur I, Zerah F, Bendib B, Cesaro P, et al. Neurogenic pain relief by repetitive transcranial magnetic cortical stimulation depends on the origin and the site of pain. J Neurol Neurosurg Psychiatry. 2004; 75(4):612-6. Publisher Full Text | DOI | PubMed
  10. Ayache SS, Ahdab R, Chalah MA, Farhat WH, Mylius V, Goujon C, et al. Analgesic effects of navigated motor cortex rTMS in patients with chronic neuropathic pain. Eur J Pain. 2016; 20(9):1413-22. DOI | PubMed
  11. Ambriz-Tututi M, Alvarado-Reynoso B, Drucker-Colín R. Analgesic effect of repetitive transcranial magnetic stimulation (rTMS) in patients with chronic low back pain. Bioelectromagnetics. 2016; 37(8):527-35. DOI | PubMed
  12. Johnson S, Summers J, Pridmore S. Changes to somatosensory detection and pain thresholds following high frequency repetitive TMS of the motor cortex in individuals suffering from chronic pain. Pain. 2006; 123(1-2):187-92. DOI | PubMed
  13. Choi GS, Chang MC. Effects of high-frequency repetitive transcranial magnetic stimulation on reducing hemiplegic shoulder pain in patients with chronic stoke: a randomized controlled trial. Int J Neurosci. 2018; 128(2):110-116. DOI | PubMed
  14. Williams JA, Imamura M, Fregni F. Updates on the use of non-invasive brain stimulation in physical and rehabilitation medicine. J Rehabil Med. 2009; 41(5):305-11. DOI | PubMed
  15. Zhou C, Tang X, You W, Wang X, Zhang X, Zhang X, Yu M. Altered Patterns of the Fractional Amplitude of Low-Frequency Fluctuation and Functional Connectivity Between Deficit and Non-Deficit Schizophrenia. Front Psychiatry. 2019; 10:680. Publisher Full Text | DOI | PubMed
  16. Far YN, Alam NR, Jalalvandi M, Hashemi H. Investigate the Frequency Spectrum of Noise Signals Affecting Synchronous Signals of fMRI and fNIRS in Resting and Moving Wrists. Frontiers in Biomedical Technologies. 2019; 6(3):139-47. DOI
  17. Zhou F, Gu L, Hong S, Liu J, Jiang J, Huang M, Zhang Y, Gong H. Altered low-frequency oscillation amplitude of resting state-fMRI in patients with discogenic low-back and leg pain. J Pain Res. 2018; 11:165-76. Publisher Full Text | DOI | PubMed
  18. Zhang B, Jung M, Tu Y, Gollub R, Lang C, et al. Identifying brain regions associated with the neuropathology of chronic low back pain: a resting-state amplitude of low-frequency fluctuation study. Br J Anaesth. 2019; 123(2):e303-11. Publisher Full Text | DOI | PubMed
  19. Wang YM, Cai XL, Zhang RT, Wang Y, Madsen KH, Sørensen TA, et al. Searchlight classification based on Amplitude of Low Frequency Fluctuation and functional connectivity in individuals with obsessive-compulsive symptoms. Cogn Neuropsychiatry. 2019; 24(5):322-34. DOI | PubMed
  20. Yang L, Yan Y, Wang Y, Hu X, Lu J, Chan P, et al. Gradual Disturbances of the Amplitude of Low-Frequency Fluctuations (ALFF) and Fractional ALFF in Alzheimer Spectrum. Front Neurosci. 2018; 12:975. Publisher Full Text | DOI | PubMed
  21. Jing B, Liu CH, Ma X, Yan HG, Zhuo ZZ, Zhang Y, Wang SH, Li HY, Wang CY. Difference in amplitude of low-frequency fluctuation between currently depressed and remitted females with major depressive disorder. Brain Res. 2013; 1540:74-83. DOI | PubMed
  22. Zang YF, He Y, Zhu CZ, Cao QJ, Sui MQ, Liang M, et al. Altered baseline brain activity in children with ADHD revealed by resting-state functional MRI. Brain Dev. 2007; 29(2):83-91. DOI | PubMed
  23. Pan ZM, Li HJ, Bao J, Jiang N, Yuan Q, Freeberg S, Zhu PW, Ye L, Ma MY, Huang X, Shao Y. Altered intrinsic brain activities in patients with acute eye pain using amplitude of low-frequency fluctuation: a resting-state fMRI study. Neuropsychiatr Dis Treat. 2018; 14:251-7. Publisher Full Text | DOI | PubMed
  24. Zhang SS, Wu W, Yang JM, Wang CH. Abnormal Spontaneous Brain Activity in Acute Low-Back Pain Revealed by Resting-State Functional MRI. Am J Phys Med Rehabil. 2017; 96(4):253-9. DOI | PubMed
  25. Wang JJ, Chen X, Sah SK, Zeng C, Li YM, Li N, et al. Amplitude of low-frequency fluctuation (ALFF) and fractional ALFF in migraine patients: a resting-state functional MRI study. Clin Radiol. 2016; 71(6):558-64. DOI | PubMed
  26. Xiang A, Yu Y, Jia X, Ma H, Liu H, Zhang Y, et al. The low-frequency BOLD signal oscillation response in the insular associated to immediate analgesia of ankle acupuncture in patients with chronic low back pain. J Pain Res. 2019; 12:841-50. Publisher Full Text | DOI | PubMed
  27. Yang S, Chang MC. Effect of Repetitive Transcranial Magnetic Stimulation on Pain Management: A Systematic Narrative Review. Front Neurol. 2020; 11:114. Publisher Full Text | DOI | PubMed
  28. Balenzuela P, Chernomoretz A, Fraiman D, Cifre I, Sitges C, Montoya P, Chialvo DR. Modular organization of brain resting state networks in chronic back pain patients. Front Neuroinform. 2010; 4:116. Publisher Full Text | DOI | PubMed
  29. Tu Y, Jung M, Gollub RL, Napadow V, Gerber J, et al. Abnormal medial prefrontal cortex functional connectivity and its association with clinical symptoms in chronic low back pain. Pain. 2019; 160(6):1308-18. Publisher Full Text | DOI | PubMed
  30. Lefaucheur JP, Drouot X, Nguyen JP. Interventional neurophysiology for pain control: duration of pain relief following repetitive transcranial magnetic stimulation of the motor cortex. Neurophysiol Clin. 2001; 31(4):247-52. DOI | PubMed
  31. Quesada C, Pommier B, Fauchon C, Bradley C, Créac’h C, Murat M, et al. New procedure of high-frequency repetitive transcranial magnetic stimulation for central neuropathic pain: a placebo-controlled randomized crossover study. Pain. 2020; 161(4):718-28. DOI | PubMed
  32. Lefaucheur JP. Cortical neurostimulation for neuropathic pain: state of the art and perspectives. Pain. 2016; 157:S81-9. DOI | PubMed
  33. Ng SK, Urquhart DM, Fitzgerald PB, Cicuttini FM, Hussain SM, Fitzgibbon BM. The Relationship Between Structural and Functional Brain Changes and Altered Emotion and Cognition in Chronic Low Back Pain Brain Changes: A Systematic Review of MRI and fMRI Studies. Clin J Pain. 2018; 34(3):237-261. DOI | PubMed
  34. Baliki MN, Petre B, Torbey S, Herrmann KM, Huang L, Schnitzer TJ, et al. Corticostriatal functional connectivity predicts transition to chronic back pain. Nat Neurosci. 2012; 15(8):1117-9. Publisher Full Text | DOI | PubMed
  35. Kobayashi Y, Kurata J, Sekiguchi M, Kokubun M, Akaishizawa T, Chiba Y, et al. Augmented cerebral activation by lumbar mechanical stimulus in chronic low back pain patients: an FMRI study. Spine (Phila Pa 1976).. 2009; 34(22):2431-6. DOI | PubMed