Skip to main content

The role of the sympathetic component of the autonomic nervous system on pain before and after third molar extraction– an observational cohort study

Abstract

Background

Assessing heart rate variability (HRV) before a standardized surgery would help to explore further the relationship between the autonomic nervous system and pain.

Methods

A single-center prospective cohort of 117 patients (55% female) scheduled for third molar extraction underwent a preoperative resting measurement of arterial pressure followed by an HRV recording, then potentiated by a Valsalva maneuver and a deep breathing challenge. Finally, pain sensitivity was assessed by hand immersion in hot water. All surgeries were conducted under local anesthesia, with or without sedation. The primary outcome was a composite pain/analgesia score (CPAS) incorporating both pain intensity and analgesic drug intake; it was adjusted to the type of anesthesia by within-subgroup ranking.

Results

The increase in heart rate in the Valsalva maneuver, and the low- to high-frequency ratio (LF/HF) in the deep breathing, were inversely correlated to preoperative heat pain, which was correlated itself to the CPAS (ρ = 0.195; p = 0.035). The only other parameter influencing CPAS was the increase in heart rate in the Valsalva maneuver, with an inverse correlation (ρ = − 0.191; p = 0.046). While age tended to impair HRV, particularly in its parasympathetic component, and while men displayed a stronger parasympathetic response than women, neither age nor sex interacted with these effects. Neither preoperative arterial pressure nor the occurrence of parental hypertension influenced the pain outcomes.

Discussion

Although the identified relationships were not particularly strong, they are consistent with an influence of the sympathetic component of the autonomic nervous system. However, they do not support the interest of HRV assessment to predict postoperative pain in current practice.

Trial registration

Not applicable.

Peer Review reports

Introduction

Both nociception and the autonomic nervous system (ANS) participate in the response of the body to a potentially harmful threat, the sympathetic response being an essential element of the behavioral reaction (“fight-or-flight”) to threat. Literature data also suggest interactions between those two systems [1]. For example, hypertension in humans is often linked to a reduced sensitivity to pain [2,3,4], and attempts to increase arterial pressure or sympathetic tone may reduce pain sensitivity [5,6,7,8]. However, it is still unknown (i) whether this analgesic effect follows an independent central or descending nervous pathway, or depends on the autonomic response; and (ii), if dependent, whether it results directly from hemodynamic changes, or from a hormonal stress response. In addition, while the stress response is mostly driven by the sympathetic response, there is also a parasympathetic response to threat [9], also possibly interacting with pain sensitivity [10,11,12].

Heart rate variability (HRV) is considered as a more reliable indicator of the ANS than the sole arterial pressure or heart rate; it has long been used to diagnose autonomic dysfunction [13]. The wide availability of low-cost non-invasive bedside devices has broadened this field of research [14, 15]. HRV outcomes are based either on descriptive statistics of the beat-to-beat intervals, or on fast Fourier transform to determine a frequency spectrum (see Table 1). These analyses help to estimate the balance between the sympathetic and parasympathetic components of the ANS, although their accuracy– especially for the sympathetic– have been challenged [16]. Nevertheless, the accuracy of HRV to assess the ANS function can be enhanced by standardized tests such as the Valsalva maneuver (VM) and the deep breathing (DB) challenge [13].

Table 1 Selected measures of heart rate variability (HRV)

Relationships between the autonomic response and pain sensitivity have mostly been studied in young healthy volunteers without disturbance to the ANS, and who ethically cannot be subjected to anything worse than transient painful stimuli. Thus, because it involves a tissue damage that, for ethical reasons, does not exist in experimental pain, postsurgical pain appears to be of the best models to study a physiological determinant of pain, as long as the surgical procedure and postoperative analgesia are homogeneous. As it is commonly practiced and often used as a model for testing analgesic drugs [17, 18], we chose to use extraction of the third molar (“wisdom tooth”) to study the effects of the preoperative status of the ANS on postoperative pain. As patients are generally young and healthy, pain outcomes are also likely to be less biased by comorbidities and concomitant medications, and this condition therefore more reliable to explore physiological aspects. HRV was tested in a resting condition, then during VM and DB. In addition, to understand the mechanisms underlying a relationship between HRV and postoperative pain (if any), the patients were also preoperatively submitted to an acute pain test by hand immersion in hot water, to assess their natural pain sensitivity [19, 20]. Finally, if we could identify easy-to-record predictors of postoperative pain, they could be used in the future to help devise personalized analgesia [21,22,23].

Materials and methods

All patients were recruited from the Dental Surgery unit of the University Hospital of Clermont-Ferrand (France). The inclusion criteria for the patients were that they had to be adult, aged 18 or over, ASA I-II, naive from previous 3rd molar extraction, and eligible for outpatient extraction of either two ipsilateral 3rd molars under local anesthesia without sedation (“two wisdom teeth” = 2WT subgroup), or the four 3rd molars under local anesthesia plus intravenous sedation (“four wisdom teeth” = 4WT subgroup). The mandibular tooth/teeth had to be impacted, i.e. completely covered by mucosa and at least partially impacted against bone; the maxillary one(s) had to be either similarly impacted, or impacted in the arch. The exclusion criteria were: weight under 50 kg; difficult access to the teeth to be extracted (e.g. limited mouth opening, temporomandibular joint dysfunction); pregnancy; high gag reflex; need for a regional anaesthesia of the mandibular nerve; contraindication to any of the drugs used in the protocol; any disease likely to impact the surgical outcomes (such as diabetes, renal insufficiency, haemostatic disturbance, autoimmune disease, or immunodeficiency); any ongoing or chronic disease likely to impact the study outcomes such as: orofacial chronic pain (e.g. due to endodontic, periodontal, gingival, osseous, or radiotherapy-induced disease), other invalidating chronic pain, or psychiatric disease; declared risky alcohol drinking (more than 2 standard drinks a day for men and more than 1 for women) or suspected alcohol use disorder; smoking more than 10 cigarettes or equivalent a day; declared or suspected use of recreational drugs; or daily ingestion of more than 4 units of caffeine-containing beverage. In addition, medications likely to impact the study outcomes were exclusion criteria; this included opioid daily intake during the past 3 months, or any intake during the past 7 days of: analgesic drugs (e.g. opioids, NSAIDs or steroids), drugs interacting with the sympathetic system (e.g. alpha/beta agonists or antagonists), drugs interacting with the parasympathetic system (e.g. muscarinic agonists or antagonists, including phenothiazines or antidepressants). Occasional intakes of weak analgesic drugs prior to 7 days before surgery were allowed.

