Background:

Novel pulse photoplethysmographic–derived indices have been proposed as tools to measure autonomic nervous system (ANS) modulation in anesthetized and awake patients, but nowadays their experimental validation is lacking. The authors aimed to investigate the ability of pulse photoplethysmographic amplitude (PPGA), ANS state (ANSS), and ANSS index (ANSSi) to measure changes of ANS modulation in response to sympathetic stimulation.

Methods:

Ten awake healthy volunteers underwent two passive head-up tilts at 45° and 90°. The heart rate variability (HRV) and systolic arterial pressure variability were analyzed in the frequency domain as a measure of ANS modulation directed to the heart and the vessels. HRV, baroreflex sensitivity, and pulse photoplethysmographic indices were measured at baseline and after tilt maneuvers. The agreement between HRV-derived indices and pulse photoplethysmographic indices was assessed using Bland–Altman plots.

Results:

PPGA, ANSS, and ANSSi changed significantly during the study protocol. Head-up tilt decreased PPGA and ANSS and increased ANNSi. There was a good agreement between ANSSi and baroreflex sensitivity explored in the high-frequency band (bias, 0.23; 95% CI, −22.7 to 23.2 normalized units) and between ANSSi and the sympathovagal modulation directed to the heart (bias, 0.96; 95% CI, −8.7 to 10.8 normalized units).

Conclusions:

In controlled experimental conditions, novel pulse plethysmographic indices seem to estimate the changes of the sympathetic outflow directed to the vessels and the sympathovagal balance modulating heart rate. These indices might be useful in the future to monitor the fluctuation of sympathetic activity in anesthetized patients.

In controlled experimental conditions, novel pulse plethysmograhic indices estimated changes of the sympathetic outflow directed to vessels and the sympathovagal balance modulating heart rate.

What We Already Know about This Topic
  • Surgical stimulation induces a stress response characterized by sympathetic activation increasing heart rate, blood pressure, catecholamine, and corticosteroid release

  • This study determined the ability of pulse photoplethysmographic amplitude, autonomic nervous system state, and autonomic nervous system state index to measure changes of autonomic nervous system modulation in response to a gravitational challenge in healthy adults

What This Article Tells Us That Is New
  • In controlled experimental conditions, novel pulse plethysmograhic indices estimated changes of the sympathetic outflow directed to vessels and the sympathovagal balance modulating heart rate

THE autonomic nervous system (ANS) modulates the heart beat-to-beat interval and peripheral muscle vascular tone in response to a multitude of stimuli. It is known that surgical stimulation induces a stress response1,2  characterized by sympathetic activation increasing heart rate, blood pressure, catecholamine, and corticosteroid release, and it has been supposed influencing postoperative outcome.3,4 

Measuring changes of ANS activity during general anesthesia is an exciting challenge for the anesthesiologists. In a daily clinical practice, titrating the depth of hypnosis and analgesia to maintain a stable level of sympathetic activity might exert interesting effects on patients’ outcome. Unfortunately, the direct measure of the sympathetic and vagal activity seems not feasible in a clinical setting. Traditionally, indirect measurement of the ANS modulation on cardiovascular system was based on the analysis of the heart rate variability (HRV) and systolic arterial pressure (SAP) variability in well-established experimental conditions.5,6  This type of analysis is still confined in a research setting because of complexity of interplay of cardiovascular signals and it does not seem to provide a clear and simple parameter to inform the anesthesiologist in real time about ANS activity. Recently, some authors stated that two indices derived from pulse photoplethysmographic amplitude (PPGA) analysis of the photoplethysmograph, namely ANS state (ANSS) and ANSS index (ANSSi), would be useful tools to measure the autonomic modulation.7,8  In fact PPGA is because of pulsatile changes in tissue volume, mainly arterial blood, and it decreases during sympathetic mediated vasoconstriction. Unfortunately, until today, the experimental validation of ANSS and ANSSi was never been provided. The aim of this study is to assess the accuracy of ANSS and ANSSi to measure changes in ANS modulation directed to the heart and the vessels in response to a sympathetic stimulus, elicited through a gravitational load, in healthy adult humans.

After the approval of local institutional ethical committee (Luigi Sacco Hospital, Milan, Italy), we studied the effects of two orthostatic challenges on HRV, SAP variability, PPGA, ANSS, and ANSSi in 10 awake adult volunteers. The study adhered to the principles of the Declaration of Helsinki for medical research involving human subjects, and written informed consent was obtained from all volunteers.

