<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Physiol.</journal-id>
<journal-title>Frontiers in Physiology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Physiol.</abbrev-journal-title>
<issn pub-type="epub">1664-042X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1090696</article-id>
<article-id pub-id-type="doi">10.3389/fphys.2023.1090696</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Physiology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Induced pain affects auricular and body biosignals: From cold stressor to deep breathing</article-title>
<alt-title alt-title-type="left-running-head">Rapalis et&#xa0;al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fphys.2023.1090696">10.3389/fphys.2023.1090696</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Rapalis</surname>
<given-names>Andrius</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1373603/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Piartli</surname>
<given-names>Povilas</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Jankauskait&#x117;</surname>
<given-names>Lina</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Marozas</surname>
<given-names>Vaidotas</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/774114/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Kaniusas</surname>
<given-names>Eugenijus</given-names>
</name>
<xref ref-type="aff" rid="aff5">
<sup>5</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/626382/overview"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Biomedical Engineering Institute</institution>, <institution>Kaunas University of Technology</institution>, <addr-line>Kaunas</addr-line>, <country>Lithuania</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Department of Electronics Engineering</institution>, <institution>Faculty of Electrical and Electronics Engineering</institution>, <institution>Kaunas University of Technology</institution>, <addr-line>Kaunas</addr-line>, <country>Lithuania</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Department of Pediatrics</institution>, <institution>Faculty of Medicine</institution>, <institution>Medical Academy</institution>, <institution>Lithuanian University of Health Sciences</institution>, <addr-line>Kaunas</addr-line>, <country>Lithuania</country>
</aff>
<aff id="aff4">
<sup>4</sup>
<institution>Institute of Physiology and Pharmacology</institution>, <institution>Medical Academy</institution>, <institution>Lithuanian University of Health Sciences Kauno Klinikos</institution>, <addr-line>Kaunas</addr-line>, <country>Lithuania</country>
</aff>
<aff id="aff5">
<sup>5</sup>
<institution>Institute of Biomedical Electronics</institution>, <institution>Faculty of Electrical Engineering and Information Technology</institution>, <institution>Vienna University of Technology (TU Wien)</institution>, <addr-line>Vienna</addr-line>, <country>Austria</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/613289/overview">John Allen</ext-link>, Coventry University, United Kingdom</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/919384/overview">Youngsun Kong</ext-link>, University of Connecticut, United States</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1372096/overview">Christopher Austelle</ext-link>, Medical University of South Carolina, United States</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Eugenijus Kaniusas, <email>eugenijus.kaniusas@tuwien.ac.at</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Computational Physiology and Medicine, a section of the journal Frontiers in Physiology</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>17</day>
<month>01</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>14</volume>
<elocation-id>1090696</elocation-id>
<history>
<date date-type="received">
<day>05</day>
<month>11</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>06</day>
<month>01</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Rapalis, Piartli, Jankauskait&#x117;, Marozas and Kaniusas.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Rapalis, Piartli, Jankauskait&#x117;, Marozas and Kaniusas</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>Pain affects every fifth adult worldwide and is a significant health problem. From a physiological perspective, pain is a protective reaction that restricts physical functions and causes responses in physiological systems. These responses are accessible for evaluation <italic>via</italic> recorded biosignals and can be favorably used as feedback in active pain therapy <italic>via</italic> auricular vagus nerve stimulation (aVNS). The aim of this study is to assess the significance of diverse parameters of biosignals with respect to their deflection from cold stressor to deep breathing and their suitability for use as biofeedback in aVNS stimulator. Seventy-eight volunteers participated in two cold pressors and one deep breathing test. Three targeted physiological parameters (<italic>RR</italic> interval of electrocardiogram, cardiac deflection magnitude <italic>Z</italic>
<sub>AC</sub> of ear impedance signal, and cardiac deflection magnitude <italic>PPG</italic>
<sub>AC</sub> of finger photoplethysmogram) and two reference parameters (systolic and diastolic blood pressures <italic>BP</italic>
<sub>S</sub> and <italic>BP</italic>
<sub>D</sub>) were derived and monitored. The results show that the cold water decreases the medians of targeted parameters (by 5.6, 9.3%, and 8.0% of <italic>RR</italic>, <italic>Z</italic>
<sub>AC</sub>, and <italic>PPG</italic>
<sub>AC</sub>, respectively) and increases the medians of reference parameters (by 7.1% and 6.1% of <italic>BP</italic>
<sub>S</sub> and <italic>BP</italic>
<sub>D</sub>, respectively), with opposite changes in deep breathing. Increasing pain level from relatively mild to moderate/strong with cold stressor varies the medians of targeted and reference parameters in the range from 0.5% to 6.0% (e.g., 2.9% for <italic>RR</italic>, <italic>Z</italic>
<sub>AC</sub> and 6.0% for <italic>BP</italic>
<sub>D</sub>). The physiological footprints of painful cold stressor and relaxing deep breathing were shown for auricular and non-auricular biosignals. The investigated targeted parameters can be used as biofeedback to close the loop in aVNS to personalize the pain therapy and increase its compliance.</p>
</abstract>
<kwd-group>
<kwd>auricular bioimpedance</kwd>
<kwd>auricular vagus nerve stimulation</kwd>
<kwd>blood pressure</kwd>
<kwd>cold pressor test</kwd>
<kwd>deep breathing</kwd>
<kwd>physiological biofeedback</kwd>
<kwd>photoplethysmography</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>Acute or chronic pain is one of the main complaints for seeking medical care. According to the International Association for the Study of Pain, pain is defined as &#x201c;unpleasant sensory and emotional experience associated with, or resembling that associated with, actual or potential tissue damage&#x201d; <xref ref-type="bibr" rid="B33">Raja&#xa0;et&#xa0;al. (2020)</xref>. Pain is a protective reaction restricting physical functions with various physiological parameters, such as heart rate, respiratory rate, and arterial blood pressure, which are potential indicators of pain intensity <xref ref-type="bibr" rid="B4">Arbour&#xa0;et&#xa0;al. (2014)</xref>; <xref ref-type="bibr" rid="B8">Cowen&#xa0;et&#xa0;al. (2015)</xref>; <xref ref-type="bibr" rid="B31">Peters and Schmidt (1991)</xref>. The specific changes can be observed by the reactivity and reflexivity of the autonomic nervous, cardiovascular, and respirator systems <xref ref-type="bibr" rid="B22">Kyle and McNeil (2014)</xref>.</p>
<p>Adults and children suffer the pain associated with different medical conditions, undergo different painful procedures, or are referred with acute pain to the emergency department <xref ref-type="bibr" rid="B9">Dahlhamer&#xa0;et&#xa0;al. (2018)</xref>; <xref ref-type="bibr" rid="B21">Keating and Smith (2011)</xref>; <xref ref-type="bibr" rid="B26">Mura&#xa0;et&#xa0;al. (2017)</xref>; <xref ref-type="bibr" rid="B30">Othow&#xa0;et&#xa0;al. (2022)</xref>. Data suggest that 7 out of 10 patients come to the emergency department due to pain <xref ref-type="bibr" rid="B39">Todd&#xa0;et&#xa0;al. (2007)</xref>. Meanwhile, chronic pain affects about 20.5% of adults in the United States <xref ref-type="bibr" rid="B41">Yong&#xa0;et&#xa0;al. (2022)</xref> and about 10%&#x2013;30% in Europe <xref ref-type="bibr" rid="B6">Breivik&#xa0;et&#xa0;al. (2006)</xref>; <xref ref-type="bibr" rid="B34">Reid&#xa0;et&#xa0;al. (2011)</xref>. The meta-analysis showed that the prevalence of chronic pain ranges between 0% and 24% globally <xref ref-type="bibr" rid="B24">Mansfield&#xa0;et&#xa0;al. (2016)</xref>.</p>
<p>Different medications, such as non-steroidal anti-inflammatory drugs, opioids, or others, are used daily against pain. Despite various pain medications and strategies, pain treatment faces many adversities <xref ref-type="bibr" rid="B13">Fishman (2007)</xref>, such as severe side effects, the use of illicit drugs, opioid crisis <xref ref-type="bibr" rid="B37">St.&#xa0;Marie and Broglio (2020)</xref>, peptic ulcers <xref ref-type="bibr" rid="B38">Tai and McAlindon (2021)</xref>, and others. All the more, pain management is a fundamental human right <xref ref-type="bibr" rid="B11">Enright and Goucke (2016)</xref>; <xref ref-type="bibr" rid="B13">Fishman (2007)</xref>.</p>
<p>Vagal nerve stimulation (VNS)&#x2014;as a pain neuromodulation technique, as reviewed in <xref ref-type="bibr" rid="B18">Kaniusas&#xa0;et&#xa0;al. (2019a)</xref>&#x2014;has been investigated in humans and animals. VNS can affect the autonomic nervous system and is an approved treatment for pharmacoresistant depression and drug-resistant epilepsy <xref ref-type="bibr" rid="B27">Nemeroff&#xa0;et&#xa0;al. (2006)</xref>; <xref ref-type="bibr" rid="B29">O&#x2019;Reardon&#xa0;et&#xa0;al. (2006)</xref>. Non-invasive transcutaneous modalities of VNS emerge <xref ref-type="bibr" rid="B7">Busch&#xa0;et&#xa0;al. (2013)</xref>; <xref ref-type="bibr" rid="B28">Nesbitt&#xa0;et&#xa0;al. (2015)</xref>, such as the electrical stimulation of the external surface of the ear innervated by the afferent auricular branch of the vagus nerve, known as auricular vagus nerve stimulation (aVNS). aVNS is performed using miniature electrodes tightly fixed inside the auricular concha. The current intensity is individually adjusted at the beginning of the aVNS session to a level without evoking pain. However, the initial personalization of the intensity of current alone does not ensure adequate vagus nerve stimulation for the relatively long treatment duration (from days to weeks). Here habituation effects, varying physiology, and deterioration of the electrode-tissue interface contribute to this uncertainty in the treatment <xref ref-type="bibr" rid="B5">Bouton (2017)</xref>; <xref ref-type="bibr" rid="B17">Kaniusas (2019)</xref>; <xref ref-type="bibr" rid="B19">Kaniusas&#xa0;et&#xa0;al. (2019b)</xref>. Therefore, aVNS can be hypothesized to avoid under- or over-stimulation, reduce side effects, and save stimulation energy when based on individual physiological biofeedback.</p>
<p>Biofeedback can be assessed using data from internal (in-the-ear) and external (outside-the-ear) sensors, i.e., auricular and non-auricular biosignals. However, it is not known which biosignals and extracted parameters help estimate the balance between the stimulated parasympathetic system and the complementary sympathetic system, which is generally predominant in chronic ailments such as pain. In the ideal case, this balance should be provided to the aVNS stimulator to avoid the disadvantages of the non-personalized aVNS. Thus, easy-to-access biosignals are of high interest which could estimate this balance in favor of the efficiency of aVNS therapy.</p>
<p>The present study proposes a cold stressor as a sympathetically driven stimulus (usually accompanied by acute pain) and deep breathing as a mainly parasympathetically driven stimulus (with relaxing effects) to manipulate the sympathovagal balance from sympathetic to parasympathetic dominance while recording a set of auricular and non-auricular biosignals. The aim of this study is to assess the significance of diverse parameters of biosignals with respect to their deflection from cold stressor to deep breathing and their suitability for use as biofeedback in aVNS stimulator.</p>
</sec>
<sec sec-type="materials|methods" id="s2">
<title>2 Materials and methods</title>
<sec id="s2-1">
<title>2.1 Study population and data acquisition</title>
<p>Seventy-eight healthy volunteers (36 women), 32.6 &#xb1; 10.7&#xa0;years old (range 20&#x2013;64&#xa0;years, with 23 men and 19 women <inline-formula id="inf1">
<mml:math id="m1">
<mml:mo>&#x3c;</mml:mo>
</mml:math>
</inline-formula>30&#xa0;years), with a height of 1.76 &#xb1; 0.1&#xa0;m, a weight of 75.0 &#xb1; 13.6&#xa0;kg, and a body mass index of 24.1 &#xb1; 3.7&#xa0;kg/m<sup>2</sup> participated in the study. All participants met the following criteria: 1) age &#x2265;18&#xa0;years; 2) no chronic pain; 3) no documented cardiovascular, respiratory, diabetes, and depression diseases; 4) no medication with b-blockers or calcium channel antagonists; and 5) no pregnancy or breastfeeding. Participants were instructed to avoid taking painkillers or anti-inflammatory drugs for at least 24&#xa0;h and activities that could affect the cardiovascular system (smoking, coffee, alcohol, physical activity, medication, etc.) for at least 4&#xa0;h before the study.</p>
<p>Data collection took place indoors at the Biomedical Engineering Institute (Kaunas, Lithuania) in a quiet and temperature-controlled (24.0&#xb0;C &#xb1; 1.0&#xb0;C) laboratory at the same time of the day (08:00&#x2013;13:00) to minimize the circadian influence. Four synchronous biosignals were recorded in the study, as illustrated in <xref ref-type="fig" rid="F1">Figure&#xa0;1</xref>: 1) a modified bipolar three-lead electrocardiogram (ECG) signal (sampling rate 2&#xa0;kHz); 2) a red wavelength finger photoplethysmogram (PPG) signal (sampling rate 1&#xa0;kHz) using a proprietary multimodal signals recording system Nautilus II (Biomedical Engineering Institute, Kaunas, Lithuania); 3) an ear impedance signal (at the frequency of 12.5&#xa0;kHz, sampling rate 1&#xa0;kHz) using the data acquisition system Biopac MP150 (Biopac Systems Inc., Aero Camino, Goleta, CA, United States); and 4) arterial blood pressure signal (sampling rate 100&#xa0;Hz) using the non-invasive arterial blood pressure monitoring system CNAP Monitor 500 (CNSystems, Graz, Austria). The subjective/perceived pain was recorded by a volunteer self-report (announced verbally and aloud to an experimenter) using the numerical rating scale NRS (range 0&#x2013;100, with 0 for no pain and 100 for unbearable pain) at least every 30&#xa0;s (or even more often based on a volunteer&#x2019;s initiative).</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Placement of all sensors and electrodes for the recording of electrocardiogram <italic>ECG</italic>(<italic>t</italic>), pulse plethysmogram <italic>PPG</italic>(<italic>t</italic>), auricular impedance <italic>Z</italic>(<italic>t</italic>), and blood pressure <italic>BP</italic>(<italic>t</italic>) signals.</p>
</caption>
<graphic xlink:href="fphys-14-1090696-g001.tif"/>
</fig>
<p>Well-known and effective pain-causing (the cold pressor test) and relaxation (deep breathing) tests were used in the study. Namely, the study protocol consisted of eight phases, as depicted in <xref ref-type="fig" rid="F2">Figure&#xa0;2A</xref>: 1) the first rest phase (Rest &#x23;1) lasting 10&#xa0;min&#xa0;in the sitting position; 2) the warm water for 1&#xa0;min&#xa0;(equalization phase), in which the participant immersed his left hand into warm water (32.0&#xb0;C &#xb1; 0.1&#xb0;C); 3) the first cold water phase (CPT &#x23;1), in which the participant immersed his left hand into cold water (7.0&#xb0;C &#xb1; 0.1&#xb0;C) for 2&#xa0;min&#xa0;or even shorter if the volunteer felt very uncomfortable and voluntarily resumed; 4) the second rest phase (Rest &#x23;2) for 5&#xa0;min&#xa0;where the participant took his left hand out from cold water and rested in the sitting position; 5) the second cold water phase (CPT &#x23;2), in which the participant immersed his left hand into a little less cold water (10.0&#xb0;C &#xb1; 0.1&#xb0;C) for 2&#xa0;min&#xa0;or even shorter if the volunteer felt very uncomfortable and voluntarily resumed; 6) the third rest phase (Rest &#x23;3) for 10&#xa0;min&#xa0;in analogy with Rest &#x23;2; 7) the deep breathing phase (DB) for 1&#xa0;min&#xa0;with the paced breathing rate 6 1/min (paced <italic>via</italic> a monitor and a bar rising/falling periodically every 10&#xa0;s); 8) the fourth rest phase (Rest &#x23;4) for 5&#xa0;min&#xa0;in analogy with Rest &#x23;2. Participants were verbally instructed to immerse their left hand (up to the middle of the forearm) in warm or cold water, indicate their subjective pain level (in cold water), and take out their hand after 1&#xa0;min&#xa0;in warm and 2&#xa0;min&#xa0;in cold water.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>
<bold>(A)</bold> The protocol of the study and <bold>(B)</bold> analyzed intervals.</p>
</caption>
<graphic xlink:href="fphys-14-1090696-g002.tif"/>
</fig>
<p>The study was conducted following the ethical principles of the Declaration of Helsinki and with ethics approval from the Kaunas Region Biomedical Research Ethics Committee (No. BE-2-24), including informed consent and voluntary participation. Personal information was removed from the collected data to ensure participants&#x2019; anonymity.</p>
</sec>
<sec id="s2-2">
<title>2.2 Signal processing and parameters extraction</title>
<p>ECG was filtered using zero-phase Butterworth high-pass and low-pass filters (cut-off frequencies 0.5 and 35&#xa0;Hz, respectively), R waves were detected using the modified Tompkins algorithm <xref ref-type="bibr" rid="B14">Hamilton and Tompkins (1986)</xref>, and <italic>RR</italic> was estimated as the time interval between the successive R peaks. The ear impedance signal reflects local changes in the blood perfusion and blood vessel size, accounting for the local changes in capacitance and resistance. The impedance signal is morphologically similar to PPG so that both <italic>PPG</italic> and <italic>Z</italic> were filtered using high-pass and low-pass zero-phase Butterworth filters (cut-off frequencies 0.5 and 10&#xa0;Hz, respectively). The associated peak and valley fiducial points in <italic>PPG</italic>, <italic>Z</italic>, and <italic>BP</italic> signals were detected in line with the detected R waves of ECG. Five parameters were extracted out of the four recorded biosignals (<xref ref-type="fig" rid="F3">Figure&#xa0;3</xref>): 1) time interval <italic>RR</italic> between R peaks of ECG; 2) cardiac deflection magnitude <italic>PPG</italic>
<sub>AC</sub> of PPG; 3) cardiac deflection magnitude <italic>Z</italic>
<sub>AC</sub> of ear impedance signal; 4) systolic blood pressure <italic>BP</italic>
<sub>S</sub>; and 5) diastolic blood pressure <italic>BP</italic>
<sub>D</sub>. Please note that the analyzed <italic>PPG</italic>
<sub>AC</sub> is mainly related to the pulsatile arterial blood, proportional to the local systolic-diastolic deflection of the blood pressure and the arterial compliance of the vascular wall <xref ref-type="bibr" rid="B16">Kaniusas (2015)</xref>.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Instructive qualitative changes of the interbeat period <italic>RR</italic>, the pulse plethysmography <italic>PPG</italic>, bioimpedance <italic>Z</italic>, and blood pressure <italic>BP</italic> from the first rest phase (Rest &#x23;1), the first cold water phase (CPT &#x23;1), and the deep breathing phase (DB) of a single participant. The maximum reported NRS value of the CPT &#x23;1 was 80.</p>
</caption>
<graphic xlink:href="fphys-14-1090696-g003.tif"/>
</fig>
<p>The entire periods of stimulation phases CPT &#x23;1, CPT &#x23;2, and DB were included in the analysis, only the last 4&#xa0;min&#xa0;of rest phases Rest &#x23;1, Rest &#x23;2, and Rest &#x23;3 were included to avoid the transient influence of the preceding phase (<xref ref-type="fig" rid="F2">Figure&#xa0;2B</xref>). The medians of evaluated parameters from CPT &#x23;1, CPT &#x23;2, and DB phases were compared with the medians of the respective Rest &#x23;1, Rest &#x23;2, and Rest &#x23;3 phases, without any averaging. The analysis was performed using different pain levels, ages, and gender. The pain level threshold was chosen at 40 points, corresponding to mild pain <xref ref-type="bibr" rid="B20">Karcioglu&#xa0;et&#xa0;al. (2018)</xref>. The age threshold was selected at 30&#xa0;years in order to end up with comparably populated groups of men and women.</p>
</sec>
<sec id="s2-3">
<title>2.3 Statistical analysis</title>
<p>The Shapiro-Wilk test was used to assess data normality. Because of the non-normal distribution, the results are summarized using boxplots with medians and quartiles. The Wilcoxon signed-rank test with the Bonferroni&#x2019;s adjustment for dependent samples was used to compute the <italic>p</italic>-value, and statistical significance was set at <italic>p</italic> <inline-formula id="inf2">
<mml:math id="m2">
<mml:mo>&#x3c;</mml:mo>
</mml:math>
</inline-formula>0.05.</p>
</sec>
</sec>
<sec sec-type="results" id="s3">
<title>3 Results</title>
<p>Out of 78 recorded data sets, two <italic>ECG</italic>, thirteen <italic>PPG</italic>, twenty-four <italic>Z</italic>, and six <italic>BP</italic> traces were eliminated from the analysis due to poor quality. Eight participants retreated earlier and did not finish the CPT &#x23;1 phase, i.e., two women (age &#x2264;30&#xa0;years), two women (age <inline-formula id="inf3">
<mml:math id="m3">
<mml:mo>&#x3e;</mml:mo>
</mml:math>
</inline-formula>30&#xa0;years), one man (age &#x2264;30&#xa0;years), and one man (age <inline-formula id="inf4">
<mml:math id="m4">
<mml:mo>&#x3e;</mml:mo>
</mml:math>
</inline-formula>30&#xa0;years). Two other participants did not finish CPT &#x23;2, i.e., one woman and one man, both aged <inline-formula id="inf5">
<mml:math id="m5">
<mml:mo>&#x3e;</mml:mo>
</mml:math>
</inline-formula>30&#xa0;years.</p>
<sec id="s3-1">
<title>3.1 General tendencies</title>
<p>
<xref ref-type="fig" rid="F3">Figure&#xa0;3</xref> illustrates the temporal courses of <italic>RR</italic>, <italic>PPG</italic>, <italic>Z</italic>, and <italic>BP</italic> during phases Rest &#x23;1, CPT &#x23;1, and DB. Compared to Rest &#x23;1, the cold stimulus CPT &#x23;1 shows reduced both <italic>RR</italic> and its variability, as well as reduced cardiac deflection magnitude <italic>PPG</italic>
<sub>AC</sub> of <italic>PPG</italic> and reduced cardiac deflection <italic>Z</italic>
<sub>AC</sub> of <italic>Z</italic>. The associated mean <italic>BP</italic> is larger during CPT &#x23;1 than during Rest &#x23;1. The subsequent DB phase contrasts CPT &#x23;1 in that <italic>PPG</italic>
<sub>AC</sub> and <italic>Z</italic>
<sub>AC</sub> increase in DB. The respiration-related variability of all four <italic>RR</italic>, <italic>PPG</italic>
<sub>AC</sub>, <italic>Z</italic>
<sub>AC</sub>, and <italic>BP</italic> dominates in DB, with the indicated respiration rate <italic>f</italic>
<sub>R</sub> (<xref ref-type="fig" rid="F3">Figure&#xa0;3</xref>).</p>
</sec>
<sec id="s3-2">
<title>3.2 Cold water versus deep breathing</title>
<p>As shown in <xref ref-type="fig" rid="F4">Figure&#xa0;4</xref>, the first cold water stimulation CPT &#x23;1 decreases the median of <italic>RR</italic> (&#x2212;5.5%) and <italic>Z</italic>
<sub>AC</sub> (&#x2212;9.8%) while increasing that of <italic>BP</italic>
<sub>S</sub> (&#x2b;12.6%) and <italic>BP</italic>
<sub>D</sub> (&#x2b;13.4%) of BP, as compared with the first rest phase Rest &#x23;1. Here the associated <italic>PPG</italic>
<sub>AC</sub> remains almost constant (&#x2b;0.9%). The second cold water stimulation CPT &#x23;2 decreases the median of <italic>RR</italic> (&#x2212;5.6%), <italic>PPG</italic>
<sub>AC</sub> (&#x2212;8.0%), and <italic>Z</italic>
<sub>AC</sub> (&#x2212;9.3%) while increasing that of <italic>BP</italic>
<sub>S</sub> (&#x2b;7.1%) and <italic>BP</italic>
<sub>D</sub> (&#x2b;6.1%), as compared with the second rest phase Rest &#x23;2. The subsequent deep breathing DB produces opposite effects: the median of <italic>RR</italic> (&#x2b;1.8%), <italic>PPG</italic>
<sub>AC</sub> (&#x2b;5.1%), and <italic>Z</italic>
<sub>AC</sub> (&#x2b;5.4%) increase, while that of <italic>BP</italic>
<sub>S</sub> (&#x2212;0.9%) and <italic>BP</italic>
<sub>D</sub> (&#x2212;5.6%) decrease, as compared with the third rest phase Rest &#x23;3. The observed changes in DB are significantly different compared to CPT &#x23;2 in all five parameters.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Relative changes <bold>(A)</bold> &#x394;<italic>RR</italic> of <italic>RR</italic>, <bold>(B)</bold> &#x394;<italic>PPG</italic>
<sub>AC</sub> of <italic>PPG</italic>
<sub>AC</sub>, <bold>(C)</bold> &#x394;<italic>Z</italic>
<sub>AC</sub> of <italic>Z</italic>
<sub>AC</sub>, <bold>(D)</bold> &#x394;<italic>BP</italic>
<sub>S</sub> of systolic BP values, and <bold>(E)</bold> &#x394;<italic>BP</italic>
<sub>D</sub> of diastolic BP values from CPT &#x23;1, CPT &#x23;2, and DB as related to the respective Rest &#x23;1, Rest &#x23;2, and Rest &#x23;3. The asterisk &#x201d;&#x2a;&#x201d; indicates significant changes (<italic>p</italic> <inline-formula id="inf6">
<mml:math id="m6">
<mml:mo>&#x3c;</mml:mo>
</mml:math>
</inline-formula>0.05) between CPT &#x23;2 and DB.</p>
</caption>
<graphic xlink:href="fphys-14-1090696-g004.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F5">Figure&#xa0;5</xref> summarizes and contrasts the observed changes for all parameters in CPT &#x23;2 (<xref ref-type="fig" rid="F5">Figure&#xa0;5A</xref>) versus DB (<xref ref-type="fig" rid="F5">Figure&#xa0;5B</xref>), with the indicated interquartile range from 25% to 75%. In line with <xref ref-type="fig" rid="F4">Figure&#xa0;4</xref>, CPT &#x23;2 reduces <italic>RR</italic>, <italic>PPG</italic>
<sub>AC</sub>, and <italic>Z</italic>
<sub>AC</sub> and increases <italic>BP</italic>
<sub>S</sub> and <italic>BP</italic>
<sub>D</sub>, whereas DB causes physiological processes with reversed tendencies, i.e., <italic>RR</italic>, <italic>PPG</italic>
<sub>AC</sub>, and <italic>Z</italic>
<sub>AC</sub> increase while <italic>BP</italic>
<sub>S</sub> and <italic>BP</italic>
<sub>D</sub> decrease.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Medians and interquartile ranges of &#x394;<italic>RR</italic>, &#x394;<italic>PPG</italic>
<sub>AC</sub>, &#x394;<italic>Z</italic>
<sub>AC</sub>, &#x394;<italic>BP</italic>
<sub>S</sub>, and &#x394;<italic>BP</italic>
<sub>D</sub> (compare <xref ref-type="fig" rid="F4">Figure&#xa0;4</xref>) during <bold>(A)</bold> CPT &#x23;2 and <bold>(B)</bold> DB.</p>
</caption>
<graphic xlink:href="fphys-14-1090696-g005.tif"/>
</fig>
</sec>
<sec id="s3-3">
<title>3.3 Pain level differences</title>
<p>
<xref ref-type="fig" rid="F6">Figure&#xa0;6</xref> illustrates the relative changes in the parameters in CPT &#x23;2 for relatively mild pain with the associated NRS &#x2264;40 (<xref ref-type="fig" rid="F6">Figure&#xa0;6A</xref>) in comparison with moderate to strong pain with NRS <inline-formula id="inf7">
<mml:math id="m7">
<mml:mo>&#x3e;</mml:mo>
</mml:math>
</inline-formula>40 (<xref ref-type="fig" rid="F6">Figure&#xa0;6B</xref>). It can be observed that the physiological changes for NRS &#x2264;40 are more closely located to the 100% reference line, i.e., to the values in Rest &#x23;2, than for NRS <inline-formula id="inf8">
<mml:math id="m8">
<mml:mo>&#x3e;</mml:mo>
</mml:math>
</inline-formula>40. Namely, the median &#x394;<italic>RR</italic> decreases by &#x2212;3.8% and &#x2212;6.7% for NRS &#x2264;40 and NRS <inline-formula id="inf9">
<mml:math id="m9">
<mml:mo>&#x3e;</mml:mo>
</mml:math>
</inline-formula>40, respectively; the associated &#x394;<italic>PPG</italic>
<sub>AC</sub> decreases by &#x2212;7.5% and &#x2212;8.0%, &#x394;<italic>Z</italic>
<sub>AC</sub> decreases by &#x2212;6.5% and &#x2212;9.4%, &#x394;<italic>BP</italic>
<sub>S</sub> increases by &#x2b;6.1% and &#x2b;8.1%, and &#x394;<italic>BP</italic>
<sub>D</sub> increases by &#x2b;3.4% and &#x2b;9.4%. When comparing NRS &#x2264;40 and NRS <inline-formula id="inf10">
<mml:math id="m10">
<mml:mo>&#x3e;</mml:mo>
</mml:math>
</inline-formula>40, statistically significant changes are observed in <italic>BP</italic>
<sub>D</sub> only.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Median and interquartile ranges of &#x394;<italic>RR</italic>, &#x394;<italic>PPG</italic>
<sub>AC</sub>, &#x394;<italic>Z</italic>
<sub>AC</sub>, &#x394;<italic>BP</italic>
<sub>S</sub>, and &#x394;<italic>BP</italic>
<sub>D</sub> (compare <xref ref-type="fig" rid="F4">Figure&#xa0;4</xref>) during CPT &#x23;2 for <bold>(A)</bold> mild pain with NRS &#x2264;40 and <bold>(B)</bold> moderate to severe pain with NRS <inline-formula id="inf11">
<mml:math id="m11">
<mml:mo>&#x3e;</mml:mo>
</mml:math>
</inline-formula>40. The distribution of maximum self-report NRS of the CPT &#x23;2 phase is presented in a bar diagram.</p>
</caption>
<graphic xlink:href="fphys-14-1090696-g006.tif"/>
</fig>
</sec>
<sec id="s3-4">
<title>3.4 Gender and age tendencies</title>
<p>The influence of gender and age is depicted in <xref ref-type="fig" rid="F7">Figure&#xa0;7</xref> considering CPT &#x23;2 and DB (compare <xref ref-type="fig" rid="F5">Figure&#xa0;5</xref>). In CPT &#x23;2, the relative values of &#x394;<italic>RR</italic> decrease by 1.5%&#x2013;8.4%, with a minor decrease for young men (<inline-formula id="inf12">
<mml:math id="m12">
<mml:mo>&#x3c;</mml:mo>
</mml:math>
</inline-formula>30 years) and the largest decrease for adult men (&#x2265;30&#xa0;years). Here &#x394;<italic>PPG</italic>
<sub>AC</sub> decreases by 8.4%&#x2013;12.3%, with almost no changes for adult women (&#x2265;30&#xa0;years). &#x394;<italic>Z</italic>
<sub>AC</sub> decreases by 3.9%&#x2013;19.7%, with little changes for adult women and maximum changes for young women (<inline-formula id="inf13">
<mml:math id="m13">
<mml:mo>&#x3c;</mml:mo>
</mml:math>
</inline-formula>30 years). &#x394;<italic>BP</italic>
<sub>S</sub> increases by 5.3%&#x2013;8.2% with minor changes for young women, whereas &#x394;<italic>BP</italic>
<sub>D</sub> increases by 5.1%&#x2013;9.7%, with the largest changes for adult women.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Medians of &#x394;<italic>RR</italic>, &#x394;<italic>PPG</italic>
<sub>AC</sub>, &#x394;<italic>Z</italic>
<sub>AC</sub>, &#x394;<italic>BP</italic>
<sub>S</sub>, and &#x394;<italic>BP</italic>
<sub>D</sub> (compare <xref ref-type="fig" rid="F4">Figure&#xa0;4</xref>) differentiated by age and gender during <bold>(A)</bold> CPT &#x23;2 and <bold>(B)</bold> DB.</p>
</caption>
<graphic xlink:href="fphys-14-1090696-g007.tif"/>
</fig>
<p>During DB, in line with <xref ref-type="fig" rid="F5">Figure&#xa0;5</xref>, the observed changes in all gender and age groups mainly follow the opposite behavior. &#x394;<italic>RR</italic> increases by 2.1%&#x2013;4.5% except for adult men (&#x2265;30&#xa0;years); &#x394;<italic>PPG</italic>
<sub>AC</sub> increases by 4.5%&#x2013;12.1% but also except for adult men; &#x394;<italic>Z</italic>
<sub>AC</sub> increases by 2.5%&#x2013;9.6% with a minor increase for young men (<inline-formula id="inf14">
<mml:math id="m14">
<mml:mo>&#x3c;</mml:mo>
</mml:math>
</inline-formula>30&#xa0;years); &#x394;<italic>BP</italic>
<sub>S</sub> increases very slightly for men (0.5%&#x2013;1.1%) but decreases more strongly for women (2.1%&#x2013;2.8%); &#x394;<italic>BP</italic>
<sub>D</sub> decreases in all cases by 1.0%&#x2013;8.5% with a minor decrease for adult men (&#x2265;30&#xa0;years) and the largest decrease for young women (<inline-formula id="inf15">
<mml:math id="m15">
<mml:mo>&#x3c;</mml:mo>
</mml:math>
</inline-formula>30&#xa0;years).</p>
</sec>
</sec>
<sec sec-type="discussion" id="s4">
<title>4 Discussion</title>
<p>The present study investigates the physiological footprints of auricular and non-auricular biosignals in response to a cold stressor and deep breathing. While a cold stressor is a sympathetically driven stimulus (accompanied by acute pain), deep breathing is a mainly parasympathetically driven stimulus (with relaxing effects). Thus, it was investigated how the opposing sympathetic and parasympathetic stimuli are reflected by the auricular biosignals, namely, its parameter <italic>Z</italic>
<sub>AC</sub>, and by parameters accessible from the auricular biosignals such as <italic>RR</italic> and <italic>PPG</italic>
<sub>AC</sub>. All these three parameters can be used as biofeedback to close the loop in aVNS, i.e., in a targeted stimulation of the parasympathetic system. The closed-loop set-up personalizes aVNS with an expected tendency to avoid over and under-stimulation of the vagus nerve/parasympathetic system. Thus, the closed-loop aVNS may minimize both the energy consumption of the aVNS stimulator and potential side effects (no over-stimulation) while optimizing and personalizing the aVNS therapy (no recruitment of pain fibers), e.g., in chronic pain. Here, the non-auricular biosignals with their parameters <italic>BP</italic>
<sub>S</sub> and <italic>BP</italic>
<sub>D</sub> serve as a necessary reference to monitor stimuli-related vital functions of the body and as an instructive substrate for their comparison with stimuli-related changes in auricular biosignals.</p>
<p>The auricular biosignal <italic>Z</italic>, namely its parameter <italic>Z</italic>
<sub>AC</sub> (<xref ref-type="fig" rid="F1">Figure&#xa0;1</xref>), decreases significantly during the sympathetic stimulus (CPT &#x23;2) as compared with the parasympathetic one (DB) (<xref ref-type="fig" rid="F3">Figures&#xa0;3</xref>&#x2013;<xref ref-type="fig" rid="F5">5</xref>), as well as decreases tendentially with increasing pain perception (<xref ref-type="fig" rid="F6">Figure&#xa0;6</xref>). This behavior indicates the potential suitability of <italic>Z</italic>
<sub>AC</sub> in assessing changes in the balance of the parasympathetic and sympathetic stimuli, or, more generally, in the balance of the parasympathetic and sympathetic systems of the human body (sympathovagal balance). On the other hand, this balance, especially its normalization from a derailed state, is usually a therapeutic target in aVNS when applied to different chronic ailments <xref ref-type="bibr" rid="B18">Kaniusas&#xa0;et&#xa0;al. (2019a)</xref>. Thus <italic>Z</italic>
<sub>AC</sub> can be hypothesized to be reasonable auricular biofeedback for the closed-loop aVNS without using any sensors external to the ear, which may obstruct the patient.</p>
<p>The parameters <italic>RR</italic> and <italic>PPG</italic>
<sub>AC</sub> also reflect sympathovagal balance. <italic>RR</italic> and <italic>PPG</italic>
<sub>AC</sub> decrease significantly during the sympathetic CPT &#x23;2 compared to the parasympathetic DB (<xref ref-type="fig" rid="F3">Figures&#xa0;3</xref>&#x2013;<xref ref-type="fig" rid="F5">5</xref>). While <italic>RR</italic> tends to decrease with increasing pain, the level of <italic>PPG</italic>
<sub>AC</sub> does not (<xref ref-type="fig" rid="F6">Figure&#xa0;6</xref>). Therefore, <italic>RR</italic> and <italic>PPG</italic>
<sub>AC</sub>, the former to a larger extent, can also be hypothesized to be reasonable auricular biofeedback for the closed-loop aVNS targeting a derailed sympathovagal balance. Please note that <italic>RR</italic> could be estimated from the period of the cardiac oscillation of the auricular <italic>Z</italic> (<xref ref-type="fig" rid="F3">Figure&#xa0;3</xref>), whereas <italic>PPG</italic>
<sub>AC</sub> from the cardiac deflection of <italic>PPG</italic> from the earlobe <xref ref-type="bibr" rid="B3">Allen (2007)</xref>. However, limitations in the precision of the estimated <italic>RR</italic> may apply in the former case due to a rather smooth waveform of <italic>Z</italic> in contrast to the spiky R peak of <italic>ECG</italic>. Likewise, limitations in <italic>PPG</italic>
<sub>AC</sub> may apply in the latter case due to a rather central connection of the ear perfusion in contrast to the peripheral perfusion of the finger (<xref ref-type="fig" rid="F1">Figure&#xa0;1</xref>).</p>
<p>The non-auricular parameters <italic>BP</italic>
<sub>S</sub> and <italic>BP</italic>
<sub>D</sub> reflect the sympathovagal balance as well. Both increase significantly during the sympathetic CPT &#x23;2 as compared with the parasympathetic DB (<xref ref-type="fig" rid="F3">Figures&#xa0;3</xref>&#x2013;<xref ref-type="fig" rid="F5">5</xref>), while this increase in CPT &#x23;2 tends to be larger for stronger pain (<xref ref-type="fig" rid="F6">Figure&#xa0;6</xref>). The level of <italic>BP</italic>
<sub>D</sub> appears to depend even stronger on the stimuli-induced sympathovagal balance with the observed changes of 11.7% (from CPT &#x23;2 to DB) in contrast to the associated changes in <italic>BP</italic>
<sub>S</sub> of 8.0% (<xref ref-type="fig" rid="F4">Figure&#xa0;4</xref>). Likewise, the sympathetically-governed vasoconstriction (governing <italic>BP</italic>
<sub>D</sub>) may be more dominant than stroke volume changes (governing <italic>BP</italic>
<sub>S</sub>) <xref ref-type="bibr" rid="B15">Kaniusas (2012)</xref>. This leads to a hypothesis that <italic>BP</italic>
<sub>S</sub> and <italic>BP</italic>
<sub>D</sub> could be used as non-auricular biofeedback for the closed-loop aVNS when external sensors are used outside the ear.</p>
<p>In terms of gender and age, the largest changes from CPT &#x23;2 to DB were shown in <italic>RR</italic> for adult women (&#x2265;30&#xa0;years), <italic>PPG</italic>
<sub>AC</sub> for young women (<inline-formula id="inf16">
<mml:math id="m16">
<mml:mo>&#x3c;</mml:mo>
</mml:math>
</inline-formula>30&#xa0;years), <italic>Z</italic>
<sub>AC</sub> for young women (<inline-formula id="inf17">
<mml:math id="m17">
<mml:mo>&#x3c;</mml:mo>
</mml:math>
</inline-formula>30&#xa0;years), <italic>BP</italic>
<sub>S</sub> for adult women (&#x2265;30&#xa0;years), <italic>BP</italic>
<sub>D</sub> for young women (<inline-formula id="inf18">
<mml:math id="m18">
<mml:mo>&#x3c;</mml:mo>
</mml:math>
</inline-formula>30&#xa0;years). In contrast, the minor changes from CPT &#x23;2 to DB were shown in <italic>RR</italic> for young men (<inline-formula id="inf19">
<mml:math id="m19">
<mml:mo>&#x3c;</mml:mo>
</mml:math>
</inline-formula>30&#xa0;years), <italic>PPG</italic>
<sub>AC</sub> for adult men (&#x2265;30&#xa0;years), <italic>Z</italic>
<sub>AC</sub> for young men (<inline-formula id="inf20">
<mml:math id="m20">
<mml:mo>&#x3c;</mml:mo>
</mml:math>
</inline-formula>30&#xa0;years), <italic>BP</italic>
<sub>S</sub> for adult men (&#x2265;30&#xa0;years), <italic>BP</italic>
<sub>D</sub> for adult men (&#x2265;30&#xa0;years). Overall, men seem to show fewer changes from CPT &#x23;2 to DB than women. This conclusion is in line with previous studies, which conclude that women are more sensitive to pain <xref ref-type="bibr" rid="B12">Fillingim&#xa0;et&#xa0;al. (2009)</xref>; <xref ref-type="bibr" rid="B25">Mogil (2012)</xref>; <xref ref-type="bibr" rid="B32">Popescu&#xa0;et&#xa0;al. (2010)</xref>, but it depends on the method of pain induction and assessment. In most cases, the study also supports the still controversial claims that older individuals are more tolerant of pain and show fewer physiological effects than younger individuals <xref ref-type="bibr" rid="B10">Edwards&#xa0;et&#xa0;al. (2003)</xref>; <xref ref-type="bibr" rid="B35">Riley&#xa0;et&#xa0;al. (2010)</xref>; <xref ref-type="bibr" rid="B36">Rittger&#xa0;et&#xa0;al. (2011)</xref>. However, these statements are very limited in their validity due to the small sample in this study.</p>
<p>A limitation of the present study is the relatively small database of recordings representing the elder part of the population which has tendentially a larger prevalence of suffering pain. Collecting and analyzing a more representative database is planned as a future research direction in the research of the aVNS stimulator. Since the warm water phase immediately preceding CPT &#x23;1 strongly affected the results in CPT &#x23;1, we focused our investigations on the comparison of CPR &#x23;2 and DB, both preceded by rest phases. Another limitation is that the order of the different phases of the protocol were not randomized, especially the order of CPT and DB. Therefore, the results may have been influenced by other factors such as expectation, adaptation, prolonged exposure.</p>
<p>Lastly, it should be noted that the recorded pain level, in contrast to nociception with its physiological encoding and processing of nociceptive stimuli, is a subjective feeling connected with the emotional experience to impeding or actual harm <xref ref-type="bibr" rid="B23">Loeser and Treede (2008)</xref> but also altering autonomic nervous system <xref ref-type="bibr" rid="B40">Woo&#xa0;et&#xa0;al. (2017)</xref>; <xref ref-type="bibr" rid="B2">Adamczyk&#xa0;et&#xa0;al. (2020)</xref>; <xref ref-type="bibr" rid="B1">Abdallah and Geha (2017)</xref>. Thus, the investigated objective characteristics of the autonomic system may be useful for a continuous and objective personalization of aVNS in chronic ailments such as pain.</p>
</sec>
<sec sec-type="conclusion" id="s5">
<title>5 Conclusion</title>
<p>The three parameters <italic>RR</italic>, <italic>PPG</italic>
<sub>AC</sub>, and <italic>Z</italic>
<sub>AC</sub> accessible from auricular biosignals reflect the artificially-induced stimuli with sympathetic or parasympathetic dominance and thus the sympathovagal balance derailed in pain and other chronic conditions. Therefore, auricular biosignals can be used as biofeedback to close the loop in auricular vagus nerve stimulation to personalize the strength and timing of the stimulation in favor of therapy, patient compliance, and resourceful energy use.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s6">
<title>Data availability statement</title>
<p>The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.</p>
</sec>
<sec id="s7">
<title>Ethics statement</title>
<p>The studies involving human participants were reviewed and approved by Kaunas Region Biomedical Research Ethics Committee (No. BE-2-24). The patients/participants provided their written informed consent to participate in this study.</p>
</sec>
<sec id="s8">
<title>Author contributions</title>
<p>AR contributed to the study design, collected data, developed algorithms, conducted data analysis, interpreted the results, and wrote and revised the manuscript. PP collected data, contributed to algorithm development, and revised the manuscript. LJ interpreted the results and revised the manuscript. VM contributed to the study design and revised the manuscript. EK initiated and designed the study, interpreted the results, and designed, wrote, and revised the manuscript. All authors have read the submitted manuscript and approved the final version.</p>
</sec>
<sec id="s9">
<title>Funding</title>
<p>This research was funded by the Research and Innovation Fund of Kaunas University of Technology (project grant No. PP2021/5) and the Research Fund of Lithuanian University of Health Sciences (2021-JV-00006). This study received funding from the European Union&#x2019;s Horizon 2020 research and innovation program under grant agreement No. 880603 (SzeleSTIM GmbH). The funder SzeleSTIM GmbH had the following involvement with the study: the study design and collection of data. All authors declare no other competing interests.</p>
</sec>
<ack>
<p>The authors acknowledge TU Wien Bibliothek for financial support through its Open Access Funding Programme. We thank all the volunteers for their participation in the study.</p>
</ack>
<sec sec-type="COI-statement" id="s10">
<title>Conflict of interest</title>
<p>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.</p>
<p>The handling editor JA declared a shared research group VascAgeNet WG3 with the author VM at the time of review.</p>
</sec>
<sec sec-type="disclaimer" id="s11">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Abdallah</surname>
<given-names>C. G.</given-names>
</name>
<name>
<surname>Geha</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Chronic pain and chronic stress: Two sides of the same coin?</article-title> <source>Chronic Stress</source> <volume>1</volume>, <fpage>2470547017704763</fpage>. <pub-id pub-id-type="doi">10.1177/2470547017704763</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Adamczyk</surname>
<given-names>A. K.</given-names>
</name>
<name>
<surname>Ligeza</surname>
<given-names>T. S.</given-names>
</name>
<name>
<surname>Wyczesany</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>The dynamics of pain reappraisal: The joint contribution of cognitive change and mental load</article-title>. <source>Cognitive, Affect. Behav. Neurosci.</source> <volume>20</volume>, <fpage>276</fpage>&#x2013;<lpage>293</lpage>. <pub-id pub-id-type="doi">10.3758/s13415-020-00768-7</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Allen</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Photoplethysmography and its application in clinical physiological measurement</article-title>. <source>Physiol. Meas.</source> <volume>28</volume>, <fpage>R1</fpage>&#x2013;<lpage>R39</lpage>. <pub-id pub-id-type="doi">10.1088/0967-3334/28/3/R01</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Arbour</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Choini&#xe8;re</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Topolovec-Vranic</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Loiselle</surname>
<given-names>C. G.</given-names>
</name>
<name>
<surname>G&#xe9;linas</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Can fluctuations in vital signs be used for pain assessment in critically ill patients with a traumatic brain injury?</article-title> <source>Pain Res. Treat.</source> <volume>2014</volume>, <fpage>175794</fpage>. <pub-id pub-id-type="doi">10.1155/2014/175794</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bouton</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Cracking the neural code, treating paralysis and the future of bioelectronic medicine</article-title>. <source>J. Intern. Med.</source> <volume>282</volume>, <fpage>37</fpage>&#x2013;<lpage>45</lpage>. <pub-id pub-id-type="doi">10.1111/joim.12610</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Breivik</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Collett</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Ventafridda</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Cohen</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Gallacher</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Survey of chronic pain in Europe: Prevalence, impact on daily life, and treatment</article-title>. <source>Eur. J. Pain</source> <volume>10</volume>, <fpage>287</fpage>&#x2013;<lpage>333</lpage>. <pub-id pub-id-type="doi">10.1016/j.ejpain.2005.06.009</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Busch</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Zeman</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Heckel</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Menne</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Ellrich</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Eichhammer</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>The effect of transcutaneous vagus nerve stimulation on pain perception-an experimental study</article-title>. <source>Brain Stimul.</source> <volume>6</volume>, <fpage>202</fpage>&#x2013;<lpage>209</lpage>. <pub-id pub-id-type="doi">10.1016/j.brs.2012.04.006</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cowen</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Stasiowska</surname>
<given-names>M. K.</given-names>
</name>
<name>
<surname>Laycock</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Bantel</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Assessing pain objectively: The use of physiological markers</article-title>. <source>Anaesthesia</source> <volume>70</volume>, <fpage>828</fpage>&#x2013;<lpage>847</lpage>. <pub-id pub-id-type="doi">10.1111/anae.13018</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dahlhamer</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Lucas</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zelaya</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Nahin</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Mackey</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>DeBar</surname>
<given-names>L.</given-names>
</name>
<etal/>
</person-group> (<year>2018</year>). <article-title>Prevalence of chronic pain and high-impact chronic pain among adults &#x2014; United States, 2016</article-title>. <source>MMWR. Morb. Mortal. Wkly. Rep.</source> <volume>67</volume>, <fpage>1001</fpage>&#x2013;<lpage>1006</lpage>. <pub-id pub-id-type="doi">10.15585/mmwr.mm6736a2</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Edwards</surname>
<given-names>R. R.</given-names>
</name>
<name>
<surname>Fillingim</surname>
<given-names>R. B.</given-names>
</name>
<name>
<surname>Ness</surname>
<given-names>T. J.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Age-related differences in endogenous pain modulation: A comparison of diffuse noxious inhibitory controls in healthy older and younger adults</article-title>. <source>Pain</source> <volume>101</volume>, <fpage>155</fpage>&#x2013;<lpage>165</lpage>. <pub-id pub-id-type="doi">10.1016/S0304-3959(02)00324-X</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Enright</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Goucke</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>The global burden of pain: The tip of the iceberg?</article-title> <source>Anesth. Analgesia</source> <volume>123</volume>, <fpage>529</fpage>&#x2013;<lpage>530</lpage>. <pub-id pub-id-type="doi">10.1213/ANE.0000000000001519</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fillingim</surname>
<given-names>R. B.</given-names>
</name>
<name>
<surname>King</surname>
<given-names>C. D.</given-names>
</name>
<name>
<surname>Ribeiro-Dasilva</surname>
<given-names>M. C.</given-names>
</name>
<name>
<surname>Rahim-Williams</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Riley</surname>
<given-names>J. L.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Sex, gender, and pain: A review of recent clinical and experimental findings</article-title>. <source>J. Pain</source> <volume>10</volume>, <fpage>447</fpage>&#x2013;<lpage>485</lpage>. <pub-id pub-id-type="doi">10.1016/j.jpain.2008.12.001</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fishman</surname>
<given-names>S. M.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Recognizing pain management as a human right: A first step</article-title>. <source>Anesth. Analgesia</source> <volume>105</volume>, <fpage>8</fpage>&#x2013;<lpage>9</lpage>. <pub-id pub-id-type="doi">10.1213/01.ane.0000267526.37663.41</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hamilton</surname>
<given-names>P. S.</given-names>
</name>
<name>
<surname>Tompkins</surname>
<given-names>W. J.</given-names>
</name>
</person-group> (<year>1986</year>). <article-title>Quantitative investigation of qrs detection rules using the mit/bih arrhythmia database</article-title>. <source>IEEE Trans. Biomed. Eng. BME-</source> <volume>33</volume>, <fpage>1157</fpage>&#x2013;<lpage>1165</lpage>. <pub-id pub-id-type="doi">10.1109/TBME.1986.325695</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Kaniusas</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2012</year>). <source>Biomedical signals and sensors I</source>. <publisher-name>Springer Berlin, Heidelberg</publisher-name>. <pub-id pub-id-type="doi">10.1007/978-3-642-24843-6</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Kaniusas</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2015</year>). <source>Biomedical signals and sensors II</source>. <publisher-name>Springer Berlin, Heidelberg</publisher-name>. <pub-id pub-id-type="doi">10.1007/978-3-662-45106-9</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Kaniusas</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2019</year>). <source>Biomedical signals and sensors III</source>. <publisher-name>Springer Cham</publisher-name>. <pub-id pub-id-type="doi">10.1007/978-3-319-74917-4</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kaniusas</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Kampusch</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Tittgemeyer</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Panetsos</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Gines</surname>
<given-names>R. F.</given-names>
</name>
<name>
<surname>Papa</surname>
<given-names>M.</given-names>
</name>
<etal/>
</person-group> (<year>2019a</year>). <article-title>Current directions in the auricular vagus nerve stimulation i &#x2013; A physiological perspective</article-title>. <source>Front. Neurosci.</source> <volume>13</volume>, <fpage>854</fpage>&#x2013;<lpage>923</lpage>. <pub-id pub-id-type="doi">10.3389/fnins.2019.00854</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kaniusas</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Kampusch</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Tittgemeyer</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Panetsos</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Gines</surname>
<given-names>R. F.</given-names>
</name>
<name>
<surname>Papa</surname>
<given-names>M.</given-names>
</name>
<etal/>
</person-group> (<year>2019b</year>). <article-title>Current directions in the auricular vagus nerve stimulation ii - an engineering perspective</article-title>. <source>Front. Neurosci.</source> <volume>13</volume>, <fpage>772</fpage>. <pub-id pub-id-type="doi">10.3389/fnins.2019.00772</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Karcioglu</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Topacoglu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Dikme</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Dikme</surname>
<given-names>O.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>A systematic review of the pain scales in adults: Which to use?</article-title> <source>Am. J. Emerg. Med.</source> <volume>36</volume>, <fpage>707</fpage>&#x2013;<lpage>714</lpage>. <pub-id pub-id-type="doi">10.1016/j.ajem.2018.01.008</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Keating</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Smith</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Acute pain in the emergency department: The challenges</article-title>. <source>Rev. Pain</source> <volume>5</volume>, <fpage>13</fpage>&#x2013;<lpage>17</lpage>. <pub-id pub-id-type="doi">10.1177/204946371100500304</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kyle</surname>
<given-names>B. N.</given-names>
</name>
<name>
<surname>McNeil</surname>
<given-names>D. W.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Autonomic arousal and experimentally induced pain: A critical review of the literature</article-title>. <source>Pain Res. Manag.</source> <volume>19</volume>, <fpage>159</fpage>&#x2013;<lpage>167</lpage>. <pub-id pub-id-type="doi">10.1155/2014/536859</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Loeser</surname>
<given-names>J. D.</given-names>
</name>
<name>
<surname>Treede</surname>
<given-names>R. D.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>The kyoto protocol of iasp basic pain terminology</article-title>. <source>Pain</source> <volume>137</volume>, <fpage>473</fpage>&#x2013;<lpage>477</lpage>. <pub-id pub-id-type="doi">10.1016/j.pain.2008.04.025</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mansfield</surname>
<given-names>K. E.</given-names>
</name>
<name>
<surname>Sim</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Jordan</surname>
<given-names>J. L.</given-names>
</name>
<name>
<surname>Jordan</surname>
<given-names>K. P.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>A systematic review and meta-analysis of the prevalence of chronic widespread pain in the general population</article-title>. <source>Pain</source> <volume>157</volume>, <fpage>55</fpage>&#x2013;<lpage>64</lpage>. <pub-id pub-id-type="doi">10.1097/j.pain.0000000000000314</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mogil</surname>
<given-names>J. S.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Sex differences in pain and pain inhibition: Multiple explanations of a controversial phenomenon</article-title>. <source>Nat. Rev. Neurosci.</source> <volume>13</volume>, <fpage>859</fpage>&#x2013;<lpage>866</lpage>. <pub-id pub-id-type="doi">10.1038/nrn3360</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mura</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Serra</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Marinangeli</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Patti</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Musu</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Piras</surname>
<given-names>I.</given-names>
</name>
<etal/>
</person-group> (<year>2017</year>). <article-title>Prospective study on prevalence, intensity, type, and therapy of acute pain in a second-level urban emergency department</article-title>. <source>J. Pain Res.</source> <volume>10</volume>, <fpage>2781</fpage>&#x2013;<lpage>2788</lpage>. <pub-id pub-id-type="doi">10.2147/JPR.S137992</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nemeroff</surname>
<given-names>C. B.</given-names>
</name>
<name>
<surname>Mayberg</surname>
<given-names>H. S.</given-names>
</name>
<name>
<surname>Krahl</surname>
<given-names>S. E.</given-names>
</name>
<name>
<surname>McNamara</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Frazer</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Henry</surname>
<given-names>T. R.</given-names>
</name>
<etal/>
</person-group> (<year>2006</year>). <article-title>Vns therapy in treatment-resistant depression: Clinical evidence and putative neurobiological mechanisms</article-title>. <source>Neuropsychopharmacology</source> <volume>31</volume>, <fpage>1345</fpage>&#x2013;<lpage>1355</lpage>. <pub-id pub-id-type="doi">10.1038/sj.npp.1301082</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nesbitt</surname>
<given-names>A. D.</given-names>
</name>
<name>
<surname>Marin</surname>
<given-names>J. C.</given-names>
</name>
<name>
<surname>Tompkins</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Ruttledge</surname>
<given-names>M. H.</given-names>
</name>
<name>
<surname>Goadsby</surname>
<given-names>P. J.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Initial use of a novel noninvasive vagus nerve stimulator for cluster headache treatment</article-title>. <source>Neurology</source> <volume>84</volume>, <fpage>1249</fpage>&#x2013;<lpage>1253</lpage>. <pub-id pub-id-type="doi">10.1212/WNL.0000000000001394</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>O&#x2019;Reardon</surname>
<given-names>J. P.</given-names>
</name>
<name>
<surname>Cristancho</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Peshek</surname>
<given-names>A. D.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Vagus nerve stimulation (vns) and treatment of depression: To the brainstem and beyond</article-title>. <source>Psychiatry (Edgmont (Pa, Townsh.</source> <volume>3</volume>, <fpage>54</fpage>&#x2013;<lpage>63</lpage>.</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Othow</surname>
<given-names>C. O.</given-names>
</name>
<name>
<surname>Ferede</surname>
<given-names>Y. A.</given-names>
</name>
<name>
<surname>Tawuye</surname>
<given-names>H. Y.</given-names>
</name>
<name>
<surname>Aytolign</surname>
<given-names>H. A.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>The magnitude and associated factors of post-operative pain among adult patients</article-title>. <source>Ann. Med. Surg.</source> <volume>81</volume>, <fpage>104406</fpage>. <pub-id pub-id-type="doi">10.1016/j.amsu.2022.104406</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Peters</surname>
<given-names>M. L.</given-names>
</name>
<name>
<surname>Schmidt</surname>
<given-names>A. J.</given-names>
</name>
</person-group> (<year>1991</year>). <article-title>Psychophysiological responses to repeated acute pain stimulation in chronic low back pain patients</article-title>. <source>J. Psychosomatic Res.</source> <volume>35</volume>, <fpage>59</fpage>&#x2013;<lpage>74</lpage>. <pub-id pub-id-type="doi">10.1016/0022-3999(91)90007-B</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Popescu</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Leresche</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Truelove</surname>
<given-names>E. L.</given-names>
</name>
<name>
<surname>Drangsholt</surname>
<given-names>M. T.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Gender differences in pain modulation by diffuse noxious inhibitory controls: A systematic review</article-title>. <source>Pain</source> <volume>150</volume>, <fpage>309</fpage>&#x2013;<lpage>318</lpage>. <pub-id pub-id-type="doi">10.1016/j.pain.2010.05.013</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Raja</surname>
<given-names>S. N.</given-names>
</name>
<name>
<surname>Carr</surname>
<given-names>D. B.</given-names>
</name>
<name>
<surname>Cohen</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Finnerup</surname>
<given-names>N. B.</given-names>
</name>
<name>
<surname>Flor</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Gibson</surname>
<given-names>S.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>The revised international association for the study of pain definition of pain: Concepts, challenges, and compromises</article-title>. <source>Pain</source> <volume>161</volume>, <fpage>1976</fpage>&#x2013;<lpage>1982</lpage>. <pub-id pub-id-type="doi">10.1097/j.pain.0000000000001939</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Reid</surname>
<given-names>K. J.</given-names>
</name>
<name>
<surname>Harker</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Bala</surname>
<given-names>M. M.</given-names>
</name>
<name>
<surname>Truyers</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Kellen</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Bekkering</surname>
<given-names>G. E.</given-names>
</name>
<etal/>
</person-group> (<year>2011</year>). <article-title>Epidemiology of chronic non-cancer pain in Europe: Narrative review of prevalence, pain treatments and pain impact</article-title>. <source>Curr. Med. Res. Opin.</source> <volume>27</volume>, <fpage>449</fpage>&#x2013;<lpage>462</lpage>. <pub-id pub-id-type="doi">10.1185/03007995.2010.545813</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Riley</surname>
<given-names>J. L.</given-names>
</name>
<name>
<surname>King</surname>
<given-names>C. D.</given-names>
</name>
<name>
<surname>Wong</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Fillingim</surname>
<given-names>R. B.</given-names>
</name>
<name>
<surname>Mauderli</surname>
<given-names>A. P.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Lack of endogenous modulation and reduced decay of prolonged heat pain in older adults</article-title>. <source>Pain</source> <volume>150</volume>, <fpage>153</fpage>&#x2013;<lpage>160</lpage>. <pub-id pub-id-type="doi">10.1016/j.pain.2010.04.020</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rittger</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Rieber</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Breithardt</surname>
<given-names>O. A.</given-names>
</name>
<name>
<surname>D&#xfc;cker</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Schmidt</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Abbara</surname>
<given-names>S.</given-names>
</name>
<etal/>
</person-group> (<year>2011</year>). <article-title>Influence of age on pain perception in acute myocardial ischemia: A possible cause for delayed treatment in elderly patients</article-title>. <source>Int. J. Cardiol.</source> <volume>149</volume>, <fpage>63</fpage>&#x2013;<lpage>67</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijcard.2009.11.046</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>St. Marie</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Broglio</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Managing pain in the setting of opioid use disorder</article-title>. <source>Pain Manag. Nurs.</source> <volume>21</volume>, <fpage>26</fpage>&#x2013;<lpage>34</lpage>. <pub-id pub-id-type="doi">10.1016/j.pmn.2019.08.003</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tai</surname>
<given-names>F. W. D.</given-names>
</name>
<name>
<surname>McAlindon</surname>
<given-names>M. E.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Non-steroidal anti-inflammatory drugs and the gastrointestinal tract</article-title>. <source>Clin. Med. J. R. Coll. Physicians Lond.</source> <volume>21</volume>, <fpage>131</fpage>&#x2013;<lpage>134</lpage>. <pub-id pub-id-type="doi">10.7861/CLINMED.2021-0039</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Todd</surname>
<given-names>K. H.</given-names>
</name>
<name>
<surname>Ducharme</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Choiniere</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Crandall</surname>
<given-names>C. S.</given-names>
</name>
<name>
<surname>Fosnocht</surname>
<given-names>D. E.</given-names>
</name>
<name>
<surname>Homel</surname>
<given-names>P.</given-names>
</name>
<etal/>
</person-group> (<year>2007</year>). <article-title>Pain in the emergency department: Results of the pain and emergency medicine initiative (pemi) multicenter study</article-title>. <source>J. Pain</source> <volume>8</volume>, <fpage>460</fpage>&#x2013;<lpage>466</lpage>. <pub-id pub-id-type="doi">10.1016/j.jpain.2006.12.005</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Woo</surname>
<given-names>C. W.</given-names>
</name>
<name>
<surname>Schmidt</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Krishnan</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Jepma</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Roy</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Lindquist</surname>
<given-names>M. A.</given-names>
</name>
<etal/>
</person-group> (<year>2017</year>). <article-title>Quantifying cerebral contributions to pain beyond nociception</article-title>. <source>Nat. Commun.</source> <volume>8</volume>, <fpage>14211</fpage>. <pub-id pub-id-type="doi">10.1038/ncomms14211</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yong</surname>
<given-names>R. J.</given-names>
</name>
<name>
<surname>Mullins</surname>
<given-names>P. M.</given-names>
</name>
<name>
<surname>Bhattacharyya</surname>
<given-names>N.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Prevalence of chronic pain among adults in the United States</article-title>. <source>Pain</source> <volume>163</volume>, <fpage>E328</fpage>&#x2013;<lpage>E332</lpage>. <pub-id pub-id-type="doi">10.1097/j.pain.0000000000002291</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>