A preoperative visit was conducted 7 to 14 days before surgery. The demographic and morphometric data were collected and the functional explorations were conducted in a quiet room specifically reserved for this purpose. Included in the demographic data noted was a parental hypertension score, based on the number of cases of hypertension amongst the parents and the grandparents, with 0 for no cases, 1 for one case, 2 for two cases, and 3 for more.

Firstly, the arterial pressure was first measured with a semi-automated device (Dinamap ProCare 300, General Electric, Boston, MA, USA) after 5-min rest in a semi-sitting position. Then, HRV was monitored by recording a continuous electrocardiographic signal (sampling frequency: 1000 Hz) transferred to a Powerlab system equipped with a LabChart application (ADInstruments, Colorado Springs, CO, USA), with respective low- and high-pass filters at 1 kHz and 0.3 Hz and a pitch at 60 Hz. The circuit was earthed via an electrode placed on the subject’s knee. An initial 5-min recording of HRV at rest was carried out on the patient lying supine. Then further 5-min recordings were conducted during two consecutive challenges with the patients back in a semi-sitting position. Firstly, they underwent a series of ten Valsalva maneuvers (VM), at the rate of one every 30 s; each VM consisted of exhaling forcibly through a manometer against a pressure of 40 mmHg for 15 s. After a resting time of 5 to 10 min (i.e. long enough to return to close to the pre-challenge arterial pressure and heart rate), the patients underwent a series of six deep breathing maneuvers (DB), at the rate of one per minute [13, 24]; each deep breath lasted 10 s, with 5 s of inspiration and 5 s of expiration. A tachogram was obtained with measurements of the normal-to-normal (NN) intervals for these three conditions (baseline, VM and DB). The calculated parameters are listed in Table 1, along with their definitions and interpretations [14, 25]. The HRV measurements were totally machine-dependent, and the data collection and management was made unaware of the postoperative outcomes.

Immediately after the HRV assessment, the patient was subjected to a heat pain challenge, in which he/she had to put the non-dominant hand in a bath of water at 47 °C (SW22, JULABO GmbH, Seelbach, D) for a targeted maximum of 5 min. Every minute, he/she was asked to rate the heat-induced pain on a visual analogue scale (VAS) ranging from 0 (no pain) to 100 (the worst possible pain) and the values were averaged to obtain the heat pain intensity (HPI). In the event of the hand being withdrawn before completion of the 5-min protocol, a value of 100 was attributed to any subsequent observations that were missing.

In both groups, avulsions were conducted under local anesthesia by buccal and lingual infiltration at the operation sites. In addition, patients in the 4WT subgroup received intravenous sedation administered by a senior anesthetist unaware of the protocol. The current sedation protocol included as a first-line a single dose of i.v. remifentanil, to which was added if necessary a single i.v. dose of either midazolam or propofol. All surgeries were performed by one of the three co-authoring surgeons (CDe, LD and YS), all trained in a standardized technique at the same dental school. Immediately after surgery, the surgeon recorded information about the conduction of anesthesia and surgery.

At discharge, the patient was given a typed prescription of analgesia for 5 days, along with an instruction for self-administration. We followed the protocol currently used in in our unit, which differed between subgroups for the molecule of nonsteroidal anti-inflammatory drug (NSAID) only. Acetaminophen was systematically prescribed (1 g first thing in the morning, at midday and last thing at night). In addition, a NSAID could be taken as a first-line rescue treatment, preferentially with food; this was ibuprofen (200 mg per tablet, 1–2 tablets depending on pain intensity, up to a maximum of 8 tablets/24 hrs) in the 2WT subgroup, and ketoprofen (50 mg per tablet, 1–2 tablets depending on pain intensity, up to a maximum of 4 tablets/24 hrs) in the 4WT subgroup. As a second-line rescue treatment, tramadol could be taken (50 mg per tablet, up to a maximum of 3 tablets/24 hrs). All patients were given a booklet to note down postoperative pain and the intake of analgesic drugs. Pain was assessed on the evening of the day of surgery, then every morning (when getting up) and evening (when going to bed) of the 2nd to the 5th postoperative day, by a pen mark on a 10-cm VAS paper with the same scale as described above. The analgesic drugs taken were noted in a table which had a number of rows corresponding to the number of hours from discharge to the end of the 5th postoperative day, and one column per drug, plus one for any other drug than those prescribed. For each intake, the number of tablets had to be noted. The patient was unaware of the results of the HRV and experimental pain challenge, and the investigators had no contact with the patient throughout the postoperative observation period.

As postoperative analgesia was semi-directed, the primary outcome for representing postoperative pain was a composite score based on both the pain intensity and analgesic drug intake. This followed the Silverman integrated approach, which aims at correcting the bias induced by analgesic treatments on pain by summing pain and analgesia [26]. We therefore called “composite pain/analgesia score” (CPAS) this primary outcome. CPAS is the sum of a total pain score and a total analgesia score, each representing respectively the amount of pain and analgesic drug intake throughout the postoperative observation period. Firstly, to calculate the total pain score, we summed (i) the averaged pain score (i.e. the mean of the nine scores measured over the five postoperative days) and (ii) the maximal pain score reported during this period. This maximal value was considered as equally important than the average pain score because it indicates the patient’s preferences of analgesic target [27]; summation was not weighted as both values covered the same range. The relevance of this summation-based outcome had been tested from personal data [4]. Then, this sum obtained for each individual, was centered and reduced by the following transformation: (ranked value within the sample) / [(mean theoretical rank) − 1], with mean theoretical rank = (sample size + 1) / 2. The rank within the sample was calculated in decreasing order. In addition, to control the selection bias as well as the bias induced by the additional sedation or the different NSAID prescribed, the ranking was conducted within each surgery subgroup (2WT/4WT). Secondly, to calculate the total analgesia score, we defined for each individual and each day of observation a parameter of analgesic drug intake. This parameter was devised based on a literature review of the reported number needed to treat (NNT, calculated vs. placebo) for 50% pain relief after a 3rd molar extraction, during the 4th -6th postoperative hours [17, 18, 28,29,30,31,32,33,34,35,36,37,38,39,40,41]. Based on the NNT for each drug dose or combination, an analgesic drug score ranging from 0.5 to 4 was given to each observation of analgesic intake (see Table 2 for details). Then, the sum of the five daily analgesia sub-scores was transformed as done for the total analgesia score. This way, the two components of the CPAS, the total pain score and the total analgesia score, had the same natural range (0 to 2), and CPAS as their sum naturally ranges from 0 to 4. As neither pain nor analgesic drug consumption follow a Gaussian distribution, transforming the data into ranks reduced the random effects of distribution and harmonized the two components in terms of shape and range. An example of calculation is shown in Fig. 1.