All subjects were healthy, did not take any drug, beverages containing caffeine or alcohol during the previous 24 h, and were asked to fast 3 h before the study protocol. Electrocardiogram and photoplethysmographic waves were recorded through a S/5 Avance monitor (GE, Finland) to a laptop computer provided with S/5 Collect software (GE), and they were sampled at 300 Hz. The photoplethysmograph probe was a standard pulse oximetry probe in endowment to the commercial S/5 Avance monitor and was positioned on the third finger of the right hand. Noninvasive continuous arterial pressures were recorded by Nexfin (Edwards Lifesciences, USA) through a cuff positioned on the third finger’s middle phalanx of the right hand, sampled at 400 Hz, and stored to a laptop computer running LabChart Pro 7 (ADInstruments, New Zealand). All signals were synchronized at the beginning of the recording and analyzed offline. During the study sequence, the healthy volunteers were asked to lie calmly, breathing in rhythm with a metronome at 18 breaths per minutes (0.3 Hz), and they were not allowed to talk.

It is well known that changes in position of patients can lead to changes in the balance of the ANS. One of the most accepted maneuver to stimulate ANS-mediated cardiovascular control is head-up tilt (HUT).9–13  Tilt causes shift of blood toward lower body and reduction of the venous return. These effects induce a compensatory baroreflex-mediated increase of heart rate and peripheral vascular resistance aiming to maintain arterial pressure near to the level preceding the challenge.13  The experiment protocol consisted of a sequence of passive postural changes, each of them lasting 10 min. The sequence was (1) baseline recording in supine position (bas), (2) HUT at 45° (HUT45), (3) recovery after HUT45, (4) HUT at 90° (HUT90), and (5) recovery after HUT90. The postural changes were gained with a tilt table (TILT TEST®, Gardhen Bilance, Italy) present in our institution.

Measurements and Extraction of Beat-to-Beat Variability Series, ANSS, and ANSSi

After detecting the QRS complex on the electrocardiographic wave and locating the R-apex using parabolic interpolation, the temporal distance between two consecutive R parabolic apexes was computed and used as an approximation of the heart period. The maximum of arterial pressure inside heart period was taken as SAP. The occurrences of QRS and SAP peaks were manually checked to avoid erroneous detections or missed beats. If isolated ectopic beats affected the heart period and SAP values, these measures were linearly interpolated using the closest values unaffected by ectopic beats. Heart period and SAP measures were performed on a beat-to-beat basis. The sequences of 300 values after 2 min from changing position were selected inside each experimental step. The mean (μ) and the variance (σ2) of heart period and SAP are expressed in ms, mmHg, ms2, and mmHg2, respectively.

The maximum pulse plethysmographic amplitude within each heart period was detected on the photoplethysmographic wave (fig. 1). Pulse-to-pulse interval was detected as interval between consecutive pulse photoplethysmographic peaks. ANSS was calculated as a mean of ANSS of 300 consecutive pulse beats, corresponding to beat-to-beat interval analyzed from electrocardiogram, and ANSSi was calculated from the same series following the equations7,8 :

Fig. 1.

Sample of recorded waves from a healthy subject. BP = noninvasive continuous blood pressure wave; ECG = electrocardiographic wave; PPGA = pulse plethysmographic amplitude; PPI = pulse-to-pulse interval on photoplethysmographic wave; Pulse Pleth = photoplethysmographic wave; RR = heart period.

Fig. 1.

Sample of recorded waves from a healthy subject. BP = noninvasive continuous blood pressure wave; ECG = electrocardiographic wave; PPGA = pulse plethysmographic amplitude; PPI = pulse-to-pulse interval on photoplethysmographic wave; Pulse Pleth = photoplethysmographic wave; RR = heart period.

Close modal
(1)
(2)

where PPI is the pulse-to-pulse interval of the photoplethysmogram, PPGA is the pulse PPGA, and ANSSmax is the largest ANSS for the subject among 300 beats during each experimental condition.

Power Spectral Analysis

The power spectrum was estimated according to an univariate parametric approach fitting the series according to an autoregressive model.9  Autoregressive spectral density was factorized into components, each of them characterized by a central frequency. A spectral component was labeled as low frequency (LF) if its central frequency was between 0.04 and 0.15 Hz, whereas it was classified as high frequency (HF) if its central frequency was between 0.15 and 0.5 Hz.6  The LF and HF powers were defined as the sum of the powers of all LF and HF spectral components, respectively. The HF spectral density of heart beat intervals (HFRR), expressed in absolute units (ms2), was used as a marker of vagal modulation directed to the heart,5  whereas the LF spectral density of systolic arterial pressure oscillations (LFSAP), expressed in absolute units (mmHg2) was used as a marker of sympathetic modulation directed to the vessels.14  The ratio of the LFRR power to the HFRR (LFRR/HFRR) was considered an indicator of the balance between sympathetic and vagal modulation directed to the heart (sympathovagal balance).10  When LFRR/HFRR increases, it indicates the predominance of sympathetic over vagal activity. The power of heart period and SAP series in LF and HF bands was calculated as well and utilized for the estimation of the baroreflex sensitivity (BRS) (equations 3 and 4).

BRS Assessment

BRS was assessed according to a spectral analysis.15  This method was grounded on the evaluation of LFRR, LFSAP, HFRR, and HF spectral density of systolic arterial pressure oscillations (HFSAP). BRS was computed as follows:

