Postoperative Pain Assessment Indices Based on Photoplethysmography Waveform Analysis

The purpose of this study was to derive parameters that might reflect postoperative pain from photoplethysmography (PPG) and verify the derived parameters in postoperative pain assessment. We obtained preoperative and postoperative PPG and 100-mm visual analog scale (VAS) from 65 surgical patients and extracted a total of 51 PPG morphology-based parameters and their normalized parameters from these PPGs obtained. Pain discrimination performances of these derived parameters were assessed by statistical analyses, including Wilcoxon signed rank test with Bonferroni correction, classification accuracy based on logistic regression, and 4-fold cross validation. After comparing these parameters derived from PPG in pre- and post-operative conditions, statistically significant difference was found in 36 of the 51 parameters. Using logistic classification, dynamic between-pulse parameters such as normalized systolic amplitude variation and normalized diastolic amplitude variation showed better pain classification performance than the static within-pulse parameters. VAS score was 0 in every pre-operation condition, but >60 VAS was observed in the post-operative condition. Systolic peak amplitude variation normalized by PPG AC amplitude showed the best performance in classifying post-operative pain, with accuracy, sensitivity, specificity, and positive predictivity values of 79.5, 74.0, 86.0, and 84.5%, respectively. These results are superior to those of the surgical pleth index (SPI, GE Healthcare, Chicago, IL, United States) at 65.9, 65.9, 66.5, and 66.5%, respectively.

Introduction

Up to 50% of surgical patients are known to experience moderate or even severe pain immediately after surgery (Boselli et al., 2013; Gerbershagen et al., 2013). While opioids are widely used as effective analgesics, opioid sensitivity is well known to vary widely among individuals. Pain sensitivity and analgesic susceptibility are major factors limiting the objectivity of post-operative pain management. Opioids can result in underdose or overdose of medication. Overdosing of analgesics is a life-threatening condition with disadvantageous effects such as prolonged sedation or prolonged respiratory insufficiency (Macintyre et al., 2011; Apfel et al., 2012). In addition, opioid analgesic overdose encompasses a range of clinical findings such as respiratory depression, miosis, stupor, hepatic injury from acetaminophen or hypoxemia, myoglobinuric renal failure, rhabdomyolysis, absent or hypoactive bowel sounds, compartment syndrome, and hypothermia (Boyer, 2012).

Therefore, the challenge for an effective analgesic strategy is to estimate the dose that can exactly counterbalance the amount of ongoing nociception with respect to individual susceptibility to analgesics and individual sensitivity to pain (Ip et al., 2009). Maintaining an appropriate level of analgesics is very important to improve the patient’s prognosis (Musizza and Ribaric, 2010). However, it is difficult to judge the adequate amount of opioid because the degree of pain experienced during surgery varies depending on the patient as well as the environmental and emotional factors of the same patient. This is because indices such as anxiety, sensitivity, and distress tolerance closely related to pain experience can vary greatly among individuals.

To quantify surgical pain, previous studies have focused on autonomic nervous system responses such as changes in heart rate (HR) and blood pressure, amount of tears, changes in pupil size, respiratory rate and volume, and patient movement (Russell, 1989; Moerman et al., 1993). Among these responses, the surgical pleth index (SPI, GE Healthcare, Chicago, IL, United States), an index based on photoplethysmography (PPG), is a promising surgical pain quantification index. SPI is a non-invasive, dimensionless score (0–100), which allows an estimate of intraoperative nociception. SPI is a model of PPG amplitude (PPGA) and heart beat interval (HBI). It is modeled as SPI = 100 × (0.67 × PPGA_norm+0.33 × HBI_norm) (Huiku et al., 2007). PPGA is known to be strongly correlated with sympathetic response by surgical stimulation while HBI is known to correlate with opioid effect (Huiku et al., 2007). The value of SPI ranges from 0 to 100, with a value of 0 indicating weaker surgical stimulus, while a value closer to 100 indicated stronger surgical stimulus. In addition, SPI 50 can be used to differentiate a patient’s pain intensity in general anesthesia (Chen et al., 2010). It has been reported that SPI can vary with pain-stimuli or anesthetic concentration in propofol-remifentanil anesthesia (Huiku et al., 2007). SPI has demonstrated excellent performance compared to various other pain estimation methods based on HR, blood pressure, or response entropy (Wennervirta et al., 2008). Several studies have reported that SPI is a potential indicator for pain assessment in unconscious patients (Ahonen et al., 2007; Huiku et al., 2007; Kallio et al., 2008). However, according to Ledowski et al. (2009), SPI has limited accuracy in assessing conscious subjects. In fact, the manufacturer recommends the use of SPI only during anesthesia; this is presumably because SPI is normalized with data from anesthesia patients only. At present, SPI are frequently used in the clinical setting. However, a complete assessment technique for pain is not available, especially for post-operative pain.