Table 2 Calculation of the analgesic drug score
Fig. 1
figure 1

Example of calculation of the composite pain/analgesia score (CPAS) for one subject of the study sample, who belonged to the “4 wisdom teeth” extraction (4WT) subgroup. How the two components of CPAS are calculated is further detailed in the Methods’ chapter

Data analysis was performed using Microsoft Office Excel 2010 (Redmond, WA, USA) and XLStat (Addinsoft, Paris, F). Descriptive statistics were conducted for each variable, and normality was checked by visualization of Q-Q plots. Numerical data were expressed by mean ± SD in the case of a normal distribution, and otherwise as quartiles. Nominal data were expressed by the number of cases and the percentage. The analyses were conducted on full datasets; if the information about HRV, pain or analgesia missed for one patient, then the case was excluded from the analysis. The primary endpoint was to test the influence of a set of hemodynamic and HRV variables on postoperative pain (assessed by the CPAS). The tested variables were the raw systolic, diastolic and mean arterial pressure, the heart rate, and the raw HRV parameters listed in Table 1, for each of the three conditions (baseline, VM and DB), as well as the percentage changes from the baseline for arterial pressure and heart rate during VM and DB. Those percentage changes were calculated by dividing the raw change (difference) by the baseline value. We did not calculate the percentage changes for the HRV parameters because the baseline value was equal to 0 for some of the observations. The main secondary endpoint was to test the influence of these explanatory variables on the intensity of pain induced by the preoperative challenge (HPI). As both the CPAS and heat-induced pain are non-Gaussian numerical variables, we used Spearman’s correlation to study the effect of numerical variables; the significance of the correlation was tested by comparing the ρ coefficient to the null value and its strength was expressed by the coefficient of determination, R2. The same calculation was done to test the relationship between the HPI and the CPAS. In addition, to study the effect of a binary variable on a numerical outcome, comparisons were carried out using either the Student’s t-test or the Mann-Whitney test depending on the distribution. To control the quality of the HRV measurement, we also conducted within-measurement principal component analyses (PCA) on Spearman’s correlation matrixes, and analyzed the effect of VM and DB for repeated measurements. For this we used either the ANOVA-RM or the Friedman’s test depending on the distribution, followed by post hoc tests if needed (respectively, the Student’s t-test for paired data or the Wilcoxon signed rank test). If multivariable analyses were needed (e.g. to test interactions), a generalized linear model fitted to the type of variables was used. For inferential analyses, the significance threshold was set at 5%. As we wished to know whether either the sympathetic or the parasympathetic, or both, had some influence, and based on the hypothesis is that, although both contributing to a balance, those two components of the ANS act independently [42], the type-I error inflation was not corrected for the separate analyses of these two components. Although there are many different outcomes in HRV analysis, most of them explore a unique domain (especially the parasympathetic one), so the results must be more read by means of internal coherence. Nevertheless, within each domain, depending on the level of multicollinearity between parameters and on the general trends (e.g. whether an effect was observed for several parameters or not), complementary Bonferroni’s corrections were conducted to check the reliability of our results. For the analysis of repeatedly measured outcomes, the preliminary trend test avoided the type-I error inflation.

We made an initial estimate of the sample size in order to allow for the analysis of 12 putative factors, resulting in an estimation of 104 + 12 = 116 subjects. We based this assumption on a rule for multivariable analysis [43, 44], although the primary endpoint was based on univariate analyses. Our objective was to include 120 patients, considering that this sample size allowed to identify a correlation of 0.3 with 1–β = 90%, and a correlation of 0.25 with 1–β = 80% (with α = 5%) [45, 46].

Results

One hundred and twenty-one patients were included in the study, of which four were unable to undergo the HRV assessment; of the 117 who completed the study, 50 (42.7%) were in the 2WT and 67 (57.2%) in the 4WT subgroup. The characteristics of these patients are presented in Table 3, along with the differences between subgroups. The only noticeable differences are a slightly younger age, a higher maximal Pederson score, and a longer duration of surgery in the 4WT subgroup.

Table 3 Characteristics of the cohort and its two subgroups

The hemodynamic and HRV data for the whole cohort are presented in Figs. 2 and 3, and Additional file 1 for details. VM increased systolic arterial pressure and heart rate, while DB decreased systolic arterial pressure and increased the heart rate; the increase in heart rate was significantly greater under VM (+ 11.6%) than under DB (+ 1.6%). Diastolic arterial pressure was unaffected. Both challenges increased the time-domain parameters of HRV, but these effects were stronger under VM for SDNN, while they were stronger under DB for NN50 and pNN50. Also, both challenges increased the frequency-domain parameters of HRV, but these effects were stronger under VM for total power, while they were stronger under DB for LF and LF/HF. The PCA conducted on all the HRV parameters showed that all the time-domain parameters, along with total power and HF (for the frequency-domain) were collinear and highly represented on the same axis, while LF/HF was highly represented respectively on the second axis. When the challenge-induced change in heart rate was added to the models, it was mostly represented on a third axis, but also on the same axis as LF/HF (Additional file 2). This supports the hypothesis that most of the parameters represent the same domain, theoretically the parasympathetic component of the ANS, while only LF/HF and change in heart rate represent another one, theoretically the sympathetic component.

Fig. 2
figure 2