(3)
(4)

and expressed in ms·mmHg−1. The prerequisites of high correlation (i.e., >0.5) and negative phase between heart period and SAP series, indicating that heart period changes lagged behind SAP variations, were tested according to the calculation of squared coherence (K2RR,SAP) and phase spectrum (Ph2RR-SAP).16  The cross-spectrum was estimated with a bivariate parametric approach fitting the heart period and SAP series according to a bivariate autoregressive model.17 

Statistics

Sample size was calculated with PS Power and Sample Size Calculator Software for Windows (Vanderbilt University, U.S.A.).18  In a previous study, ANSSi increased from 23.3 (SD, 13.7) to 78.1 (SD, 7.4) in response to painful stimulus.7  We needed to study seven subjects to detect a mean ANSSi difference of 30 (SD, 15) between baseline and HUT90 with a power of 0.8 and α error of 0.05. However, we planned to enroll 10 subjects in the eventuality of artifacts precluding HRV analysis. For continuous variables, we used mean (SD) or median (interquartile range) when not normally distributed. The normality distribution was checked using the Kolmogorov–Smirnov test. We performed the one-way ANOVA, preceded by Levene’s test, for repeated measures and post hoc Bonferroni test for multiple comparison. The correlation between HRV variables and pulse photopletysmographic–derived indices were checked with the Spearman correlation coefficient. With the hypothesis assessed by some authors7,8  that pulse photoplethysmographic–derived indices measure ANS activity, we compared HRV variables and pulse photoplethysmographic variables through Bland–Altman plot.19  Because the studied variables have different measuring units, we normalized each data series to a common scale to test their agreement on measuring ANS activity. Normalization consisted of 0 to 1 normalization following the formula and multiplying the obtained value per 100, where X was a value of the series, and Xmin and Xmax, respectively, were the lowest and the highest value of the series. We proceeded to the normalization of our measures to limit the excursion of a set of values within a certain predefined range. However, normalization—as it is a mathematical procedure—cannot correct error in measuring; moreover, our measures were determined using interval level of measurement. Because we recorded data at five time points for all patients, we performed the Bland–Altman analysis adjusted for a random effect model to estimate within-subject variance, in which each subject has a different intercept and slope over the observation period.20  In this analysis, we entered the “subject” and the “time point” (without interaction term) into the model as random effects. As fixed effects, we entered “position” (supine or upright) and the “mean of values” of 10 subjects at each time point. Statistical analysis was carried out using GraphPad Prism 5 (GraphPad Software Inc., USA) and R 3.0 with lme4 package for the mixed model analysis. For all test, P value less than 0.05 was considered as significant.

All subjects completed the study protocol, and none of them experienced symptoms of hypotension or syncope. Results of univariate analysis of heart period and SAP series from 10 healthy volunteers (median age, 27.5 yr; range, 25 to 34 yr; M/F = 4/6) are summarized in table 1 and fig. 2. Mean heart rate increased during head-up position respect to the supine position, whereas SAP remained unchanged. The indices derived from the analysis of heart period and SAP in the frequency domain (HFRR, LFRR/HFRR, and LFSAP) showed a significant variability during the study protocol. LFSAP increased significantly from baseline to HUT45 and HUT90 and LFRR/HFRR increased significantly at HUT 90, whereas HFRR decreased significantly from baseline to HUT90.

Table 1.

Heart Rate and Systolic Arterial Pressure Variability Characteristics of the Studied Subjects

Heart Rate and Systolic Arterial Pressure Variability Characteristics of the Studied Subjects
Heart Rate and Systolic Arterial Pressure Variability Characteristics of the Studied Subjects
Fig. 2.

The median (10th to 90th percentiles) of baroreflex sensitivity in the high-frequency band (αHF, A), pulse photoplethysmographic amplitude (PPGA, B), autonomic nervous system state (ANSS, C), and autonomic nervous system state index (ANSSi, D) during the study protocol. Data during recovery period are not shown. * A significant difference between boxes (P < 0.05) assessed with ANOVA for repeated measurements at five time points and Bonferroni post hoc test. Bas = baseline; HUT45 = head-up tilt at 45°; HUT90 = head-up tilt at 90°.

Fig. 2.

The median (10th to 90th percentiles) of baroreflex sensitivity in the high-frequency band (αHF, A), pulse photoplethysmographic amplitude (PPGA, B), autonomic nervous system state (ANSS, C), and autonomic nervous system state index (ANSSi, D) during the study protocol. Data during recovery period are not shown. * A significant difference between boxes (P < 0.05) assessed with ANOVA for repeated measurements at five time points and Bonferroni post hoc test. Bas = baseline; HUT45 = head-up tilt at 45°; HUT90 = head-up tilt at 90°.

Close modal

