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.

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.

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).

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.

Fig. 1.

Flow chart of included studies.

Fig. 1.

Flow chart of included studies.

Close modal
Table 1.

Study characteristics of included articles

StudyN (F)Age, years(µ±SD or range)Repeated measures (time points)Duration between time pointsStress typeContextHRV measures significantly associated with stress and reported direction of relationship
Carnevali et al. [31] (2018) 42 (22) 23.7±1.1 (F)23.0±1.0 (M) 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 
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 
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 
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 
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 
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 
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 
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 
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 
StudyN (F)Age, years(µ±SD or range)Repeated measures (time points)Duration between time pointsStress typeContextHRV measures significantly associated with stress and reported direction of relationship
Carnevali et al. [31] (2018) 42 (22) 23.7±1.1 (F)23.0±1.0 (M) 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 
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 
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 
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 
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 
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 
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 
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 
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).

Table 2.

Psychometric evaluation instruments used in each study, the measures derived, and the findings

StudyEvaluation tool usedDescriptionContext of the studyFindings
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 
StudyEvaluation tool usedDescriptionContext of the studyFindings
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).

Table 3.

HRV metrics analyzed in the included studies

HRV measureDescriptionPhysiological 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 measureDescriptionPhysiological 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 

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.

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.

An ethics statement is not applicable because this study is based exclusively on the published literature.

The authors have no conflicts of interest to declare.

There were no funding sources.

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.

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.