Summary of the effects of the two challenges, namely the Valsalva maneuver (VM) and the deep breathing challenge (DB) on the hemodynamic parameters of the study sample. The variables displayed on the graph are the percentage of change from the baseline values; the black squares represent the mean value, and the whiskers represent the 95% confidence interval (CI) limits. When the null value is out of CI limits, this corresponds to a significant difference (p < 0.05) compared to baseline; when the respective CI limits of the variable for each challenge do not overlap, this corresponds to a significant difference between the two challenges. Those significances are without Bonferroni’s correction. The aggregated data are detailed in the Supplementary Table 1

Fig. 3
figure 3

Summary of the effects of the two challenges, namely the Valsalva maneuver (VM) and the deep breathing challenge (DB) on the heart rate variability (HRV) parameters of the study sample, in comparison to the baseline values (BL). The variables displayed on the graph are box plots of the raw values. The interpretation of the HRV parameters is detailed in Table 1. The symbols ‘#’ and ‘$’ signal a significant difference (p < 0.05), respectively with baseline and between VM and DB (see Methods’ chapter for detailed description of the statistics). The aggregated data are detailed in the Supplementary Table 1

The two subgroups (2WT/4WT) behaved similarly (Additional file 3), but heart rate was higher at the baseline in the 4WT subgroup (p = 0.041), and increased less in this subgroup under both challenges (p = 0.001 and 0.003 respectively), in such a way that the raw values became similar to those of the 2WT subgroup. Also, the slight decrease in diastolic arterial pressure under VM was only observed in the 2WT subgroup.

The effects of age on the hemodynamic and HRV parameters are detailed in Additional file 4. Age was positively correlated to both diastolic and mean arterial pressure; it was negatively correlated to the change in systolic arterial pressure during VM (i.e. VM-induced hypertension decreased with age) and positively correlated to the change in mean arterial pressure during DB (i.e. DB-induced hypotension increased with age). Age did not influence heart rate or its changes during VM or DB. In the time-domain, age was negatively correlated to all the baseline parameters for HRV, as well as to all those parameters measured during DB; none of these correlations were observed for measurements during VM. In the frequency-domain, similar trends were observed for total power and HF at baseline and during DB. Conversely, LF/HF tended to increase with age, and this correlation was significant during VM.

The effects of sex on the hemodynamic and HRV parameters are detailed on Additional file 5. For most of the measurements of arterial pressure, higher values were observed in men than in women, and the DB-induced change (i.e. decrease) in mean arterial pressure was greater in men. In the time-domain analysis of HRV, SDNN and RMSSD were higher in male, but this only under VM, and in the frequency-domain, higher values in male under VM were also observed for total power and HF. The LF/HF ratio was unaffected by sex.

The intensity of pain during the hot water challenge (HPI) was similar for the 2WT and 4WT subgroups, and was not influenced by age or sex. Additional file 6 shows the results of the correlation analyses of HPI against the hemodynamic and HRV parameters. Only the increase of heart rate during VM and LF/HF under DB were found to significantly influence HPI, with an inverse correlation. Significance was still present after Bonferroni’s correction for two tests. Those effects are displayed in Fig. 4. Neither age nor sex interacted with those two significant correlations.

Fig. 4
figure 4

Scatterplots showing the relationships between hemodynamic or heart rate variability (HRV) parameters that have been found significant (p < 0.05) in the correlation analyses without Bonferroni’s correction. The dependent variables are either heat pain intensity (experimental pain) or the composite pain/analgesia score (CPAS) (see Methods’ chapter for details). The parameters tested are either the increase in heart rate as observed under Valsalva maneuver (VM) or the LF/HF ratio (HRV parameter, see Table 1 for interpretation) under deep breathing challenge (DB). The values shown for each relationship scatterplot are the Spearman’s correlation coefficient (ρ), the p value for the comparison of ρ to the null value (without Bonferroni’s correction), and the coefficient of determination R2

The CPAS was similar for the 2WT and 4WT subgroups, as expected due to the within-subgroup ranking. Neither age nor sex significantly influenced the CPAS, although it was slightly higher in females (2.4 [1.4–2.9]) than in males (1.8 [1.0–2.6]). There was a significant positive correlation between HPI and CPAS (ρ = 0.195; p = 0.035; R2 = 3.8%); correlation was weaker with the total pain score and the analgesia pain score taken separately (ρ = 0.185 and 0.165 respectively). Neither age nor sex interacted with this relationship (tested by ANCOVA). Additional file 7 shows the results of the correlation analyses of CPAS against the hemodynamic and HRV parameters. Only the increase of heart rate during VM influenced the CPAS, with an inverse correlation; however, significance was no longer present after Bonferroni’s correction for two tests. This effect is displayed in Fig. 4. Neither age nor sex interacted with this relationship (tested by ANCOVA).

As only one patient in the cohort reported chronic hypertension, we did not take this factor into consideration in our analyses. We did study the influence of parental hypertension score, which was unrelated to either the subgroup, the sex or the age, and was also unrelated to the HPI and to CPAS.

Discussion

In this exploratory study conducted in a healthy young population exposed to a standardized surgical painful condition, we were able to identify a relationship between the sympathetic component of the ANS and pain, while no effect of the parasympathetic component could be shown. As the observed effects (inverse correlations, i.e. the sympathetic response may protect form pain) were significant only for the increase in heart rate in the Valsalva maneuver, we must consider that this relationship was mild, compared to the multiple putative determinants of pain. A limitation of the study lies on the existence of two surgery subgroups, including difference in the anesthetic protocol and the NSAID molecule; nevertheless, the other typological differences between those subgroups were mild, and furthermore this bias was corrected as the primary outcome was ranked within each subgroup.