All measurement of BRS satisfied the prerequisite of squared coherence (K2RR,SAP) and phase spectrum (Ph2RR-SAP).16  The baroreflex sensitivity in the HF band (αHF) decreased significantly from baseline (14.9 ms·mmHg−1 [11.7 to 20.4 ms·mmHg−1]) to HUT45 (7.6 ms·mmHg−1 [5.8 to 11.6 ms·mmHg−1], P < 0.05) and HUT90 (4.8 ms·mmHg−1 [2.4 to 6.3 ms·mmHg−1], P < 0.05).

All pulse photoplethysmographic–derived indices (PPGA, ANSS, and ANSSi) changed significantly during the study protocol (one-way ANOVA, P < 0.0001). PPGA decreased from baseline (4.4% [1.9 to 11%]) to HUT45 (2.9% [1.7 to 3.7%], P < 0.05) and HUT90 (1.8% [0.9 to 2.6%], P < 0.05) (fig. 2B). ANSS decreased from baseline (3.7% [2 to 7.4%], P < 0.05) to HUT45 (2.6% [1.5 to 3.2%], P < 0.05) and HUT90 (1.4% [0.6 to 2%], P < 0.05) (fig. 2C). ANSSi increased from baseline (64.9 [37 to 80.3]) to HUT45 (78.3 [73.7 to 83], P < 0.05) and HUT90 (87 [84.6 to 89.9], P < 0.05) (fig. 2D). There were not significant differences of PPGA, ANSS, and ANSSi between HUT45 and HUT90. Additional data are shown in the appendices 1 and 2.

There were significant correlations between crude values of αHF and pulse photoplethysmographic indices and between the sympathovagal balance (LFRR/HFRR) and ANSSi (table 2). After normalization of the series and the application of the linear mixed model (correcting for within-subjects variance), the resulting intercepts of ANS variables of all subjects are compared through Bland–Altman plots (fig. 3). There was good agreement between (1) αHF and ANSSi (bias, 0.23 normalized units [NU]; 95% CI, −22.7 to 23.2 NU) and (2) LF/HF and ANSSi (bias, 0.96 NU; 95% CI, −8.7 to 10.8 NU) in measuring ANS modulation during the study protocol.

Table 2.

Correlations between Heart Rate Variability Parameters and Pulse Photoplethysmographic Indices

Correlations between Heart Rate Variability Parameters and Pulse Photoplethysmographic Indices
Correlations between Heart Rate Variability Parameters and Pulse Photoplethysmographic Indices
Fig. 3.

Bland–Altman plot of agreement between couple of tested variables measuring autonomic nervous system modulation, corrected for within-subject variance with a linear mixed model. Thick lines represent the bias, and dotted lines represent ±1.96 SD. (A) αHF versus PPGA, bias = −0.04 (95% CI, −44.3 to 44.2); (B) αHF versus ANSS, bias = 0.09 (95% CI, −34.8 to 34.9); (C) αHF versus ANSSi, bias = 0.23 (95% CI, −22.7 to 23.2); (D) LF/HF versus ANSSi, bias = 0.96 (95% CI, −8.7 to 10.8). All data are given in normalized units. ANSS = autonomic nervous system state; ANSSi = autonomic nervous system state index; LF/HF = sympathovagal balance; PPGA = pulse photo-plethysmographic amplitude; αHF = baroreflex sensitivity in high-frequency band.

Fig. 3.

Bland–Altman plot of agreement between couple of tested variables measuring autonomic nervous system modulation, corrected for within-subject variance with a linear mixed model. Thick lines represent the bias, and dotted lines represent ±1.96 SD. (A) αHF versus PPGA, bias = −0.04 (95% CI, −44.3 to 44.2); (B) αHF versus ANSS, bias = 0.09 (95% CI, −34.8 to 34.9); (C) αHF versus ANSSi, bias = 0.23 (95% CI, −22.7 to 23.2); (D) LF/HF versus ANSSi, bias = 0.96 (95% CI, −8.7 to 10.8). All data are given in normalized units. ANSS = autonomic nervous system state; ANSSi = autonomic nervous system state index; LF/HF = sympathovagal balance; PPGA = pulse photo-plethysmographic amplitude; αHF = baroreflex sensitivity in high-frequency band.

Close modal