In this study, we proposed new indicators in post-operative pain assessment by deriving characteristic features not presented in previous studies through PPG waveform analysis. We then verified the pain correlation of the derived candidate parameters. In particular, in the previous studies, pain caused by specific stimuli was evaluated, such as incision during surgery (Ahonen et al., 2007; Huiku et al., 2007; Kallio et al., 2008; Boselli et al., 2013), but this may include changes due to various factors such as the effect of an electrosurgical unit other than the pain itself. Moreover, in a clinical setting, it is important to distinguish pain from non-pain because the specified dose is repeatedly administered depending on the presence or absence of pain without adjusting the dose according to the pain intensity. Therefore, in order to observe the pure changes due to pain, in this study, two conditions were defined rather than the various pain intensities: preoperative condition as a state without pain and post-operative condition as a state of pain.

Materials and Methods

Ethics and Dataset

This study was approved by the Asan Medical Center Institutional Review Board (approval number: 2016-0477) and registered on the international clinical trials registry platform (¹KCT0002080). A total of 81 subjects were enrolled. All experiments were approved by the Human Research Ethics Committee of Asan Medical Center and all methods were performed in accordance with the relevant guidelines and regulations. All participants gave written informed consent for their participation in the study. Finally, data from 65 subjects (28 males, 37 females, mean age of 50.5 ± 10.5 years) were used for analysis, excluding data from 16 subjects: data from 2 subjects had recording errors, 1 subject had severe arrhythmia; and data from 13 subjects were recorded without SPI. Patient’s information and type of surgery are described in Table 1.

TABLE 1

TABLE 1. Patient characteristics and type of surgery.

Experimental Protocol

Photoplethysmography data were recorded at pre-operation and post-operation. All surgeries included skin incision, and were performed under general anesthesia to reduce surgical pain. Pre-operative data were recorded in the surgery waiting room immediately before surgery. First, the 100 mm visual analog scale (VAS) was recorded and the PPG and SPI were simultaneously recorded for 6 min. Post-operative data were recorded after transfer of the patient to the recovery room immediately after the operation. At this time, PPG and SPI were recorded for 6 min after first recording VAS. Opioid was not used for 6 min immediately after surgery, and it was administrated after signal measurement was completed.

The VAS consists of a straight line with endpoints defining the extreme limits of the signal. For pain intensity, the scale is most commonly anchored by “no pain” (score of 0) and “pain as bad as it could be” or “worst imaginable pain” [score of 100 (100 mm scale)] (Haefeli and Elfering, 2006). PPG was recorded on the index finger of the opposite arm with intravenous therapy at a sampling frequency of 300 Hz for 6 min in each period. SPI was simultaneously recorded every 10 s. A GE Datex-Ohmeda S/5 Anesthesia Monitor, disposable PPG sensor, and a portable computer running the S/5 Collect software (GE Healthcare, Chicago, IL, United States) were used for monitoring and recording, respectively. Every recording is performed at room temperature.

Preprocessing

The obtained PPG was filtered using a finite impulse response (FIR) bandpass filter at 0.5–10 Hz. Systolic and diastolic peak points were detected using the adaptive threshold method (Shin et al., 2009). The detected peak position was verified by experienced researchers and manually corrected for false detected or missed peaks. In the case of post-operative measurements, motion artifact may also occur, which can distort the data, even rendering it impossible to analyze. Therefore, the section where motion noise occurred was manually excluded from the analysis. Matlab 2016b (Mathworks, Inc., MA, United States) was used for all procedures of peak detection and signal analysis.

PPG Parameters

Basic Parameters

Pain affects the autonomic nervous system, leading to the vasoconstriction of sympathetic nerves (Huiku et al., 2007). Most parameters of this study were derived based on this clinical fact. We derived 23 parameters from PPG morphology analysis. These parameters were based on the spatial characteristic, amplitude, temporal characteristic, and the length or interval of PPG. Each parameter for each condition was calculated by averaging the values obtained for each pulse for 6 min. Graphical representation of derived parameters is shown in Figure 1. Detailed descriptions are summarized in Table 2.

FIGURE 1

FIGURE 1. Graphical representation of extracted parameters from PPG morphology analysis.

TABLE 2

TABLE 2. Parameters extracted from PPG morphology analysis.

Normalized Parameters

Amplitudes of PPG have arbitrary units (a.u.) because they are affected by various factors such as skin color of the subject, human body composition of the measurement site, anatomical structure, nail condition, ambient light, and sensor wearing condition at the time of measurement (Allen, 2007). Therefore, amplitude-based characteristics should not be quantified, and normalization for arbitrary units was considered. The parameters of the baseline-systolic amplitude of PPG (ACA_bl) and the adjacent diastolic-systolic amplitude (ACA_dia) had arbitrary units that required normalization. Area-related parameters related to ACA_bl and ACA_dia also had arbitrary units. Calculating the ratio between parameters including arbitrary units could give a dimensionless parameter instead of arbitrary units. Therefore, individual or measured conditional characteristics related to amplitude were removed. Time unit parameters are affected by HR, which has an inter-individual difference. Therefore, the inter-individual difference due to HR needs to be corrected. By calculating the ratio between time unit parameters, it is possible to eliminate inter-individual differences. Therefore, time unit parameters were normalized to adjacent diastolic peak interval or systolic peak-diastolic peak interval in the present study. Consequently, 28 normalized parameters including 27 non-dimensional parameters were derived by normalization. These normalized parameters are listed in Table 3.

