Abstract
The utility of heart rate variability (HRV) for characterizing psychological stress is primarily impacted by methodological considerations such as study populations, experienced versus induced stress, and method of stress assessment. Here, we review studies on the associations between HRV and psychological stress, examining the nature of stress, ways stress was assessed, and HRV metrics used. The review was performed according to the PRISMA guidelines on select databases. Studies that examined the HRV-stress relationship via repeated measurements and validated psychometric instruments were included (n = 15). Participant numbers and ages ranged between 10 and 403 subjects and 18 and 60 years, respectively. Both experimental (n = 9) and real-life stress (n = 6) have been explored. While RMSSD was the most reported HRV metric (n = 10) significantly associated with stress, other metrics, including LF/HF (n = 7) and HF power (n = 6) were also reported. Various linear and nonlinear HRV metrics have been utilized, with nonlinear metrics used less often. The most frequently used psychometric instrument was the State-Trait Anxiety Inventory (n = 10), though various other instruments have been reported. In conclusion, HRV is a valid measure of the psychological stress response. Standard stress induction and assessment protocols combined with validated HRV measures in different domains will improve the validity of findings.
Introduction
With advances in wearable technologies, it has become much easier to monitor the general population’s physiological health and well-being objectively. Heart rate variability (HRV), referring to variations in the time intervals between consecutive heart beats [1, 2], has emerged as a non-invasive tool to estimate psychological states such as fatigue, stress, anxiety, and burnout and is considered an indirect indicator of general mental well-being [3, 4]. HRV refers to the fluctuations in time between consecutive heartbeat cycles. HRV biofeedback, widely adopted in various therapeutic interventions, uses real-time electronic feedback of the moment-to-moment changes in HRV to achieve therapeutic goals [5, 6]. Several studies in psychology and neuroscience have confirmed that HRV measurements reflect the continuous interplay between the sympathetic and parasympathetic branches of the autonomic nervous system that maintain homeostasis of physiological arousal [7‒9]. A strict periodicity of reduced HRV is not a sign of good health which often reflects the experience of physiological distress and can be associated with several pathological conditions [10]. As such, novel wearable devices increasingly incorporate non-invasive physiological monitoring methods to capture real-time HRV-based stress changes. Studies comparing the quality of HRV measurements acquired from conventional ECG and those obtained from wearable devices have revealed that HRV obtained from trackers and wearable devices resulted in a small acceptable error, the method is valid, more practical, and cost-effective in tracking stress and well-being [11]. Given the vast interest toward using HRV metrics in unobtrusive stress tracking, studies that report HRV and its associations with stress assessed with stress-screening tools warrant a thorough examination.
Any physical or psychological stimuli can invoke a stress response due to the disruption of homeostasis. Or in other words, an environmental demand (stressor) translates both into psychological (self-appraisal or perception) and biological responses (stress response) in the body [12]. Response to these stressors is shown to be mediated by a complex interplay of nervous, endocrine, and immune mechanisms, which in turn produce objective physiological changes such as altered HRV [13]. The adaptability of an individual to these stressors and their physiological responses are influenced by several factors including the nature of stressors, intensity, frequency, duration of exposure to the stressful stimuli, and other health conditions [14].
Over the last two decades, several original research and review articles have been published examining the relationship between HRV and different types of stressful situations. Heightened occupational stress is associated with lowered HRV, specifically with reduced parasympathetic activation [15‒17]. Low parasympathetic activity, characterized by a decrease in high-frequency power and an increase in low-frequency power, was reported in a review as the most common factor associated with changes in stress [18]. As a hallmark of depression, blunted stress reactivity manifested as reduced high-frequency HRV has also been reported [19]. Neuroimaging studies on cerebral blood flow have shown that threat and safety-related acute stress significantly modulates HRV and supports the idea that the vagus nerve serves as a structural and functional link between the brain and the heart [20].
Psychological stress occurs when a stressful condition causes negative affective states (e.g., feelings of anxiety and depression) in an individual where they perceive that the environment’s demand exceeds their capacity [12] and invoke physiological changes in the body. Studies have shown that psychological stressors such as public speaking or interpersonal conflict have been shown to evoke increases in blood pressure and heart rate [21]. While previous reviews support the notion that changes in HRV is a valuable physiological indicator of psychological stress, it is unclear whether psychometrically validated instruments underpin stress assessments in all the considered studies. If quantified using psychometrically validated instruments, the perceived subjective component of stress [12] is more likely to reduce errors and biases associated with self-reports [22]. Second, findings from studies based on different stress contexts or demands such as physiological (e.g., pain, hunger), psychosocial, cognitive load, anticipatory etc., (e.g., mathematical tasks, social exclusion, achievement, or competitive situations) and different protocols of stress inducement were all congregated in the synthesis and meta-analysis of reviews as though the autonomic responses experienced in all these scenarios are homogenous [18]. This aspect requires some dismantling to investigate how responses vary with contexts and to explore any specific associations between the contexts and certain HRV metrics. Such metrics could be further explored within those contexts [23]. Third, studies do not always include a comparison with a stress-assessed at baseline or restful condition, which is also important to investigate and address, again based on the above-mentioned reasons that emphasize the need for quantifying and validating the intra-individual perceived stress changes and statistically comparing them with their baseline. This scoping review was carried out to address these gaps by (i) reviewing HRV-based studies with a well-designed within-subjects design that associated HRV with psychological stress states, (ii) including studies that used a validated psychometric instrument to elicit stress levels to systematically compare stress and HRV responses between stressful states, and (iii) including studies that provide a comparison of stress state measures with a no-stress baseline or restful condition. In the review process, as secondary objectives, we also explored whether included studies reported the mediating and/or moderating effects of underlying chronic stress, short-term HRV analysis windows, and the impact of breathing on HRV.
Methods
Search Strategy
A detailed stepwise search strategy was developed in collaboration with a research librarian at Flinders University. Two primary key terms (psychological stress and heart rate variability) were developed as “concepts ”based on the research question. Detailed alternate terms and synonyms for each key concept were also identified as free terms and used to search articles’ titles, abstracts, keywords, and key concept fields in the following electronic databases: Medline, PsycINFO, Scopus, and Web of Science. Detailed search strategies are provided in the online supplement (online suppl. Table. S3; for all online suppl. material, see https://doi.org/10.1159/000530376).
Eligibility Criteria
Studies that (i) examined the relationship between HRV and psychological stress levels, (ii) with measures of HRV acquired in two or more time points that are established or hypothesized within the protocol as different stress states, (iii) used validation instruments to quantify mental or psychological stress, (iv) included healthy adult population, and (v) written in English were included. Conference articles, any form of review articles, case studies, letters to the editors, dissertations, and books were excluded.
Studies that purposefully sampled clinical populations with a diagnosed condition (e.g., insomnia, dysphoria, clinically depressed states, metabolic syndrome, irritable bowel syndrome, neuroticism, psychosis, panic disorder, eating disorders, patients who underwent surgery, hypertensive women, borderlines) were excluded. Thus, “healthy adult population” here refers to populations recruited from the general community without using any sampling criteria based on clinical status, and not based on whether the “healthy” status was ascertained with a psychometric questionnaire.
Data Extraction
A total of 15 journal articles evaluating HRV in response to psychological stress and meeting the eligibility criteria mentioned above were included in this review. Study characteristics and findings were extracted from every publication and listed. During the process of data acquisition and integration, this list was completed with additional variables that appeared relevant and eventually covered the following information: author, publication year, sample size, gender, age range, type of stress or experimental manipulation, duration of analyzed task period, the time between repeated measures, the context of the study, outcome measures, and key findings. A mixed methods appraisal tool was used to appraise the quality of the studies and is included in the online supplement (online suppl. Table S2).
Results
Included Studies
Using the search criteria described in the methods section, electronic databases were searched for HRV research in psychological stress. The initial search fetched 8,297 studies as of December 2020; these were reduced to 5,638 once 2,659 duplicates were removed, to 90 after screening for titles and abstracts, and then to 15 after reviewing the full texts for study methods (Fig. 1). A detailed list of the studies included in the present review is shown in Table 1.
Study characteristics of included articles
. | Study . | N (F) . | Age, years(µ±SD or range) . | Repeated measures (time points) . | Duration between time points . | Stress type . | Context . | HRV measures significantly associated with stress and reported direction of relationship . |
---|---|---|---|---|---|---|---|---|
1 | Carnevali et al. [31] (2018) | 42 (22) | 23.7±1.1 (F)23.0±1.0 (M) | 3 | 12.8 months between T0 and T120.7 months between T1 and T2 | Real life, chronic | Depressive symptoms and rumination over a period of 3 years | (↓) RMSSD |
2 | Brugnera et al. [42] (2018) | 60 (51.7%) | 25.6±3.8 | 5 (baseline, three tasks and recovery) | 5–7 min; entire procedure was finished in an hour | Experimentally induced,acute | Mental stress before, during and after MIST, Stroop and Speech tasks | MIST: (↓) SDNN, RMSSD, TP, HFpow, LFpow, SD1, SD2, D2(↑) SI, SampEn, α1Speech: (↓) RMSSD, Hfnu, SD1, SampEn(↑) Lfnu, LF/HF, α1Stroop: (↓)Hfnu, SampEn(↑) LFpow, Lfnu, LF/HFSD2, α1 |
3 | Cervantes Blasquez et al. [28] (2009) | 10 (6) swimmers | 47 | 2 (training, competition) | 7 days | Real-life, anticipatory, acute | State anxiety before entering a swimming competition | (↓) RMSSD, HF, Hfnu, SD1(↑) LF/HF |
4 | Delaney [36] (2000) | 30 (16) | 30.9±3.9 (F)34.4±8.7 (M) | 2 (baseline, task) | Not mentioned clearly; stress situation was followed by 10 min of baseline rest period | Experimental, acute | Emotional state and muscle tension during a Stroop test and a mental arithmetic component | (↓) SDNN, RMSSD, PNN50, TP, HFpow, LFpow, Hfnu, LFnu(↑) LF/HF |
5 | Dimitriev et al. [27] (2016) | 96 (81) students | 20.5±0.1 | 2 (rest day, before an exam) | Not clearly mentioned; article states as “two different days” | Real-life, Acute | State anxiety just before an academic examination | (↓) ApEn, SampEn, SD1, SD1/SD2, CCM, LLE.(↑) α1 |
6 | Kanthak et al. [29] (2017) | 403 (66.6%) | 42.2±11.2 | 3 (Procedure, recovery, post-stress rest) | 20 min | Real life, acute | Emotional stress during a blood sampling procedure | (↓) RMSSD with emotional exhaustion component of burnout |
7 | Logan et al. [35] (2020) | 85 (100%) | 28.8±9.8 | 2 (before and after stress task) | Sequential TSST | Experimental, acute | Psychological distress and anxiety associated with TSST | No significant changes in HF-HRV in the pre and post TSST |
8 | Lucini et al. [26] (2002) | 30*university students | 22.0±1.0 | 2 (before an exam, restful holiday) | 3 months | Real life, acute | State anxiety before an academic exam | (↓) HFnu, alpha index(↑) LFnu, LF/HF, SAP, DAP |
9 | Mohammadi et al. [34] (2019) | 33 (10) | 18–30 | 3 (rest, stress, and recovery) | Sequential TSST | Experimental, acute | Mental stress before, during and after a TSST | (↓) SpeEn (in both genders)(↑) LF/HF, SD1 (only in women) |
10 | Pereira et al. [47] (2017) | 14 (5) | 20–26 | 5 (baseline, silence, reading presentation, and counting phases of TSST) | Sequential TSST | Experimental, acute | Mental stress before and during a TSST | (↓) AVNN, SDNN, RMSSD and pNN20, LF, HF(↑) LF/HF, SD1, SD2, SampEn, α1 |
11 | Rajčáni et al. [30] (2016) | 27(17) | 18–26 | Laboratory: 6 (adapt, prepare, speech 1 and 2, relax 1 and 2)Real-life setting: 2 (stressful day and a relax day) | Laboratory setting: sequential PSSTReal-life setting: variable number of days with the entire study spanning a period of 1 year | Experimental, acuteReal-life, perceived | Laboratory setting: mental stress before, during, and after a PSSTReal-life setting: subjectively experienced stress | Laboratory setting – (↓) SDNN, RMSSD, HF and LF/HFReal-life setting: Similar as laboratory setting but not statistically significant |
12 | Schubert et al. [41] (2009) | 50 (28) | 30.3±4.7 | 2 (baseline and speech task) | Three-minute preparation time between baseline and speech | Experimental, acute | Short-term acute stress before and during a speech taskChronic stress was also measured before the experiment | Acute stress: (↓) D2(↑) SDRR, LF, and HFChronic stress: (↓) D2 |
13 | Spangler [50] (2015) | 84 (40)Undergraduate students | 19.6±1.8 | 3 (baseline, tasks and recovery for each stressor) | A three-minute adapted “vanilla” baseline was presented before each stressor | Experimental, acute | Cardiac vagal influences during a mental arithmetic PASAT, verbal fluency test COWAT, and a speech preparation task | (↓) ln(HF) in all three stressor tasks |
14 | Spellenberg et al. [46] (2020) | 24 (12) | 20–35 | 9 (before (T1, T2), during (T3–T5) and after (T6–T9) stress test | T3, T4, and T5 are sequential and other time points with a 5 min duration between each | Experimental, acute | Psychological stress before, during, and after a TSST | (↓) RR, SDNN, RMSSD, LF, HF, P1V%, P2V%, P1Vτ%, P2Vτ%(↑) LF%, P0V%, P0Vτ% |
15 | Traina et al. [33] (2011) | 13 (6) | 30–60 | 2 (baseline, stress task) | Sequential | Experimental, acute | Anxiety and worry before and during a metal arithmetic test | (↓) Hfnu(↑) LFnu, LnPLF (ms2), LF/HF |
. | Study . | N (F) . | Age, years(µ±SD or range) . | Repeated measures (time points) . | Duration between time points . | Stress type . | Context . | HRV measures significantly associated with stress and reported direction of relationship . |
---|---|---|---|---|---|---|---|---|
1 | Carnevali et al. [31] (2018) | 42 (22) | 23.7±1.1 (F)23.0±1.0 (M) | 3 | 12.8 months between T0 and T120.7 months between T1 and T2 | Real life, chronic | Depressive symptoms and rumination over a period of 3 years | (↓) RMSSD |
2 | Brugnera et al. [42] (2018) | 60 (51.7%) | 25.6±3.8 | 5 (baseline, three tasks and recovery) | 5–7 min; entire procedure was finished in an hour | Experimentally induced,acute | Mental stress before, during and after MIST, Stroop and Speech tasks | MIST: (↓) SDNN, RMSSD, TP, HFpow, LFpow, SD1, SD2, D2(↑) SI, SampEn, α1Speech: (↓) RMSSD, Hfnu, SD1, SampEn(↑) Lfnu, LF/HF, α1Stroop: (↓)Hfnu, SampEn(↑) LFpow, Lfnu, LF/HFSD2, α1 |
3 | Cervantes Blasquez et al. [28] (2009) | 10 (6) swimmers | 47 | 2 (training, competition) | 7 days | Real-life, anticipatory, acute | State anxiety before entering a swimming competition | (↓) RMSSD, HF, Hfnu, SD1(↑) LF/HF |
4 | Delaney [36] (2000) | 30 (16) | 30.9±3.9 (F)34.4±8.7 (M) | 2 (baseline, task) | Not mentioned clearly; stress situation was followed by 10 min of baseline rest period | Experimental, acute | Emotional state and muscle tension during a Stroop test and a mental arithmetic component | (↓) SDNN, RMSSD, PNN50, TP, HFpow, LFpow, Hfnu, LFnu(↑) LF/HF |
5 | Dimitriev et al. [27] (2016) | 96 (81) students | 20.5±0.1 | 2 (rest day, before an exam) | Not clearly mentioned; article states as “two different days” | Real-life, Acute | State anxiety just before an academic examination | (↓) ApEn, SampEn, SD1, SD1/SD2, CCM, LLE.(↑) α1 |
6 | Kanthak et al. [29] (2017) | 403 (66.6%) | 42.2±11.2 | 3 (Procedure, recovery, post-stress rest) | 20 min | Real life, acute | Emotional stress during a blood sampling procedure | (↓) RMSSD with emotional exhaustion component of burnout |
7 | Logan et al. [35] (2020) | 85 (100%) | 28.8±9.8 | 2 (before and after stress task) | Sequential TSST | Experimental, acute | Psychological distress and anxiety associated with TSST | No significant changes in HF-HRV in the pre and post TSST |
8 | Lucini et al. [26] (2002) | 30*university students | 22.0±1.0 | 2 (before an exam, restful holiday) | 3 months | Real life, acute | State anxiety before an academic exam | (↓) HFnu, alpha index(↑) LFnu, LF/HF, SAP, DAP |
9 | Mohammadi et al. [34] (2019) | 33 (10) | 18–30 | 3 (rest, stress, and recovery) | Sequential TSST | Experimental, acute | Mental stress before, during and after a TSST | (↓) SpeEn (in both genders)(↑) LF/HF, SD1 (only in women) |
10 | Pereira et al. [47] (2017) | 14 (5) | 20–26 | 5 (baseline, silence, reading presentation, and counting phases of TSST) | Sequential TSST | Experimental, acute | Mental stress before and during a TSST | (↓) AVNN, SDNN, RMSSD and pNN20, LF, HF(↑) LF/HF, SD1, SD2, SampEn, α1 |
11 | Rajčáni et al. [30] (2016) | 27(17) | 18–26 | Laboratory: 6 (adapt, prepare, speech 1 and 2, relax 1 and 2)Real-life setting: 2 (stressful day and a relax day) | Laboratory setting: sequential PSSTReal-life setting: variable number of days with the entire study spanning a period of 1 year | Experimental, acuteReal-life, perceived | Laboratory setting: mental stress before, during, and after a PSSTReal-life setting: subjectively experienced stress | Laboratory setting – (↓) SDNN, RMSSD, HF and LF/HFReal-life setting: Similar as laboratory setting but not statistically significant |
12 | Schubert et al. [41] (2009) | 50 (28) | 30.3±4.7 | 2 (baseline and speech task) | Three-minute preparation time between baseline and speech | Experimental, acute | Short-term acute stress before and during a speech taskChronic stress was also measured before the experiment | Acute stress: (↓) D2(↑) SDRR, LF, and HFChronic stress: (↓) D2 |
13 | Spangler [50] (2015) | 84 (40)Undergraduate students | 19.6±1.8 | 3 (baseline, tasks and recovery for each stressor) | A three-minute adapted “vanilla” baseline was presented before each stressor | Experimental, acute | Cardiac vagal influences during a mental arithmetic PASAT, verbal fluency test COWAT, and a speech preparation task | (↓) ln(HF) in all three stressor tasks |
14 | Spellenberg et al. [46] (2020) | 24 (12) | 20–35 | 9 (before (T1, T2), during (T3–T5) and after (T6–T9) stress test | T3, T4, and T5 are sequential and other time points with a 5 min duration between each | Experimental, acute | Psychological stress before, during, and after a TSST | (↓) RR, SDNN, RMSSD, LF, HF, P1V%, P2V%, P1Vτ%, P2Vτ%(↑) LF%, P0V%, P0Vτ% |
15 | Traina et al. [33] (2011) | 13 (6) | 30–60 | 2 (baseline, stress task) | Sequential | Experimental, acute | Anxiety and worry before and during a metal arithmetic test | (↓) Hfnu(↑) LFnu, LnPLF (ms2), LF/HF |
α1 and α2,DFA, detrended fluctuation analysis measures; ApEn, approximate entropy; AVNN, average value of NN (RR) intervals; D2, correlation dimension; HF, high-frequency power (0.15–0.4 Hz); LF, low-frequency power (0.04–0.15 Hz); LF/HF, frequency domain ratio between LF and HF power; LLE, largest Lyapunov exponent; N (F), number of participants; F, female; pNN20, proportion of consecutive NN intervals that differ by more than 20 ms; pNN50, proportion of consecutive NN intervals that differ by more than 50 ms; RMSSD, square root of the mean of the sum of the squares of differences between adjacent NN interval; SampEn, sample entropy; SD1, the standard deviation of instantaneous beat-to-beat RR interval variability; SD2, the standard deviation of continuous long-term RR interval variability; SDNN, standard deviation of NN intervals; SpeEn, spectral entropy; *, the number of male and female participants were not mentioned clearly; µ ± SD, mean ± standard deviation; TSST, Trier Social Stress Test; MIST, Montreal Imaging Stress Task.
Study Characteristics
Sample sizes ranged between 10 (Cervantes Blasquez 2009) and 403 (Kanthak 2017) participants, with a mean of 67 participants. All but one study included both genders, ages ranging between 18 and 60. The nature of stress, whether it was real-life situational stress or experimentally induced stress, was also elicited from the studies. Ten studies involved a controlled setting where stress was induced and quantified using HRV, both at baseline and while participants performed cognitive control-based timed tasks such Stroop task, mental arithmetic, speech preparation and delivery, Trier Social Stress Test (TSST) etc. To collect HRV, four studies used single-lead ECG signals, while Dimitriev et al. (2016) acquired 12 lead ECGs. Wearable devices (n = 7) such as Polar, Pulse, Firstbeat, Holter (n = 1) etc., were also used as RR acquisition modes. The study characteristics are summarized in Table 1.
Assessment of Stress
The most common stress assessment method used in the reviewed studies (n = 10) is the State-Trait Anxiety Inventory (STAI) which captures state anxiety [24]. Different versions of the STAI have been used within these studies, complemented by other instruments like the visual analog scale, and the perceived stress scale has also been used (Table 2). Some studies use specific symptom-focused instruments like the Maslach Burnout Inventory (MBI), Centre for Epidemiological Studies Depression Scale (CES-D), and Ruminative Response Scale (RRS), among others, to quantify psychological stress. One study (Schubert 2009) measured chronic stress using the hassles frequency subscale of the Combined Hassles and Uplifts Scale (CHUS).
Psychometric evaluation instruments used in each study, the measures derived, and the findings
Study . | Evaluation tool used . | Description . | Context of the study . | Findings . |
---|---|---|---|---|
Carnevali et al. [31] (2018) | Center for Epidemiological Studies Depression Scale (CES-D)Ruminative Response Scale (RRS) | A 20-item self-report scale designed to measure depressive symptomatology during the past weekDesigned to measure depressive rumination by how often people engage in responses to depressed mood that are self-focused, symptom focused, and focused on the possible consequences, and causes of one’s mood | Depressive symptoms and rumination over a period of 3 years | A significant time effect was observed in the CES-D score but not in the RRS score (CES-D score decreased from time 0 to time 1 and again increased at time 2) |
Brugnera et al. [42] (2018) | Stress rating questionnaire (SRQ)Task engagement questionnaire (TEQ) | A self-report 5-item questionnaire developed to assess change in stress awareness where current stress is rated on a 7-point Likert scale. An internally adapted version of the original measure was used in this studyA single item 10-point Likert scale to rate stress felt at that moment, where zero indicates the lowest level of stress and 10 being the highest level of stress | Mental stress before, during, and after MIST, Stroop and speech tasks | SRQ and TEQ scores increased from rest during all stress tasks and decreased during recovery, but no differences were found between MIST, speech and Stroop tasks |
Cervantes Blasquez et al. [28] (2009) | Competitive State Anxiety Inventory-2 | Designed to measure preperformance cognitive anxiety, somatic anxiety, and self-confidence of participants entering a competition. Comprises 27 items, with nine items in each subscale | Comparing competitive state anxiety between a baseline training condition (TC) and a competition condition (CC) | Precompetitive somatic anxiety measured by CSAI-2 showed a significant increase from the training to competitive condition |
Delaney [36] (2000) | Visual analog scales (VAS) | Self-evaluation of physical tension and emotional state at that moment measured by a simple analog scale from 1 to 10 | Emotional state and muscle tension before and during an experimental state consisting of a Stroop test and a mental arithmetic component | The experimental state was associated with a significant increase in VAS scores as compared to baseline |
Dimitriev et al. [27] (2016) | Russian version of Spielberger’s State-Trait Anxiety Inventory (STAI) subscale | Evaluates the current state of anxiety and intensity of emotional experiences by measuring subjective feelings of apprehension, tension, nervousness, and worry. Scores were classified as low (0–30), moderate (31–45), and high (≥46) | State anxiety just before an academic examination as compared to a rest day | Most participants showed an increase in SA during the exam session with the changes either from low to moderate, moderate to high, or low to high |
Kanthak et al. [29] (2017) | Maslach Burnout Inventory (MBI)The German version of patient health questionnaire PHQ-9 | Burnout was assessed with the German version of the Maslach Burnout Inventory-General Survey which consists of 16 items forming three subscales (EE, cynicism, reduced personal efficacy). The items are rated on a 7-point Likert scale. A weighted MBI total score and the three MBI sub scores were considered as continuous variables and analyzed separatelyDepression was measured with the German version (PHQ-9-D; 37) of the patient health questionnaire. The PHQ-9 consists of nine items, each scored on a 4-point rank, which quantify the frequency, over the last 2 weeks, of each of the nine diagnostic criteria for a depressive disorder of the Diagnostic and Statistical Manual of Mental Disorders | Investigating ANS modulations in an emotionally arousing situation in relation to burnout symptomatology as well as to examine potential overlap between burnout and depressivesymptomatology. A seated resting condition, a blood sampling procedure, and a recovery period were compared | EE component of burnout was significantly associated with ANS modulations during the blood sampling and the seated resting periodDepression symptoms were negatively associated with ANS modulation during rest and recovery |
Logan et al. [35] (2020) | Subjective Units of Distress Scale (SUDS)Spielberger State-Trait Anxiety Inventory (STAI) | Peak subjective distress is measured using the SUDS which is a one-item scale from 0 to 10This 20-item Likert-type scale contains response options ranging from 0 to 4 to assess the intensity of an individual’s momentary feelings. The total score is derived from the sum of the items with higher scores indicating greater anxiety | Acute psychological distress operationalized as peak subjective distress and self-reported state anxiety is compared pre and post stress tasks based on the TSST protocol | Both SUDS and STAI scores were significantly higher post-test as compared to pre-test conditions |
Lucini et al. [26] (2002) | A battery of questionnaires with self-rated scales focuses on appraisal, coping, and health | Subjective appraisal of stress, tiredness, or activation by a global scoring index (0–30), coping by a graphic questionnaire (scores 0–10) and somatic complaints by a global scoring index (scores 0–50) | University students, a model of mild real-life stress, were examined shortly before a university examination, and a second time 3 months afterward, during holiday | Significant differences in stress and symptoms scores were found between the stress day and the control day |
Mohammadi et al. [34] (2019) | Emotional Visual Analog Scale (EVAS) | Subjective appraisal of stress perception was marked by subjects as a point on a horizontal line ranging from 0 (feeling good without distress) to 10 (feeling highly and unbearable distress) | Evaluate the persistent effect of acute stress by comparing stress before, during, and after a TSST protocol | TSST protocol increased the negative mood in all participants based on the EVAS score during the test. The EVAS score returned to the baseline state 20 min after recovery |
Pereira et al. [47] (2017) | Portuguese version of the State-Trait Anxiety Inventory (STAI) | Assesses psychological stress based on 20 questions, and the scale uses a 4-point Likert-type response scale anchored at 1 (not at all) to 4 (very much) | Mental stress before and after a TSST protocol | STAI scores were significantly higher after the experiments compared to the baseline condition |
Rajčáni et al. [30] (2016) | Slovak version of STAIPerceived stress scale (PSS) | Laboratory setting: trait anxiety assessed using STAI which is a 20-item Likert-type scale contains response options ranging from 0 to 4 to assess momentary stressReal-life setting: PSS is a subjective perceived level of stress evaluated on a 10-point scale | Laboratory setting: mental stress before and after a PSSTReal-life setting: compare subjective stress between a subjectively selected stressful day and a relaxed day | Laboratory setting: significantly higher state anxiety both before and after public speech based on STAI.Real-life setting: PSS scores were significantly different between stressful and the relaxed day |
Schubert et al. [41] (2009) | The hassles frequency subscale of the Combined Hassles and Uplifts Scale (CHUS) | This 53-item subscale measures chronic stress by scoring based on perceived number of events over the past month | Examine effect of chronic stress on short-term stress reactivity before and during a laboratory-based speech task | A nonlinear HRV measure of acute short-term stress, D2, was significantly correlated with chronic stress |
Spangler [50] (2015) | Adult Temperament Questionnaire (ATQ)Block’s Ego Resiliency ScaleSelf-report emotion questionnaire (SEQ) | Attentional control, inhibitory control, and activation control scales from the short form of the ATQ were used to measure EC. Items consisted of 19 statements that were rated on a 7-point Likert scale Resiliency was measured using Block’s Ego Resiliency Scale consisting of 14 items each on 4-point Likert scale Experienced emotions were assessed with a SEQ where subjects rated 14 emotional adjectives on a 9-point Likert scale | To examine if resiliency and EC differ as self-regulatory traits in the context of stress reactivity and recovery assessed based on a mental arithmetic PASAT, verbal fluency test COWAT, and a speech preparation task | Neither EC nor resiliency moderated relations representing autonomic control of acute short-term stress reactivitySEQ based frustration significantly increased from baseline while contentment decreased after the exposure to acute short-term stressors |
Spellenberg et al. [46] (2020) | Visual analog scale (VAS) | Perceived stress quantified using a visual analog scale (VAS) score ranging from zero (feeling no stress at all) to a 100 (feeling maximally stressed) | Investigating acute stress and accompanying alterations of cardiac autonomic regulation by assessing psychological stress before, during, and after a TSST protocol | Across the procedure, changes in the VAS scores were significantly different, increasing from baseline and maximum after the arithmetic task and decreasing gradually thereafter |
Traina et al. [33] (2011) | Y-form of the State-Trait Anxiety Inventory (STAI) | Composed of two questionnaires, one to evaluate how the subject feels in the actual state (A-state) and how the subject feels habitually in everyday life (A-trait) | Anxiety and worry before and after a metal arithmetic task. A-Trait was administered before, and A-state was administered before and after the test | A-State score was not significantly different between baseline ad during the mental arithmetic task |
Study . | Evaluation tool used . | Description . | Context of the study . | Findings . |
---|---|---|---|---|
Carnevali et al. [31] (2018) | Center for Epidemiological Studies Depression Scale (CES-D)Ruminative Response Scale (RRS) | A 20-item self-report scale designed to measure depressive symptomatology during the past weekDesigned to measure depressive rumination by how often people engage in responses to depressed mood that are self-focused, symptom focused, and focused on the possible consequences, and causes of one’s mood | Depressive symptoms and rumination over a period of 3 years | A significant time effect was observed in the CES-D score but not in the RRS score (CES-D score decreased from time 0 to time 1 and again increased at time 2) |
Brugnera et al. [42] (2018) | Stress rating questionnaire (SRQ)Task engagement questionnaire (TEQ) | A self-report 5-item questionnaire developed to assess change in stress awareness where current stress is rated on a 7-point Likert scale. An internally adapted version of the original measure was used in this studyA single item 10-point Likert scale to rate stress felt at that moment, where zero indicates the lowest level of stress and 10 being the highest level of stress | Mental stress before, during, and after MIST, Stroop and speech tasks | SRQ and TEQ scores increased from rest during all stress tasks and decreased during recovery, but no differences were found between MIST, speech and Stroop tasks |
Cervantes Blasquez et al. [28] (2009) | Competitive State Anxiety Inventory-2 | Designed to measure preperformance cognitive anxiety, somatic anxiety, and self-confidence of participants entering a competition. Comprises 27 items, with nine items in each subscale | Comparing competitive state anxiety between a baseline training condition (TC) and a competition condition (CC) | Precompetitive somatic anxiety measured by CSAI-2 showed a significant increase from the training to competitive condition |
Delaney [36] (2000) | Visual analog scales (VAS) | Self-evaluation of physical tension and emotional state at that moment measured by a simple analog scale from 1 to 10 | Emotional state and muscle tension before and during an experimental state consisting of a Stroop test and a mental arithmetic component | The experimental state was associated with a significant increase in VAS scores as compared to baseline |
Dimitriev et al. [27] (2016) | Russian version of Spielberger’s State-Trait Anxiety Inventory (STAI) subscale | Evaluates the current state of anxiety and intensity of emotional experiences by measuring subjective feelings of apprehension, tension, nervousness, and worry. Scores were classified as low (0–30), moderate (31–45), and high (≥46) | State anxiety just before an academic examination as compared to a rest day | Most participants showed an increase in SA during the exam session with the changes either from low to moderate, moderate to high, or low to high |
Kanthak et al. [29] (2017) | Maslach Burnout Inventory (MBI)The German version of patient health questionnaire PHQ-9 | Burnout was assessed with the German version of the Maslach Burnout Inventory-General Survey which consists of 16 items forming three subscales (EE, cynicism, reduced personal efficacy). The items are rated on a 7-point Likert scale. A weighted MBI total score and the three MBI sub scores were considered as continuous variables and analyzed separatelyDepression was measured with the German version (PHQ-9-D; 37) of the patient health questionnaire. The PHQ-9 consists of nine items, each scored on a 4-point rank, which quantify the frequency, over the last 2 weeks, of each of the nine diagnostic criteria for a depressive disorder of the Diagnostic and Statistical Manual of Mental Disorders | Investigating ANS modulations in an emotionally arousing situation in relation to burnout symptomatology as well as to examine potential overlap between burnout and depressivesymptomatology. A seated resting condition, a blood sampling procedure, and a recovery period were compared | EE component of burnout was significantly associated with ANS modulations during the blood sampling and the seated resting periodDepression symptoms were negatively associated with ANS modulation during rest and recovery |
Logan et al. [35] (2020) | Subjective Units of Distress Scale (SUDS)Spielberger State-Trait Anxiety Inventory (STAI) | Peak subjective distress is measured using the SUDS which is a one-item scale from 0 to 10This 20-item Likert-type scale contains response options ranging from 0 to 4 to assess the intensity of an individual’s momentary feelings. The total score is derived from the sum of the items with higher scores indicating greater anxiety | Acute psychological distress operationalized as peak subjective distress and self-reported state anxiety is compared pre and post stress tasks based on the TSST protocol | Both SUDS and STAI scores were significantly higher post-test as compared to pre-test conditions |
Lucini et al. [26] (2002) | A battery of questionnaires with self-rated scales focuses on appraisal, coping, and health | Subjective appraisal of stress, tiredness, or activation by a global scoring index (0–30), coping by a graphic questionnaire (scores 0–10) and somatic complaints by a global scoring index (scores 0–50) | University students, a model of mild real-life stress, were examined shortly before a university examination, and a second time 3 months afterward, during holiday | Significant differences in stress and symptoms scores were found between the stress day and the control day |
Mohammadi et al. [34] (2019) | Emotional Visual Analog Scale (EVAS) | Subjective appraisal of stress perception was marked by subjects as a point on a horizontal line ranging from 0 (feeling good without distress) to 10 (feeling highly and unbearable distress) | Evaluate the persistent effect of acute stress by comparing stress before, during, and after a TSST protocol | TSST protocol increased the negative mood in all participants based on the EVAS score during the test. The EVAS score returned to the baseline state 20 min after recovery |
Pereira et al. [47] (2017) | Portuguese version of the State-Trait Anxiety Inventory (STAI) | Assesses psychological stress based on 20 questions, and the scale uses a 4-point Likert-type response scale anchored at 1 (not at all) to 4 (very much) | Mental stress before and after a TSST protocol | STAI scores were significantly higher after the experiments compared to the baseline condition |
Rajčáni et al. [30] (2016) | Slovak version of STAIPerceived stress scale (PSS) | Laboratory setting: trait anxiety assessed using STAI which is a 20-item Likert-type scale contains response options ranging from 0 to 4 to assess momentary stressReal-life setting: PSS is a subjective perceived level of stress evaluated on a 10-point scale | Laboratory setting: mental stress before and after a PSSTReal-life setting: compare subjective stress between a subjectively selected stressful day and a relaxed day | Laboratory setting: significantly higher state anxiety both before and after public speech based on STAI.Real-life setting: PSS scores were significantly different between stressful and the relaxed day |
Schubert et al. [41] (2009) | The hassles frequency subscale of the Combined Hassles and Uplifts Scale (CHUS) | This 53-item subscale measures chronic stress by scoring based on perceived number of events over the past month | Examine effect of chronic stress on short-term stress reactivity before and during a laboratory-based speech task | A nonlinear HRV measure of acute short-term stress, D2, was significantly correlated with chronic stress |
Spangler [50] (2015) | Adult Temperament Questionnaire (ATQ)Block’s Ego Resiliency ScaleSelf-report emotion questionnaire (SEQ) | Attentional control, inhibitory control, and activation control scales from the short form of the ATQ were used to measure EC. Items consisted of 19 statements that were rated on a 7-point Likert scale Resiliency was measured using Block’s Ego Resiliency Scale consisting of 14 items each on 4-point Likert scale Experienced emotions were assessed with a SEQ where subjects rated 14 emotional adjectives on a 9-point Likert scale | To examine if resiliency and EC differ as self-regulatory traits in the context of stress reactivity and recovery assessed based on a mental arithmetic PASAT, verbal fluency test COWAT, and a speech preparation task | Neither EC nor resiliency moderated relations representing autonomic control of acute short-term stress reactivitySEQ based frustration significantly increased from baseline while contentment decreased after the exposure to acute short-term stressors |
Spellenberg et al. [46] (2020) | Visual analog scale (VAS) | Perceived stress quantified using a visual analog scale (VAS) score ranging from zero (feeling no stress at all) to a 100 (feeling maximally stressed) | Investigating acute stress and accompanying alterations of cardiac autonomic regulation by assessing psychological stress before, during, and after a TSST protocol | Across the procedure, changes in the VAS scores were significantly different, increasing from baseline and maximum after the arithmetic task and decreasing gradually thereafter |
Traina et al. [33] (2011) | Y-form of the State-Trait Anxiety Inventory (STAI) | Composed of two questionnaires, one to evaluate how the subject feels in the actual state (A-state) and how the subject feels habitually in everyday life (A-trait) | Anxiety and worry before and after a metal arithmetic task. A-Trait was administered before, and A-state was administered before and after the test | A-State score was not significantly different between baseline ad during the mental arithmetic task |
MIST, Montreal Imaging Stress Task; EC, effortful control.
Measures of HRV
In the reviewed studies, a wide range of linear and nonlinear HRV metrics recommended by the Task Force of the European Society of Cardiology [25] were used. The number of metrics used in individual studies varied from one parameter to 11 different parameters. The most used time-domain parameter was the RMSSD (n = 10), and the second most common time-domain parameter was the SDNN (n = 8). The most used frequency-domain parameters were the LF/HF ratio (n = 7) and HF power (n = 6). Overall, nonlinear metrics have been applied less frequently in stress assessment. However, within included studies based on nonlinear analysis, SD1 (n = 5), SampEn, and α1 (n = 3) have been shown to be associated with stress. The online supplementary material (online suppl. Table S1) shows a detailed summary table showing all HRV metrics analyzed in each study and highlighting the metrics reported to be significantly associated with stress in the context of each study.
Nature of Stress
Anticipatory, acute stress in real-life situations such as before sitting for a university exam [26, 27] and before swimming [28] was investigated to evaluate stress. Reviewed studies also addressed other real-life stressors, such as chronic stress due to depression and rumination, emotional exhaustion, and burnout. Experimental acute stressors were the most common form employed in nine reviewed studies (Table 1). Acute task-based stress was induced in these studies by making participants perform mental arithmetic, speech, verbal fluency, imaging, and Stroop test tasks.
Quantification and Coding of the Relationship between HRV and Stress
In the included studies, a range of linear time-domain, frequency-domain, and nonlinear metrics derived from the RR interval time series have been used to quantify HRV. The definition of these metrics, the studies in which they have been included, and the physiological changes that these metrics have shown to elicit are summarized in Table 3. To integrate the findings on HRV changes associated with stress, we coded baseline-to-task changes for every HRV measure as either significantly increasing or directly related (↑), decreasing or inversely related (↓) as shown in Table 1. Similarly, different psychometric evaluation instruments were used, and measures were derived in each study (Table 2).
HRV metrics analyzed in the included studies
HRV measure . | Description . | Physiological significance . |
---|---|---|
SDNN | The standard deviation of NN intervals | Represents the global autonomic function of the heart, reflecting joint sympathetic and parasympathetic modulation of heart rate |
RMSSD | The square root of the mean of the sum of the squares of differences between adjacent NN interval | Represents predominantly parasympathetic modulation of heart rate |
pNN50 | The proportion of consecutive NN intervals that differ by more than 50 ms | Represents short-term changes in heart rate variability, reflecting the rapid alterations characteristic of the parasympathetic nervous system |
pNN20 | The proportion of consecutive NN intervals that differ by more than 20 ms | Represents short-term changes in heart rate variability, reflecting the rapid alterations characteristic of the parasympathetic nervous system |
HF | High-frequency power (0.15–0.4 Hz) | Reflects respiratory sinus arrhythmia, mediated by changes in vagal tone. A decrease in HF is associated with vagal withdrawal |
HFn | High-frequency power (0.15–0.4 Hz), normalized to HF and LF power | Often interpreted as a marker of sympathovagal balance, where an increase is viewed as vagal predominance, note that HFn+LFn = 1 |
LF | Low-frequency power (0.04–0.15 Hz) | Reflects a mix of the sympathetic and vagal influences on heart rate variability |
LFn | Low-frequency power (0.04–0.15 Hz), normalized to HF and LF power | Often interpreted as a marker of sympathovagal balance, where an increase is viewed as sympathetic predominance, note that HFn+LFn = 1 |
VLF | Very low-frequency power (0.0033–0.04 Hz) | Represent long-term regulation mechanisms, thermoregulation, and hormonal mechanisms |
LF/HF | Frequency-domain ratio between LF and HF power | Conceptualized as a global index of sympathovagal balance, but controversial |
SD1, SD2 | SD1 and SD2 can be derived from Poincaré plots that plot RR intervals against the succeeding RR intervals or, mathematically RMSSD and SDNN. | SD1 is strongly influenced by parasympathetic activity, and SD2 reflects overall variability |
D2 | The correlation dimension (D2) gives information on the complexity of HRV. | The interpretation of D2 is challenging |
α1 and α2 | Detrended fluctuation analysis (DFA) measures short- and long-range correlations | α1 is modulated by sympathetic activity [51] |
APEn, SampEn | ApEn and SampEn measure the randomness (i.e., unpredictability) and irregularity of the NN interval series | Lower values suggest a reduced HRV complexity |
SpeEn | SpeEn is a measure of the randomness (i.e., unpredictability) and irregularity of the power density spectrum | Lower values suggest a reduced HRV complexity |
LLE | The largest Lyapunov exponent quantifies the trajectory of a nonlinear dynamic system. LLE >0 indicates a chaotic system | HRV of sinus rhythm has characteristics of chaos-like determinism, which reflects in a positive Lyapunov exponent, and changes in LLE are believed to be associated with parasympathetic activation |
Guzik’s index of asymmetry (GI) | An index used to analyze the asymmetry of HRV along the line of identity in the Poincare plot | Asymmetry of heart rate accelerations and decelerations is a normal physiological phenomenon related to respiratory sinus arrhythmia [52] |
Complex correlation measure (CCM) | An extension of the Poincaré plot analysis, measures beat-to-beat dynamics | A decrease in CCM indicates increased regularity or decreased variability and is associated with parasympathetic activity |
Symbolic dynamics 0V%, 1V%, 2LV% | Symbolic dynamics allows studying the nonlinear dynamics of a system using coarse-grained symbol patterns (“words”) | 0V% and 2LV% have been associated with sympathetic activity [53] |
Baevsky’s Stress Index (SI) | An index used to quantify interbeat interval shapes and distribution during sympathetic activation and is based on a statistical analysis of histogram of RR interval distribution | Reflects stress-related sympathetic activation levels |
HRV measure . | Description . | Physiological significance . |
---|---|---|
SDNN | The standard deviation of NN intervals | Represents the global autonomic function of the heart, reflecting joint sympathetic and parasympathetic modulation of heart rate |
RMSSD | The square root of the mean of the sum of the squares of differences between adjacent NN interval | Represents predominantly parasympathetic modulation of heart rate |
pNN50 | The proportion of consecutive NN intervals that differ by more than 50 ms | Represents short-term changes in heart rate variability, reflecting the rapid alterations characteristic of the parasympathetic nervous system |
pNN20 | The proportion of consecutive NN intervals that differ by more than 20 ms | Represents short-term changes in heart rate variability, reflecting the rapid alterations characteristic of the parasympathetic nervous system |
HF | High-frequency power (0.15–0.4 Hz) | Reflects respiratory sinus arrhythmia, mediated by changes in vagal tone. A decrease in HF is associated with vagal withdrawal |
HFn | High-frequency power (0.15–0.4 Hz), normalized to HF and LF power | Often interpreted as a marker of sympathovagal balance, where an increase is viewed as vagal predominance, note that HFn+LFn = 1 |
LF | Low-frequency power (0.04–0.15 Hz) | Reflects a mix of the sympathetic and vagal influences on heart rate variability |
LFn | Low-frequency power (0.04–0.15 Hz), normalized to HF and LF power | Often interpreted as a marker of sympathovagal balance, where an increase is viewed as sympathetic predominance, note that HFn+LFn = 1 |
VLF | Very low-frequency power (0.0033–0.04 Hz) | Represent long-term regulation mechanisms, thermoregulation, and hormonal mechanisms |
LF/HF | Frequency-domain ratio between LF and HF power | Conceptualized as a global index of sympathovagal balance, but controversial |
SD1, SD2 | SD1 and SD2 can be derived from Poincaré plots that plot RR intervals against the succeeding RR intervals or, mathematically RMSSD and SDNN. | SD1 is strongly influenced by parasympathetic activity, and SD2 reflects overall variability |
D2 | The correlation dimension (D2) gives information on the complexity of HRV. | The interpretation of D2 is challenging |
α1 and α2 | Detrended fluctuation analysis (DFA) measures short- and long-range correlations | α1 is modulated by sympathetic activity [51] |
APEn, SampEn | ApEn and SampEn measure the randomness (i.e., unpredictability) and irregularity of the NN interval series | Lower values suggest a reduced HRV complexity |
SpeEn | SpeEn is a measure of the randomness (i.e., unpredictability) and irregularity of the power density spectrum | Lower values suggest a reduced HRV complexity |
LLE | The largest Lyapunov exponent quantifies the trajectory of a nonlinear dynamic system. LLE >0 indicates a chaotic system | HRV of sinus rhythm has characteristics of chaos-like determinism, which reflects in a positive Lyapunov exponent, and changes in LLE are believed to be associated with parasympathetic activation |
Guzik’s index of asymmetry (GI) | An index used to analyze the asymmetry of HRV along the line of identity in the Poincare plot | Asymmetry of heart rate accelerations and decelerations is a normal physiological phenomenon related to respiratory sinus arrhythmia [52] |
Complex correlation measure (CCM) | An extension of the Poincaré plot analysis, measures beat-to-beat dynamics | A decrease in CCM indicates increased regularity or decreased variability and is associated with parasympathetic activity |
Symbolic dynamics 0V%, 1V%, 2LV% | Symbolic dynamics allows studying the nonlinear dynamics of a system using coarse-grained symbol patterns (“words”) | 0V% and 2LV% have been associated with sympathetic activity [53] |
Baevsky’s Stress Index (SI) | An index used to quantify interbeat interval shapes and distribution during sympathetic activation and is based on a statistical analysis of histogram of RR interval distribution | Reflects stress-related sympathetic activation levels |
Discussion
This study was conducted to review the available literature for studies investigating various HRV metrics in healthy adult populations under different psychological stress states, including a baseline or a restful state. The review thus presents a synthesized summary derived from the included studies on different categories of HRV metrics, associated stress types, stress elicitation methods used, and the psychometric validation instruments used to evaluate stress.
Real-Life Psychological Stressors
Common real-life stressors of moderate intensity, such as preparation for a major university examination, may alter basal hemodynamic indices of resting autonomic cardiac control. Such stressors were associated with frequency domain and nonlinear HRV metrics. Lucini et al. [26] evaluated psychological involvement in students with exam-related anticipatory stress using a battery of questionnaires providing self-rated scales that focus on the appraisal of stress, coping, and health. Anticipatory stress was characterized by significantly higher values of LFn. In a similar study, analysis of state anxiety and stress in a cohort of students just before a university exam and their nonlinear HRV revealed that anxiety is associated with alterations in the complexity of HRV captured as an increase in the short-term fractal exponent α1[27]. Competitive sports, where adaptation to training loads and anticipation to compete causes anxiety, are another context of real-life anticipatory stress. Precompetitive anxiety before a swimming competition was studied [28] using HRV metrics and a Competitive State Anxiety Inventory-2 (CSAI-2). High precompetitive anxiety elicited a shift toward sympathetic predominance due to parasympathetic withdrawal. This study also concludes that RMSSD seems to be the most valid indicator of emotional state in pre-competitive situations among time-domain parameters.
In addition to anticipatory stress, other aspects, such as chronic fatigue, emotional exhaustion, and burnout, were also involved in the reviewed studies. During a blood-sampling procedure, emotional stress was examined and correlated with vagally mediated HRV in a large population-based sample of the Dresden Burnout Study [29]. The emotional exhaustion component of burnout, assessed using a Maslach Burnout Inventory (MBI), was associated with reduced vagal cardiac control during an emotionally arousing situation, reflecting an individual’s capacity to respond flexibly and adapt to changing environmental needs. Rajcani et al. 2016, investigated psychophysiological changes in the stress reaction, both in an experimental setting using PSST based on a simulated speech task and in a real-life setting by comparing a stressful and a relaxed day. Significantly diminished SDNN, LF/HF, RMSSD, and high-frequency power, all pointing to reduced HRV during a stressful state, have been observed in both settings. However, the effects in naturalistic settings were weaker but analogous to results based on experimental protocol [30]. The reduced HRV correlated with significant anxiety and perceived stress differences based on STAI and PSS scores. Observations in real-life settings have also been used by Carnevali et al. to evaluate the interplay between HRV and depressive and rumination symptoms causing chronic stress in young, healthy adults. In this study, RMSSD was negatively correlated with both rumination and depressive symptoms at each time point, implying that autonomic dysfunction, predominately low vagal tone, characterizes individuals with higher rumination traits and is prospectively linked to the generation of depressive symptoms in a non-clinical setting [31].
Experimental Acute Stressors
Experimental acute task-based stress induced by performing mental arithmetic, speech task, verbal fluency components of the TSST, Montreal Imaging Stress Task, Stroop test etc., have been used to investigate HRV-based physiological changes in nine of the reviewed studies. TSST is a standardized laboratory social stressor that induces robust and reliable increases in psychological, physiological, and neuroendocrine measures of stress. It is a helpful alternative to physical stressors and reproduces the more naturalistic psychological stress of performance in the presence of an evaluative audience [32]. In this review, Torino et al. employed only the mental arithmetic phase of TSST to study the degree of stress-associated changes in sympathetic activity and vagal tone using LF and HF components of HRV [33]. Mohammadi et al. observed stress before, during, and after TSST and noticed that HRV metrics change during stress tasks, but the changes persisted during recovery, i.e., after removing stress. This persistence of reduced HR complexity after stress is speculated to represent lower adaptability and a functional restriction of the participating cardiovascular elements. Gender differences have also been shown in this study with some significant correlations, such as an increase in LF/HF, and SD1 observed only in women [34].
Logan et al. [35] used only a single HRV metric, the HF-HRV to study autonomic function where TSST-induced stress. To quantify the effect of TSST, baseline measures of HF-HRV and state anxiety using STAI were compared with measurements obtained after the test, and no significant effects were observed. It must be noted that in this analysis, HRV before and after were only considered and HRV during stress was not included in the study. In a similar approach toward using HF-HRV, Spangler et al. used HF-HRV to assess cardiac vagal influences before, during, and after stress tasks. They demonstrated an inverse relationship between stress during tasks and HF power. The commonality seen in the analysis and findings of all these laboratory stressor-based studies is that induced stress causes the ANS to shift toward sympathetic predominance because of parasympathetic withdrawal, which provokes a characteristic defense-arousal reaction in the physiological system. A similar result was validated by Delaney et al [36]. In addition to demonstrating stress responses in a task-performing group, an intervention effect was also shown based on comparisons with a control group.
The Interplay between Chronic and Acute Stress
Schubert et al. performed an interesting study that demonstrated the influence of baseline chronic stress on the association between acute stress and HRV metrics. Short-term stressor reactivity was assessed with a speech task and various HRV metrics. In addition to other conventional measures, D2, a nonlinear HR complexity measure, significantly decreased with acute stress. The higher the HR D2, the more degrees of freedom of the cardiac pacemaker and, therefore, the greater range of possible adaptive responses [37]. This HR complexity measure was significantly inversely correlated with baseline chronic stress levels and further reduced during acute stress, suggesting that long-term chronic stress levels mediate or may be influenced by HRV responses to short-term stressors. Studies have also demonstrated the role of laboratory-based acute HRV metrics in predicting autonomic modulation of ecological emotional stress [38].
The literature shows that people with high trait anxiety tend to exhibit increased HR and diminished HRV [39, 40]. Two studies in this review have investigated the interplay of these trait features in HRV stress associations. Effortful control, a trait measure of adaptability to stress, was used by Spangler et al. (2009) to investigate the extent to which such trait-based self-regulatory constructs influence autonomic control of cardiac responses to acute stress. Resting cardiac vagal control was significantly correlated with effortful control, further influencing the vagal recovery after the verbal fluency stress task. The study by Medica-Torino et al. [33] had a similar approach in testing the influence of chronic stress, assessed using trait anxiety measure at baseline, on associations between state anxiety-based acute stress and HRV during experimental stressors. Significant associations were found between the level of state anxiety and worry associated with stress tasks and the degree of sympathetic activation and vagal tone reduction; however, no significant moderating role of trait anxiety was observed.
Due to a reduction in vagal tone, acute stress is generally associated with a typical increase in heart rate, which has been observed and reported in five reviewed studies [28, 35, 36, 41, 42]. Elevated blood pressure due to an increase in sympathetic drive is another physiological phenomenon observed during acute stress conditions which were reported in one of the reviewed studies [26].
Limitations and Future Perspectives
As discussed above, variations in HRV due to acute stress, along with the underlying non-clinical yet persistent trait anxiety, depression, burnout, and chronic stressors to which individuals might adapt differently, might make the interpretations of observed HRV changes more complex. There is an enormous scope of well-designed studies that untangle the various layers of these stress responses and mediating factors which is a recommendation for future researchers.
It has been argued that acute psychological stress acts on cardiac autonomic regulation in a way that may lead to non-stationarities in the interbeat interval series, making linear methods, primarily the frequency domain-based metrics, unreliable compared to nonlinear approaches [43]. Despite the robustness and reliability of nonlinear metrics in revealing useful additional information about HRV characteristics in different applications and patient groups [44], nonlinear methods are not very popular in HRV studies investigating psychological stress. A robust nonlinear approach termed binary symbolic dynamics analysis [45], which considers signal non-stationarity, has been used by one study to detect acute stress-associated changes in HRV and is thought to be more informative than linear methods [46]. It is recommended in this review that more studies be designed to compare and prove the usefulness of several nonlinear metrics and their analysis windows, which inherently can overcome the limitations of non-stationarities.
In addition to non-stationarity, fine-grained stress assessments are precluded by longer analysis windows, especially in real-time unobtrusive measurement using sensors, wearables etc. Short-term analysis (≤5 min) thus would be an advantage for stress studies due to the rapid physiological response time. Discrimination of short-term stress conditions using short-term HRV metrics was explored by Pereira et al. [47]. They used the TSST protocol and quantified HRV at five time points (baseline, silence, reading, presentation, and counting phases), allowing five short-term time windows for evaluation. In this study, several linear and nonlinear metrics demonstrated reduced HRV during stress, with much higher reduction strength in time-domain metrics. In particular, AVNN, RMSSD, and SDNN measures using windows as short as 50 s have been shown as the most discriminating HRV metrics distinguishing stress and non-stress states [47]. More studies demonstrating the discriminating power of such short-term HRV in real-life stress detection are recommended for future research on wearable smart device-based health trackers. Brugnera et al. performed a similar protocol in experimental stress but on five-minute analysis windows. Their results indicate that stress led to a sympathetic shift in sympathovagal balance and reduced the complexity of the cardiac signal. Montreal Imaging Stress Task induced the strongest cardiovascular response compared to the other speech and Stroop tasks and was associated with a specific profile of cardiovascular activity, characterized by sympathovagal co-inhibition, i.e., decreases in both HF and LF power [42].
In this review, spectral metrics such as HF power, LF power, and LF/HF have been claimed to be stress indicators in studies (n = 3) involving a speech task. Spectral analysis-based HRV changes, especially during the speech-based task, are argued in the literature to be interpreted with caution, as the respiratory changes produced by speech markedly alter variability and the spectral component of HRV without necessarily involving respective changes in autonomic activation [48, 49]. Respiratory patterns observed during speech tasks have been shown far from being sinusoidal, highly erratic, and associated with markedly broadband characteristics in the HR power spectrum, with considerable power present in both the LF and HF bands. In this review, we observed that the analysis was controlled for respiratory frequency in only one study [41, 42]. In contrast, two studies acknowledged respiratory influences on HRV during speech tasks [47] and mentioned not controlling for breathing and non-stationarities as their study limitation [46]. Future studies are recommended to rigorously account for and control analysis for the influence of respiratory patterns in speech-based stress tasks to improve the validity of findings. As with many other reviews, this paper included only published journal articles, and a possible publication bias may have affected our findings. Also, the length of recordings obtained and analyzed varied between the included papers, and this might have influenced our findings and interpretations. However, due to the heterogeneity, we did not further group studies based on these data and compare the findings.
Conclusion
In this review, we focused on identifying specific HRV metrics and their association with change in psychological stress under different stress states in studies using repeated sampling protocol in healthy populations. It was believed that these associations would allow us to make inferences about the metrics most suited to act as physiological markers to monitor and track the impacts of psychological stress in different stress contexts, especially real-life stressors. The review does highlight the impact of psychological stress on sympathovagal balance which is closely reflected by HRV metrics, thus evidencing that HRV can be used as an informative marker of the physiological effects of psychological stressors in healthy adult populations. The traditional time-domain RMSSD is still the most explored and reported measure. Most of the reviewed studies were performed under laboratory conditions instead of natural working life settings. Few novel studies have considered both scenarios, explored nonlinear novel metrics, and the influence of underlying chronic stressors and trait characteristics on the effects observed under acute stress conditions. However, of the few, little has been validated with large samples. In addition to short-term laboratory measurements based on conventional metrics, acute real-life stress contexts, novel HRV metrics, long-term HR monitoring, repeated sampling, and accounting for underlying chronic stressors are imperative for validating HRV metrics. Such validated metrics can accurately assess and predict stress and recovery patterns, and such validation approaches warrant future research and large-scale study designs. In addition, HRV measurements in real-life ambulatory settings often involve known confounders such as physical exertion and other unidentified confounding factors that can rarely be controlled entirely but must be carefully considered in interpreting study results. These effects can be reduced to some extent by cross-validation of findings using consistent, validated subjective screening instruments with objective HRV quantification methods in large-scale study settings. In addition, this review identified the diversity of instruments used in studies assessing psychological stress and the heterogeneity in methodological approaches of quantifying HRV, which makes the interpretation of measurements inconclusive in the context of this review. In summary, utilization of standard stress theories/models, uniform validated stress indicators, and standardized methods would improve the comparability of results in future studies. More unified HRV assessment and analysis methods utilizing both conventional and contemporary metrics and longitudinal real-life study settings are needed.
Statement of Ethics
An ethics statement is not applicable because this study is based exclusively on the published literature.
Conflict of Interest Statement
The authors have no conflicts of interest to declare.
Funding Sources
There were no funding sources.
Author Contributions
Sarah Immanuel: conception and design of the work, search criteria formulation, data acquisition, and analysis and interpretation of data for the work. Meseret Teferra: data extraction and screening and drafting of the work. Niranjan Bidargaddi: study design and drafting and critical revision of the work. Mathias Baumert: drafting, providing context, data synthesis, and critical revision for intellectual content.
Data Availability Statement
All data generated or analyzed during this study are included in this article and its supplementary materials. Further inquiries can be directed to the corresponding author.