In 10 awake healthy subjects, we found that during gravitational sympathetic stimulation, as expected, SAP power in the LF band increased, baroreflex sensitivity assessed in the HF band decreased, and sympathovagal balance increased toward a prevalence of sympathetic modulation over vagal activity. Moreover, the HF power of HRV was reduced during tilt. These findings suggest a significant increase of sympathetic modulation directed to the vessels and a reduction of vagal modulation directed to the heart. In the orthostatic position, the blood shift to the lower body determines the decrease in central blood volume and pressure unloading cardiopulmonary and carotid baroreceptors, thus leading to tachycardia, splanchnic, and cutaneous vasoconstriction.21–23  We found that pulse photoplethysmograph–derived indices changed significantly during tilt, thus suggesting that they reflect mainly the sympathetic-mediated vasoconstriction. It is well demonstrated that different types of sympathetic stimuli (painful, orthostatic, and lower body negative pressure) increase the firing rates of sympathetic muscular fibers resulting in vessels constriction.11,24–26  These effects are the principal determinants of PPGA changes because of sympathetic stimulation. The pulse PPGA is because of pulsatile changes of the arteriolar blood volume into the tissue. The blood volume pulsations (ΔV) are related to both the systemic intravascular pulse pressure (ΔP) and the distensibility of the vascular wall (D) according to the relationship: .27  The distensibility is influenced by intravascular volume status and sympathetic activity directed to the vessels, whereas vagal fibers are absent in the peripheral vascular bed. As a consequence, indices derived from the pulsatile PPGA are affected by a wide variety of stimuli exciting the sympathetic branch of the ANS. It is well known that nociception induces changes of ANS modulation toward a sympathetic activation.24,28  Surgical Pleth Index, a PPGA-derived index, was proposed as a measure of nociception–antinociception balance at first29,30 ; afterward, it was demonstrated being affected by several confounding factors such as atropine administration, level of sedation, spinal anesthesia, intravascular volume status, and patient’s position.31–34  All these factors alter the ANS without affecting the nociception. We previously demonstrated that the changes of ANS modulation correlate with the changes of Surgical Pleth Index during general anesthesia.35  We believe that this point is the focus of the misleading about pulse photoplethysmographic–derived indices as a measure of nociception–antinociception: In the present study, we provide the demonstration that pulse photoplethysmographic indices seem to reflect both sympathetic-mediated changes in vascular tone and sympathovagal efference to the heart.

It is well known that passive HUT determines a pooling of blood toward the lower body accordingly to gravitational attraction and leads to a reduction of the venous return. These changes induce a reflected sympathetic activity and increase of muscle sympathetic nerve’s firing rate aiming to squeeze the blood toward the upper body to maintain arterial pressure near to the level preceding the challenge.9–13  Many human studies reported that increasing the angle tilt results in a greater muscle sympathetic nerve activity.12,36–38  Because the sympathetic activation does not discriminate between upper and lower body, its effect is “massive” causing tachycardia, generalized vasoconstriction, and reducing blood flow even in the forearm.26  We found that a greater intensity of sympathetic stimulation (at higher degree of tilt angle) did not changed significantly PPGA, ANSSi, and ANSS respect to HUT45. The reason of significance dearth between the two tilt angles might be due to a high activation of balancing sympathetic reflexes in response to the orthostatic challenge in healthy young subject since from lower degree angles. Furthermore, small differences between two different degrees of tilt could be detected with a large sample.

HRV analysis is a well-documented method to measure the sympathetic and vagal modulation directed to the heart and the vessels.9–11  Conversely, some authors, in absence of adequate experimental validation, postulated that novel pulse photoplethysmograph–derived indices represent an adequate tool to measure ANS activity.7,8  With this study, we provide an experimental evidence that pulse photoplethysmographic indices, especially ANSSi, increased in response to increases of sympathetic activity. To quantify the agreement between HRV and pulse photoplethysmograph–derived indices on measuring ANS modulation, we considered the ANS activity as an independent variable. Because the ANS activity is a concept that encompasses a dynamic sympathetic and vagal outflow (that is difficult to quantify as absolute size) and the analyzed variables have quite different measuring units, we normalize each series to better compare the variables with each others. We found that pulse photoplethysmograph–derived indices have a good agreement with the more validated HRV-derived indices as a measure of ANS modulation.

Methods for the HRV analysis in the frequency domain may be generally classified as nonparametric (i.e., fast Fourier transform) and parametric (i.e., autoregressive).6  Although in most instances, both methods provide comparable results, in this study, we choose to analyze the heart rate and the systolic pressure variability with an autoregressive model. The advantages of this model are an easy identification of the central frequency of each component and an accurate estimation of frequency spectrum even on a small number of samples in which the signal is steady.6  Moreover, autoregressive analysis provides a more reliable estimation of the HF component than the fast Fourier transform does.39,40  Because the vagal efference detected in the HF band is synchronous with respiratory activity, the healthy subjects in our studies breathed at a fixed respiratory rate of 0.3 Hz. This allows us to better measure the HF spectral component, because its central frequency was about in the middle of the HF band (0.15 to 0.5 Hz).

Our study has few limitations. We analyzed, from a physiologic point of view, the effects of sympathetic activity on pulse photoplethysmographic signals without administration of hypnotics or opiates. It has been demonstrated that general anesthesia reduces significantly ANS modulation on cardiovascular system.41  Surgical stimuli can increase sympathetic outflow in these circumstances, but they introduce new variables in this research field: the intensity of surgical stimulation and the sickness of the studied patient. We strongly believe that a physiologic correlation between ANS activity and photoplethysmographic waves should be first assessed in a stable experimental condition in which all confounders should be reduced as much as possible. Furthermore, because the linear regressions showed highly skewed data, we applied a nonparametric correlation that ignores the impact of the nonindependence of observations that is created with the repeated measurements for each individual. This approach allows us to consider the crude association between the HRV parameters and the photoplethysmographic indices, but confounds within-subject variance and between-subject variance. These indices are only intended to serve as a gross measure of the association between the two measurements. The use of a mixed model (i.e., with subjects having random effects) would be a large improvement in these circumstances. However, under most situations, this model would also be subject to extreme outliers.