The reliability of the explorations of ANS we performed, as well as the meaning of our observations, must be discussed under the light of the current knowledge. While it is currently admitted that HRV may reflect the ANS [14, 25], its ability to explore its sympathetic component is a matter of debate [16]. Also, the frequency-domain analyses are considered to be a better fit for short recordings [14]. HRV measurements in resting conditions reflect the natural balance between the two components of ANS, and normative data obtained in very similar conditions in young adults fell into similar ranges to those of our observations [47]. However, and probably because the ANS balance was in the physiological range in this healthy population, those resting conditions were not favorable to observe an effect of any of the two components of ANS. The VM and DB challenges have been designed to elicit the ANS activity and therefore to potentiate the putative effects of its components [13]. Firstly, as we observed an increase of all the HRV parameters under the VM and DB challenges, we can consider that ANS was generally functional in our cohort. However, while both VM and DB evoked significant increases of all the HRV parameters– of which most are supposed to be parasympathetic-driven– they acted differently on the main hemodynamic parameters (Fig. 2). According to physiological knowledge, VM induces a multiphase response with an initial direct effect of intrathoracic overpressure eliciting a sympathetic response (with tachycardia), then delayed mixed effects due to both a rebound improvement of circulation and a parasympathetic response [25]. As our recording averaged the events occurring during all the phases and beyond, we could not discriminate the successive responses, but we observed that VM stimulated both components of the ANS, with a lesser sympathetic response (as expressed by tachycardia). Contrary to the orthostatic test (which we did not perform), VM does not increase diastolic arterial pressure, probably because the intrathoracic overpressure stresses the heart cavities and because the sympathetic response is more hormonal that neural [25]. How DB influences ANS is more controversial, although it might stimulate preferentially the cardiac afferent innervation through extracardiac baroreceptors. While a close relationship between the parasympathetic control and DB-induced HRV has been demonstrated [48, 49], a predominantly sympathetic response was suggested by a study on young healthy volunteers, in which DB increased sensibly more LF/HF than HF [50]. Such between-studies discrepancies can be explained by differences in the technique used and in the intrathoracic pressure regimens. As, in our DB test, HF and LF/HF both increased in a similar range (by 150 and 143%, respectively), an activation of both components is possible, but a recent review pointed out that LF/HF may not reflect properly the sympathetic component of ANS, especially if respiratory frequency is out of its physiological range [16]. Finally, we observed that HPI was able to predict CPAS, confirming that preoperative pain sensitivity partly explains postoperative pain [21]. Nevertheless, HPI explained only 4.3% of the whole variance of CPAS, highlighting the multiple other factors influencing postoperative pain (e.g. psychology, surgery), and this may explain why VM-induced tachycardia predicted HPI but not CPAS.

Our results suggest that pain sensitivity is related to the subject’s capacity for a sympathetic response. Nevertheless, a recent systematic review investigated the ability of ECG-derived assessments of ANS function to predict pain in various clinical conditions (i.e. experimental, postoperative, and persistent pain), and highlighted a general trend in showing antinociceptive effects of the parasympathetic, but not of the sympathetic tone [12]. However, while such trend was also noted in predicting surgical acute pain– which we consider as particularly relevant to study pain mechanisms– the results of those studies were not so strongly conclusive. One study conducted in 30 patients undergoing open carpal surgery (moderate postoperative pain) reported that the parasympathetic component (HF at rest) predicted a lesser pain [10]; LF/HF was not tested. A British team worked on a HRV-based polynomial model to predict postoperative pain in small samples of patients (respectively 17 and 29) undergoing varicose vein surgery, and concluded that most of the autonomic contribution in the model was parasympathetic; nevertheless, the role of LF (non-specific) seemed to overcome this of HF [51, 52]. In 30 patients undergoing surgery for epilepsy, postoperative pain was also negatively correlated with HF, but positively correlated with LF/HF [53]; of note, HRV was recorded in patients entering the operating room, i.e. in stressful conditions. A study conducted in neonates reported a predictive value of a marker of parasympathetic activity in postoperative pain, but here again, the sympathetic component was not explored [54]. Finally, the preoperative (not HRV-based) cardiovagal baroreflex sensitivity has been studied in 55 patients before cardiothoracic surgery, and– contrary to the previous studies– this marker of parasympathetic tone was found as favoring postoperative pain [55]; heterogeneity in surgery, anesthesia and postoperative analgesia could have influenced those results. With a sensibly bigger and homogeneous sample of patients, the results of the current study contradict those previous reports of a parasympathetic influence on pain, as none of the related parameters– even under potentiation of VM and DB and without Bonferroni’s correction– was show to be influential.

Neither age nor sex were found to influence the pain outcomes in our cohort (despite a small trend for sex), and none of these biometric factors interacted with the effects of HRV on pain, while both had some influence on the hemodynamic and HRV parameters. While pain sensitivity is known to decrease with age [56, 57], and women are more sensitive to pain than men [58, 59], none of these effects were observed in our cohort. This could be explained by the particularities of our sample or of our pain models, or by insufficient statistical power (age variation was limited in our sample). On the other hand, we observed that age correlated with arterial pressure (a phenomenon currently explained by vascular stiffness) [60], and tended to impair many HRV parameters, also consistently with published data [47, 49, 61]. Of note, such age-induced impairment was not observed for LF/HF, also consistent with previous reports [47], and suggesting that ageing has a greater effect on the parasympathetic component of HRV. This could also explain why the age effect on HRV was not converted into an effect on pain in our study. In terms of sex, we observed that our male patients had a higher arterial pressure, as well as higher values for several HRV parameters, but this was only during VM and was not observed for LF/HF. This was also in accordance with reports in young resting subjects [47, 60]. Whereas this sex effect we observed on arterial pressure cannot currently be explained [60], this effect on HRV is explained by a higher vagal tone in young males, while the sympathetic response seems unaffected by sex [62, 63]. This could support the hypothesis that sex effects on pain– if any– are not transmitted by the ANS.

Conclusions

Within a standardized model of postoperative pain, we identified a relationship between the sympathetic component of ANS, and pain perception. Although the effects we observed were quite mild, this contradicts previous reports of a preferential predictive role of parasympathetic activity. We also observed a gap between experimental pain– as measured in healthy preoperative conditions– and postoperative pain, probably because the surgery itself adds many other physiological and psychological factors to the strictly neurosensory ones. Furthermore, the expression of pain itself as a relevant marker of the nociceptive processes is also matter of debate, as it is highly dependent from individual prior pain/life experiences. Anyway, a simple preoperative assessment of VM-induced tachycardia could be interesting for prediction of postoperative pain, but probably with a small added value besides the already known psychosocial predictors [64,65,66]. Finally, attempts at modifying the ANS balance by medical or psychological interventions to prevent pain deserve further attention in future research.

Data availability

No datasets were generated or analysed during the current study.