TABLE 3

TABLE 3. List of normalized parameters.

Logistic Classification

Logistic classification has been used in binary classification problems. Logistic function is a monotonic, continuous function between 0 and 1. Mathematical expression of logistic function is represented in Eq. 1. Sigmoid function is a representative example of logistic function (L = 1, k = 1, x₀ = 0). In logistic classification, output is a probability. In our study, preoperative data were set as class 0 while post-operative data were set as class 1, and we used a 0.5 as a classification criterion of pain or non-pain assessment. This means that we classified the result is a pain if f(x) > 0.5, and non-pain if f(x) < 0.5.

$$\begin{array}{c}f(x)=\frac{1}{1+e^{-k(x-x_0)}}\left(1\right)\end{array}$$

Validation

Validation was obtained using two methods: statistical significance test and accuracy evaluation of pain classification. Statistical tests were performed to determine if the derived parameters had significant changes due to pain. In statistical analysis, proposed parameters extracted from data obtained in pre- and post-operation were analyzed for statistical significance. To determine inter-subject variance, we calculated the coefficient of variance (CV). For testing normality, the Kolmogorov–Smirnov test was used. Because all data do not have normality, we used the Wilcoxon signed rank test to verify the significant difference between parameters derived with PPG measured on pre- and post-operation. In multiple comparisons, the Bonferroni correction was used. P < 0.05 was considered statistically significant.

The accuracy of pain classification was evaluated by applying logistic classification. For the random partitioning of samples into a training set and a test set, we used 4-fold cross validation. Classification performance was evaluated based on accuracy (AC), sensitivity (SE), specificity (SP), and positive predictivity value (PPV). AC indicates the proportion of correctly classified numbers in the total case. SE is the probability that a classification result will indicate pain among those with pain. SE is the fraction of those without pain who have a negative classification result. PPV is the probability of actual pain among those classified as pain.

Results

Statistical Validation

In our experiment, we confirmed that 0 VAS was observed in the pre-operation condition, while >60 VAS was observed in the post-operative condition. Table 4 shows the mean, SD, CV, and P-value for each parameter before and after surgery. The results showed that SPI was significantly increased after the operation. Among the derived parameters, area-based indices such as A, A_sys, A_dia, TriA, TriA_sys, and TriA_dia were significantly decreased in the post-operative condition. Amplitude-based indices, ACA_dia and ACA_bl were also significantly decreased in the post-operation. For temporal indices, L_sys was significantly increased in post-operation while L_dia was not significantly changed. Slope-based parameters, RS, L_RS, RS_max, FS, and L_FS were decreased significantly in post-operation. Pulse widths, PW₃₀, PW₅₀, and PW₇₀ showed no significant change. However, PW₉₀ were increased significantly in post-operation. Normalized area-based indices, A_sys/A, A_dia/A, A_sys/A_dia, TriA_sys/TriA, TriA_diai/TriA, and TriA_dia/TriA, were found to have significant differences among the derived dimensionless parameters. The area of systolic phase was relatively increased, while the area of diastolic phase was decreased in post-operation. Similarly, A_sys/ACA_bl showed a significant increase in the post-operation, whereas A_dia/ACA_bl was not significantly changed. All of the RS/FS and RS/ACA_bl calculated as slope ratios were significantly decreased in the post-operative condition. However, the ratio of maximum rising slope to average rising slope, RS_max/RS, was not significantly changed. This might be due to an increase in the systolic interval caused by postoperative pain, consistent with the results of parameters derived from area ratio. In addition, parameters (A/ACA_bl)/L_sys, (A/ACA_bl)/L_dia, (A_sys/ACA_bl)/L_sys, and (A_sys/ACA_bl)/L_dia showed significant differences between pre- and post-operation. Results showed that a total of 15 parameters (L_dia, PPI_sys, PPI_dia, ACV_sys, ACV_dia, PW₃₀, PW₅₀, PW₇₀, A/ACA_bl, A_dia/ACA_bl, RS_max/RS, (A/ACA_bl)/L_sys, (A/ACA_bl)/PPI_dia, (A_dia/ACA_bl)/L_dia, and (A_dia/ACA_bl)/PPI_dia) were not significantly changed after the operation. For parameters representing dynamic inter-pulse variation, a significant increase was observed in the parameters of ACV_sys/ACA_dia, ACV_sys/ACA_bl, ACV_dia/ACA_dia, and ACV_dia/ACA_bl. This suggests that the fluctuation of amplitude is increased with pain.

TABLE 4

TABLE 4. Statistical changes of derived parameters according to pain.

Figure 2 shows a boxplot of the top 10 parameters with high accuracy. The lighter box in the figure shows preoperative results, while the darker box shows post-operative results. The boxplot in Figure 2 graphically depicts the groups of values of each parameter through their quartiles. Lines extending vertically from these boxes indicate variability outside the upper and lower quartiles. The cross mark (“×”) represents an outlier. Crosses on the right side refer to the outliers of more than 3/2 times the size of the upper quartile. Crosses on the left-side refer to the outliers of less than 3/2 times the size of the lower quartile. The spacing between different parts of the box indicates the degree of dispersion and skewness in the data.