1.
Malik
M
,
Bigger
JT
,
Camm
AJ
,
Kleiger
RE
,
Malliani
A
,
Moss
AJ
.
Heart rate variability: standards of measurement, physiological interpretation, and clinical use
.
Eur Heart J
.
1996
;
17
(
3
):
354
81
.
2.
McCraty
R
,
Shaffer
F
.
Heart rate variability: new perspectives on physiological mechanisms, assessment of self-regulatory capacity, and health risk
.
Glob Adv Health Med
.
2015
;
4
(
1
):
46
61
.
3.
Lennartsson
AK
,
Jonsdottir
I
,
Sjors
A
.
Low heart rate variability in patients with clinical burnout
.
Int J Psychophysiol
.
2016 Dec
110
171
8
.
4.
Gonzalez
K
,
Sasangohar
F
,
Mehta
RK
,
Lawley
M
,
Erraguntla
M
.
Measuring fatigue through Heart Rate Variability and activity recognition: a scoping literature review of machine learning techniques
.
Proceedings of the human factors and ergonomics society annual meeting
Los Angeles, CA
SAGE Publications Sage CA
2017
. p.
1748
52
.
5.
Gevirtz
R
.
The promise of heart rate variability biofeedback: evidence-based applications
.
Biofeedback (Online)
.
2013
;
41
(
3
):
110
20
.
6.
Moss
D
,
Shaffer
F
.
The application of heart rate variability biofeedback to medical and mental health disorders
.
Biofeedback
.
2017
;
45
(
1
):
2
8
.
7.
Hughson
R
,
Yamamoto
Y
,
McCullough
R
,
Sutton
J
,
Reeves
J
.
Sympathetic and parasympathetic indicators of heart rate control at altitude studied by spectral analysis
.
J Appl Physiol
.
1994
;
77
(
6
):
2537
42
.
8.
Chiu
HW
,
Kao
T
.
A mathematical model for autonomic control of heart rate variation
.
IEEE Eng Med Biol Mag
.
2001 Mar–Apr
20
2
69
76
.
9.
Guzzetti
S
,
Borroni
E
,
Garbelli
PE
,
Ceriani
E
,
Della Bella
P
,
Montano
N
.
Symbolic dynamics of heart rate variability: a probe to investigate cardiac autonomic modulation
.
Circulation
.
2005 Jul 26
112
4
465
70
.
10.
Stauss
HM
.
Heart rate variability
.
Am J Physiol Regul Integr Comp Physiol
.
2003 Nov
285
5
R927
31
.
11.
Georgiou
K
,
Larentzakis
AV
,
Khamis
NN
,
Alsuhaibani
GI
,
Alaska
YA
,
Giallafos
EJ
.
Can wearable devices accurately measure heart rate variability? A systematic review
.
Folia Med
.
2018
;
60
(
1
):
7
20
.
12.
Cohen
S
,
Kessler
RC
,
Gordon
LU
.
Strategies for measuring stress in studies of psychiatric and physical disorders
.
Measuring Stress A Guide Health Soc Sci
.
1995
;
28
:
3
26
.
13.
Muscatell
KA
,
Eisenberger
NI
.
A social neuroscience perspective on stress and health
.
Soc Personal Psychol Compass
.
2012
;
6
(
12
):
890
904
.
14.
Chu
B
,
Marwaha
K
,
Sanvictores
T
,
Ayers
D
Physiology, stress reaction
StatPearls Publishing
2021
. StatPearls [Internet].
15.
Jarvelin-Pasanen
S
,
Sinikallio
S
,
Tarvainen
MP
.
Heart rate variability and occupational stress-systematic review
.
Ind Health
.
2018 Nov 21
56
6
500
11
.
16.
Reijmerink
I
,
van der Laan
M
,
Cnossen
F
.
Heart rate variability as a measure of mental stress in surgery: a systematic review
.
Int Arch Occup Environ Health
.
2020
1
17
.
17.
Thielmann
B
,
Pohl
R
,
Böckelmann
I
.
Heart rate variability as a strain indicator for psychological stress for emergency physicians during work and alert intervention: a systematic review
.
J Occup Med Toxicol
.
2021
;
16
(
1
):
24
9
.
18.
Kim
HG
,
Cheon
EJ
,
Bai
DS
,
Lee
YH
,
Koo
BH
.
Stress and heart rate variability: a meta-analysis and review of the literature
.
Psychiatry Investig
.
2018 Mar
15
3
235
45
.
19.
Schiweck
C
,
Piette
D
,
Berckmans
D
,
Claes
S
,
Vrieze
E
.
Heart rate and high frequency heart rate variability during stress as biomarker for clinical depression. A systematic review
.
Psychol Med
.
2019
;
49
(
2
):
200
11
.
20.
Thayer
JF
,
Åhs
F
,
Fredrikson
M
,
Sollers
JJ
III
,
Wager
TD
.
A meta-analysis of heart rate variability and neuroimaging studies: implications for heart rate variability as a marker of stress and health
.
Neurosci Biobehav Rev
.
2012
;
36
(
2
):
747
56
.
21.
Kamarck
TW
,
Lovallo
WR
.
Cardiovascular reactivity to psychological challenge: conceptual and measurement considerations
.
Psychosom Med
.
2003
;
65
(
1
):
9
21
.
22.
Ebner
K
,
Singewald
N
.
Individual differences in stress susceptibility and stress inhibitory mechanisms
.
Curr Opin Behav Sci
.
2017
;
14
:
54
64
.
23.
Pattyn
N
,
Migeotte
PF
,
Neyt
X
,
den Nest
A
,
Cluydts
R
.
Comparing real-life and laboratory-induced stress reactivity on cardio-respiratory parameters: differentiation of a tonic and a phasic component
.
Physiol Behav
.
2010 Sep 1
101
2
218
23
.
24.
Spielberger
CD
The Corsini encyclopedia of psychology State‐Trait anxiety inventory
.
2010
. p.
1
.
25.
Electrophysiology TFotESoCtNASoP
.
Heart rate variability: standards of measurement, physiological interpretation, and clinical use
.
Circulation
.
1996
;
93
(
5
):
1043
65
.
26.
Lucini
D
,
Norbiato
G
,
Clerici
M
,
Pagani
M
.
Hemodynamic and autonomic adjustments to real life stress conditions in humans
.
Hypertension
.
2002 Jan
39
1
184
8
.
27.
Dimitriev
DA
,
Saperova
EV
,
Dimitriev
AD
.
State anxiety and nonlinear dynamics of heart rate variability in students
.
PLos One
.
2016
;
11
(
1
):
e0146131
.
28.
Cervantes Blásquez
JC
,
Rodas Font
G
,
Capdevila Ortís
L
.
Heart-rate variability and precompetitive anxiety in swimmers
.
Psicothema
.
2009
;
21
(
4
):
531
6
.
29.
Kanthak
MK
,
Stalder
T
,
Hill
LK
,
Thayer
JF
,
Penz
M
,
Kirschbaum
C
.
Autonomic dysregulation in burnout and depression: evidence for the central role of exhaustion
.
Scand J Work Environ Health
.
2017 Sep 1
43
5
475
84
.
30.
Rajcani
J
,
Solarikova
P
,
Turonova
D
,
Brezina
I
.
Heart rate variability in psychosocial stress: comparison between laboratory and real-life setting
.
Activitas Nervosa Superior Rediviva
.
2016
;
58
:
77
82
.
31.
Carnevali
L
,
Thayer
JF
,
Brosschot
JF
,
Ottaviani
C
.
Heart rate variability mediates the link between rumination and depressive symptoms: a longitudinal study
.
Int J Psychophysiol
.
2018
;
131
:
131
8
.
32.
Birkett
MA
.
The Trier Social Stress Test protocol for inducing psychological stress
.
J Vis Exp
.
2011 Oct 19
56
3238
.
33.
Medica-Torino
EM
.
Effects of anxiety due to mental stress on heart rate variability in healthy subjects
.
Minerva Psichiatr
.
2011
;
52
:
227
31
.
34.
Mohammadi
A
,
Emamgoli
A
,
Shirinkalam
M
,
Meftahi
GH
,
Yagoobi
K
,
Hatef
B
.
The persistent effect of acute psychosocial stress on heart rate variability
.
Egypt Heart J
.
2019 Sep 11
71
1
18
.
35.
Logan
JG
,
Teachman
BA
,
Liu
X
,
Farber
CR
,
Liu
Z
,
Annex
BH
.
Acute psychological stress, autonomic function, and arterial stiffness among women
.
Int J Psychophysiol
.
2020
;
155
:
219
26
.
36.
Delaney
JP
,
Brodie
DA
.
Effects of short-term psychological stress on the time and frequency domains of heart-rate variability
.
Percept Mot Skills
.
2000 Oct
91
2
515
24
.
37.
Nahshoni
E
,
Adler
E
,
Laniado
S
,
Keren
G
.
Fractal organization of the pointwise correlation dimension of the heart rate
.
Med Hypotheses
.
1998 Nov
51
5
367
76
.
38.
Dikecligil
GN
,
Mujica-Parodi
LR
.
Ambulatory and challenge-associated heart rate variability measures predict cardiac responses to real-world acute emotional stress
.
Biol Psychiatry
.
2010 Jun 15
67
12
1185
90
.
39.
Brosschot
JF
,
Van Dijk
E
,
Thayer
JF
.
Daily worry is related to low heart rate variability during waking and the subsequent nocturnal sleep period
.
Int J Psychophysiol
.
2007
;
63
(
1
):
39
47
.
40.
Miu
AC
,
Heilman
RM
,
Miclea
M
.
Reduced heart rate variability and vagal tone in anxiety: trait versus state, and the effects of autogenic training
.
Auton Neurosci
.
2009 Jan 28
145
1–2
99
103
.
41.
Schubert
C
,
Lambertz
M
,
Nelesen
R
,
Bardwell
W
,
Choi
JB
,
Dimsdale
J
.
Effects of stress on heart rate complexity: a comparison between short-term and chronic stress
.
Biol Psychol
.
2009
;
80
(
3
):
325
32
.
42.
Brugnera
A
,
Zarbo
C
,
Tarvainen
MP
,
Marchettini
P
,
Adorni
R
,
Compare
A
.
Heart rate variability during acute psychosocial stress: a randomized cross-over trial of verbal and non-verbal laboratory stressors
.
Int J Psychophysiol
.
2018 May
127
17
25
.
43.
Kuusela
T
.
Methodological aspects of heart rate variability analysis
.
Heart rate variability (HRV) signal analysis: clinical applications
.
2013
. p.
10
42
.
44.
Voss
A
,
Schulz
S
,
Schroeder
R
,
Baumert
M
,
Caminal
P
.
Methods derived from nonlinear dynamics for analysing heart rate variability
.
Philos Trans A Math Phys Eng Sci
.
2009
;
367
(
1887
):
277
96
.
45.
Voss
A
,
Kurths
J
,
Kleiner
HJ
,
Witt
A
,
Wessel
N
.
Improved analysis of heart rate variability by methods of nonlinear dynamics
.
J Electrocardiol
.
1995
28
Suppl l
81
8
.
46.
Spellenberg
C
,
Heusser
P
,
Bussing
A
,
Savelsbergh
A
,
Cysarz
D
.
Binary symbolic dynamics analysis to detect stress-associated changes of nonstationary heart rate variability
.
Sci Rep
.
2020 Sep 22
10
1
15440
.
47.
Pereira
T
,
Almeida
PR
,
Cunha
JPS
,
Aguiar
A
.
Heart rate variability metrics for fine-grained stress level assessment
.
Comput Methods Programs Biomed
.
2017 Sep
148
71
80
.
48.
Bernardi
L
,
Wdowczyk-Szulc
J
,
Valenti
C
,
Castoldi
S
,
Passino
C
,
Spadacini
G
.
Effects of controlled breathing, mental activity and mental stress with or without verbalization on heart rate variability
.
J Am Coll Cardiol
.
2000
;
35
(
6
):
1462
9
.
49.
Beda
A
,
Jandre
FC
,
Phillips
DI
,
Giannella-Neto
A
,
Simpson
DM
.
Heart-rate and blood-pressure variability during psychophysiological tasks involving speech: influence of respiration
.
Psychophysiology
.
2007 Sep
44
5
767
78
.
50.
Spangler
DP
,
Friedman
BH
.
Effortful control and resiliency exhibit different patterns of cardiac autonomic control
.
Int J Psychophysiol
.
2015 May
96
2
95
103
.
51.
Tulppo
MP
,
Mäkikallio
TH
,
Seppänen
T
,
Shoemaker
K
,
Tutungi
E
,
Hughson
RL
.
Effects of pharmacological adrenergic and vagal modulation on fractal heart rate dynamics
.
Clin Physiol
.
2001
;
21
(
5
):
515
23
.
52.
Klintworth
A
,
Ajtay
Z
,
Paljunite
A
,
Szabados
S
,
Hejjel
L
.
Heart rate asymmetry follows the inspiration/expiration ratio in healthy volunteers
.
Physiol Meas
.
2012
;
33
(
10
):
1717
31
.
53.
Baumert
M
,
Lambert
GW
,
Dawood
T
,
Lambert
EA
,
Esler
MD
,
McGrane
M
.
Short-term heart rate variability and cardiac norepinephrine spillover in patients with depression and panic disorder
.
Am J Physiol Heart Circ Physiol
.
2009
297
2
H674
9
.