Carmen Schiweck
Ali Gholamrezaei3
Maxim Hellyn
Thomas Vaessen- 1Department for Psychiatry, Psychosomatic Medicine and Psychotherapy, University Hospital, Goethe University, Frankfurt, Germany
- 2Psychiatry Research Group, Department of Neurosciences, KU Leuven, Leuven, Belgium
- 3Faculty of Medicine and Health, Pain Management Research Institute, The University of Sydney, Sydney, NSW, Australia
- 4University Psychiatric Centre, KU Leuven, Leuven, Belgium
Past research links depression and blunted cardiac vagal reactivity to chronic stress. Yet, to our knowledge no experiment investigates heart rate (variability) responses to a repeated laboratory stressor in patients with depression. Repeated exposure may provide valuable information on stress reactivity in depression. Fifty-nine women (30 inpatients diagnosed with depression and 29 matched controls) underwent two consecutive runs of a mental arithmetic stress paradigm consisting of one baseline and two exposures to control, stress, and recovery phases of 5 min each, in a case-control design. Subjective stress and electrocardiography were recorded. Variance of heart rate (HR) and root mean square of successive RR interval differences (RMSSD) were analyzed using linear mixed models. Overall, physiological parameters (HR and RMSSD) and subjective stress showed a strong group effect (all p < 0.001). In both groups, subjective stress and HR increased in response to stress, but the subjective stress levels of patients with depression did not return to baseline levels after the first stressor and for the remainder of the experiment (all p < 0.004 compared to baseline). Patients’ HR reactivity responded oppositely: while HR recovered after the first stress exposure, no reactivity was observed in response to the second exposure. These findings may suggest that the often-reported blunted HR/HRV response to stressors results from exhaustion rather than an incapacity to react to stress. The altered HR reactivity could indicate allostatic (over-) load in depression.
Introduction
Major Depressive Disorder (MDD) is currently diagnosed in an estimated 300 million people worldwide (1), conferring profound disability for those affected and an incremental economic burden for society (2). That the experience of psychological stress is among the most critical factors preceding the development of a first depressive episode has long been established (3, 4) and it is now also clear that stress acts not only as a catalyst, but also as an important factor for the maintenance of depression (5, 6). This auto-conserving cycle of chronic stress exposure and depressive episodes does bear deep and long-lasting consequences for underlying psychological and biological systems (7). Next to a generally overactive hypothalamus-pituitary-adrenal (HPA) axis (8, 9) and altered immunological reactivity (10), the Autonomic Nervous System (ANS) also shows significantly altered activity in both patients with clinical depression (11–14) and healthy participants with depressive symptoms (15–17).
We recently reviewed the literature to investigate physiological reactivity to self-relevant and active stressors in patients with clinical MDD (18). The majority of studies showed an overall blunted response during stress or impaired recovery in MDD, both for heart rate (HR) and heart rate variability (HRV). Contemporary theories of depression pathogenesis posit that this blunted stress reactivity is the effect of allostatic overload (19–21). Allostatic overload can be defined as dysregulation or wear and tear of the body, reflecting a state of physiological exhaustion, brought about by chronic and recurrent exposure to stressors. The mechanisms leading to such exhaustion include the failure to return to homeostasis after the stressor is eliminated, an insufficient response to cope with a stressor and the failure to physiologically and/or psychologically adapt to repeated stressors [also see Guidi et al. (22) for a detailed overview of these postulates]. The first and second of these factors are relevant for both the clinical assessment of allostatic load (20, 23) and for the clinical/biological presentation and treatment of MDD.
Regarding adaptation to repeated stressors, in an interesting recent article, the group of Ginty has investigated the cardiovascular adaptation capacity to repeated stressors in healthy controls using an identical stressor within one session (24). The authors have shown that greater perceived life stress was associated with blunted HR responses and less adaptation to the repeated stressor. However, they did not account for depressive symptoms, nor did they assess a patient group. As a result, it is still unclear whether the cardiovascular blunting of the stress response in patients with depression is due to chronic stress itself rather than depression and if, or how, the magnitude of the acute response is modulated by repeated stress experience in patients with depression. The quick adaptability of HR/HRV responses may provide deeper insights in the acute mechanisms linked to allostatic overload. Understanding the exact nature of the presumed exhaustion (i.e., gradual blunting) may not only advance the understanding of stress-related disorders and the further development of the allostatic load model but can also open new treatment avenues. To clarify whether and how repeated exposure to a stressor is altered in depression, and how this is affected by chronic stress we here set out to explore the effect of an identical repeated stressor on HR and HRV in patients with depression in response to an active stress task.
Several factors need to be taken into consideration in order to assess the above questions. Firstly, patients with MDD commonly encounter recurrent stressors of the same nature (3). Yet, to our knowledge, all current laboratory studies used one single, or two different stressors when investigating stress reactivity in MDD, showing (mostly) blunted stress reactivity to active stressors and (mostly) hyperactivity to passive stressors [for an overview see Schwerdtfeger and Rosenkaimer (17)]. Therefore, for the first time, we here explore the effect of an identical, repeated stressor on cardiac parameters in patients with depression in response to an active stress task within the same session.
Secondly, it is well known that variables such as sex (25, 26), age (27), smoking status (28), antidepressant medication (29) and Body Mass Index (BMI) (30) influence HR/HRV and potentially physiological stress reactivity and should be taken into account to draw firm conclusions. To address these potentially confounding factors, we here opted to focus on women, and match patients with MDD and controls for age and where possible BMI, and restricted smoking immediately prior to the experiment. We chose women, since depression is more prevalent in women and it has been suggested that blunting of the stress response may be particularly relevant (31).
In addition, the effect of chronic stress is highly relevant for stress reactivity analyses in patients with depression, since stress/stressful life events often precede the onset of depressive symptoms (14, 32) and chronic stress per se may be linked to blunted stress reactivity: Ginty et al. (33) showed that in healthy individuals, higher perceived stress (compared to objective stress exposure) was linked to blunted cardiovascular stress reactivity (33). However, the effect of perceived stress is often not examined in patients with depression. To disentangle the effects of stress and depression diagnosis we will therefore attempt to assess the moderating role of chronic stress levels on the effect of depression diagnosis. Lastly, recent studies have successfully used HR/HRV as biomarker for MDD (34–36). In order to use HR reactivity during stress as stratification marker for new treatment options, this would require a reasonable sensitivity and specificity for depression. Therefore, we will here attempt to replicate stress reactivity as marker to distinguish patients with MDD.
Thus, we here explore HR and HRV reactivity to an identical, repeated stressor in a carefully controlled design to test whether (a) stress reactivity/recovery to a repeated stress task is altered in depression (exploratory analysis); (b) higher levels of (perceived) chronic stress are associated with a blunted stress response; (c) whether this is different in healthy individuals compared to patients; and lastly, (d) whether HR/HRV reactivity to- and recovery from stress can be used to distinguish patients and controls. We hypothesize that (a.1.) congruent with the literature, HR reactivity to the first stressor is blunted and (a.2.) that recovery from stress is impaired in patients with MDD compared to controls. Additionally, we explore whether the second stress exposure differs from the first exposure with regard to magnitude (e.g., increased blunting, increased reactivity, or no qualitative change); for (b) we hypothesize that participants with depression who have higher levels of chronic stress show blunted stress reactivity; (c) in contrast, we expect that healthy controls with high levels of chronic stress show increased reactivity; and lastly that (d) combined HR and root mean square of successive difference (RMSSD) reactivity during stress and/or recovery can be used to distinguish patients and controls.
Materials and Methods
This study was performed according to Good Clinical Practice and approved by the local ethics committee (EC number: S59441, UZ Leuven Medical Ethics Committee, Belgium). Participants gave written informed consent prior to study related activities. The study recruitment and tests ran from 2016 to 2018.
Study Design
The study consisted of a screening visit during which the diagnostic interviews and questionnaires were administered, and a test session held between 12.00 and 18.30 during which participants performed two runs of the 5 min repeated mental arithmetic stress task separated by 15 min of break. During the entire task, participants were wearing an ECG. See Figure 1 for an overview of the study design and below for detailed procedures.
Figure 1. Participants underwent extensive psychological screening and completed a computerized, repeated version of the mental arithmetic stress task including social evaluative components. This procedure was administered twice in one session. For a complete description see Supplementary Information A. Design: “images: Flaticon.com.” This figure has been designed using resources from Flaticon.com. Stress icons created by Freepik – Flaticon (https://www.flaticon.com/free-icons/stress). Psychology icons created by iconixar – Flaticon (https://www.flaticon.com/free-icons/psychology). Calm icons created by Becris – Flaticon (https://www.flaticon.com/free-icons/calm).
Participant Selection and Procedure
Thirty women with MDD (56% undergoing antidepressant pharmacological therapy) and 29 women without mental health diagnosis matched for age (tolerance ± up to 5 years) and where possible BMI, underwent an extensive diagnostic interview (MINI PLUS Version 6.0.0, Dutch version) following the Diagnostic and Statistical Manual of Mental Disorder 4th edition (DSM–IV). Patients were currently treated for depression as in- or outpatients in the UPC KU Leuven Psychiatric Hospital. Controls had to be free of any lifetime psychiatric disorder. All participants had to be free of cardiovascular disease based on self-report and were excluded if they were taking medication with significant effects on cardiovascular function (e.g., antihistamines, beta-blockers, and tricyclic antidepressants) at the exception of SSRIS, SNRIS, atypical antidepressants or benzodiazepines (as on-demand medication), and contraceptive medication. Demographic data including menstrual cycle phase was recorded. One control participant and 3 patients had to be excluded from analysis [due to prohibited medication on the testing day (antihistamine) n = 1, arrhythmias n = 1 discovered during analysis of data, drop-out before stress task n = 2], leaving a total sample of 28 controls and 27 patients.
Questionnaires
The 10-item version of the self-rated Perceived Stress Scale (PSS) was used to assess levels of chronic stress (37). It is a well-validated scale commonly used to assess stress perception during the past 4 weeks and has good psychometric properties (38). The Dutch self-rated version of the Rumination Response scale (RRS) was used to assess whether patients had higher levels of rumination than controls. The Dutch RRS has shown good reliability and satisfactory validity (39). The General Practice Physical Activity Questionnaire was used to assess activity levels of participants yielding four different categories of physical activity from inactive to very active [National Collaborating Centre for Nursing and Supportive Care (United Kingdom) (40)]. The self-rated, Depression Anxiety and Stress Scale (DASS) 21-item version was used to assess general symptoms of stress, depression and anxiety in the past week (41, 42). It has good internal consistency and validity was also good in Dutch participants (43). In addition, to assess depressive symptoms, the Hamilton Rating Scale for Depression (HRSD) 17-item version was used to assess depression severity in participants (44, 45). The HRSD is the gold standard to assess clinician-rated, clinical symptoms of depression and has generally good validity and reliability (46).
Repeated Modified Mental Arithmetic Task
Detailed information on the study design can be found in Supplementary Information 1. In short, the repeated modified mental arithmetic task (RMMAT) was computed in Affect version 4.0 (47). It is a modified version of the mental arithmetic/Montreal Stress Imaging task, and includes a baseline, a control task (i.e., a mental arithmetic without stressful elements), a mental arithmetic (stress), and a recovery phase. In this modified, computerized version, several elements for social evaluation, time pressure and feedback by the researcher are integrated to induce stress during the stress phases. After the first exposure to the RMMAT and the first recovery period, the researcher asked the participants unexpectedly to perform the RMMAT a second time. The participants were not aware in advance that this second administration would take place. This was done on purpose to be able to assess the reactivity to repeated stressors without confounding the effect by anticipatory stress.
Procedure
Within 7 days after the screening visit (including the diagnostic interview and questionnaires), participants underwent the RMMAT. If more than 7 days had elapsed between screening and the stress task, questionnaires were re-administered. To avoid bias, participants were instructed not to smoke at least 1 h prior to the start of the experiment, were asked not to consume caffeine and restrain from sport on the day of, and from heavy exercise the day prior to the experiment.
Upon arrival, participants were placed in a comfortable chair in front of a computer, instructed to find a comfortable position and move as little as feasible. Electrodes for cardiovascular measurement and a respiratory belt were attached as described below and tested for proper functioning. Next, the experimenter explained the setup of the experiment, omitting to mention the second exposure to the stress task. After the first stress task the experimenter re-entered the room and informed participants that unexpectedly, they had to perform the task again. At baseline, after the stress and during the recovery phases of both stress exposure experiments, participants were asked to rate their stress level on a 10-point Likert scale (0 = not at all to 10 = extremely stressed). Total task duration was approximately 60 min and the time of assessment was between 12.00 and 18.30 (i.e., after the diurnal normal increase in HR and before the diurnal drop in the evening). The choice to use the timeslot between 12.00 and 18.30 was based on our pilot study on diurnal variations in HR in patients with MDD and healthy controls, where a prominent increase in basal HR occurred during the morning, whereas the most significant decrease occurred later in the evening. The period-between 12.00 and 18.30 yielded a plateau ideal for behavioral analyses (48). Participants were in a seated position in a room with a controlled temperature of 22°C. After completion, participants were debriefed, thanked and discharged.
Electrocardiography and Respiratory Recording
Cardiac electrical activity (ECG) and respiratory movement were recorded using a 10-channel isolated bioamplifier (NeXus-10 MKII, Mind Media BV, Netherlands). Three disposable electrodes (Kendall™, H66LG electrodes) were placed in a lead II configuration on the chest to record ECG at sampling frequency of 1,024 Hz. Simultaneously, respiratory (chest) movement was recorded by placing a strain gauge belt over the abdomen, with the central part just above the navel and a sampling rate of 32 Hz was used, as instructed by the manufacturer. Correct functioning of the recordings was checked using the BioTrace+ Software (Mind Media BV, Netherlands).
Data Analysis and Cleaning
The ECG signal was processed and R-wave peaks were detected using Kubios software version 2.0 (49). A visual inspection of raw data for possible movement artifacts and ectopic beats was performed per segment (i.e., each 5-min interval). Where applicable, an inbuilt artifact correction based on cubic spline interpolation was applied. Segments with more than 5% of beats in need for correction were excluded from analysis (one segment). Data was excluded for one patient showing a high number of ectopic beats (arrhythmia), and two segments of data were lost for two participants due to equipment failure. The respiratory movement signal was processed using custom-written algorithms in MATLAB software (R2018b, Mathworks, Inc., Natick, MA, United States), see Gholamrezaei et al. (50), and breathing rate was measured for each segment.
Heart Rate and Heart Rate Variability
Heart rate variability analysis was performed using the Kubios software on segments with duration of 5 min according to the guidelines (51). HR, RMSSD, and high-frequency HRV (HF-HRV) were extracted for each segment. RMSSD is a time-domain HRV index reflecting beat-to-beat variation in HR, and HF-HRV is a frequency-domain HRV index reflecting variability of HR at spontaneous breathing frequency (i.e., respiratory sinus arrhythmia or RSA). Both indexes are mainly derived by cardiac parasympathetic (vagal) activity as is shown by studies using vagolytic (e.g., high dose atropine) (52–55) or vagotonic (e.g., low dose atropine) agents (56–58) to inhibit or increase cardiac parasympathetic stimulation, respectively (59). We decided to use RMSSD as the primary HRV index in our analysis for specific reasons: Evidence has shown that RMSSD is less influenced by breathing behavior, particularly the breathing frequency, compared to HF-HRV (54, 55, 60). Controlled breathing is recommended for a reliable assessment of cardiac vagal activity when using frequency-domain measures (54, 55). Therefore, RMSSD is a more suitable measure of cardiac vagal activity when controlled breathing is not possible/feasible, as was the case in the current study. Notwithstanding, to provide more comparable data with previous studies, we decided to also analyze HF-HRV while taking breathing frequency into account. Results of HF-HRV and for completeness, low frequency (LF-HRV) are provided in the Supplementary Information D.
Statistical Analysis
Statistical analyses were performed in R (61) and matching of participants and controls was performed using SPSS version 25 (62). The power analysis on simulated data can be found in Supplementary Information B. We conducted a repeated-measures mixed model with random intercept, a phase (baseline, stress, recovery) × group (MDD/control) interaction and PSS × group interaction as independent variables, and either subjective stress rating, or HR/HRV as dependent variables. Assumptions of mixed models were checked by visual inspection of q–q-plots and histograms. Due to slight heterogeneity of residuals, HR and HRV indexes were transformed using the logarithmic (log) transformation for all statistical analyses. The initially planned interaction of PSS × phase to test the moderating effect was collinear with the group × phase interaction and was therefore not explored in the overall model. Nevertheless, to investigate this effect we conducted a separate model within the groups (either MDD or HC) with a PSS × phase. Significant time × group interactions were followed-up by planned contrasts adjusted for multiple testing using Benjamini–Hochberg correction between baseline and stress, between stress and recovery and between baseline and recovery of each run. Since no difference occurred between the baseline and the control condition (see Supplementary Information C), and the baseline phase was more comparable to the recovery phases, the control condition was not considered in the model.
To investigate classification accuracy of physiological parameters, we used binary logistic regression with leave one out cross validation (LOOCV) from the caret package (50).
Results
Demographic Data and Baseline Differences
Patients and controls were female, between 18 and 60 years old, matched for age (using fuzzy matching, tolerance up to ±5 years) with a similar BMI. Menstrual phases (contraception, post-menopausal, follicular, or luteal phase) did not differ significantly between groups (X2 = 3.54, df = 3, p = 0.315). The mean HRSD score for patients was 24.63 indicating a severe form of depression (51). Controls were more likely to be non-smokers (X2 = 13.55, df = 2, p = 0.001) and were significantly more active than patients (W = 634, p < 0.001). Additionally, patients had significantly higher depression severity, rumination- and perceived stress levels in the last month (all p < 0.001, see Table 1). None of the control population and around half (15 of 27) of the patients were taking SSRIS/SNRIs/SARIS (n = 11) or atypical antidepressants (n = 4). Four patients were exclusively treated with benzodiazepines.
Table 1. Demographic data and statistical comparison.
Effect of Covariates on Baseline Levels of Heart Rate/Root Mean Square of Successive Difference
While antidepressant medication intake was limited to SSRIs/SNRIs/SARIs and atypical antidepressants and smoking immediately prior to the experiment (at least for 1 h) was prohibited, these parameters may constitute an important confounder. We therefore tested whether the demographic variables that were different between controls and patients (i.e., medication, smoking, and physical activity levels) contributed to HR/HRV group differences at baseline. No main or moderating effect occurred for smoking (all p > 0.231) or physical activity (all p > 0.765). Antidepressant medication had no main effect on HR at baseline (t = −0.33, p = 0.748), but significantly decreased RMSSD (t = −2.46, p = 0.021).
Subjective Stress in Response to the Stress Task
As a manipulation check, we tested whether subjective stress induction was achieved. A significant group by phase interaction occurred (F = 3.11, df = 4, p = 0.016). Post hoc contrasts of phase by group showed the expected increase of SS in healthy controls, during both runs (compared to baseline t = 7.57, p.adj < 0.001 and t = 6.64, p.adj < 0.001), returning to approximate baseline levels during recovery (t = −1.48, p.adj = 0.154 and t = −0.37, p.adj = 0.712). Congruent with the control’s stress response, in patients, SS increased in response to the first stressor (t = 6.21, p < 0.001). However, while there was a subtle recovery after stress (t = −0.89, p.adj = 0.030), SS levels did not return to baseline during the recovery period (t = −3.95, p < 0.001) and stayed significantly elevated throughout the remainder of the experiment (all p < 0.004 compared to baseline), indicating incomplete subjective stress recovery (Figure 2 and Supplementary Table 1).
Figure 2. Subjective stress ratings per phase. The black (bottom) line represents control participants, the red (top) line represents patients. BL, baseline; S1, first stressor; R1, first recovery phase; S2, second stressor; R2, second recovery phase. Summary statistics indicate mean ± standard errors. All change scores were significantly different for the control group. Patients and healthy control individuals showed significant reactivity to the first stressor. In healthy individuals, stress levels recovered after both stress exposures, while patients did not return to baseline.
Heart Rate Reactivity to Stress
Mean levels of HR and RMSSD reactivity can be found in Supplementary Table 2. The final model included a group × phase interaction (and respective main effects), a group × PSS interaction (and respective main effects) and antidepressant medication. The PSS × group interaction (F = 10.41, p = 0.002) and the group × phase interaction (F = 2.99, p = 0.019), and antidepressant medication (F = 5.78, p = 0.020) were significant. Patients with MDD had overall higher HR throughout the experiment (Figure 3A). Post hoc contrasts revealed that in controls, HR increased significantly during the first (t = −5.11, p.adj < 0.001) and second stressor compared to baseline (t = −4.01, p.adj = 0.001) and returned to approximate baseline levels (or slightly lower levels) during first (t = 0.89, p.adj = 0.421) and second recovery phase (t = 2.02, p.adj = 0.080) compared to baseline (see Table 2 for exact statistics).
Table 2. Results of planned within group contrasts for split linear mixed models for HR reactivity.
Figure 3. Physiological reactivity to stress. (A) Heart rate (HR): significant overall reactivity and a group × timepoint interaction occurred. Controls show significant reactivity and recovery, while patients show an increase of HR to the first stressor and a decrease during recovery 1, but no longer react to the second stressor (neither compared to BL nor compared to recovery 1). (B) Root mean square of successive differences (RMSSD). Controls show significant response after the first stressor (recovery 1), the second stressor and second recovery, MDD patients show no fluctuation at any time point. Heart rate is displayed in beats per minute; MDD, Major Depressive Disorder; HC, healthy controls.
In patients with MDD, HR significantly increased during the first stressor (t = −2.90, p.adj = 0.011 compared to baseline) and returned to baseline levels during recovery (t = 0.43, p.adj = 0.667). However, no significant reaction to the second stressor was observed, neither compared to baseline, nor the first recovery period (t = −0.67, p.adj = 0.534 and t = −1.09, p.adj = 0.346, respectively). During the second recovery phase, a non-significant reduction of HR occurred compared to the second stressor (t = 2.16, p.adj = 0.062) (Figure 3A). Exploration of the significant group × PSS interaction showed an inverse U shape at baseline, where high stress levels were associated with higher HR in healthy controls, but with lower HR in patients with depression (see Supplementary Figure 1). For an overview of planned contrasts per phase for patients and controls, see Table 2. Of note, in an additional exploratory model within the MDD group, the main effect of PSS (F = 19.16, p = 0.002), but not the phase × PSS interaction (F = 1.22, p = 0.306) was significant, indicating that chronic stress influenced overall HR levels but did not (additionally) alter the reactivity in patients with MDD. No effect of PSS (F = 1.26, p = 0.271) or the phase × PSS interaction (F = 0.370, p = 0.830) occurred for healthy controls. It should however be noted that healthy controls had in general much lower stress levels than patients with MDD and an exact disentanglement of group and stress effects was therefore not possible.
Root Mean Square of Successive Difference Reactivity to Stress
For RMSSD reactivity, group (F = 13.98, p < 0.001), phase (F = 6.36, p < 0.001) and a group × PSS interaction (F = 4.26, p = 0.043) were significant. Antidepressant medication use was not a significant predictor (F = 0.35, p = 0.557), but since RMSSD is known to be influenced by this measure, was kept in the model. Post hoc contrasts by group revealed that controls showed significantly lower RMSSD levels during the second stress task when compared to the first recovery phase (t = 3.18, p.adj = 0.010) and significant recovery after the second stress exposure (t = −4.81, p < 0.001). Reactivity of the first stressor compared to baseline was not significant (t = −1.92, p.adj = 0.168). Patients with depression showed no significant reactivity at any time point (all p.adj > 0.170). The data is displayed in Figure 3B and Supplementary Table 2.
Heart Rate During Recovery as Diagnostic Marker for Depression
To assess if, and which phase of HR reactivity could best be used to distinguish depressed patients from controls, we performed binary logistic regression with LOOCV. The best distinction of patients and controls was achieved using HR with RMSSD during recovery. Supplementary Table 3 shows that, as a unique predictor, HR outperformed HRV indexes for classification. During recovery, combining these indexes led to a slightly better classification than for HR alone. Predictive accuracy was lowest during stress (65.45 and 63.63%) and highest during recovery from stress (first recovery 75.47% and second recovery 72.22%, respectively). Receiver operating characteristic (ROC) curves comparing various conditions can be found in Figure 4.
Figure 4. Receiver operating characteristic curves demonstrating sensitivity and specificity of the classification. Classification based on heart rate (HR) during various conditions. BL, baseline; STRESS, stress 1 and 2; and RECOVERY, recovery 1 and 2.
Discussion
The present study aimed to explore the effect of repeated exposure to an identical stressor in depression. Results showed marked elevations in HR, lower vagally mediated HRV at rest and insufficient subjective stress recovery in women with depression. The results confirm previous findings showing blunted cardiac vagal modulation during stress in patients with MDD, e.g., Schiweck et al. (18) and Duschek et al. (63). In addition, we extend the literature by showing that repeated stress exposure provoked strongly reduced HR responses and poor subjective recovery in patients with depression, as opposed to healthy controls. Lastly, this pattern of cardiac reactivity differentiated between individuals with MDD and healthy volunteers with reasonably good accuracy and may have value as a complementary assessment tool for depression.
Blunted Heart Rate Responses in Major Depressive Disorder: Adaptive Habituation or System Exhaustion?
As expected from the literature (64–66), our data clearly show overall higher subjective stress ratings in patients with MDD compared with healthy controls, both over the past month (as measured by the PSS) and during all phases of the experimental manipulation (as measured with the VAS). We here also replicate the finding that patients with depression show higher basal levels of HR and lower levels of HRV than controls (29), thereby confirming the representability of our sample. Most intriguing, however, is the finding that patients could not subjectively recover from stress and showed blunted reactivity to the second stressor.
Physiological habituation, conceived as significantly attenuated responses to a second stress exposure, can be influenced by many factors. In their review, Hughes et al. (67) conclude that healthy controls often show habituation of HR to stressors comparable to the here used MAT, but that this process can be disrupted by the presence and observation of the participants by the experimenter (60, 61, 68, 69). Since we designed the task to contain psychosocial stress elements including observation by the experimenter, it is not surprising that healthy controls do not show (significant) adaptation to the stressor. In our opinion, the psychosocial stress elements are crucial to test relevant stressors since psychosocial stress is highly important for depression pathogenesis. Opposed to subjective stress reactivity, HR levels in patients with depression did fully recover from the first stress exposure but did not show any reactivity to the second stressor. This may suggest a rapid habituation in MDD patients, despite the above-mentioned psychosocial stress elements, or alternatively may point to a dysfunctional or exhausted stress response. The first scenario would suggest that the observed physiological non-responsivity is an adaptive strategy to reduce energy expenditure. From this perspective, while essential to mount sufficient stress responses in the face of acute danger, a fast adaptation to a similar stressor would be very beneficial for energy resources. Yet, adaptive habituation would fail to account for the chronic high stress levels, increased HR and consistently blunted HR variability that have been observed in MDD (18, 70) and is therefore an unlikely hypothesis.
The physiological response could thus perhaps find a more integrative explanation when conceived as the consequence of allostatic load, as introduced by McEwen and Stellar (19, 21, 22, 71, 72). The authors propose that an initially adaptive mechanism to mobilize energy resources (allostasis) can lead to considerable wear and tear if it persists long term, as is the case with chronic stress (allostatic load). In line with these observations is a recent study: In a large population of healthy young adults, high levels of stress also led to blunted HR responses and additionally led to lower habituation over time (24), although the authors did not account for depressive symptomatology. In our sample, the group of healthy people was too small to investigate this differential effect in detail. However, we did observe a qualitatively different effect of chronic stress on baseline HR in healthy controls compared to patients with depression: While chronic stress increased HR in the controls, it decreased HR in patients with depression. It is therefore essential to account for this relationship in studies involving patients with MDD and would be highly informative to assess the effect of depressive symptoms in a large healthy working population.
The here observed discrepancy of subjective and physiological responses to stress in depression has also been found by others. In their study, Söder et al. found a higher discrepancy of recorded HR and reported subjective stress in patients with depression compared to controls during a noise stressor (73). This discrepancy of higher self-reported stress and relatively lower HR reactivity may contribute to the maintenance of chronic stress levels and deserves further investigation.
The repeated stress paradigm as tested in this study, albeit an artificial exposure, still seems highly relevant for the daily experience of patients with clinical depression, who often report a continuous feeling of high subjective stress and suffer from repeated minor stressors with which they seem unable to cope adequately. This inability to cope with repeated stressors might be somehow related to the decoupling of the biological stress response, caused by exhaustion. A next step to test this hypothesis would be to test patients with depression in daily life, carefully recording stressors, subjective responses, and physiological parameters using wearable devices.
Blunted Cardiac Vagal Modulation to Psychosocial Stress
In response to the first stressor, patients with MDD showed increased subjective stress levels of approximately equal magnitude as controls, and a corresponding increase in HR. However, RMSSD showed nearly no variation throughout all phases of the task within the MDD population. While healthy controls also did not show vagal withdrawal to the first stressor, this observation could possibly be attributed to anticipatory stress levels at baseline. This is supported by the apparent recovery after the first stressor and significant reactivity to and recovery from the second stressor.
Ample research has documented blunted HR and HRV responses to stress in both participants with depressive symptoms (70) and patients with clinical depression (18). While we did not observe altered HR reactivity in patients with depression in response to the first stressor, the response to the second stress exposure was blunted, although the blunting could perhaps better be described as exhausted stress response (see section “Blunted Heart Rate Responses in Major Depressive Disorder: Adaptive Habituation or System Exhaustion?”). On the other hand, the association between depression and blunted cardiac vagal reactivity was confirmed and may point to consistently exhausted cardiac vagal modulation.
Blunted reactivity to stress in depression is not a novel concept in stress research: several reviews have investigated the matter with a focus on HPA axis reactivity. The most recent meta-analysis reveals that active, self-relevant stressors provoke similar, though sex-specific changes in cortisol release: whereas female participants with MDD show blunted reactivity to stress, male participants with MDD show increased responses (74). This corroborates findings with regard to the here described blunted cardiac vagal modulation in response to stress in our sample of female participants with MDD and highlights the importance of a careful study design accounting for relevant variables such as sex. The biological background behind this blunted reactivity alongside elevated HR to the first stressor may point toward primarily sympathetic activation rather than parasympathetic withdrawal during stress. This could be highly relevant for current therapies such as vagus nerve stimulation for depression, since restoring the balance between parasympathetic and sympathetic nervous systems should evaluate both sympathetic and parasympathetic output to quantify efficacy. However, as nerve activity was not assessed directly in this study, it is impossible to draw conclusions on this hypothesis.
Decoupled Subjective and Physiological Stress Recovery: A Hallmark of Depression?
Despite the lack of (complete) subjective recovery after the first stress exposure, the patient group did show significant physiological recovery. This is interesting, as these results may inform about a possible process contributing to the maintenance of subjective stress levels: a failure to subjectively recover from minor stressors, as observed here, may lead to continuous maintenance of perceived stress and ultimately, increase the impact estimation of a stressor (66). This altered stress perception may also be the cause for the change in physiological parameters at baseline: e.g., chronic stress may increase basal HR levels and reduce parasympathetic activity.
The “why’s” and “how’s” of these physiological changes and their link to subjective stress can find a possible explanation in the literature of “perseverative cognition”: Brosschot et al. proposed that particularly rumination (both retrospective and prospective), as frequently encountered in depression, can negatively impact the consequences of psychological stress (75). As such, repetitive or “perseverative” thoughts about stress, and not necessarily the stressor itself, would lead to persistent physiological changes such as a chronically increased HR, blood pressure and altered cortisol levels (69). While the present study was not designed to directly test hypotheses related to perseverative cognition, it is conceivable that it impacts these physiological parameters. Particularly patients with depression did have high rumination scores. Yet, in this study, perceived stress rather than rumination scores increased model performance for the group with depression considerably, while rumination and perceived stress did not improve model fit for controls (data not shown). Therefore, preservative cognition may not be the most relevant factor for physiological changes in depression when accounting for chronic stress levels.
Classification of Patients and Clinical Relevance: A Marker for Depression?
In the present study, the most accurate classification of patients could be achieved based on HR, with slightly improved classification when combined with RMSSD. Differences were largest during recovery from stress and lower accuracy was achieved during stress. Others have recently published results on classifying individuals with depression based on cardiovascular parameters, with similar accuracy as here reported. To illustrate, Kuang et al. found high classification accuracy (86.4%) with RMSSD throughout a series of physiological ANS stimuli (36). In a stress task, Byun et al. also found the clearest distinction with a HRV feature (SDNN) during stress recovery (34). Of note, although HR or HF-HRV were not the most important features reported by Byun et al., we here specifically focused on these parameters, given their biological relevance and relative robustness regarding interpretability.
The usefulness of an objective biomarker is obvious in a time where remote monitoring for the clinic becomes tangible. While this additional data does not replace diagnoses or contact with a clinician it may help making informed decisions for possible treatment augmentation, in patients with partial response to therapy, or for secondary prevention in recovered patients. These findings could also be particularly interesting in the light of cardiovascular health. This is certainly important in light of the increased risk for cardiovascular disease observed in depression (76). While we here did not have longitudinal data to ascertain whether HR could be used as a state marker, others have shown that in remitted patients, cardiovascular parameters normalized toward the control group (77). More recently, Hartmann et al. showed that in their brief longitudinal study, indeed, HF-HRV (significantly) and RMSSD (non-significantly) increased in patients with depression in the course of a 2-week antidepressant treatment regimen, emphasizing the need for further longitudinal studies to assess the feasibility of cardiovascular measurements as state markers for depression (78).
Limitations
Firstly, the generalizability of this study is limited to female participants. Due to the well-documented sex differences in autonomic activity and sample sizes necessary to control for this effect, this was a necessary choice. However, menstrual status is known to have important effects on autonomic (re)activity. While recorded in this study, and an approximately equal proportion of patients and controls in each menstrual phase, analyses were underpowered to include this factor as covariate. Larger studies should include menstrual cycle as covariate and explore whether our findings are also applicable to men. It should also be mentioned that subjective stress experience and physiological reactivity were not always congruent. This well-known apparent incongruence (79) can arise from confounding factors such as age (80), sex (81) or the nature of the stressor. Other, external factors in larger samples may show different results. For instance, it has been shown that early life adversity can have profound effects on (physiological) stress reactivity (82), and future studies should account for this variable. In addition, habituation might better be tested in a stress paradigm without social stressors. However, this was not yet clear when this study was conceived in 2015, and social stress may actually be crucial to test ecologically relevant habituation (i.e., social compared to non-social stressors). A frequently reported measure for HRV analyses is LF-HRV. LF-HRV has previously been interpreted as reflecting cardiac sympathetic innervation. Yet, in resting conditions (as is the case here for parts of our experiments, the LF band reflects baroreflex activity and not cardiac sympathetic innervation (83). Interpretation of LF-HRV thus is complex and is not suitable to investigate our study aims. However, for the sake of completeness, it is reported in Supplementary Section. Finally, this is a cross-sectional study and therefore by definition, our results are correlational. Large-scale, longitudinal studies are needed to further investigate the potential of HR and HRV during stress as state markers for depression and its usefulness for the clinic.
Conclusion and Future Research
We here show that stress reactivity is functional in patients with depression but appears exhausted in face of recurrent stress. The specific design of this study enabled us, for the first time, to experimentally demonstrate physiological (HR/HRV) exhaustion of the stress response in depression, lending support to Mc Ewen’s conception of “allostatic load” as destructive consequence of chronic stress, and in particular its importance for depression. Importantly, our results also show that re-exposure to a stressor of the same nature can be used to study exhaustion. Future studies should therefore replicate these findings and extend them to male subjects.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, doi: 10.17632/6jzfpzmks6.1.
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethische Commissie Onderzoek UZ/KU Leuven (EC Onderzoek). The patients/participants provided their written informed consent to participate in this study.
Author Contributions
CS: conceptualization, data curation, formal analysis, investigation, methodology, project administration, software, visualization, writing—original draft, and writing—review and editing. AG: formal analysis, methodology, software, and writing-review and editing. MH: data curation, investigation, and writing—review and editing. TV: methodology, and writing—review and editing. EV: conceptualization, methodology, project administration, funding acquisition, supervision, roles/writing—original draft, and writing—review and editing. SC: conceptualization, funding acquisition, methodology, project administration, resources, validation, roles/writing—original draft, and writing—review and editing. All authors contributed to the article and approved the submitted version.
Funding
This work was sponsored by the Agency for Innovation by Science and Technology (IWT) with grant number 140745.
Conflict of Interest
CS is a senior clinical researcher from the Fund for Scientific Research Flanders (FWO Vlaanderen).
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank everyone who participated in this study. We would like to thank Mathijs Franssen and Jeroen Clarysse for providing technical support with affect. This data was presented as a poster at two meetings with published abstracts (71, 72). R packages used include ggpubr package in R (73), caret, MumIn, psych, lme4, car, Hmisc, gplots, HLMDiag, DHARMa, pROC, and broman.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpsyt.2022.869608/full#supplementary-material
References
1. WHO.Depression and other Common Mental Disorders: Global Health Estimates. Geneva: World Health Organization (2017).
2. Greenberg PE, Fournier A-A, Sisitsky T, Pike CT, Kessler RC. The economic burden of adults with major depressive disorder in the United States (2005 and 2010). J Clin Psychiatry. (2015) 76:155–62. doi: 10.4088/jcp.14m09298
3. Hammen C. Stress and Depression. Annu Rev Clin Psychol. (2004) 1:293–319. doi: 10.1146/annurev.clinpsy.1.102803.143938
4. Mazure CM. Life Stressors as Risk Factors in Depression. Clin Psychol Sci Pract. (1998) 5:291–313. doi: 10.1111/j.1468-2850.1998.tb00151.x
5. Hammen C. Generation of stress in the course of unipolar depression. J Abnorm Psychol. (1991) 100:555–61. doi: 10.1037//0021-843x.100.4.555
6. Liu RT, Alloy LB. Stress generation in depression: a systematic review of the empirical literature and recommendations for future study. Clincal Psychol Rev. (2010) 30:582–93. doi: 10.1016/j.cpr.2010.04.010
7. Brotman DJ, Golden SH, Wittstein IS. The cardiovascular toll of stress. Lancet. (2007) 370:1089–100. doi: 10.1016/S0140-6736(07)61305-1
8. Dinan TG. Glucocorticoids and the genesis of depressive illness a psychobiological model. Br J Psychiatry. (1994) 164:365–71.
9. Pariante CM, Lightman SL. The HPA axis in major depression: classical theories and new developments. Trends Neurosci. (2008) 31:464–8. doi: 10.1016/j.tins.2008.06.006
10. Wohleb ES, Franklin T, Iwata M, Duman RS. Integrating neuroimmune systems in the neurobiology of depression. Nat Rev Neurosci. (2016) 17:497.
11. Brindle RC, Ginty AT, Conklin SM. Is the association between depression and blunted cardiovascular stress reactions mediated by perceptions of stress? Int J Psychophysiol. (2013) 90:66–72. doi: 10.1016/j.ijpsycho.2013.06.003
12. Bylsma LM, Salomon K, Taylor-Clift A, Morris BH, Rottenberg J. Respiratory sinus arrhythmia reactivity in current and remitted major depressive disorder. Psychosom Med. (2014) 76:66–73. doi: 10.1097/psy.0000000000000019
13. Hu MX, Lamers F, de Geus EJ, Penninx BW. Differential Autonomic Nervous System Reactivity in Depression and Anxiety During Stress Depending on Type of Stressor. Psychosom Med. (2016) 78:562–72.
14. Lloyd C. Life events and depressive disorder reviewed. II. Events as precipitating factors. Arch Gen Psychiatry. (1980) 37:541–8. doi: 10.1001/archpsyc.1980.01780180055005
15. Carroll D, Phillips AC, Hunt K, Der G. Symptoms of depression and cardiovascular reactions to acute psychological stress: evidence from a population study. Biol Psychol. (2007) 75:68–74. doi: 10.1016/j.biopsycho.2006.12.002
16. de Rooij SR, Schene AH, Phillips DI, Roseboom TJ. Depression and anxiety: associations with biological and perceived stress reactivity to a psychological stress protocol in a middle-aged population. Psychoneuroendocrinology. (2010) 35:866–77. doi: 10.1016/j.psyneuen.2009.11.011
17. Schwerdtfeger A, Rosenkaimer AK. Depressive symptoms and attenuated physiological reactivity to laboratory stressors. Biol Psychol. (2011) 87:430–8. doi: 10.1016/j.biopsycho.2011.05.009
18. Schiweck C, Gholamrezaei A, Vaessen T, Vrieze E, Claes S. Exhausted cardiac vagal adaptation to a recurrent laboratory stressor in women with major depressive disorder. 32nd European College of Neuropsychopharmacology (ECNP). Copenhagen: ECNP (2019).
20. Fava GA, McEwen BS, Guidi J, Gostoli S, Offidani E, Sonino N. Clinical characterization of allostatic overload. Psychoneuroendocrinology. (2019) 108:94–101. doi: 10.1016/j.psyneuen.2019.05.028
21. McEwen BS. Stress, adaptation, and disease: allostasis and allostatic load. Ann NY Acad Sci. (1998) 840:33–44. doi: 10.1111/j.1749-6632.1998.tb09546.x
22. Guidi J, Lucente M, Sonino N, Fava GA. Allostatic load and its impact on health: a systematic review. Psychother Psychosom. (2021) 90:11–27.
23. Fava GA, Guidi J, Semprini F, Tomba E, Sonino N. Clinical assessment of allostatic load and clinimetric criteria. Psychother Psychosom. (2010) 79:280–4.
24. Tyra AT, Soto SM, Young DA, Ginty AT. Frequency and perceptions of life stress are associated with reduced cardiovascular stress-response adaptation. Int J Psychophysiol Off J Int Organ Psychophysiol. (2020) 157:51–60. doi: 10.1016/j.ijpsycho.2020.08.006
25. Smetana P, Malik M. Sex differences in cardiac autonomic regulation and in repolarisation electrocardiography. Pflügers Arch Eur J Physiol. (2013) 465:699–717. doi: 10.1007/s00424-013-1228-x
26. Whited MC, Larkin KT. Sex differences in cardiovascular reactivity: influence of the gender role relevance of social tasks. J Psychophysiol. (2009) 23:77–84. doi: 10.1027/0269-8803.23.2.77
27. Kudielka BM, Buske-Kirschbaum A, Hellhammer DH, Kirschbaum C. Differential heart rate reactivity and recovery after psychosocial stress (TSST) in healthy children, younger adults, and elderly adults: the impact of age and gender. Int J Behav Med. (2004) 11:116–21.
28. Evans BE, Greaves-Lord K, Euser AS, Tulen JHM, Franken IHA, Huizink AC. Alcohol and tobacco use and heart rate reactivity to a psychosocial stressor in an adolescent population. Drug Alcohol Depend. (2012) 126:296–303. doi: 10.1016/j.drugalcdep.2012.05.031
29. Licht CMM, de Geus EJC, Zitman FG, Hoogendijk WJG, van Dyck R, Penninx BWJH. Association Between Major Depressive Disorder and Heart Rate Variability in the Netherlands Study of Depression and Anxiety (NESDA). Arch Gen Psychiatry. (2008) 65:1358–67. doi: 10.1001/archpsyc.65.12.1358
30. Phillips AC. Blunted cardiovascular reactivity relates to depression, obesity, and self-reported health. Biol Psychol. (2011) 86:106–13. doi: 10.1016/j.biopsycho.2010.03.016
31. Burke HM, Davis MC, Otte C, Mohr DC. Depression and cortisol responses to psychological stress: a meta-analysis. Psychoneuroendocrinology. (2005) 30:846–56. doi: 10.1016/j.psyneuen.2005.02.010
32. Cristóbal-Narváez P, Haro JM, Koyanagi A. Perceived stress and depression in 45 low- and middle-income countries. J Affect Disord. (2020) 274:799–805. doi: 10.1016/j.jad.2020.04.020
33. Ginty AT, Conklin SM. High perceived stress in relation to life events is associated with blunted cardiac reactivity. Biol Psychol. (2011) 86:383–5. doi: 10.1016/j.biopsycho.2011.01.002
34. Byun S, Kim AY, Jang EH, Kim S, Choi KW, Yu HY, et al. Detection of major depressive disorder from linear and nonlinear heart rate variability features during mental task protocol. Comput Biol Med. (2019) 112:103381. doi: 10.1016/j.compbiomed.2019.103381
35. Gaetz M, Iverson GL, Rzempoluck EJ, Remick R, McLean P, Linden W. Self-Organizing Neural Network Analyses of Cardiac Data in Depression. Neuropsychobiology. (2004) 49:30–7. doi: 10.1159/000075336
36. Kuang D, Yang R, Chen X, Lao G, Wu F, Huang X, et al. Depression recognition according to heart rate variability using Bayesian Networks. J Psychiatr Res. (2017) 95:282–7. doi: 10.1016/j.jpsychires.2017.09.012
37. Cohen S, Kamarck T, Mermelstein R. Perceived stress scale. In: S Cohen editor. Measuring Stress: A Guide for Health and Social Scientists Kindle Edition. Oxford: Oxford University Press (1994). 10 p.
38. Cohen S. Perceived stress in a probability sample of the United States. In: Spacapan S, Oskamp S editors. The Social Psychology of Health. Thousand Oaks, CA: Sage Publications, Inc (1988). p. 31–67.
39. Raes F, Hermans D, Eelen P. Kort instrumenteel De Nederlandstalige versie van de Ruminative Response Scale (RRS-NL) en de Rumination on Sadness Scale (RSS-NL). Gedragstherapie. (2003) 36:97–104.
40. National Collaborating Centre for Nursing and Supportive Care (UK).The General Practice Physical Activity Questionnaire (GPPAQ), Irritable Bowel Syndrome in Adults: Diagnosis and Management of Irritable Bowel Syndrome in Primary Care [Internet]. London: Royal College of Nursing (2008).
41. de Beurs E, Van Dyck R, Marquenie LA, Lange A, Blonk RWB. De DASS: een vragenlijst voor het meten van depressie, angst en stress. Gedragstherapie. (2001) 34:35–54.
42. Lovibond PF, Lovibond SH. The structure of negative emotional states: comparison of the depression anxiety stress scales (DASS) with the beck depression and anxiety inventories. Behav Res Ther. (1995) 33:335–43. doi: 10.1016/0005-7967(94)00075-u
43. Wardenaar KJ, Wanders RBK, Jeronimus BF, de Jonge P. The Psychometric Properties of an Internet-Administered Version of the Depression Anxiety and Stress Scales (DASS) in a Sample of Dutch Adults. J Psychopathol Behav Assess. (2018) 40:318–33. doi: 10.1007/s10862-017-9626-6
44. Hamilton M. The Hamilton rating scale for depression. In: Norman S, Thomas AB editors. Assessment of Depression. Berlin: Springer (1986). p. 143–52. doi: 10.1007/978-3-642-70486-4_14
46. Trajković G, Starčević V, Latas M, Leštarević M, Ille T, Bukumirić Z, et al. Reliability of the Hamilton Rating Scale for Depression: a meta-analysis over a period of 49years. Psychiatry Res. (2011) 189:1–9. doi: 10.1016/j.psychres.2010.12.007
47. Spruyt A, Clarysse J, Vansteenwegen D, Baeyens F, Hermans D. Affect 4.0: a free software package for implementing psychological and psychophysiological experiments. Exp Psychol. (2009) 57:36–45. doi: 10.1027/1618-3169/a000005
48. Schiweck C, Lutin E, De Raedt W, Cools O, Coppens V, Morrens M, et al. Twenty-Four-Hour Heart Rate Is a Trait but Not State Marker for Depression in a Pilot Randomized Controlled Trial With a Single Infusion of Ketamine. Front Psychiatry. (2021) 12:696170. doi: 10.3389/fpsyt.2021.696170
49. Tarvainen MP, Niskanen J-P, Lipponen JA, Ranta-Aho PO, Karjalainen PA. Kubios HRV—a software for advanced heart rate variability analysis. 4th European Conference of the International Federation for Medical and Biological Engineering. Berlin: Springer (2009). p. 1022–5.
50. Gholamrezaei A, Van Diest I, Aziz Q, Vlaeyen JWS, Van Oudenhove L. Influence of inspiratory threshold load on cardiovascular responses to controlled breathing at 0.1 Hz. Psychophysiology. (2019) 56:e13447.
51. Task Force of the European Society of Cardiology the North American Society of Pacing Electrophysiology. Heart Rate Variability. Circulation. (1996) 93:1043–65. doi: 10.1161/01.CIR.93.5.1043
52. Fontolliet T, Pichot V, Bringard A, Fagoni N, Adami A, Tam E, et al. Testing the vagal withdrawal hypothesis during light exercise under autonomic blockade: a heart rate variability study. J Appl Physiol. (2018) 125:1804–11. doi: 10.1152/japplphysiol.00619.2018
53. Ng J, Sundaram S, Kadish AH, Goldberger JJ. Autonomic effects on the spectral analysis of heart rate variability after exercise. Am J Physiol Circ Physiol. (2009) 297:H1421–8. doi: 10.1152/ajpheart.00217.2009
54. Penttilä J, Kuusela T, Scheinin H. Analysis of rapid heart rate variability in the assessment of anticholinergic drug effects in humans. Eur J Clin Pharmacol. (2005) 61:559–65. doi: 10.1007/s00228-005-0953-2
55. Penttilä J, Helminen A, Jartti T, Kuusela T, Huikuri HV, Tulppo MP, et al. Time domain, geometrical and frequency domain analysis of cardiac vagal outflow: effects of various respiratory patterns. Clin Physiol. (2001) 21:365–76. doi: 10.1046/j.1365-2281.2001.00337.x
56. Cepeda FX, Lapointe M, Tan CO, Taylor JA. Inconsistent relation of nonlinear heart rate variability indices to increasing vagal tone in healthy humans. Auton Neurosci. (2018) 213:1–7. doi: 10.1016/j.autneu.2018.04.007
57. Cho SK, Hwang GS, Kim YK, Huh IY, Hahm KD, Han SM. Low-dose atropine amplifies cardiac vagal modulation and increases dynamic baroreflex function in humans. Aut Neurosci. (2005) 118:108–15.
58. Picard G, Tan CO, Zafonte R, Taylor JA. Incongruous changes in heart period and heart rate variability with vagotonic atropine: implications for rehabilitation medicine. PMR. (2009) 1:820–6. doi: 10.1016/j.pmrj.2009.07.017
59. Shaffer F, Ginsberg JP. An Overview of Heart Rate Variability Metrics and Norms. Front Public Heal. (2017) 5:258. doi: 10.3389/fpubh.2017.00258
60. Alvares GA, Quintana DS, Hickie IB, Guastella AJ. Autonomic nervous system dysfunction in psychiatric disorders and the impact of psychotropic medications: a systematic review and meta-analysis. J Psychiatry Neurosci. (2016) 41:89–104. doi: 10.1503/jpn.140217
61. R Core Team.R: A Language and Environment for Statistical Computing. Vienna: R Core Team (2013).
63. Duschek S, Bair A, Hoffmann A, Marksteiner J, Montoro CI, Reyes del Paso GA. Cardiovascular variability and reactivity in major depressive disorder. J Psychophysiol. (2021) 2021:a000277. doi: 10.1027/0269-8803/a000277
64. Bylsma LM, Taylor-Clift A, Rottenberg J. Emotional reactivity to daily events in major and minor depression. J Abnorm Psychol. (2011) 120:155.
65. Fava M, Rosenbaum JF, McCarthy M, Pava JA, Steingard R, Fox R. Correlations between perceived stress and depressive symptoms among depressive outpatients. Stress Med. (1992) 8:73–6. doi: 10.1002/smi.2460080203
66. Felsten G. Stress reactivity and vulnerability to depressed mood in college students. Pers Individ Dif. (2004) 36:789–800.
67. Hughes BM, Lü W, Howard S. Cardiovascular stress-response adaptation: conceptual basis, empirical findings, and implications for disease processes. Int J Psychophysiol. (2018) 131:4–12. doi: 10.1016/j.ijpsycho.2018.02.003
68. Kelsey RM, Blascovich J, Leitten CL, Schneider TR, Tomaka J, Wiens S. Cardiovascular reactivity and adaptation to recurrent psychological stress: the moderating effects of evaluative observation. Psychophysiology. (2000) 37:748–56. doi: 10.1111/1469-8986.3760748
69. Pieper S, Brosschot JF. Prolonged stress-related cardiovascular activation: is there any? Ann Behav Med. (2005) 30:91–103. doi: 10.1207/s15324796abm3002_1
70. Hamilton JL, Alloy LB. Atypical reactivity of heart rate variability to stress and depression across development: systematic review of the literature and directions for future research. Clincal Psychol Rev. (2016) 50:67–79. doi: 10.1016/j.cpr.2016.09.003
71. McEwen BS. Physiology and neurobiology of stress and adaptation: central role of the brain. Physiol Rev. (2007) 87:873–904. doi: 10.1152/physrev.00041.2006
72. McEwen BS, Stellar E. Stress and the individual: mechanisms leading to disease. Arch Intern Med. (1993) 153:2093–101. doi: 10.1001/archinte.153.18.2093
73. Söder E, Clamor A, Kempkensteffen J, Moritz S, Lincoln TM. Stress levels in psychosis: do body and mind diverge? Biol Psychol. (2018) 138:156–64. doi: 10.1016/j.biopsycho.2018.08.016
74. Zorn JV, Schür RR, Boks MP, Kahn RS, Joëls M, Vinkers CH. Cortisol stress reactivity across psychiatric disorders: a systematic review and meta-analysis. Psychoneuroendocrinology. (2017) 77:25–36. doi: 10.1016/j.psyneuen.2016.11.036
75. Brosschot JF, Gerin W, Thayer JF. The perseverative cognition hypothesis: a review of worry, prolonged stress-related physiological activation, and health. J Psychosom Res. (2006) 60:113–24. doi: 10.1016/j.jpsychores.2005.06.074
76. Grippo AJ, Johnson AK. Stress, depression and cardiovascular dysregulation: a review of neurobiological mechanisms and the integration of research from preclinical disease models. Stress. (2009) 12:1–21. doi: 10.1080/10253890802046281
77. Salomon K, Bylsma LM, White KE, Panaite V, Rottenberg J. Is blunted cardiovascular reactivity in depression mood-state dependent? A comparison of major depressive disorder remitted depression and healthy controls. Int J Psychophysiol. (2013) 90:50–7. doi: 10.1016/j.ijpsycho.2013.05.018
78. Hartmann R, Schmidt FM, Sander C, Hegerl U. Heart Rate Variability as Indicator of Clinical State in Depression. Front Psychiatry. (2019) 9:735. doi: 10.3389/fpsyt.2018.00735
79. Campbell J, Ehlert U. Acute psychosocial stress: does the emotional stress response correspond with physiological responses? Psychoneuroendocrinology. (2012) 37:1111–34. doi: 10.1016/j.psyneuen.2011.12.010
80. Mather M, Ready R, Dixon-Gordon K, Martins B. Differential Coherence In Physiological And Subjective Stress Reactivity In Older And Younger Adults. Innov Aging. (2018) 2:482. doi: 10.1093/geroni/igy023.1799
81. Ordaz S, Luna B. Sex differences in physiological reactivity to acute psychosocial stress in adolescence. Psychoneuroendocrinology. (2012) 37:1135–57. doi: 10.1016/j.psyneuen.2012.01.002
82. Lovallo WR, Farag NH, Sorocco KH, Cohoon AJ, Vincent AS. Lifetime adversity leads to blunted stress axis reactivity: studies from the Oklahoma Family Health Patterns Project. Biol Psychiatry. (2012) 71:344–9. doi: 10.1016/j.biopsych.2011.10.018
Keywords: stress reactivity, heart rate variability, allostatic load, depression, biomarker, repeated stress
Citation: Schiweck C, Gholamrezaei A, Hellyn M, Vaessen T, Vrieze E and Claes S (2022) Exhausted Heart Rate Responses to Repeated Psychological Stress in Women With Major Depressive Disorder. Front. Psychiatry 13:869608. doi: 10.3389/fpsyt.2022.869608
Received: 04 February 2022; Accepted: 29 March 2022;
Published: 18 April 2022.
Edited by:
Maria R. Dauvermann, University of Birmingham, United KingdomReviewed by:
Agorastos Agorastos, Aristotle University of Thessaloniki, GreeceFolajinmi Oluwasina, University of Alberta, Canada
Copyright © 2022 Schiweck, Gholamrezaei, Hellyn, Vaessen, Vrieze and Claes. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Carmen Schiweck, carmen.schiweck@kgu.de