FIGURE 2

FIGURE 2. Boxplot of the top 10 accuracy parameters and SPI before and after surgery. The lighter box shows pre-operative results while the darker box shows post-operative results. Cross mark (“×”) represents the outlier. Crosses on the right side refer to outliers of more than 1.5 times the interquartile range above the third quartile. Crosses on the left-side refer to outliers of less than 3/2 times the size of the lower quartile. The spacings between different parts of the box indicate the degree of dispersion (spread) and skewness in the data. Every parameter shows a significant difference by surgical stress (^∗∗∗P < 0.001, ^∗∗P < 0.01).

Classification Performance

Thirteen out of the 51 indices showed better performance than SPI in post-operative pain classification accuracy. Among these, the four indicators normalizing the difference of the adjacent pulse amplitudes, ACV_sys/ACA_dia, ACV_sys/ACA_bl, ACV_dia/ACA_dia, and ACV_dia/ACA_bl, showed remarkably high performance. The remaining top 10 accuracy of classification was higher in the order of pulse amplitude from the adjacent diastolic peak (ACA_dia), normalized diastolic waveform area (A_dia/A), pulse amplitude from baseline (ACA_bl), normalized systolic waveform area (A_sys/A), average slope of systolic interval (RS), and maximum value of systolic slope (RS_max). Table 5 shows the classification of the performance of parameters with statistically significant differences.

TABLE 5

TABLE 5. Classification performance based on logistic classification.

Discussion

Possibility of PPG as a Post-operative Pain Indicator

Photoplethysmography is an optical technique used to determine the blood volume changes by measuring the amount of light that is transmitted, reflected, or scattered after the light has been irradiated to the peripheral site of the body. PPG is mainly measured at the fingertips, tips of toes, or earlobes. PPG has a larger or smaller value depending on the amount of blood ejected from the heart. It generally defines the upward convex pole as the systolic peak and the downward convex pole as the diastolic peak (Allen, 2007). Previous studies using PPG waveform analysis for surgical pain assessment have found that the PPG amplitude and systolic peak interval are significantly reduced according to the skin incision during propofol-remifentanil anesthesia (Rantanen et al., 2006). In another study it was found that the maximum slope between the diastolic pole and the systolic pole is significantly reduced according to the skin incision during fentanyl, propofol, and sevoflurane anesthesia (Seitsonen et al., 2005). Huiku et al. (2007) extracted normalized PPGA (PPGA_norm), pulse transit time (PTT), systolic blood pressure (SBP), normalized heart beat interval (HBI_norm), RE, and sympatho-vagal ratio (SVR) from PPG or other signals. After analyzing pain relevance with these parameters, they found that PPG_norm, PTT, and SBP were highly correlated with pain stimuli (Huiku et al., 2007). Moreover, studies using state and RE, HR, and PPGA (Rantanen et al., 2006) or using HBI, SD1, SD2, SD1/SD2 of HRV Poincaré analysis results, RE, difference of RE and state entropy, PPGA, PPG notch amplitude, and PPG notch y position (Rantanen et al., 2006) have shown that SPI is superior to other indices.