In conclusion, in these healthy adults, we found that PPGA, ANSS, and ANSSi changed significantly in response to a sympathetic gravitational challenge. We believe that the results of this study provide a more extensive understanding of the physiologic basis of ANSSi: It seems to be a good and easy to keep surrogate to measure changes of the sympathovagal balance and sympathetic outflow directed to the vessels. HRV analysis provides a more comprehensive measure of ANS cardiovascular modulation taking in account each autonomic efference: (1) the vagal activity directed to the heart (HFRR), (2) the sympathetic outflow directed to the vessels (LFSAP), and (3) their interaction with the baroreflex control (αHF and baroreflex sensitivity in the LF band [αLF]). However, these measures are difficult to obtain in a clinical setting lacking automated, real-time calculation of these variables. Furthermore HRV analysis in the frequency domain requires the stationarity of the signals6  often difficult to keep in a daily clinical practice. In this scenario, ANSSi might be a valuable tool displaying on anesthesia monitor the changes of global cardiovascular sympathetic modulation.

We thank Alberto Porta, M.Eng. Ph.D. (Dipartimento di Scienze Biomediche per la Salute, University of Milan, Milano, Italy), for providing us with the software performing autoregressive heart rate variability, arterial pressure variability, and baroreflex sensitivity analysis.

Support was provided solely from institutional and/or departmental sources.

The authors declare no competing interests.

