Intro
Approximately 1.18 billion people regularly smoke tobacco ( Dai et al., 2022 ), and the prevalence of smoking prior to surgery exceeds 20% ( Howard et al., 2021 ). Nicotine, the main addictive substance in tobacco, has been shown to induce hyperalgesia during withdrawal ( Ditre et al., 2018 ; LaRowe and Ditre, 2020 ). Surgical patients are typically asked to quit smoking due to the elevated risk of postoperative complications associated with tobacco use, including delayed wound healing, tissue infections, respiratory dysfunction, and/or sepsis ( Aspera-Werz et al., 2022 ). Numerous clinical studies have reported that abstinent smokers experience greater postoperative pain and require higher doses of analgesics compared with nonsmokers ( Weingarten et al., 2011 ; Chiang et al., 2016 ; Shen et al., 2017 ; Iida et al., 2022 ).
Nicotine withdrawal-induced hyperalgesia may be caused by changes in the nervous system. Several animal studies have examined the spinal mechanisms underlying this phenomenon. For example, Lina et al. demonstrated that interferon regulatory factor 8 plays a crucial role in nicotine withdrawal-induced hyperalgesia by enhancing microglia activation and the expression of spinal P2X purinoceptor 4 and brain-derived neurotrophic factor in mice ( Guo et al., 2020 ). Zhang et al. found that nicotine withdrawal alters basal levels of glutamate decarboxylase 67 and 65, opioid receptors, endorphins, and γ-aminobutyric acid, contributing to heightened pain sensitivity in an animal model ( Zhang et al., 2019 ). Although animal studies (including our work in raphe nuclei) implicate supraspinal mechanisms in nicotine withdrawal-induced hyperalgesia ( Schmidt et al., 2001 ; Shen et al., 2021 ), direct neuroimaging evidence in humans remains scarce.
Over the past decade, resting-state functional magnetic resonance imaging (rs-fMRI) has been widely used to investigate the brain dysfunction related to nicotine dependence and cessation ( Beaver et al., 2011 ; Lesage et al., 2020 ). A recent meta-analysis which encompassed a total of 35 rs-fMRI studies revealed that individuals with tobacco use disorder exhibited increased intrinsic function in the right cerebellum crus2, left superior frontal gyrus, left inferior parietal gyrus, and left supplementary motor area and decreased intrinsic function in the right gyrus rectus, right superior/middle frontal gyrus, and left inferior frontal gyrus ( Ma et al., 2025 ). Another rs-fMRI study found that compared with satiated state, abstinent smokers showed altered functional connectivity (FC) between striatal and cingulo-insular network ( Sweitzer et al., 2016 ). However, few studies have focused on comparing the difference of brain function between abstinent smokers and nonsmokers using rs-fMRI, specifically exploring the altered brain activity associated with nicotine withdrawal-induced hyperalgesia. In this study, we used rs-fMRI to compare regional brain activity/connectivity—assessed via fractional amplitude of low-frequency fluctuation (fALFF; 0.01–0.1 Hz), regional homogeneity (ReHo), and rs-FC—between abstinent smokers and nonsmokers. To explore the relationship between altered regional brain activity/connectivity and hyperalgesia, we measured preoperative pain thresholds and the postoperative analgesic requirements in all patients. We hypothesized that nicotine abstinence would lead to specific alterations in regional brain activity/connectivity, which could be associated with lower pain thresholds and increased analgesic use in abstinent smokers. We explored the relationship between smoking behaviors and pain thresholds or postoperative analgesic requirements, as well as the role of different brain regions. These findings may suggest that altered brain activity contributes to nicotine withdrawal-induced hyperalgesia and could offer novel insights into managing postoperative pain in this population.
Results
Sixty right-handed male patients scheduled for surgery due to liver masses were enrolled in the study between April 1, 2017, and November 8, 2017. Among them, 30 were nonsmokers (nonsmoking group), and 30 were smokers who had been asked to quit smoking upon hospitalization (abstinent smoker group). There were no significant differences in demographic and clinical characteristics between the groups, except for white blood cell count ( Table 1 ). Plasma and urine cotinine concentrations were significantly higher in the abstinent smoker group compared with the nonsmoker group ( Table 1 ), confirming the validity of group classification based on smoking history. Details on smoking behavior, nicotine dependence severity, and withdrawal symptoms in the abstinent smoker group are shown in Table 2 .
Patient characteristics of the abstinent smoker and nonsmoker groups
ASA, American Society of Anesthesiologists. Data are presented as means (standard deviations) or numbers. Group differences for continuous variables were evaluated using unpaired t tests or the Mann–Whitney U test. Group differences for categorical variables were tested using the chi-square test. p < 0.05 was considered statistically significant.
p < 0.05.
Psychological variables, pain sensitivity, and smoking behaviors in the abstinent smoker and nonsmoker groups
SAI, State Anxiety Inventory; TAI, Trait Anxiety Inventory; BDI, Beck Depression Inventory; PSQ, Pain Sensitivity Questionnaire; FTND, Fagerstrom Test for Nicotine Dependence; MNWS, Minnesota Nicotine Withdrawal Scale; QSU, Questionnaire of Smoking Urges. Data are presented as means (standard deviations). Group differences for continuous variables were evaluated with unpaired t tests or the Mann–Whitney U test. Pain thresholds were evaluated using an electrical stimulator a and von Frey filaments b .
Repeated analysis of variance (repeated ANOVA) was used for analysis, and post hoc comparisons were performed with post hoc Bonferroni’s test to compare group differences at each time point.
Represents U values calculated from Mann–Whitney U test.
Represents F values of interaction effect calculated from repeated ANOVA.
p < 0.05.
p < 0.01.
p < 0.001.
As shown in Table 2 , compared with nonsmokers, abstinent smokers exhibited higher PSQ scores ( t = −3.28; p = 0.002; 95% CI, 8.72–36.08), lower EPT ( t = 2.30; p = 0.03; 95% CI, −0.90 to −0.06), and lower MPT ( U = 184; p < 0.001; 95% CI, −27.19 to −1.38), indicating higher pain sensitivity and lower pain threshold in abstinent smokers before surgery. In addition, repeated-measure ANOVA revealed a significant group × time interaction for postoperative analgesic requirements ( F [3,174] = 9.37; p < 0.001; Fig. 2 D ). Post hoc tests indicated that abstinent smokers required significantly more analgesics than nonsmokers at both 24 h (95% CI, 0.19–0.51; p < 0.001, Bonferroni’s corrected) and 48 h (95% CI, 0.07–0.40; p = 0.0012, Bonferroni’s corrected) postoperatively. These results suggest that, compared with nonsmokers, abstinent smokers had lower pain thresholds and greater postoperative analgesic requirements.
Brain regions with altered regional activity and FC in abstinent smokers. A , The voxel-wise fALFF results suggested that bilateral medial vmPFC activity was significantly lower in the abstinent smoker group than in the nonsmoker group. B , The voxel-wise ReHo results suggested that left CAL activity was significantly higher in the abstinent smoker group than in the nonsmoker group. C , The ROI-based rs-FC analysis revealed that rs-FC of vmPFC to bilateral MTG and PCUN was lower in abstinent smokers.
After applying exclusion criteria related to head movement, 58 participants (29 abstinent smokers and 29 nonsmokers) remained eligible for rs-fMRI analysis. Group differences in fALFF and ReHo are illustrated in Figure 2 and summarized in Table 3 . Compared with nonsmokers, abstinent smokers showed significantly lower fALFF in the bilateral ventromedial prefrontal cortex (vmPFC; Fig. 2 A ), with a total cluster size of 13 voxels. In contrast, ReHo was significantly higher in the left calcarine (CAL; Fig. 2 B ), with a cluster of 34 voxels, in abstinent smokers relative to nonsmokers. These two brain regions were subsequently defined as ROIs for further rs-FC analysis. Compared with nonsmokers, abstinent smokers exhibited significantly lower rs-FC between vmPFC and both the bilateral middle temporal gyrus (MTG) and the precuneus (PCUN; Fig. 2 C , D ; Table 4 ). No other group differences reached statistical significance. These findings suggest that abstinent smokers and nonsmokers differ in preoperative regional brain activity and FC patterns.
Brain regions with significant differences in fMRI measures (fALFF, ReHo, and rs-FC) between abstinent smokers and nonsmokers
fALFF, fractional amplitude of low-frequency fluctuation; ReHo, regional homogeneity; rs-FC, resting-state functional connectivity; MNI, Montreal Neurological Institute. Results are Gaussian random field-corrected [voxel-level, p < 0.001; cluster-level p < 0.05 (two-tailed)].
Correlation between regional fMRI activity and pain thresholds and postoperative analgesic requirements at 48 h postsurgery in the abstinent smokers
fMRI, functional magnetic resonance imaging, fALFF, fractional amplitude of low-frequency fluctuation; ReHo, regional homogeneity; rs-FC, resting-state functional connectivity; MNI, Montreal Neurological Institute. Results are Gaussian random field-corrected [voxel-level, p < 0.001; cluster-level, p < 0.05 (two-tailed)].
Whole-brain correlation analysis between functional fMRI maps with pain thresholds ( Table 4 ) showed that in the abstinent smoker group, EPT was significantly positively correlated with fALFF values ( Fig. 3 A ) in the right CAL, with a cluster size of 20 voxels. For ReHo ( Fig. 3 B ), significant positive correlations were found between EPT and ReHo in the left CAL (cluster size of 43 voxels) and the left posterior cingulate cortex (PCC; cluster size of 36 voxels). No significant correlations were found between EPT and rs-FC values. However, as shown in Figure 3 C , MPT was significantly positively correlated with rs-FC from the vmPFC to the right dorsolateral prefrontal cortex (DLPFC). No significant correlations were found between MPT and either fALFF or ReHo values. In the nonsmoking group, no significant correlations were observed between brain activity/FC and either EPT or MPT.
Association between pain thresholds and whole-brain rs-fMRI. A , EPT was positively correlated with fALFF value in the right CAL cortex (CAL). B , EPT were positively correlated with the ReHo value in the left CAL and the left PCC. C , MPT was positively correlated with rs-FC of the bilateral vmPFC to the right DLPFC. D , The scatterplots showed the correlation between duration of abstinence and ReHo values in the left CAL and left PCC and between duration of abstinence and EPT. The shaded regions represent 95% credibility interval. E , The mediation effect of left CAL and PCC in association with duration of abstinence and EPT. A1, A2, effect of duration of abstinence on mediator. B1, B2, effect of mediator on EPT. C, total effect of duration of abstinence on EPT. C’, direct effect of duration of abstinence on EPT.
Next, we explored the correlations between smoking behaviors (i.e., cigarettes/day, smoking duration, abstinence duration, FTND, MNWS, QSU, and plasma/urine cotinine levels) and the brain regions associated with pain thresholds. Duration of abstinence was positively correlated with ReHo values in both the left CAL and left PCC ( Fig. 3 D ). Notably, duration of abstinence was also significantly positively correlated with EPT ( r = 0.6211; p = 0.0003; Fig. 3 D ). Based on these findings, we hypothesized that the lower pain thresholds observed in abstinent smokers may be mediated by altered regional activity in the left CAL and left PCC. To test this hypothesis, a mediation analysis was conducted, with duration of abstinence as independent variable, EPT as dependent variable, ReHo values in the left CAL and left PCC as mediators, and age as a covariate. The results indicated that both the left CAL and left PCC fully mediated the relationship between duration of abstinence and EPT in abstinent smokers ( Fig. 3 E ).
The result of whole-brain correlation analysis between postoperative analgesic requirements at 48 h and functional fMRI maps are shown in Figure 4, A and B , and Table 4 . In the abstinent smoker group, significant negative correlations were found between postoperative analgesic requirements and fALFF values in the left anterior cingulate cortex (ACC) and bilateral vmPFC ( Fig. 4 A ). Additionally, a significant positive correlation was found between postoperative analgesic requirements and rs-FC between the vmPFC with the bilateral medial/middle orbitofrontal cortex (OFC; Fig. 4 B ; Table 4 ). No significant correlations were detected for ReHo values. Similarly, no significant correlations were found between postoperative analgesic requirements and regional fMRI activity in the nonsmoking group.
Association between postoperative analgesic requirements and whole-brain rs-fMRI. A , Postoperative analgesic requirements at 48 h was negatively correlated with fALFF values in the left ACC and bilateral vmPFC. B , Postoperative analgesic requirements at 48 h was positively correlated with rs-FC of vmPFC to the bilateral medial and middle orbital frontal cortex (OFC). C , The scatterplots showed the correlation between MNWS with fALFF value in the right vmPFC, and left vmPFC and ACC, and between MNWS and postoperative analgesic requirements at 48 h. The shaded regions represent 95% credibility interval. D , The mediation effect of vmPFC and ACC in association with MNWS score and postoperative analgesic requirements at 48 h. A1, A2, effect of duration of abstinence on mediator. B1, B2, effect of mediator on EPT. C, total effect of duration of abstinence on EPT. C’, direct effect of duration of abstinence on EPT.
To further explore the relationship between smoking behavior and the brain regions associated with postoperative analgesic requirements in abstinent smokers, we conducted a correlation analysis. The results showed that the MNWS score was negatively correlated with the fALFF value in left ACC and bilateral vmPFC ( Fig. 4 C ). Moreover, the MNWS score was positively correlated with postoperative analgesic requirements ( Fig. 4 C ). Based on these findings, we hypothesized that increased postoperative analgesic requirements in abstinent smokers may be mediated by altered regional activity in the left ACC and bilateral vmPFC. To test this hypothesis, we conducted a mediation analysis, with the MNWS score as independent variable, postoperative analgesic requirements at 48 h as dependent variable, fALFF values in the left ACC and bilateral vmPFC as mediators, and age as a covariate. The analysis revealed that the left ACC and bilateral vmPFC fully mediated the relationship between MNWS scores and postoperative analgesic requirements in abstinent smokers ( Fig. 4 D ).
Conclusion
This study aimed at exploring the underlying mechanisms of nicotine withdrawal-induced hyperalgesia using fMRI in surgical patients. Clinical analysis showed that abstinent smokers exhibited lower pain thresholds and higher postoperative analgesic requirements than nonsmokers. These differences may be attributed to functional abnormalities of specific brain regions, including vmPFC, CAL, PCC, ACC, as well as rs-FC from vmPFC to MTG, PCUN, DLPFC, and medial/middle OFC. Further studies exploring the neurobiological mechanisms underlying the effects of nicotine abstinence in these brain regions are warranted.
Discussion
In this study, abstinent smokers exhibited lower pain thresholds, increased postoperative analgesic requirements, and altered neural activity—specifically lower fALFF in the vmPFC, higher ReHo in the left CAL, and decreased FC between vmPFC and bilateral MTG/PCUN. These findings are consistent with previous reports ( Garrison et al., 2017 ; Kim et al., 2017 ; Shen et al., 2017 ; Ditre et al., 2018 ; Wang et al., 2018 ; Zhao et al., 2018 ; Wei et al., 2023 ; Shepherd et al., 2024 ). Furthermore, abstinent smokers—unlike nonsmokers—showed a positive correlation between EPT and fALFF in the right CAL and between EPT and ReHo in the left CAL and PCC. Rs-FC of vmPFC was also positively correlated with MPT in right DLPFC. Pain thresholds (EPT) correlated positively with abstinence duration, an association fully mediated by CAL and PCC activity. Meanwhile, abstinent smokers showed a negative correlation between postoperative analgesic requirements and fALFF in vmPFC and left ACC and positively with rs-FC from vmPFC to bilateral medial/middle OFC. Withdrawal symptoms positively predicted postoperative analgesic requirements and mediated by vmPFC and ACC activity. While prior studies have described structural or functional neural changes in nicotine dependence and addiction ( Sutherland et al., 2016 ; Garrison et al., 2017 ; Sutherland and Stein, 2018 ; Wang et al., 2018 ; Zhao et al., 2018 ; Zhang and Volkow, 2019 ), this is the first study to link fMRI-based neural alterations to hyperalgesia in abstinent smokers accessed by both pain thresholds and postoperative analgesic requirements.
The default mode network (DMN)—a key resting-state network implicated in pain processing—typically encompasses three components: a midline core network, a medial temporal subsystem (MTL-DMN), and a dorsal medial prefrontal subsystem (dMPFC-DMN). The midline core network consists of the medial prefrontal cortex (mPFC) and PCC/PCUN ( Zhang and Volkow, 2019 ). We observed reduced activity and rs-FC in midline core network of DMN in abstinent smokers. Mediation analyses indicated that PCC and vmPFC activity contributed to elevated pain sensitivity and increased postoperative analgesic requirements, aligning with reports of DMN abnormalities in pain modulation ( Kucyi et al., 2014 ; Hsiao et al., 2017 ), nicotine addiction, withdrawal, and relapse ( Menossi et al., 2013 ; Sutherland et al., 2016 ; Zhang and Volkow, 2019 ; Yip et al., 2022 ; Zhang et al., 2024 ). The current study underscored the role of DMN dysfunction in nicotine withdrawal-induced hyperalgesia and suggested its potential as a therapeutic target for managing postoperative pain in abstinent smokers. Interestingly, different DMN components mediated heightened pain sensitivity before surgery and increased postoperative analgesic requirements, indicating distinct mechanisms. Indeed, the pain threshold measured in our study reflects acute, transient pain—elicited by the activation of nociceptive transducers without tissue injury—whereas postoperative pain results from substantial tissue damage ( Loeser and Melzack, 1999 ; Yarnitsky et al., 2013 ).
VmPFC, as the anterior midline core network in DMN, plays an important role in both nicotine withdrawal symptoms ( Cole et al., 2020 ; Huang et al., 2023 ) and pain modulation ( Ong et al., 2019 ; Pan et al., 2023 ), respectively. Short-term abstinence affects prefrontal activity related to cue reactivity ( Goldstein and Volkow, 2011 ), and neuromodulation of the DLPFC and vmPFC can alleviate withdrawal-related cognitive deficits ( Fischell et al., 2020 ). Agonism of nicotinic acetylcholine receptors improves cognition by modulating mPFC activity ( Sutherland et al., 2015 ). Animal studies further support that nicotine withdrawal alters mPFC structure and function, leading to withdrawal symptoms such as impaired cognitive memory, anxiety, and hyperalgesia ( Huang et al., 2022 , 2023 ; Shao et al., 2023 ; Gozen et al., 2024 ). The mPFC also regulates various types of pain ( Cheriyan and Sheets, 2018 ; Ong et al., 2019 ; Fang et al., 2023 ; Mosch et al., 2023 ; Lv et al., 2024 ), and its activity may even predict spontaneous pain ( Shirvalkar et al., 2023 ).
We also found reduced rs-FC between vmPFC and PCUN—a midline core network in DMN—in abstinent smokers. Notably, analysis of data from the “1000 Functional Connectomes” project revealed that the strongest FC hubs in the brain are concentrated in the midline core regions of the DMN, specifically encompassing the mPFC and PCC/PCUN ( Zhang and Volkow, 2019 ). The decrease of FC strength between mPFC and PCC/PCUN were associated with some of diseases such as schizophrenia ( Zong et al., 2019 ). Altered PCUN activity and connectivity are reported in chronic pain conditions ( Yu et al., 2023 ; Chen et al., 2025 ; Lai et al., 2025 ; Wu et al., 2025 ), and PCC/PCUN—an important hub in DMN ( Fransson and Marrelec, 2008 )—is linked to neuropathic pain ( J. Zhu et al., 2023 ) and headache ( Biggs et al., 2020 ; Ou et al., 2024 ). Our results found a decreased rs-FC between vmPFC and PCUN, and PCC mediated the association between duration of abstinence and pain thresholds in abstinent smokers, inferring that nicotine deprivation disrupts vmPFC—PCC/PCUN connectivity, contributing to hyperalgesia.
Beyond the DMN core, CAL and ACC also mediated the relationships between abstinence and pain measures. Traditionally, the CAL has been regarded as a principal structure within the primary visual cortex, primarily implicated in visual processing. However, according to the recent research, emerging evidence suggests that several components of the visual cortex, especially CAL, showed functional alterations or structural changes in various pain disorders, including endometriosis-related chronic pelvic pain ( Maulitz et al., 2024 ), primary dysmenorrhea ( Jin et al., 2024 ), brachial plexus avulsion-induced neuropathic pain ( J. Zhu et al., 2023 ), new persistent daily headaches ( Qiu et al., 2023 ), acute pain after total knee arthroplasty ( Kang et al., 2023 ), and chronic migraine ( Yuan et al., 2022 ). Furthermore, aberrant spontaneous activity in the CAL might underlie the mechanism of acupuncture analgesia in chronic stable angina pectoris ( Lan et al., 2022 ). Beyond somatic pain, patients suffering from chronic visceral pain (e.g., irritable bowel syndrome) showed higher spontaneous activity in CAL at resting-state, as well as heightened activation during pain anticipation ( Chen et al., 2021 ). As the discussion in Lan et al.'s paper ( Lan et al., 2022 ), these findings imply that the CAL may be involved not only in pain perception and pain affective processing but also in hypervigilance to pain anticipation. Another recent interesting finding was that the FC from CAL to the right nucleus accumbens (one of reward-related brain regions) was increased after repetitive transcranial magnetic stimulation in nicotine-dependent participants, indicating CAL may be involved in nicotine dependence and withdrawal ( Wang et al., 2024 ). ACC is similarly implicated in nicotine addiction ( Sutherland and Stein, 2018 ; Keeley et al., 2020 ) and withdrawal symptoms ( Fischell et al., 2020 ), especially hyperalgesia ( Long et al., 2023 ; Chen et al., 2024 ). Short-term nicotine deprivation alters dorsal ACC glutamate concentration and concomitant cingulate-cortical FC ( Abulseoud et al., 2020 ). These studies, combined with our research, suggest that CAL and ACC may play a certain role in nicotine withdrawal-induced hyperalgesia. More definitive evidence and its neural molecular mechanisms should be further studied in future research.
Notably, electrical (EPT) and mechanical (MPT) pain thresholds correlated with distinct neural markers: EPT with CAL and PCC activity and MPT with vmPFC–DLPFC connectivity. In modern physiology, three subtypes of sensory nerve fibers including Aβ, Aδ, and C fibers have been identified ( Lin et al., 2020 ). The low-intensity electrical stimuli primarily activate Aβ afferent fibers ( Yague et al., 2010 ), while the mechanical stimuli primarily activate Aδ and C fibers ( Lee et al., 2011 ). No such correlations occurred in nonsmokers, suggesting abstinence-specific neural dysfunction in CAL, ACC, and fronto-cingulate pathways underlying hyperalgesia.
Several limitations should be noted. First, only male patients were included. Nicotine withdrawal symptoms vary between the sexes ( Sussman et al., 1998 ; Allen et al., 2018 ; Wen et al., 2022 ), with women generally reporting greater pain sensitivity than men ( Bartley and Fillingim, 2013 ). Further studies should include females. Second, this study focused solely on fALFF, ReHo, and rs-FC. Future studies should investigate other neural dynamics, particularly those underlying connectivity among canonical brain networks. Third, as an observational study, mediation analyses suggest but cannot confirm causality—fMRI changes may mediate pain behaviors, but reverse or third-factor causation remains possible. Intervention studies are needed to establish causal links. Finally, given the limited sample size of the current study, it would be valuable to conduct future research with a larger cohort to validate the results.
Significance
Abstinent smokers experienced heightened pain and require more analgesics after surgery, yet the underlying neural mechanisms remain poorly understood. This prospective cohort study identified altered regional brain activity associated with reduced pain thresholds and increased postoperative analgesic requirements in abstinent smokers. We found specific brain regions that were functionally altered and correlated with pain-related outcomes, which mediated the relationship between abstinence and pain-related behaviors. These findings provided novel insights into the supraspinal mechanisms of nicotine withdrawal-induced hyperalgesia and point to potential therapeutic targets for improving postoperative pain management in abstinent smokers.
Materials|Methods
This observational study was conducted at Shanghai Eastern Hepatobiliary Surgery Hospital, China. The study protocol was approved by the Institutional Ethics Committee (Shanghai Eastern Hepatobiliary Surgery Hospital; approval number: EHBHKY2017-03-007) and registered with the Chinese Clinical Trial Registry (ChiCTR-ROC-17011004, https://www.chictr.org.cn ) on March 28, 2017. Written informed consent was obtained from all participants before enrolment. All procedures were performed in accordance with the Declaration of Helsinki. Male patients scheduled for partial hepatectomy due to liver masses were assigned to either the abstinent smoker group or nonsmoker group on the day before surgery.
Patients were eligible for inclusion if they met the following criteria: (1) male; (2) aged 18–60 years; (3) right-handed; (4) no history of substance abuse other than nicotine; (5) no history of brain lesions or trauma; (6) no communication barriers and capable of completing all assessments; and (7) classified as American Society of Anesthesiologists Grade I or II and Child–Pugh Class A for liver function. Exclusion criteria included (1) contraindications to MRI (e.g., nonremovable metallic implants, claustrophobia); (2) history of peripheral or neurological disorders; (3) chronic pain conditions; (4) long-term use of neuroactive medications (e.g., sedatives, opioids, antidepressants); (5) body mass index >28 or <18 kg/m 2 ; (6) inability to complete psychophysiological assessments; (7) contraindications to patient-controlled analgesia (PCA); or (8) inability to complete the 48 h postoperative follow-up, including due to severe postoperative complications. Patients who did not meet the predefined criteria for abstinent smokers or nonsmokers were also excluded.
The study workflow is illustrated in Figure 1 . Patients were assigned to either the nonsmoker or abstinent smoker group based on their smoking status. In brief, nonsmokers were defined as individuals who had never smoked or had smoked <100 cigarettes in their lifetime. Evidence from our previous work showed that smokers who have abstained for <1 month exhibit lower pain thresholds ( Shen et al., 2017 ). Consistently, Hooten et al. reported that current smokers required higher opioid doses at admission to complete a pain rehabilitation program than former smokers who had abstained for at least 1 month and nonsmokers ( Hooten et al., 2008 ). As part of standard preoperative care, all smokers scheduled for elective surgery were asked to quit smoking to minimize postoperative complications. Consequently, abstinent smokers in this study were defined as patients who had regularly smoked at least 10 cigarettes per day over the past year and had quit within 1 month before screening. Individuals abstinent for greater than 1 month were excluded. We also excluded patients who used pharmacological cessation aids—such as nicotine replacement therapy or varenicline—during the perioperative period, because our earlier research suggests that nicotine replacement can alleviate nicotine withdrawal-induced hyperalgesia ( C. Zhu et al., 2023 ).
Study flow chart. MRI, magnetic resonance imaging; PCA, patient-controlled analgesia.
A team of three independent researchers collected all experimental data 1 d before surgery. The first researcher recorded baseline information and conducted all psychophysiological assessments. The second researcher, blinded to group assignment, performed the pain threshold assessment. The third researcher, also blinded to group assignment, collected the fMRI data with the assistance of technicians. On the following day, all patients underwent hepatectomy under standardized general anesthesia. A standardized postoperative analgesic protocol was applied to all patients to allow for comparison of analgesic consumption. The anesthesiologists responsible for perioperative care were blinded to group assignments.
All patients completed a series of psychophysiological assessments 1 d before surgery, based on their current feelings. Specifically, state and trait anxiety were evaluated using the State–Trait Anxiety Inventory, and depressive symptoms were measured using the Beck Depression Inventory. Pain sensitivity was evaluated with the Pain Sensitivity Questionnaire (PSQ), which reflects an individual's natural disposition to perceive pain. Additionally, abstinent smokers completed three smoking-related questionnaires. The Fagerstrom Test for Nicotine Dependence (FTND) was used to evaluate the severity of nicotine dependence based on preabstinence smoking behavior. Current withdrawal symptoms and nicotine cravings were assessed using the Minnesota Nicotine Withdrawal Scale (MNWS) and the Questionnaire of Smoking Urges (QSU), respectively ( Cho et al., 2024 ).
Our previous study ( Shen et al., 2017 ) demonstrated that abstinent smokers exhibited lower pain thresholds when assessed using electrical stimulation. Based on this finding, we employed an electrical stimulator to measure the electrical pain threshold (EPT) 1 d before surgery. To further validate the increased pain sensitivity in abstinent smokers, we assessed mechanical pain threshold (MPT) using von Frey filaments on the same day, as this method is convenient, clinically applicable, and well tolerated. Both assessments were noninvasive, time-efficient, and suitable for preoperative evaluation. Specifically, electrical stimuli were delivered using a constant current stimulator (model DS7A; Digitimer) through a pair of electrodes placed on the index finger of the left hand (interelectrode distance, 1 cm; pulse width, 2 ms). The method of limits (in steps of 0.1 mA) was applied to determine the stimulus intensity at which the patient started to feel pain. This procedure was repeated three times, and the mean stimulus intensity was defined as the EPT for each patient. Mechanical stimuli were delivered to the dorsum of the left hand using Semmes–Weinstein monofilaments (Touch-Test 20 Piece Kit; North Coast Medical), with evaluator sizes ranging from 1.65 to 6.65 and target forces from 0.008 to 300 g ( Keizer et al., 2007 ). The filament was pressed against the skin at a 90° angle until it bowed and then held in place for 1.5 s. The method of limits was used to determine the evaluator size at which pain was first reported (increasing sequence) and subsequently lost (decreasing sequence). This process was repeated three times, and the geometric mean of the evaluator size values was defined as the MPT for each patient.
We measured the cotinine levels in both plasma and urine to verify the presumed difference between the two groups using previously reported methods ( Shen et al., 2017 ). Specifically, 1 d before surgery, blood and urine specimens were obtained, labeled, centrifuged at room temperature, and promptly stored at −80°C. After study completion, all samples were transferred to the pharmacy department for analysis. Cotinine levels in the serum and urine were determined using ultrapressure liquid chromatography–tandem mass spectrometry (1290 Infinity II series LC system, G6470 triple quadrupole mass spectrometer; Agilent Technologies), following an optimized protocol established in a previous study ( Wei et al., 2014 ). The limit of quantification of the method was 0.1 ng/ml. The accuracy and precision of the method were within 8% of the relevant concentration range.
All patients underwent hepatectomy with standardized general anesthesia, as reported previously ( Shen et al., 2017 ; Wei et al., 2023 ). In brief, tracheal intubation was implemented after anesthesia induction with intravenous midazolam (0.03 mg/kg), target-controlled infusion of propofol (target plasma concentration of 5.0 µg/ml achieved using the Paedfusor model), sufentanil (0.3 µg/kg), and cis -atracurium (0.2 mg/kg). Anesthesia was maintained with propofol (target plasma concentration of 3.0–5.0 µg/ml), remifentanil (0.1–0.2 µg/kg/min), and cis -atracurium (1.5 µg/kg/min) to maintain the bispectral index in the range of 40–60 and the heart rate and blood pressure to within 20% of preoperative values. Mechanical ventilation was adjusted to maintain the end-tidal partial pressure of carbon dioxide to between 35 and 45 mm Hg, with a tidal volume of 8–10 ml/kg and a respiratory rate of 10–14 breaths/min.
Postoperative analgesia was administered in a manner consistent with our previous report ( Shen et al., 2017 ). Following anesthesia, patients were transferred to a postanesthesia care unit for ongoing care. All patients received sufentanil (1 µg/ml; Humanwell Pharmaceutical) via a PCA device. The PCA device was programmed to deliver 2 ml of analgesics as an intravenous bolus, with a lockout time of 10 min and a background infusion rate of 2 ml/h. The maximum allowable sufentanil dosage was capped at 10 µg/h. If patients exhibited poor response to sufentanil, supplemental rescue analgesics were used to relieve postoperative pain. The analgesics administered in this study were injectable morphine, meperidine, and oxycodone. The doses of these opiate analgesics were converted into morphine-equivalent doses based on equianalgesic dose ratios.
The fMRI data were collected using a 3.0 Telsa MRI system (Discovery MR 750; General Electric Healthcare) with an eight-channel receive–only radiofrequency head coil array 1 d before surgery. A standard birdcage head coil, along with restraining foam pads, was used to minimize head motion and scanner noise. A T1-weighted structural image was acquired using a three-dimensional spoiled gradient recalled echo sequence with the following parameters: repetition time, 6.896 ms; echo time, 2.99 ms; flip angle, 8°; in-plane matrix size, 256 × 256; 176 slices; field of view, 256 × 256 mm 2 ; in-plane resolution, 1 × 1 mm 2 ; and slice thickness, 1 mm. The rs-fMRI images were acquired using an echoplanar imaging sequence using the following parameters: repetition time, 2,000 ms; echo time, 30 ms; flip angle, 90°; field of view, 64 × 64 mm 2 ; data matrix, 64 × 64; in-plane resolution, 3 × 3 mm 2 ; slice thickness, 3.5 mm; slice spacing, 0.5 mm; and 256 volumes. The rs-fMRI scan took ∼10 min. During the rs-fMRI data acquisition, patients were instructed to keep their eyes closed, remain still, and not think about anything systematically.
Patient characteristics, including age, height, weight, and smoking-related behaviors, were obtained via standardized questionnaires. Medical information such as blood test results, anesthesia records, and surgical data was retrieved from electronic medical records. Postoperative analgesic requirements at 1, 6, 24, and 48 h after surgery were recorded by a trained research assistant who was blinded to group assignment. To account for variations in the body size, all analgesic doses were normalized by dividing the total dose by the patient's body weight.
The rs-fMRI data were analyzed using Data Processing and Analysis for (Resting-State) Brain Imaging (DPABI; version V8.2_240510, https://rfmri.org/DPABI ) toolbox ( Yan et al., 2016 ), implemented in MATLAB 2023b (MathWorks). The first 10 volumes were discarded to minimize the effects of spin saturation and participant adaptation to the scanning environment. The remaining 246 volumes were slice-time corrected and realigned for head motion using a six-parameter rigid body transformation. Patients with head movements exceeding 2.0 mm in translation or 2.0° in rotation were excluded. Individual structural images were coregistered to the mean functional image and segmented into gray matter, white matter, and cerebrospinal fluid. Nuisance covariates—including white matter signal, cerebrospinal fluid signal, and head motion scrubbing parameters (Friston 24)—were regressed out from the fMRI data. Please note that, time points with a framewise displacement (Jenkinson FD) exceeding 0.2 were classified as “bad.” The scrubbing regressor was defined to include the identified bad time points, along with one preceding and two subsequent time points. The processed fMRI data were then normalized to standard space using DARTEL (resampling voxel size, 3 × 3 × 3 mm 3 ) and linearly detrended to remove drift. The resulting data were used for subsequent analyses.
The DPABI toolbox was used to compute fALFF, ReHo, and rs-FC from the preprocessed rs-fMRI data. For fALFF calculation, data were spatially smoothed using a 4 mm full-width at half-maximum (FWHM) Gaussian kernel. The relative amplitude of low-frequency fluctuations within the 0.01–0.1 Hz range was then calculated. Standardized fALFF values were obtained by normalizing each voxel's fALFF to the global mean and were used for subsequent statistical analysis. For ReHo, Kendall's coefficient of concordance (KCC) was used to measure the local synchronization of the time series between each voxel and its 26 neighboring voxels, after bandpass filtering (0.01–0.1 Hz). The resulting ReHo values were standardized by dividing each voxel's KCC by the mean KCC of the whole brain and then spatially smoothed with a 4 mm FWHM Gaussian kernel.
Brain regions showing significant group differences in fALFF and ReHo were defined as regions of interest (ROIs) for ROI-based rs-FC analysis. The preprocessed data were smoothed with a 4 mm FWHM Gaussian kernel and bandpass filtered (0.01–0.1 Hz). Pearson's correlation coefficients were calculated between the mean time series of each ROI and those of all other brain voxels. These correlation maps were then transformed to Fisher's Z scores (zFC maps) for further statistical analyses.
The statistical analysis plan was determined before participant enrolment. Descriptive statistics are presented as means (standard deviations), medians (interquartile ranges), or frequencies (proportions), depending on the data type and distribution. For continuous variables, normality was confirmed using the Kolmogorov–Smirnov test. Group comparisons of continuous variables were conducted using Student's t test or the Mann–Whitney U test, as appropriate. Repeated measures were compared using repeated-measure analysis of variance (ANOVA), with post hoc tests performed when a significant interaction effect was detected. Categorical variables were compared using the χ 2 test. To explore the associations between smoking behaviors (i.e., cigarettes/day, smoking duration, abstinence duration, FTND, MNWS, QSU, and cotinine levels) and pain thresholds or postoperative analgesic requirements, we performed Pearson's or Spearman’s correlation analyses in abstinent smokers, selected based on data distribution. The level of statistical significance was set at p < 0.05. All statistical analyses were performed using SPSS 22.0 (SPSS).
To examine group differences in fALFF, ReHo, and FC, two-sample t tests were performed using DPABI in MATLAB 2023b, with age and head motion (mean framewise displacement) included as covariates. Correlation analyses were used to explore associations between whole-brain regional fMRI maps and pain thresholds or postoperative analgesic requirements in both abstinent smokers and nonsmokers. Multiple comparisons were conducted using Gaussian random field correction [voxel-level, p < 0.001; cluster-level, p < 0.05 (two-tailed)] for the fALFF, ReHo, rs-FC, and correlation analyses. Furthermore, mediation analyses tested whether pain-associated brain regions mediated the effects of smoking behaviors on pain thresholds and postoperative analgesic requirements.
Text is read by the "Ask this paper" AI Q&A widget below.
Extraction quality varies by source — PMC NXML preserves structure
cleanly, OA-HTML may include some navigation residue, and OA-PDF can
have broken hyphenation. The publisher copy
(via DOI)
is the canonical version.