Previous studies have revealed the relationship between PPG waveform characteristics and pain. However, further analytical methods are needed for pain assessment using PPG. The waveform of PPG can provide a range of hemodynamic information. For example, the perfusion index (the ratio of AC amplitude and DC level of PPG) is sensitive to proximal sympathectomy (Klodell et al., 2005), peripheral perfusion in critically ill patients (Lima et al., 2002), and neonatal left heart obstruction (Granelli and Östman-Smith, 2007). Stiffness index (defined as the amplitude ratio of the systolic peak and notch point) and reflection index (defined as systolic pole-to-notch position interval proposed by Millasseau et al. (2006) can be used to indirectly estimate the stiffness of the vessel. It has been reported that the ratio of the inflection points of the second derivative photoplethysmography (SDPTG) reflects the vascular stiffness due to aging (Takazawa et al., 1998). Thus, PPG morphology analysis reveals the possibility of extracting pain parameters with high significance, besides the parameters used in existing PPG-pain studies.

Effectiveness of PPG-Based Parameters in Post-operative Pain Assessment

Similar to previous studies, A, A_sys, and A_dia were significantly decreased after the operation (Seitsonen et al., 2005). Considering that the area-related parameter of the waveform reflects changes in blood volume at the measurement point, the reduction in area might be due to the contraction of blood vessels through sympathetic activation by pain. RS_max, RS, L_RS, FS, and L_FS associated with the up and down slopes of ACA_dia, ACA_bl, and PPG representing PPG amplitudes were also significantly reduced after the operation, which is consistent with the results of previous studies (Seitsonen et al., 2005; Rantanen et al., 2006; Kallio et al., 2008). These results can be explained by increased blood vessel resistance shown in pulse area change or decreased maximum blood volume due to vasoconstriction. From a similar viewpoint, significant increase in L_sys indicates that contracted vessels might have increased the vascular resistance due to vasoconstriction of the sympathetic nerves (Jaryal et al., 2009; Hickey et al., 2016). TriA, TriA_sys, and TriA_dia in the area of a triangle where the waveform of the high frequency component such as harmonics was removed showed similar patterns as those of A, A_sys, and A_dia, respectively. This might indicate the lack of significant correlation between pain and the high frequency component of PPG. HBI derived from systolic and diastolic extremes known to be pain-related was increased after the operation but in our study the increase was not statistically significant. While this result is not common, the previous research also showed that HBI has relatively low sensitivity and specificity compared with combined feature of PPG (Seitsonen et al., 2005; Rantanen et al., 2006).

In the case of PW, PW₃₀, PW₅₀, and PW₇₀ showed no significant change between pre- and post-operation. However, PW₉₀ was significantly increased in post-operation. Changes in PW₉₀ might be due to the narrowing of the waveform at the upper part of the amplitude when pain occurred. This can be interpreted in terms of increasing systolic interval length. Physiologically, the increase of systolic length increased the severity of vascular flow resistance (Angius et al., 2012; Wang et al., 2015).

Normalization

The normalization in this study is divided into two types: that for the spatial-domain of PPG such as amplitude and that for the temporal domain of the PPG waveform such as interval or length. First, on the spatial domain, since the amplitude has an arbitrary value as mentioned above, it has an inter-individual difference as well as an intra-individual difference for each measurement. Thus, we derived ratio-based indices that could compensate for the arbitrary change of the amplitude. Secondly, on the temporal dimension, since the basal heart has individual difference, it is impossible to determine whether or not the pain is based on a specific HR. Moreover, every temporal value such as the systolic and diastolic lengths also depends on the beat-to-beat interval. However, the percentage of systolic interval during HR may be uniform for every person compared with the HR.

As a result of normalization, we confirmed that CVs were decreased in most of the indices. In the spatial domain, CV of A_sys and A_dia was decreased by dividing with A from 62.1 to 15.8% and from 70.0 to 6.1%, respectively. In the temporal domain, the CV of L_dia was decreased from 20.1 to 5.2%. Interestingly, the CV of the L_sys was increased by normalization (7.5–16.8%). This result could suggest that systolic length is less effected by the pulse length, and this concurs with the result of previous research that investigates systolic length change according to the exercises (Franculescu et al., 1982). In ACV_sys and ACV_dia, CVs were decreased by normalization from 65.9 to 41.9% and from 67.5 to 41.8%, respectively. These results suggest that the normalization could reduce the intra-individual and inter-individual differences. However, for a more generalized result, normalization with a massive dataset is needed.

Dynamic Parameters

The results of pain classification using these proposed parameters showed that parameters derived from pulse-to-pulse analysis and normalized parameters were superior to those derived from single pulse analysis. For example, the classification accuracy of parameters such as ACV_sys/ACA_bl, ACV_sys/ACA_dia, ACV_dia/ACA_bl, and ACV_dia/ACA_dia reflecting the dynamic changes between pulses was about 80%, which was more than 10% higher than the accuracy when single parameters were used without normalization. These dynamic parameters can resemble blood pressure variability (BPV). BPV is known to reflect blood pressure variations in very short-term analysis (beat-to-beat), which influences the central and reflex autonomic modulation in short term analysis (<24 h) (Conway et al., 1984; Mancia et al., 1986; Parati et al., 1995) and the elastic properties of arteries (reduced arterial compliance) (Parati, 2005; Kotsis et al., 2011; Fukui et al., 2013). However, the relationship between BPV and pain has not yet been reported. Therefore, it is more reasonable to conclude that dynamic parameters of amplitude variation can distinguish the pain better than static parameters, since stroke volume variability (SVV), a dynamic parameter, can reflect volume responses more effectively than the static parameters HR and BP in hemodynamic analysis.

Collinearity Between Features

The aim of this study is to analyze the many features of PPG, and to identify the indicators that correlate with pain. Thus, a total of 51 indicators were derived, which include interdependent indicators. Indicator collinearity does not need to be treated as important when single indicators are used to classify pain. However, as in the case of the SPI, if collinearity is used to generate a combined indicator with some indicators, it may be necessary to evaluate the collinearity between the two indicators, since this may lead to errors in accuracy evaluation. The index collinearity is usually assessed using the condition number, variance inflation factor (VIF) or correlation coefficient. In this study, the collinearity between each index was calculated from the correlation coefficient, and is shown in Figure 3. In the figure, blue is the positive correlation coefficient and red is the negative correlation coefficient. The color concentration reflects the magnitude of the correlation coefficient. The results of this collinearity evaluation can be used in developing combined indicators using single indicators.

FIGURE 3

FIGURE 3. Collinearity test result between 51 features. Blue and red represents positive and negative correlation, respectively. The darker the color, the higher the correlation. 1: A, 2: A_sys, 3: A_dia, 4: TriA, 5: TriA_sys, 6: TriA_dia, 7: L_sys, 8: L_dia, 9: PPI_sys, 10: PPI_dia, 11: ACA_dia, 12: ACA_bl, 13: ACV_sys, 14: ACV_dia, 15: RS, 16: L_RS, 17: RS_max, 18: FS, 19: L_FS, 20: PW₃₀, 21: PW₅₀, 22: PW₇₀, 23: PW₉₀, 24: A_sys/A, 25: A_dia/A, 26: A_sys/A_dia, 27: TriA_sys/TriA, 28: TriA_dia/TriA, 29: TriA_sys/TriA_dia, 30: A/ACA_bl, 31: A_sys/ACA_bl, 32: A_dia/ACA_bl, 33: L_sys/PPI_dia, 34: L_dia/PPI_dia, 35: L_sys/L_dia, 36: ACV_sys/ACA_dia, 37: ACV_sys/ACA_bl, 38: ACV_dia/ACA_dia, 39: ACV_dia/ACA_bl, 40: RS_max/RS, 41: RS/FS, 42: RS/ACA_bl, 43: (A/ACA_bl)/L_sys, 44: (A/ACA_bl)/L_dia, 45: (A/ACA_bl)/PPI_dia, 46: (A_sys/ACA_bl)/L_sys, 47: (A_sys/ACA_bl)/L_dia, 48: (A_sys/ACA_bl)/PPI_dia, 49: (A_dia/ACA_bl)/L_sys, 50: (A_dia/ACA_bl)/L_dia, 51: (A_dia/ACA_bl)/PPI_dia

Limitations

The aim of our study is to evaluate individual predictors rather than regression by combining multiple predictors. If we develop a combination feature that combines individual features in the future, we should perform a multicollinearity test between the features found in this study. Also, in order for the study to be a useful index in clinical practice, it is necessary to derive a pain index discriminant formula that can express the pain score as a normalized value. Finally, this study does not include an assessment of the pain situation suppressed by analgesics. Therefore, it is necessary to verify whether the post-operative pain candidate index reflects real-time pain change following analgesic administration in order to be used clinically in post-operative pain evaluation. Moreover, a composite test for a specific age, disease, etc., is required to improve reliability.

Conclusion

Results of this study demonstrated that dynamic parameters obtained by inter-pulse analysis were more sensitive to pain than static parameters obtained from single pulse analysis. Considering that existing surgical pain assessment studies are performed mainly on the above mentioned static indicators, dynamic indicators found in this study can be used to develop more sophisticated pain assessment algorithms in the future. However, further studies are needed before these new pain parameters identified in this study can be applied to clinical practice, including evaluation in general anesthesia. First, this study was based on a statistical analysis of total signals for 6 min obtained before and after surgery without observing changes in pain in real-time. Therefore, it is necessary to determine whether real-time pain change estimation is possible. In addition, this study only dealt with single parameters obtained from PPG waveforms. Further studies with multivariate models using multiple indicators are needed for more sophisticated pain assessment. Finally, since the pain assessment technique must be provided as a formula that can quantify the degree of pain, further research is needed to determine the normalization technique to correct inter-individual differences.

Author Contributions

B-MC and G-JN collected the data. HS designed the data analysis scheme. YLY, HSS, and HS performed the data analysis. B-MC and HS interpreted the results. All authors contributed to the writing of the manuscript, provided critical revisions, and approved the final version.

Funding

This research was supported by a grant (NRF-2015R1C1A1A02036535) from the Basic Science Research Program through the National Research Foundation (NRF) of Korea, funded by the Ministry of Science, ICT and Future Planning, Republic of Korea and by the Bio & Medical Technology Development Program of the NRF funded by the Korean government, MSIP (NRF-2016M3A9F1941328).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Footnotes

1. ^ http://cris.nih.go.kr

References

Ahonen, J., Jokela, R., Uutela, K., and Huiku, M. (2007). Surgical stress index reflects surgical stress in gynaecological laparoscopic day-case surgery. Br. J. Anaesth. 98, 456–461. doi: 10.1093/bja/aem035

PubMed Abstract | CrossRef Full Text | Google Scholar

Allen, J. (2007). Photoplethysmography and its application in clinical physiological measurement. Physiol. Meas. 28, R1–R39. doi: 10.1088/0967-3334/28/3/R01

PubMed Abstract | CrossRef Full Text | Google Scholar

Angius, G., Barcellona, D., Cauli, E., Meloni, L., and Raffo, L. (2012). “Myocardial infarction and antiphospholipid syndrome: a first study on finger PPG waveforms effects,” in Proceedings of the Computing in Cardiology (CinC), 2012 (Kraków: IEEE), 517–520.

Google Scholar

Apfel, C. C., Heidrich, F. M., Jukar-Rao, S., Jalota, L., Hornuss, C., Whelan, R. P., et al. (2012). Evidence-based analysis of risk factors for postoperative nausea and vomiting. Br. J. Anaesth. 109, 742–753. doi: 10.1093/bja/aes276

PubMed Abstract | CrossRef Full Text | Google Scholar

Boselli, E., Bouvet, L., Bégou, G., Dabouz, R., Davidson, J., Deloste, J.-Y., et al. (2013). Prediction of immediate postoperative pain using the analgesia/nociception index: a prospective observational study. Br. J. Anaesth. 112, 715–721. doi: 10.1093/bja/aet407

PubMed Abstract | CrossRef Full Text | Google Scholar

Boyer, E. W. (2012). Management of opioid analgesic overdose. N. Engl. J. Med. 367, 146–155. doi: 10.1056/NEJMra1202561

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X., Thee, C., Gruenewald, M., Wnent, J., Illies, C., Hoecker, J., et al. (2010). Comparison of surgical stress index-guided analgesia with standard clinical practice during routine general anesthesia: a pilot study. Anesthesiology 112, 1175–1183. doi: 10.1097/ALN.0b013e3181d3d641

PubMed Abstract | CrossRef Full Text | Google Scholar

Conway, J., Boon, N., Davies, C., Jones, J. V., and Sleight, P. (1984). Neural and humoral mechanisms involved in blood pressure variability. J. Hypertens. 2, 203–208. doi: 10.1097/00004872-198404000-00013

CrossRef Full Text | Google Scholar

Franculescu, V., Sallam, E., and Otoiu, M. (1982). Systolic time intervals at rest and during exercise in sportsmen. Ann. Clin. Res. 14, 160–164.

Google Scholar

Fukui, M., Ushigome, E., Tanaka, M., Hamaguchi, M., Tanaka, T., Atsuta, H., et al. (2013). Home blood pressure variability on one occasion is a novel factor associated with arterial stiffness in patients with type 2 diabetes. Hypertens. Res. 36, 219–225. doi: 10.1038/hr.2012.177

PubMed Abstract | CrossRef Full Text | Google Scholar

Gerbershagen, H. J., Aduckathil, S., Van Wijck, A. J., Peelen, L. M., Kalkman, C. J., and Meissner, W. (2013). Pain intensity on the first day after surgery: a prospective cohort study comparing 179 surgical procedures. Anesthesiology 118, 934–944. doi: 10.1097/ALN.0b013e31828866b3

PubMed Abstract | CrossRef Full Text | Google Scholar

Granelli, A., and Östman-Smith, I. (2007). Noninvasive peripheral perfusion index as a possible tool for screening for critical left heart obstruction. Acta Paediatr. 96, 1455–1459. doi: 10.1111/j.1651-2227.2007.00439.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Haefeli, M., and Elfering, A. (2006). Pain assessment. Eur. Spine J. 15(Suppl. 1), S17–S24. doi: 10.1007/s00586-005-1044-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hickey, M., Phillips, J. P., and Kyriacou, P. A. (2016). Investigation of peripheral photoplethysmographic morphology changes induced during a hand-elevation study. J. Clin. Monit. Comput. 30, 727–736. doi: 10.1007/s10877-015-9761-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Huiku, M., Uutela, K., Van Gils, M., Korhonen, I., Kymäläinen, M., Meriläinen, P., et al. (2007). Assessment of surgical stress during general anaesthesia. Br. J. Anaesth. 98, 447–455. doi: 10.1093/bja/aem004

PubMed Abstract | CrossRef Full Text | Google Scholar

Ip, H. Y., Abrishami, A., Peng, P. W., Wong, J., and Chung, F. (2009). Predictors of postoperative pain and analgesic consumption: a qualitative systematic review. Anesthesiology 111, 657–677. doi: 10.1097/ALN.0b013e3181aae87a

PubMed Abstract | CrossRef Full Text | Google Scholar

Jaryal, A. K., Selvaraj, N., Santhosh, J., Anand, S., and Deepak, K. K. (2009). Monitoring of cardiovascular reactivity to cold stress using digital volume pulse characteristics in health and diabetes. J. Clin. Monit. Comput. 23, 123–130. doi: 10.1007/s10877-009-9174-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Kallio, H., Lindberg, L., Majander, A., Uutela, K., Niskanen, M., and Paloheimo, M. (2008). Measurement of surgical stress in anaesthetized children. Br. J. Anaesth. 101, 383–389. doi: 10.1093/bja/aen204

PubMed Abstract | CrossRef Full Text | Google Scholar

Klodell, C. T., Lobato, E. B., Willert, J. L., and Gravenstein, N. (2005). Oximetry-derived perfusion index for intraoperative identification of successful thoracic sympathectomy. Ann. Thorac. Surg. 80, 467–470. doi: 10.1016/j.athoracsur.2005.02.075

PubMed Abstract | CrossRef Full Text | Google Scholar

Kotsis, V., Stabouli, S., Karafillis, I., Papakatsika, S., Rizos, Z., Miyakis, S., et al. (2011). Arterial stiffness and 24h ambulatory blood pressure monitoring in young healthy volunteers: the early vascular ageing Aristotle University Thessaloniki Study (EVA–ARIS Study). Atherosclerosis 219, 194–199. doi: 10.1016/j.atherosclerosis.2011.07.111

PubMed Abstract | CrossRef Full Text | Google Scholar

Ledowski, T., Ang, B., Schmarbeck, T., and Rhodes, J. (2009). Monitoring of sympathetic tone to assess postoperative pain: skin conductance vs surgical stress index. Anaesthesia 64, 727–731. doi: 10.1111/j.1365-2044.2008.05834.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Lima, A. P., Beelen, P., and Bakker, J. (2002). Use of a peripheral perfusion index derived from the pulse oximetry signal as a noninvasive indicator of perfusion. Crit. Care Med. 30, 1210–1213. doi: 10.1097/00003246-200206000-00006

PubMed Abstract | CrossRef Full Text | Google Scholar

Macintyre, P. E., Loadsman, J. A., and Scott, D. A. (2011). Opioids, ventilation and acute pain management. Anaesth. Intensive Care 39, 545–558.

Google Scholar

Mancia, G., Parati, G., Pomidossi, G., Casadei, R., Di Rienzo, M., and Zanchetti, A. (1986). Arterial baroreflexes and blood pressure and heart rate variabilities in humans. Hypertension 8, 147–153. doi: 10.1161/01.HYP.8.2.147

CrossRef Full Text | Google Scholar

Millasseau, S. C., Ritter, J. M., Takazawa, K., and Chowienczyk, P. J. (2006). Contour analysis of the photoplethysmographic pulse measured at the finger. J. Hypertens. 24, 1449–1456. doi: 10.1097/01.hjh.0000239277.05068.87

PubMed Abstract | CrossRef Full Text | Google Scholar

Moerman, N., Bonke, B., and Oosting, J. (1993). Awareness and recall during general anesthesia: facts and feelings. Anesthesiology 79, 454–454. doi: 10.1097/00000542-199309000-00007

PubMed Abstract | CrossRef Full Text | Google Scholar

Musizza, B., and Ribaric, S. (2010). Monitoring the depth of anaesthesia. Sensors 10, 10896–10935. doi: 10.3390/s101210896

PubMed Abstract | CrossRef Full Text | Google Scholar

Parati, G. (2005). Blood pressure variability: its measurement and significance in hypertension. J. Hypertens. 23, S19–S25. doi: 10.1097/01.hjh.0000165624.79933.d3

CrossRef Full Text | Google Scholar

Parati, G., Saul, J. P., Di Rienzo, M., and Mancia, G. (1995). Spectral analysis of blood pressure and heart rate variability in evaluating cardiovascular regulation. Hypertension 25, 1276–1286. doi: 10.1161/01.HYP.25.6.1276

CrossRef Full Text | Google Scholar

Rantanen, M., Yli-Hankala, A., Van Gils, M., Yppärilä-Wolters, H., Takala, P., Huiku, M., et al. (2006). Novel multiparameter approach for measurement of nociception at skin incision during general anaesthesia. Br. J. Anaesth. 96, 367–376. doi: 10.1093/bja/ael005

PubMed Abstract | CrossRef Full Text | Google Scholar

Russell, I. F. (1989). Conscious awareness during general anaesthesia: relevance of autonomic signs and isolated arm movements as guides to depth of anaesthesia. Baillières Clin. Anaesthesiol. 3, 511–532. doi: 10.1016/S0950-3501(89)80016-9

CrossRef Full Text | Google Scholar

Seitsonen, E., Korhonen, I., Van Gils, M., Huiku, M., Lötjönen, J., Korttila, K., et al. (2005). EEG spectral entropy, heart rate, photoplethysmography and motor responses to skin incision during sevoflurane anaesthesia. Acta Anaesthesiol. Scand. 49, 284–292. doi: 10.1111/j.1399-6576.2005.00654.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Shin, H. S., Lee, C., and Lee, M. (2009). Adaptive threshold method for the peak detection of photoplethysmographic waveform. Comput. Biol. Med. 39, 1145–1152. doi: 10.1016/j.compbiomed.2009.10.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Takazawa, K., Tanaka, N., Fujita, M., Matsuoka, O., Saiki, T., Aikawa, M., et al. (1998). Assessment of vasoactive agents and vascular aging by the second derivative of photoplethysmogram waveform. Hypertension 32, 365–370. doi: 10.1161/01.HYP.32.2.365

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, A., Yang, L., Liu, C., Cui, J., Li, Y., Yang, X., et al. (2015). Athletic differences in the characteristics of the photoplethysmographic pulse shape: effect of maximal oxygen uptake and maximal muscular voluntary contraction. Biomed Res. Int. 2015:752570. doi: 10.1155/2015/752570

PubMed Abstract | CrossRef Full Text | Google Scholar

Wennervirta, J., Hynynen, M., Koivusalo, A. M., Uutela, K., Huiku, M., and Vakkuri, A. (2008). Surgical stress index as a measure of nociception/antinociception balance during general anesthesia. Acta Anaesthesiol. Scand. 52, 1038–1045. doi: 10.1111/j.1399-6576.2008.01687.x

PubMed Abstract | CrossRef Full Text | Google Scholar