Abbreviations

ANOVA-RM:

Analysis-of-variance for repeated measures

ANS:

Autonomic nervous system

ASA:

American Association of Anesthesiologists

CPAS:

Composite pain/analgesia score

DB:

Deep breathing (challenge)

HPI:

Heat pain intensity

HRV:

Heart rate variability

NN:

Normal-to-normal (intervals)

NNT:

Number-needed-to-treat

NSAID:

Nonsteroidal anti-inflammatory drug

PCA:

Principal component analysis

VAS:

Visual analogue scale

VM:

Valsalva maneuver

2WT/4WT:

“Two” and “four wisdom teeth” subgroup, respectively

References

  1. Suarez-Roca H, Klinger RY, Podgoreanu MV, Ji RR, Sigurdsson MI, Waldron N, et al. Contribution of baroreceptor function to pain perception and perioperative outcomes. Anesthesiology. 2019;130:634–50.

    Article  PubMed  Google Scholar 

  2. France CR. Decreased pain perception and risk for hypertension: considering a common physiological mechanism. Psychophysiology. 1999;36:683–92.

    Article  PubMed  CAS  Google Scholar 

  3. Hagen K, Zwart JA, Holmen J, Svebak S, Bovim G, Stovner LJ. Does hypertension protect against chronic musculoskeletal complaints? The Nord-Trøndelag Health Study. Arch Intern Med. 2005;165:916–22.

    Article  PubMed  Google Scholar 

  4. Deschaumes C, Devoize L, Sudrat Y, Baudet-Pommel M, Dualé C, Dallel R. The relationship between resting arterial blood pressure and oral postsurgical pain. Clin Oral Investig. 2015;19:1299–305.

    Article  PubMed  Google Scholar 

  5. D’Antono B, Ditto B, Sita A, Miller SB. Cardiopulmonary baroreflex stimulation and blood pressure-related hypoalgesia. Biol Psychol. 2000;53:217–31.

    Article  PubMed  Google Scholar 

  6. Ring C, van Veldhuijzen JJ, McIntyre D, Kavussanu M. Hypervolemic hyperalgesia in healthy young adults. J Behav Med. 2007;30:449–54.

    Article  PubMed  Google Scholar 

  7. Ditto B, Lewkowski MD, Rainville P, Duncan GH. Effects of cardiopulmonary baroreceptor activation on pain may be moderated by risk for hypertension. Biol Psychol. 2009;82:195–7.

    Article  PubMed  Google Scholar 

  8. Devoize L, Dualé C, Dubray C, Dallel R. Impact of sympathetic activation on pain threshold in human subjects. Physiol Behav. 2017;177:1–3.

    Article  PubMed  CAS  Google Scholar 

  9. Blanc JJ, Alboni P, Benditt DG. Vasovagal syncope in humans and protective reactions in animals. Europace. 2015;17:345–9.

    Article  PubMed  Google Scholar 

  10. Nielsen R, Nikolajsen L, Kroner K, Molgaard H, Vase L, Jensen TS, et al. Pre-operative baroreflex sensitivity and efferent cardiac parasympathetic activity are correlated with post-operative pain. Acta Anaesthesiol Scand. 2015;59:475–85.

    Article  PubMed  CAS  Google Scholar 

  11. Forte G, Troisi G, Pazzaglia M, Pascalis V, Casagrande M. Heart rate variability and pain: a systematic review. Brain Sci. 2022;12.

  12. Wegeberg AM, Sejersgaard-Jacobsen TH, Brock C, Drewes AM. Prediction of pain using electrocardiographic-derived autonomic measures: a systematic review. Eur J Pain. 2024;28:199–213.

    Article  PubMed  Google Scholar 

  13. Ewing DJ, Martyn CN, Young RJ, Clarke BF. The value of cardiovascular autonomic function tests: 10 years experience in diabetes. Diabetes Care. 1985;8:491–8.

    Article  PubMed  CAS  Google Scholar 

  14. ESC-NASPE. Heart rate variability. Standards of measurement, physiological interpretation, and clinical use. Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. Eur Heart J. 1996;17:354–81.

    Article  Google Scholar 

  15. Mazzeo AT, La Monaca E, Di Leo R, Vita G, Santamaria LB. Heart rate variability: a diagnostic and prognostic tool in anesthesia and intensive care. Acta Anaesthesiol Scand. 2011;55:797–811.

    Article  PubMed  Google Scholar 

  16. Hayano J, Yuda E. Pitfalls of assessment of autonomic function by heart rate variability. J Physiol Anthropol. 2019;38:3.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Moore RA, Straube S, Paine J, Derry S, McQuay HJ. Minimum efficacy criteria for comparisons between treatments using individual patient meta-analysis of acute pain trials: examples of etoricoxib, Paracetamol, ibuprofen, and ibuprofen/paracetamol combinations after third molar extraction. Pain. 2011;152:982–9.

    Article  PubMed  CAS  Google Scholar 

  18. Moore PA, Hersh EV. Combining ibuprofen and acetaminophen for acute pain management after third-molar extractions: translating clinical research to dental practice. J Am Dent Assoc. 2013;144:898–908.

    Article  PubMed  CAS  Google Scholar 

  19. Tousignant-Laflamme Y, Rainville P, Marchand S. Establishing a link between heart rate and pain in healthy subjects: a gender effect. J Pain. 2005;6:341–7.

    Article  PubMed  Google Scholar 

  20. Streff A, Kuehl LK, Michaux G, Anton F. Differential physiological effects during tonic painful hand immersion tests using hot and ice water. Eur J Pain. 2010;14:266–72.

    Article  PubMed  Google Scholar 

  21. Abrishami A, Chan J, Chung F, Wong J. Preoperative pain sensitivity and its correlation with postoperative pain and analgesic consumption: a qualitative systematic review. Anesthesiology. 2011;114:445–57.

    Article  PubMed  Google Scholar 

  22. Dualé C, Leray V, Giron F, Boulliau S, Macian N, Ruscheweyh R, et al. The added value of sensitivity to nonnoxious stimuli to predict an individual’s sensitivity to pain. Pain Physician. 2021;24:E783–94.

    Article  PubMed  Google Scholar 

  23. Luedi MM, Schober P, Hammoud B, Andereggen L, Hoenemann C, Doll D. Preoperative pressure pain threshold is associated with postoperative pain in short-stay anorectal surgery: a prospective observational study. Anesth Analg. 2021;132:656–62.

    Article  PubMed  Google Scholar 

  24. Keet SW, Bulte CS, Boer C, Bouwman RA. Reproducibility of non-standardised autonomic function testing in the pre-operative assessment screening clinic. Anaesthesia. 2011;66:10–4.

    Article  PubMed  CAS  Google Scholar 

  25. Ramos G, Carrasco S, Medina V. Time-frequency analysis of the heart rate variability during the Valsalva manoeuvre. J Med Eng Technol. 2000;24:73–82.

    Article  PubMed  CAS  Google Scholar 

  26. Dai F, Silverman DG, Chelly JE, Li J, Belfer I, Qin L. Integration of pain score and morphine consumption in analgesic clinical studies. J Pain. 2013;14:767–77.

    Article  PubMed  CAS  Google Scholar 

  27. Carvalho B, Hilton G, Wen L, Weiniger CF. Prospective longitudinal cohort questionnaire assessment of labouring women’s preference both pre- and post-delivery for either reduced pain intensity for a longer duration or greater pain intensity for a shorter duration. Br J Anaesth. 2014;113:468–73.

    Article  PubMed  CAS  Google Scholar 

  28. Moore RA, McQuay HJ. Single-patient data meta-analysis of 3453 postoperative patients: oral tramadol versus placebo, codeine and combination analgesics. Pain. 1997;69:287–94.

    Article  PubMed  CAS  Google Scholar 

  29. Edwards JE, McQuay HJ, Moore RA. Combination analgesic efficacy: individual patient data meta-analysis of single-dose oral tramadol plus acetaminophen in acute postoperative pain. J Pain Symptom Manage. 2002;23:121–30.

    Article  PubMed  CAS  Google Scholar 

  30. Van Aken H, Thys L, Veekman L, Buerkle H. Assessing analgesia in single and repeated administrations of propacetamol for postoperative pain: comparison with morphine after dental surgery. Anesth Analg. 2004;98:159–65.

    Article  PubMed  Google Scholar 

  31. Moller PL, Juhl GI, Payen-Champenois C, Skoglund LA. Intravenous acetaminophen (Paracetamol): comparable analgesic efficacy, but better local safety than its prodrug, propacetamol, for postoperative pain after third molar surgery. Anesth Analg. 2005;101:90–6.

    Article  PubMed  CAS  Google Scholar 

  32. Moller PL, Sindet-Pedersen S, Petersen CT, Juhl GI, Dillenschneider A, Skoglund LA. Onset of acetaminophen analgesia: comparison of oral and intravenous routes after third molar surgery. Br J Anaesth. 2005;94:642–8.

    Article  PubMed  CAS  Google Scholar 

  33. Juhl GI, Norholt SE, Tonnesen E, Hiesse-Provost O, Jensen TS. Analgesic efficacy and safety of intravenous Paracetamol (acetaminophen) administered as a 2 g starting dose following third molar surgery. Eur J Pain. 2006;10:371–7.

    Article  PubMed  CAS  Google Scholar 

  34. Derry C, Derry S, Moore RA, McQuay HJ. Single dose oral ibuprofen for acute postoperative pain in adults. Cochrane Database Syst Rev. 2009;CD001548.

  35. Akural EI, Jarvimaki V, Lansineva A, Niinimaa A, Alahuhta S. Effects of combination treatment with ketoprofen 100 mg + acetaminophen 1000 mg on postoperative dental pain: a single-dose, 10-hour, randomized, double-blind, active- and placebo-controlled clinical trial. Clin Ther. 2009;31:560–8.

    Article  PubMed  CAS  Google Scholar 

  36. Barden J, Derry S, McQuay HJ, Moore RA. Single dose oral ketoprofen and dexketoprofen for acute postoperative pain in adults. Cochrane Database Syst Rev. 2009;CD007355.

  37. McNicol ED, Tzortzopoulou A, Cepeda MS, Francia MB, Farhat T, Schumann R. Single-dose intravenous Paracetamol or propacetamol for prevention or treatment of postoperative pain: a systematic review and meta-analysis. Br J Anaesth. 2011;106:764–75.

    Article  PubMed  CAS  Google Scholar 

  38. Derry CJ, Derry S, Moore RA. Single dose oral ibuprofen plus Paracetamol (acetaminophen) for acute postoperative pain. Cochrane Database Syst Rev. 2013;CD010210.

  39. Bailey E, Worthington HV, van Wijk A, Yates JM, Coulthard P, Afzal Z. Ibuprofen and/or Paracetamol (acetaminophen) for pain relief after surgical removal of lower wisdom teeth. Cochrane Database Syst Rev. 2013;12:CD004624.

    Google Scholar 

  40. Moore RA, Derry S, Wiffen PJ, Straube S, Aldington DJ. Overview review: comparative efficacy of oral Ibuprofen and Paracetamol (acetaminophen) across acute and chronic pain conditions. Eur J Pain. 2015;19:1213–23.

    Article  PubMed  CAS  Google Scholar 

  41. Moore RA, Derry S, Aldington D, Wiffen PJ. Single dose oral analgesics for acute postoperative pain in adults - an overview of Cochrane reviews. Cochrane Database Syst Rev. 2015;9:CD008659.

    Google Scholar 

  42. Bender R, Lange S. Adjusting for multiple testing–when and how? J Clin Epidemiol. 2001;54:343–9.

    Article  PubMed  CAS  Google Scholar 

  43. Harris RJ. A primer of multivariate statistics. 2nd ed. New York: Academic; 1985.

    Google Scholar 

  44. Green SB. How many subjects does it take to do a regression analysis. Multivar Behav Res. 1991;26:499–510.

    Article  CAS  Google Scholar 

  45. Bujang MA, Baharum N. Sample size guideline for correlation analysis. World J Soc Sci Res. 2016;3:37–46.

    Article  Google Scholar 

  46. Kohn MA, Senyak J. UFCS-CTSI. Sample size calculators for designing clinical research– Correlation sample size. Updated 2025. Available at: https://sample-size.net/correlation-sample-size

  47. Kim GM, Woo JM. Determinants for heart rate variability in a normal Korean population. J Korean Med Sci. 2011;26:1293–8.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Fouad FM, Tarazi RC, Ferrario CM, Fighaly S, Alicandri C. Assessment of parasympathetic control of heart rate by a noninvasive method. Am J Physiol. 1984;246:H838–42.

    PubMed  CAS  Google Scholar 

  49. Shields RWJr. Heart rate variability with deep breathing as a clinical test of cardiovagal function. Cleve Clin J Med. 2009;76(Suppl 2):S37–40.

    Article  Google Scholar 

  50. Palak K, Furgala A, Ciesielczyk K, Szygula Z, Thor PJ. The changes of heart rate variability in response to deep breathing in professional swimmers. Folia Med Cracov. 2013;53:43–52.

    PubMed  Google Scholar 

  51. Adjei T, Von Rosenberg W, Goverdovsky V, Powezka K, Jaffer U, Mandic DP. Pain prediction from ECG in vascular surgery. IEEE J Transl Eng Health Med. 2017;5:2800310.

    Article  PubMed  Google Scholar 

  52. Powezka K, Adjei T, Von Rosenberg W, Normahani P, Goverdovsky V, Standfield NJ, et al. A pilot study of preoperative heart rate variability predicting pain during local anesthetic varicose vein surgery. J Vasc Surg Venous Lymphat Disord. 2019;7:382–6.

    Article  PubMed  Google Scholar 

  53. Caton L, Bolzon M, Boschiero D, Thayer JF, Gidron Y. Pre-surgical heart-rate variability strongly predicts less post-operative pain in patients with epilepsy. J Psychosom Res. 2021;145:110421.

    Article  PubMed  Google Scholar 

  54. Verweij LM, Kivits JTS, Weber F. The performance of the heart rate variability-derived newborn infant parasympathetic evaluation index as a measure of early postoperative pain and discomfort in infants-A prospective observational study. Paediatr Anaesth. 2021;31:787–93.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Suarez-Roca H, Mamoun N, Watkins LL, Bortsov AV, Mathew JP. Higher cardiovagal baroreflex sensitivity predicts increased pain outcomes after cardiothoracic surgery. J Pain. 2024;25:187–201.

    Article  PubMed  Google Scholar 

  56. Neziri AY, Scaramozzino P, Andersen OK, Dickenson AH, Arendt-Nielsen L, Curatolo M. Reference values of mechanical and thermal pain tests in a pain-free population. Eur J Pain. 2011;15:376–83.

    Article  PubMed  Google Scholar 

  57. Lue YJ, Wang HH, Cheng KI, Chen CH, Lu YM. Thermal pain tolerance and pain rating in normal subjects: gender and age effects. Eur J Pain. 2018;22:1035–42.

    Article  PubMed  Google Scholar 

  58. Racine M, Tousignant-Laflamme Y, Kloda LA, Dion D, Dupuis G, Choiniere M. A systematic literature review of 10 years of research on sex/gender and experimental pain perception - part 1: are there really differences between women and men? Pain. 2012;153:602–18.

    Article  PubMed  Google Scholar 

  59. Racine M, Tousignant-Laflamme Y, Kloda LA, Dion D, Dupuis G, Choiniere M. A systematic literature review of 10 years of research on sex/gender and pain perception - part 2: do biopsychosocial factors alter pain sensitivity differently in women and men? Pain. 2012;153:619–35.

    Article  PubMed  Google Scholar 

  60. Wills AK, Lawlor DA, Matthews FE, Sayer AA, Bakra E, Ben-Shlomo Y, et al. Life course trajectories of systolic blood pressure using longitudinal data from eight UK cohorts. PLoS Med. 2011;8:e1000440.

    Article  PubMed  PubMed Central  Google Scholar 

  61. O’Brien IA, O’Hare P, Corrall RJ. Heart rate variability in healthy subjects: effect of age and the derivation of normal ranges for tests of autonomic function. Br Heart J. 1986;55:348–54.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Koenig J, Rash JA, Campbell TS, Thayer JF, Kaess M. A meta-analysis on sex differences in resting-state vagal activity in children and adolescents. Front Physiol. 2017;8:582.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Geovanini GR, Vasques ER, de Oliveira Alvim R, Mill JG, Andreão RV, Vasques B, et al. Age and sex differences in heart rate variability and vagal specific patterns– baependi heart study. Glob Heart. 2020;15:71.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Sipilä RM, Haasio L, Meretoja TJ, Ripatti S, Estlander AM, Kalso EA. Does expecting more pain make it more intense? Factors associated with the first week pain trajectories after breast cancer surgery. Pain. 2017;158:922–30.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Pinto PR, Vieira A, Pereira D, Almeida A. Predictors of acute postsurgical pain after inguinal hernioplasty. J Pain. 2017;18:947–55.

    Article  PubMed  Google Scholar 

  66. Schreiber KL, Zinboonyahgoon N, Xu X, Spivey T, King T, Dominici L, et al. Preoperative psychosocial and psychophysical phenotypes as predictors of acute pain outcomes after breast surgery. J Pain. 2019;20:540–56.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would thank the staff of the dental surgery unit and the clinical research and innovation office of the university hospital of Clermont-Ferrand (France).

Funding

The devices to explore HRV and pain sensitivity were provided by the INSERM U1107 Neuro-Dol research unit (Clermont-Ferrand, France).

Author information

Authors and Affiliations

Authors

Contributions

CDe, LD and YS recruited the patients, conducted the explorations, performed the surgeries and collected the data. CDe, LD, BP, RD and CDu designed the study. CDe coordinated the research. BP and CDu conducted the analyses. CDe and CDu drafted the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Christian Dualé.

Ethics declarations

Ethics approval and consent to participate

This study has been approved by the appropriate research ethics committee (Comité de Protection des Personnes Sud-Est VI, ref. n° AU841). All participants gave their written informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Deschaumes, C., Devoize, L., Sudrat, Y. et al. The role of the sympathetic component of the autonomic nervous system on pain before and after third molar extraction– an observational cohort study. BMC Anesthesiol 25, 117 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12871-025-02949-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12871-025-02949-8

Keywords