1.
Roizen
MF
,
Horrigan
RW
,
Frazer
BM
:
Anesthetic doses blocking adrenergic (stress) and cardiovascular responses to incision—MAC BAR.
Anesthesiology
1981
;
54
:
390
8
2.
Guyton
A
,
Hall
J
:
The autonomic nervous system and the adrenal medulla
in
Textbook of Medical Physiology
. Edited by
Hall
J
.
Philadelphia
,
Saunders
,
2006
, pp
pp 748
60
3.
Parker
SD
,
Breslow
MJ
,
Frank
SM
,
Rosenfeld
BA
,
Norris
EJ
,
Christopherson
R
,
Rock
P
,
Gottlieb
SO
,
Raff
H
,
Perler
BA
:
Catecholamine and cortisol responses to lower extremity revascularization: Correlation with outcome variables. Perioperative Ischemia Randomized Anesthesia Trial Study Group.
Crit Care Med
1995
;
23
:
1954
61
4.
Monk
TG
,
Ding
Y
,
White
PF
:
Total intravenous anesthesia: Effects of opioid versus hypnotic supplementation on autonomic responses and recovery.
Anesth Analg
1992
;
75
:
798
804
5.
Akselrod
S
,
Gordon
D
,
Ubel
FA
,
Shannon
DC
,
Berger
AC
,
Cohen
RJ
:
Power spectrum analysis of heart rate fluctuation: A quantitative probe of beat-to-beat cardiovascular control.
Science
1981
;
213
:
220
2
6.
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.
Circulation
1996
;
93
:
1043
65
7.
Hamunen
K
,
Kontinen
V
,
Hakala
E
,
Talke
P
,
Paloheimo
M
,
Kalso
E
:
Effect of pain on autonomic nervous system indices derived from photoplethysmography in healthy volunteers.
Br J Anaesth
2012
;
108
:
838
44
8.
Paloheimo
MP
,
Sahanne
S
,
Uutela
KH
:
Autonomic nervous system state: The effect of general anaesthesia and bilateral tonsillectomy after unilateral infiltration of lidocaine.
Br J Anaesth
2010
;
104
:
587
95
9.
Pagani
M
,
Lombardi
F
,
Guzzetti
S
,
Rimoldi
O
,
Furlan
R
,
Pizzinelli
P
,
Sandrone
G
,
Malfatto
G
,
Dell’Orto
S
,
Piccaluga
E
,
Turiel
M
,
Baselli
G
,
Cerutti
S
,
Malliani
A
:
Power spectral analysis of heart rate and arterial pressure variabilities as a marker of sympatho-vagal interaction in man and conscious dog.
Circ Res
1986
;
59
:
178
93
10.
Montano
N
,
Ruscone
TG
,
Porta
A
,
Lombardi
F
,
Pagani
M
,
Malliani
A
:
Power spectrum analysis of heart rate variability to assess the changes in sympathovagal balance during graded orthostatic tilt.
Circulation
1994
;
90
:
1826
31
11.
Furlan
R
,
Porta
A
,
Costa
F
,
Tank
J
,
Baker
L
,
Schiavi
R
,
Robertson
D
,
Malliani
A
,
Mosqueda-Garcia
R
:
Oscillatory patterns in sympathetic neural discharge and cardiovascular variables during orthostatic stimulus.
Circulation
2000
;
101
:
886
92
12.
Cooke
WH
,
Hoag
JB
,
Crossman
AA
,
Kuusela
TA
,
Tahvanainen
KU
,
Eckberg
DL
:
Human responses to upright tilt: A window on central autonomic integration.
J Physiol
1999
;
517
(
pt 2
):
617
28
13.
Robertson
D
,
Diedrich
A
,
Chapleau
MW
:
Editorial on arterial baroreflex issue.
Auton Neurosci
2012
;
172
:
1
3
14.
Pagani
M
,
Montano
N
,
Porta
A
,
Malliani
A
,
Abboud
FM
,
Birkett
C
,
Somers
VK
:
Relationship between spectral components of cardiovascular variabilities and direct measures of muscle sympathetic nerve activity in humans.
Circulation
1997
;
95
:
1441
8
15.
Pagani
M
,
Somers
V
,
Furlan
R
,
Dell’Orto
S
,
Conway
J
,
Baselli
G
,
Cerutti
S
,
Sleight
P
,
Malliani
A
:
Changes in autonomic regulation induced by physical training in mild hypertension.
Hypertension
1988
;
12
:
600
10
16.
de Boer
RW
,
Karemaker
JM
,
Strackee
J
:
Relationships between short-term blood-pressure fluctuations and heart-rate variability in resting subjects. II: A simple model.
Med Biol Eng Comput
1985
;
23
:
359
64
17.
Porta
A
,
Baselli
G
,
Rimoldi
O
,
Malliani
A
,
Pagani
M
:
Assessing baroreflex gain from spontaneous variability in conscious dogs: Role of causality and respiration.
Am J Physiol Heart Circ Physiol
2000
;
279
:
H2558
67
19.
Bland
JM
,
Altman
DG
:
Measuring agreement in method comparison studies.
Stat Methods Med Res
1999
;
8
:
135
60
20.
Myles
PS
,
Cui
J
:
Using the Bland-Altman method to measure agreement with repeated measures.
Br J Anaesth
2007
;
99
:
309
11
21.
Rowell
L
:
Reflex control during orthostasis
in
Human Cardiovascular Control
. Edited by
Rowell
LB
.
New York
,
Oxford University Press
,
1993
, pp
pp 37
80
22.
Jacobsen
TN
,
Morgan
BJ
,
Scherrer
U
,
Vissing
SF
,
Lange
RA
,
Johnson
N
,
Ring
WS
,
Rahko
PS
,
Hanson
P
,
Victor
RG
:
Relative contributions of cardiopulmonary and sinoaortic baroreflexes in causing sympathetic activation in the human skeletal muscle circulation during orthostatic stress.
Circ Res
1993
;
73
:
367
78
23.
Taylor
JA
,
Halliwill
JR
,
Brown
TE
,
Hayano
J
,
Eckberg
DL
:
‘Non-hypotensive’ hypovolaemia reduces ascending aortic dimensions in humans.
J Physiol
1995
;
483
(
pt 1
):
289
98
24.
Burton
AR
,
Birznieks
I
,
Bolton
PS
,
Henderson
LA
,
Macefield
VG
:
Effects of deep and superficial experimentally induced acute pain on muscle sympathetic nerve activity in human subjects.
J Physiol
2009
;
587
(
pt 1
):
183
93
25.
Cooke
WH
,
Rickards
CA
,
Ryan
KL
,
Kuusela
TA
,
Convertino
VA
:
Muscle sympathetic nerve activity during intense lower body negative pressure to presyncope in humans.
J Physiol
2009
;
587
(
pt 20
):
4987
99
26.
Fu
Q
,
Iwase
S
,
Niimi
Y
,
Kamiya
A
,
Kawanokuchi
J
,
Cui
J
,
Mano
T
,
Suzumura
A
:
Effects of lower body positive pressure on muscle sympathetic nerve activity response to head-up tilt.
Am J Physiol Regul Integr Comp Physiol
2001
;
281
:
R778
85
27.
Burton
AC
:
Relation of structure to function of the tissues of the wall of blood vessels.
Physiol Rev
1954
;
34
:
619
42
28.
Jeanne
M
,
Logier
R
,
De Jonckheere
J
,
Tavernier
B
:
Heart rate variability during total intravenous anesthesia: Effects of nociception and analgesia.
Auton Neurosci
2009
;
147
:
91
6
29.
Huiku
M
,
Uutela
K
,
van Gils
M
,
Korhonen
I
,
Kymäläinen
M
,
Meriläinen
P
,
Paloheimo
M
,
Rantanen
M
,
Takala
P
,
Viertiö-Oja
H
,
Yli-Hankala
A
:
Assessment of surgical stress during general anaesthesia.
Br J Anaesth
2007
;
98
:
447
55
30.
Rantanen
M
,
Yli-Hankala
A
,
van Gils
M
,
Yppärilä-Wolters
H
,
Takala
P
,
Huiku
M
,
Kymäläinen
M
,
Seitsonen
E
,
Korhonen
I
:
Novel multiparameter approach for measurement of nociception at skin incision during general anaesthesia.
Br J Anaesth
2006
;
96
:
367
76
31.
Bonhomme
V
,
Uutela
K
,
Hans
G
,
Maquoi
I
,
Born
JD
,
Brichant
JF
,
Lamy
M
,
Hans
P
:
Comparison of the surgical Pleth Index™ with haemodynamic variables to assess nociception-anti-nociception balance during general anaesthesia.
Br J Anaesth
2011
;
106
:
101
11
32.
Ilies
C
,
Gruenewald
M
,
Ludwigs
J
,
Thee
C
,
Höcker
J
,
Hanss
R
,
Steinfath
M
,
Bein
B
:
Evaluation of the surgical stress index during spinal and general anaesthesia.
Br J Anaesth
2010
;
105
:
533
7
33.
Ilies
C
,
Ludwigs
J
,
Gruenewald
M
,
Thee
C
,
Hanf
J
,
Hanss
R
,
Steinfath
M
,
Bein
B
:
The effect of posture and anaesthetic technique on the surgical pleth index.
Anaesthesia
2012
;
67
:
508
13
34.
Höcker
J
,
Broch
O
,
Gräsner
JT
,
Gruenewald
M
,
Ilies
C
,
Steinfath
M
,
Bein
B
:
Surgical stress index in response to pacemaker stimulation or atropine.
Br J Anaesth
2010
;
105
:
150
4
35.
Colombo
R
,
Raimondi
F
,
Corona
A
,
Rivetti
I
,
Pagani
F
,
Porta
VD
,
Guzzetti
S
:
Comparison of the Surgical Pleth Index with autonomic nervous system modulation on cardiac activity during general anaesthesia: A randomised cross-over study.
Eur J Anaesthesiol
2014
;
31
:
76
84
36.
Fu
Q
,
Okazaki
K
,
Shibata
S
,
Shook
RP
,
VanGunday
TB
,
Galbreath
MM
,
Reelick
MF
,
Levine
BD
:
Menstrual cycle effects on sympathetic neural responses to upright tilt.
J Physiol
2009
;
587
(
pt 9
):
2019
31
37.
Iwase
S
,
Mano
T
,
Watanabe
T
,
Saito
M
,
Kobayashi
F
:
Age-related changes of sympathetic outflow to muscles in humans.
J Gerontol
1991
;
46
:
M1
5
38.
Kamiya
A
,
Iwase
S
,
Sugiyama
Y
,
Mano
T
,
Sudoh
M
:
Vasomotor sympathetic nerve activity in men during bed rest and on orthostasis after bed rest.
Aviat Space Environ Med
2000
;
71
:
142
9
39.
Badilini
F
,
Maison-Blanche
P
,
Coumel
P
:
Heart rate variability in passive tilt test: Comparative evaluation of autoregressive and FFT spectral analyses.
Pacing Clin Electrophysiol
1998
;
21
:
1122
32
40.
Pichon
A
,
Roulaud
M
,
Antoine-Jonville
S
,
de Bisschop
C
,
Denjean
A
:
Spectral analysis of heart rate variability: Interchangeability between autoregressive analysis and fast Fourier transform.
J Electrocardiol
2006
;
39
:
31
7
41.
Ogawa
Y
,
Iwasaki
K
,
Shibata
S
,
Kato
J
,
Ogawa
S
,
Oi
Y
:
Different effects on circulatory control during volatile induction and maintenance of anesthesia and total intravenous anesthesia: Autonomic nervous activity and arterial cardiac baroreflex function evaluated by blood pressure and heart rate variability analysis.
J Clin Anesth
2006
;
18
:
87
95
Appendix 1.

Heart Rate and Systolic Arterial Pressure Variability Characteristics of the Studied Subjects

Heart Rate and Systolic Arterial Pressure Variability Characteristics of the Studied Subjects
Heart Rate and Systolic Arterial Pressure Variability Characteristics of the Studied Subjects

Appendix 2. Hemodynamic and Pulse Photoplethysmographic Variables during the Study Protocol

ANSS = autonomic nervous system state; ANNSi = autonomic nervous system state index; Bas = baseline; HFRR = high-frequency spectral density of heart period; HUT45 = head-up tilt at 45°; HUT90 = head-up tilt at 90°; LFRR/HFRR = low frequency to high frequency ratio of heart period’s spectral density; LFSAP = low-frequency spectral density of systolic arterial pressure; PPGA = pulse photoplethysmographic amplitude; R1 = recovery after HUT45; R2 = recovery after HUT90; RR = heart period; SAP = systolic arterial pressure; αHF = baroreflex sensitivity in high-frequency band. * P < 0.05 between columns (post hoc Bonferroni test).