Skip to main content

Heart rate variability and heart rate asymmetry in adolescents with major depressive disorder during nocturnal sleep period

Abstract

Background

Although reduced heart rate variability (HRV) has been observed in adolescents with major depressive disorder (MDD), substantial between-study heterogeneity and conflicting outcomes exist. Moreover, few studies have investigated heart rate asymmetry (HRA) features despite the high sensitivity of nonlinear indices to heart rate fluctuations. This study aimed to investigate the variations in HRV measures, especially the nonlinear features of HRA, among adolescents with MDD during the nocturnal sleep period.

Methods

Adolescents with MDD and healthy controls completed the clinical assessment of depressive symptom severity and sleep quality followed by a three-night sleep electrocardiogram (ECG) monitoring. Traditional time-domain and frequency-domain HRV measures, nonlinear HRA measures, and the prevalence of different HRA forms and HRA compensation were calculated.

Results

A total of 61 participants with 154 nocturnal ECG time series were available for analysis. Vagally-mediated HRV measures, such as RMSSD, PNN50, and HF, as well as C1d were statistically lower in clinically depressed adolescents compared with healthy controls, whereas C2d was significantly higher. A substantial decrease in the prevalence of short-term HRA, long-term HRA, and the corresponding compensation effect were also observed. In contrast to the medium to large effect sizes observed in traditional HRV indices, nonlinear HRA features showed extremely large effect sizes in discriminating MDD (C1d: Cohen’s d= − 1.38; C2d: Cohen’s d = 1.11), and exhibited a statistical correlation with the severity of depression (C1d: rho = − 0.269; C2d: rho = 0.243). Moreover, there were no significant differences in the distributions of nocturnal HRA measures collected over various nights.

Conclusion

Adolescents with MDD suffered a significant decrease in vagal tone compared to healthy controls, and the features focusing on the directionality of heart rate variations may provide further information on cardiac autonomic activity associated with depression.

Peer Review reports

Introduction

Adolescence is the lifetime peak period for the onset of major depressive disorder (MDD), with an increasing 12-month prevalence of 4–20% worldwide [1, 2]. The burden of depression is considerable since it is associated with impaired school performance, interpersonal difficulties later in life, increased risk of other mental or somatic disorders, and high suicide risk [3,4,5,6]. Over one-third of adolescent-onset depression had recurrence above 18 years of age [7]. However, MDD is often underdiagnosed with only 50% of adolescents diagnosed before reaching adulthood [8]. Adolescence represents a vulnerable age period to a probable depression-induced neuro-dysregulation and cardiac autonomic dysfunction [9,10,11]. Considering the current diagnosis of MDD is primarily dependent on the subjective assessments of clinicians, introducing an appropriate and objective diagnostic marker is critical for early detection and timely treatment.

MDD has been suggested as an independent risk factor for the onset and progression of cardiovascular diseases [12]. The autonomic nervous system plays a major role in the regulation of cardiovascular function [13], and one plausible mechanism for the correlation between MDD and cardiovascular diseases is autonomic dysfunction. Heart rate variability (HRV), a particularly sensitive measure of autonomic activity, characterizes the variations in heartbeats and reflects the dynamic interaction and balance between acceleratory sympathetic nervous system and decelerator parasympathetic nervous system [14]. HRV measurements have time and cost advantages while also being noninvasive and easily applicable. Consequently, HRV has been widely studied to reveal links between mental processes and cardiac autonomic regulation [15].

The methods for HRV analysis can be roughly divided into two categories, including the measures considering the magnitude of variations in beat-to-beat RR intervals and the measures taking into account the directionality of the variations (shortening or prolongation of RR intervals). The former approaches range from simple time-domain methods, such as the standard deviation of RR intervals (SDNN) and the root mean square of successive differences (RMSSD), to frequency-domain methods, such as the high-frequency power component (HF, 0.15–0.4 Hz) and the low-frequency power component (LF, 0.04–0.15 Hz) yielded by analyzing the Fourier spectrum, and to the complexity like symbolic dynamics index [16,17,18]. As suggested by a recent meta-analysis on MDD in adolescents [19], a considerably reduced vagal activity has been identified using HRV measures, but substantial between-study heterogeneity and conflicting outcomes exist as well. Measures with large effect sizes, such as HF and RMSSD, have shown high levels of heterogeneity, whereas the inconsistent directions of changes in the features such as SDNN and LF have been attained across studies [15, 17, 20]. Meanwhile, most studies have focused on short-term resting state HRV lasting for 5 min. Given the instability and unpredictability of short-term data, the researchers using long-term recordings are demanded to verify the reliability and validity of existing results.

Approaches to heart rate asymmetry (HRA), aiming to separately explore the behavior of decelerations and accelerations of heart rhythm, have been attracting increasing attention. The potential mechanisms involved in HRA include respiratory patterns and various physiological reflexes, such as the baroreflex regulating the cardiovascular system [21, 22]. A significant reduction in HRA features has been observed in subjects with attention-deficit/hyperactivity disorder, heart failure, vasovagal syndrome, gastric cancer, or a history of diabetes [23,24,25,26], but the influence of MDD on HRA features in adolescents has only been examined in one study [27]. Additionally, a compensatory phenomenon has been demonstrated in HRA analysis, that in healthy subjects, the contributions of heart rate decelerations are greater than those of accelerations for short-term HRV, whereas the contributions of accelerations are greater than those of decelerations for long-term HRV [28, 29]. However, no studies related to HRA compensation have been performed in adolescent major depression. In view of the high sensitivity of nonlinear indices to heart rate fluctuations in adults with MDD, the performance of HRA-related measures is promising [30]. Moreover, additional research on the correlation between depression severity and measures is required to confirm the utility of autonomous regulation of the heart rate as a possible biomarker, due to the unsatisfactory outcomes of traditional HRV measures [15, 31, 32].

In consideration of the strong link between sleep and depression [33] and the good long-term stability of sleep data not being affected by temporary emotional states, nocturnal electrocardiogram (ECG) signal was chosen as the target data for this study. We hypothesize that in adolescents with MDD, reduced HRV can be consistently identified in long-term nocturnal ECG recordings, and the expression of HRA differs from that of healthy controls with a close connection to depression severity and without altering over time. Accordingly, the objectives of this study were to investigate the following: (i) variations in traditional HRV measures derived from time domain and frequency domain, nonlinear indices of HRA, and the prevalence of HRA phenomena among clinically depressed adolescents; (ii) associations between measures and depressive symptom severity. Furthermore, the consistency of findings in HRV and HRA measures collected over various nights was analyzed.

Methods

Participants

We recruited 28 adolescent patients (13 male, average age: 15.8 ± 1.9 years) with unipolar depressive disorder from the inpatients admitted to the Affiliated Mental Health Center, Zhejiang University School of Medicine, and 33 age-matched healthy subjects (14 male, average age: 15.4 ± 0.5 years) from a high school.

The diagnosis of MDD was classified by thorough clinical investigation based on a structured diagnostic interview by a certified psychiatrist according to the Diagnostic and Statistical Manual of Mental Disorders (DSM-5) [34]. The exclusion criteria for both MDD and control groups were the following: (1) obesity with a body mass index (BMI) over 30; (2) a history of cardiovascular, respiratory, neurological, or metabolic diseases; (3) mental retardation, substance dependence or abuse, bipolar disorder, obsessive-compulsive disorder, post-traumatic stress disorder, and eating disorder defined by DSM-5; (4) presence of psychotic symptoms or suicidal tendencies. The control participants as well as their first-degree relatives had never been treated for any mental disorders. This study was approved by Hangzhou Seventh People’s Hospital Ethics Committee (No. 2024033). All participants and their parents were carefully informed about the study protocol and written informed consent was obtained prior to the examination.

Clinical measures

Trained clinical raters with appropriate inter-rater reliability (Spearman’s rho = 0.83) conducted a semi-structured interview with all adolescent participants [35] and completed the Hamilton Depression Rating Scale (HDRS-17) to determine the severity of depressive symptoms [36]. Participants also completed a self-administrated sleep assessment using the Pittsburgh Sleep Quality Index (PSQI) [37], which estimated total sleep time, sleep quality, sleep medication use, sleep disturbance, daytime dysfunction, sleep latency, and symptoms affecting sleep over the last month. These clinical assessments took place before the measurements of ECG signals.

Sleep ECG recordings and preprocessing

A three-night sleep ECG monitoring in the ward or at home was taken for each participant. Participants were instructed not to smoke or take caffeine-containing soft drinks for at least 12 h before the examination. Data collection was performed overnight in the supine position using CP-025 (Beneware, Hangzhou, China), a single-lead wearable ECG recorder with a sampling rate of 200 Hz. The participants received written guidelines for using the devices.

The raw data of nocturnal ECG were transferred offline to a PC for data storage, on-screen monitoring, and preliminary automatic identification of R-peaks with the use of HolterSystem software (Beneware, Hangzhou, China). Necessary corrections of misidentified R-peaks caused by artifacts and supplements of missing R-peaks were manually operated following a visual inspection of all ECG signals, with the assistance of tools provided by the software. As a result of the significant differences in sleep durations among participants with a median recording period of 7.7 h (IQR 7.3–9.2 h), a specific and continuous ECG recording, which starts at least 2 h after recording and lasts for 4 h, was selected for each night and the corresponding RR sequences were extracted for further HRV analysis using Python (Python Foundation, USA). The overnight recordings with large data loss or without continuous 4-h data were excluded from the analysis.

Heart rate variability analysis

Nocturnal RR recordings were analyzed using traditional HRV measures derived from the time domain and frequency domain, as well as nonlinear indices of HRA.

Traditional time domain and frequency domain analysis

Concerning the time domain analysis, the mean value of RR intervals (RR mean, ms), SDNN (ms), RMSSD (ms), and the percentage of the number of successive heartbeats differing > 50 ms (PNN50) (%) were evaluated.

Concerning the frequency domain analysis, the time series of RR intervals were resampled using cubic spline interpolation with the frequency of 4 Hz, and features of HF, LF, total power (TP), as well as the ratio of LF to HF (LF/HF) were consequently obtained using spectral analysis of fast Fourier transform.

Considering that short RR recordings of a 5-minute duration can be advantageous to examine autonomic modulation for traditional HRV measures with relatively clear physiological implications [38], the 4-hour long RR intervals of each night were divided into consecutive 5-minute segments without overlapping. Each night’s profile of RR intervals and HRV measures were determined by averaging the individual values of all these segments [39].

Heart rate asymmetry analysis

The Poincaré plot of RR intervals, a scatter plot of the current RR interval plotted against the preceding RR interval, was applied to describe HRA illustrating the behavior of decelerations and accelerations of heart rhythm. Short-term and long-term HRAs of the 4-hour RR intervals were quantified using the contributions of heart rate decelerations to short-term and long-term HRVs described by C1d and C2d [28, 40], respectively. C1d and C2d were defined as follows:

$$ {\rm{C}}{{\rm{1}}_{\rm{d}}}{\rm{ = SD}}{{\rm{1}}_{\rm{d}}}^{\rm{2}}{\rm{/SD}}{{\rm{1}}^{\rm{2}}}\,{\rm{and}}\,{\rm{C}}{{\rm{2}}_{\rm{d}}}{\rm{ = }}\,{\rm{SD}}{{\rm{2}}_{\rm{d}}}^{\rm{2}}{\rm{/SD}}{{\rm{2}}^{\rm{2}}}{\rm{,}} $$

SD12 and SD22 represent the variance of the dispersion of points in the Poincaré plot of RR intervals across the identity line and along the identity line, respectively, whereas SDid2 (i = 1, 2) represents the part of SDi2 (i = 1, 2) related to heart rate decelerations.

The phenomenon of short-term HRA is present when C1d > 0.5, and long-term HRA exists when C2d < 0.5. Heart rate asymmetry compensation occurs when both conditions are satisfied, that is, C1d > 0.5 and C2d < 0.5.

Statistical analysis

The group differences between participants with and without MDD were compared using the Mann-Whitney U test in consideration of their non-normal distribution. Specially, the differences in the prevalence of different HRA forms, HRA compensation, and sex between the two groups were examined with Fisher’s exact test. Correlation analysis was subsequently employed to investigate the associations between measures and depression severity measured by HDRS scores, using Spearman’s rank correlation test. Besides, the Wilcoxon signed-rank test was also implemented to contrast the distributions of nocturnal HRV measures collected over various nights to examine the stability of measures.

The effect sizes between the MDD and the control group were calculated by Cohen’s d. Effect sizes of 0.2, 0.5, and 0.8 were categorized as having minimal, moderate, and substantial effects, respectively [41]. All statistical analyses were performed using SPSS 25.0 (SPSS Inc., Chicago, USA) and STATA 17.0 (StataCorp, Texas, USA). All tests were two-tailed with statistical significance defined as P < 0.05.

Results

Group characteristics

A total of 61 participants with 154 nocturnal RR intervals time series were available for analysis in the study. Patients with MDD and healthy controls did not differ in sex, age, weight, height, or BMI. The patient group showed moderate to severe depression (mean HDRS score = 24.57) and poor sleep quality (mean PSQI score = 11.00). Group characteristics are provided in Table 1.

Table 1 Demographic characteristics of the MDD and control groups

Between-group comparison of HRV/HRA measures

Between-group tests on measures of heart rate and automatic nervous system activity are summarized in Table 2. HRV measures (RMSSD, PNN50, HF, LF, TP) and C1d were statistically lower in clinically depressed adolescents compared with healthy controls, whereas C2d was significantly higher. C1d yielded the largest effect size (Cohen’s d = − 1.38), followed by C2d, TP, HF, LF, PNN50, and RMSSD with effect sizes of 1.11, − 0.72, − 0.69, − 0.63, − 0.60, and − 0.17, respectively. No significant differences were observed between the MDD group and control group in terms of mean RR interval, SDNN, and LF/HF.

Table 2 Comparison between the patient group and control group

Prevalence of HRA phenomena

The prevalence of different HRA forms and HRA compensation for the two groups is presented in Fig. 1. In adolescents with MDD, short-term asymmetry was revealed in 29.5% of recordings (23/78), long-term asymmetry in 64.1% (50/78), and both types of HRA coexisted in 28.2% (22/78). The control group had significantly higher prevalence rates: 85.5% for short-term asymmetry (65/76 recordings, P < 0.001), 94.7% for long-term asymmetry (72/76 recordings, P < 0.001), and 85.5% for HRA compensation (65/76 recordings, P < 0.001). All recordings from healthy participants with short-term asymmetry also exhibited HRA compensation. Figure 2 shows two typical Poincaré plots of RR intervals from clinically depressed and healthy participants. The distribution of RR intervals in healthy subjects was more dispersed than that in patients.

Fig. 1
figure 1

Prevalence of different types of HRA and comparison between the patient (n = 78) and control (n = 76) groups. *** P < 0.001. HRA, heart rate asymmetry

Fig. 2
figure 2

Two typical Poincaré plots of RR intervals attained from healthy (left panel) and clinically depressed (right panel) participants with the scatter density plots in the plane (RR(i), RR(i + 1)) (A and B, respectively) and the corresponding distributions of ΔRR = RR(i + 1)– RR(i) (C and D, respectively)

Correlation between HRV/HRA measures and depression severity

HRV measures of RMSSD, PNN50, HF, LF, and TP, and HRA measures of C1d and C2d, which can differentiate patients from healthy adolescents, were further investigated concerning the relationship with depressive symptoms. Correlations using data from clinically depressed adolescents revealed a significant negative correlation of HDRS with RMSSD (rho = − 0.255, P = 0.024), PNN50 (rho = − 0.231, P = 0.041), HF (rho = − 0.236, P = 0.037), TP (rho = − 0.245, P = 0.032), and C1d (rho = − 0.269, P = 0.017), and a positive correlation of HDRS scores with C2d (rho = 0.243, P = 0.032). However, no significant interaction was found between LF and depression severity.

Impact of data collection time

The distributions of nocturnal HRV and HRA measures collected over various nights were contrasted (Table 3). A total of 51 participants completed two nights of high-quality ECG recordings, and 40 subjects achieved three nights. The measures of mean RR intervals, HF, LF, and TP were highly variable and sensitive to the timing of data collection.

Table 3 Comparison of the distributions of HRV and HRA measures collected over various nights

Discussion

In this study, autonomic activity in adolescents with MDD was evaluated using traditional HRV and nonlinear HRA analysis. Firstly, adolescents suffering from MDD showed a reduction in traditional time-domain and frequency-domain HRV measures using 4-hour nocturnal ECG recordings. Secondly, the phenomena of both short-term and long-term HRA in depressed adolescents as well as the corresponding compensation effect were significantly diminished. Thirdly, in contrast to the medium to large effect sizes yielded in traditional HRV indices, nonlinear HRA features revealed a strong ability to distinguish between MDD and healthy adolescents with extremely large effect sizes, showing statistical correlations with depression severity. Notably, HRA measures remained consistent across multiple nights of monitoring.

Comparison of HRV measures among clinically depressed adolescents

Instead of using single 5-minute RR intervals, the traditional HRV measures derived from time domain and frequency domain in this study were calculated using multiple consecutive 5-minute RR segments. Most findings aligned with previous studies employing short-term resting state ECG recordings. All vagally-mediated HRV measures including RMSSD, PNN50, and HF were statistically reduced and correlated with depression severity. There were no substantial associations found between MDD and LF. The recent meta-review showed no differences in LF with a high level of heterogeneity between clinically depressed and healthy adolescents [19], as well as the present study demonstrated little interaction between LF and depressive symptom severity. No alterations occurred in SDNN either, which are modulated by both sympathetic and parasympathetic activity and highly correlated with LF [16, 42]. Overall, this confirms the hypothesis proposed in the previous review that autonomic nervous system dysfunction in adolescents with MDD is predominantly driven by reduced vagal activity, rather than increased sympathetic activity or joint sympathetic and parasympathetic functioning [19, 32].

The neurovisceral integration model and the polyvagal theory have supported the idea of linking the cardiac vagal control indexed by HRV and psychological processes [43, 44]. The former model connecting the heart to the prefrontal cortex assumes that the higher the vagal tone, the better executive cognitive performance, as well as better emotional and health regulation, whereas the latter based on the phylogeny of the neural circuits emphasizes the impact of vagal activity on social functioning. Accordingly, the decreased vagally-mediated HRV measures in this study were matched to the clinical manifestations of MDD in adolescents, including unstable peer and family relationships, poor school performance, and diminished emotion differentiation [45,46,47,48].

With regard to nonlinear HRA measures, a considerable decrease appeared in the contributions of heartbeat decelerations to short-term HRVs (C1d) and an increase to long-term HRVs (C2d) in adolescents with MDD. This implied that physiological processes make an effort to maintain a balance between short-term and long-term heart rate decelerations in response to depression, but insufficiently. Compared to the medium to large effect sizes observed in linear HRV measures, both C1d and C2d showed extremely large effect sizes in identifying depression from healthy controls with a significant connection to the severity of depression. The dynamic nonlinear indices regarding heart rate deceleration or acceleration capacities and patterns therefore may play a major role in the diagnosis and assessment of depression.

In addition, this study inspected the stability of HRV measures. Large variations with data collecting time were revealed in heart rate and spectrum features of LF, HF, and TP. This may indicate that time-domain HRV and nonlinear HRA features are more adept at capturing the depression-induced disturbances that exist in a stable manner in physiological signals, and longer lengths of RR intervals are demanded for frequency-domain features to minimize the influence of temporary autonomic nervous system activity. Interestingly, the instability of short-term HF may also contribute to the high level of between-study heterogeneity in meta-review since no trustworthy sources of heterogeneity were discovered [19]. Nevertheless, Vloet et al. [49] did not find any correlations between depression and HRV in adolescents using 24-hour ECG during everyday activity. In view of the fact that adolescents have higher cardiovascular autonomic reactivity than adults [50, 51], the data involving a relatively steady psychological and physiological state is recommended to avoid depression-related alterations being obscured by transient but intense autonomic activations.

Prevalence of HRA phenomena in clinically depressed adolescents

The prevalence of short-term and long-term asymmetry as well as its compensation were all statistically decreased in clinically depressed adolescents, with the decline in short-term asymmetry and HRA compensation being particularly notable. As discovered by Sibrecht et al. [52] and Zalas et al. [29], the prevalence of various HRA features exceeds 85% in healthy volunteers including children and adults. Depression-induced heavy loss in short-term asymmetry (56%) and HRA compensation (57.3%) offered important insights into the underlying physiological mechanism.

Although the real physiological origin is still uncertain, all factors regulating momentary heart rate may cause HRA, involving the factors naturally asymmetric such as the difference of delay in the sinus node’s response to vagal tone changes and the sympathetic response, and the factors contributing to HRA such as hormones, cytokines, and emotional activity. From the perspective of cardiac autonomic regulation, some researchers have ascribed the decelerations-based part of variance and short-term HRA to the parasympathetic branch [28, 53], which was consistent with the findings in traditional HRV measures. Correspondingly, HRA compensation can be recognized as an effect of the antagonistic function of sympathetic and parasympathetic systems. However, the pathophysiological mechanism of depression-related alterations in HRA remains not fully understood. The evidence like the dysregulation of the hypothalamus-pituitary-adrenal axis, the depletion of brain serotonin, and the immune system response observed in depression might also be considered.

Limitations

Some limitations should be mentioned in this study. Firstly, several cardiovascular risk factors such as physical activity and inflammation were not carefully monitored. We just put limitations on smoking and caffeine-containing beverages before the data collection. Secondly, the lack of detailed labeling of sleep process, including wakefulness, sleep stages, and movement time makes it impossible to build a relationship among sleep, HRV, and depression. Thirdly, HRA measures were not thoroughly investigated. Although the scatter plots and the corresponding distributions of RR intervals showed a great difference between healthy and clinically depressed participants, only features of C1d and C2d and the related HRA phenomena were analyzed in this study. Furthermore, HRV analysis could benefit from methodological improvements, such as exploring algorithms like the Lomb-Scargle Periodogram, which could provide more accurate power spectral density estimation without the need for regular RR intervals or resampling.

Implications for future research

Regarding technical considerations, more novel asymmetry measures using variance-, symbol-, graph-based, and entropy approaches are needed to promote the accuracy of depression monitoring and advance our knowledge of cardiac nerve electrophysiology. Meanwhile, the application of machine learning algorithms using HRV features, especially HRA measures, as input is anticipated to enhance the efficacy of depression assessment and be utilized in clinical practice for diagnostic support.

From a clinical perspective, reductions in vagal activity have been verified, but the exploration of the impact of antidepressant treatment on HRV measures is demanded to validate the capacity of HRV as a psychophysiological biomarker for prognosis and prediction of therapy response in depression. Moreover, as understanding of the pathogenesis of depression deepens, the therapy affecting autonomic activity may become an option for addressing depression among adolescents.

Conclusion

Adolescents with MDD suffered a considerable decrease in vagal tone, along with reductions in time-domain and frequency-domain HRV measures and a heavy loss of short-term HRA and HRA compensation, compared to healthy controls. In contrast to the medium to large effect sizes observed in traditional HRV indices and their high sensitivity to the timing of data collection, nonlinear HRA features showed extremely large effect sizes in discriminating MDD with stable performance, and they also had a significant correlation with the severity of depression. The features focusing on the directionality of heart rate variations may provide further information on cardiac autonomic activity associated with depression. These findings contribute to better assessing the usefulness of HRV in adolescent depression.

Data availability

The datasets used and analyzed during the current study are available from the corresponding author (Shulin Chen) on reasonable request.

Abbreviations

BMI:

Body mass index

C1d :

Contributions of heart rate decelerations to short-term heart rate variability

C2d :

Contributions of heart rate decelerations to long-term heart rate variability

CI:

Confidence intervals

ECG:

Electrocardiogram

HDRS:

Hamilton Depression Rating Scale

HF:

High-frequency power

HRA:

Heart rate asymmetry

HRV:

Heart rate variability

LF/LF:

Ratio of low frequency to high frequency

LF:

Low-frequency power

MDD:

Major depressive disorder

PNN50:

Percentage of the number of successive heartbeats differing > 50 ms

PSQI:

The Pittsburgh Sleep Quality Index

RMSSD:

Root mean square of successive differences

SDNN:

Standard deviation of RR intervals

TP:

Total power

References

  1. Institute of Psychology, Chinese Academy of Sciences. The report on China’s national mental health development (2021–2022) (in Chinese). 2023. https://www.pishu.com.cn/skwx_ps/databasedetail?SiteID=14%26contentId=14414530%26contentType=literature%26type=%25E6%258A%25A5%25E5%2591%258A%26subLibID=.

  2. Substance Abuse and Mental Health Services Administration. Key substance use and mental health indicators in the United States: Results from the 2021 National Survey on Drug Use and Health (HHS Publication No. PEP22-07-01-005, NSDUH Series H-57). Center for Behavioral Health Statistics and Quality, Substance Abuse and Mental Health Services Administration., 2022. https://www.samhsa.gov/data/report/2021-nsduh-annual-national-report

  3. Boersma-van Dam E, Hale B, Koot H, Meeus W, Branje S. Adolescents’ and best friend’s depressive symptoms and conflict management: intraindividual and interpersonal processes over time. J Clin Child Adolesc Psychol. 2019;48(2):203–17.

    Article  PubMed  Google Scholar 

  4. Lee DS, Cederbaum JA, Davis JP, Hurlburt MS, Mennen FE. Maternal and adolescent depressive symptoms and family conflict: an autoregressive cross-lagged examination of competing models in multi-stressed mothers and adolescents. Fam Process. 2023;62(1):254–71.

    Article  PubMed  Google Scholar 

  5. Serra G, De Crescenzo F, Maisto F, Galante JR, Iannoni ME, Trasolini M, et al. Suicidal behavior in juvenile bipolar disorder and major depressive disorder patients: systematic review and meta-analysis. J Affect Disord. 2022;311:572–81.

    Article  PubMed  Google Scholar 

  6. Steffen A, Nuebel J, Jacobi F, Baetzing J, Holstiege J. Mental and somatic comorbidity of depression: a comprehensive cross-sectional analysis of 202 diagnosis groups using German nationwide ambulatory claims data. BMC Psychiatry. 2020;20(1):142.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Kiss E, Baji I, Kellner A, Mayer L, Kapornai K. Long-term follow-up of childhood-onset depression - comorbidity, suicidal behavior and prognosis in adulthood. Psychiatr Hung. 2020;35(1):58–67.

    PubMed  Google Scholar 

  8. Grp GPS, Zuckerbrot RA, Cheung A, Jensen PS, Stein REK, Laraque D. Guidelines for adolescent depression in primary care (GLAD-PC): part I. practice preparation, identification, assessment, and initial management. Pediatrics. 2018;141(3):e20174081.

    Article  Google Scholar 

  9. Thayer JF, Sollers JJ III, Labiner D, Weinand M, Herring AM, Lane RD, et al. Age-related differences in prefrontal control of heart rate in humans: a pharmacological blockade study. Int J Psychophysiol. 2009;72(1):81–8.

    Article  PubMed  Google Scholar 

  10. Tonhajzerova I, Ondrejka I, Javorka K, Turianikova Z, Farsky I, Javorka M. Cardiac autonomic regulation is impaired in girls with major depression. Prog Neuropsychopharmacol Biol Psychiatry. 2010;34(4):613–8.

    Article  PubMed  Google Scholar 

  11. Yang TT, Simmons AN, Matthews SC, Taperta SF, Bischoff-Grethe A, Frank GKW, et al. Increased amygdala activation is related to heart rate during emotion processing in adolescent subjects. Neurosci Lett. 2007;428(2–3):109–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Zambrano J, Celano CM, Januzzi JL, Massey CN, Chung WJ, Millstein RA, et al. Psychiatric and psychological interventions for depression in patients with heart disease: a scoping review. J Am Heart Assoc. 2020;9(22):e018686.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Malpas SC. Sympathetic nervous system overactivity and its role in the development of cardiovascular disease. Physiol Rev. 2010;90(2):513–57.

    Article  CAS  PubMed  Google Scholar 

  14. Thayer JF, Yamamoto SS, Brosschot JF. The relationship of autonomic imbalance, heart rate variability and cardiovascular disease risk factors. Int J Cardiol. 2010;141(2):122–31.

    Article  PubMed  Google Scholar 

  15. Park S, Lee J, Kim J, Suh S, Lee MS. Changes in heart rate variability in first-episode drug-naive adolescents with major depressive disorder: a 12-week prospective study. J Affect Disord. 2018;238:250–5.

    Article  PubMed  Google Scholar 

  16. Shaffer F, Ginsberg JP. An overview of heart rate variability metrics and norms. Front Public Health. 2017;5:258.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Blom EH, Serlachius E, Chesney MA, Olsson EMG. Adolescent girls with emotional disorders have a lower end-tidal CO2 and increased respiratory rate compared with healthy controls. Psychophysiology. 2014;51(5):412–8.

    Article  Google Scholar 

  18. Mestanikova A, Mestanik M, Ondrejka I, Hrtanek I, Cesnekova D, Jurko A, et al. Complex cardiac vagal regulation to mental and physiological stress in adolescent major depression. J Affect Disord. 2019;249:234–41.

    Article  PubMed  Google Scholar 

  19. Chen W, Zhong Q, Chen H, Chen S. Heart rate variability in children and adolescents with major depressive disorder: a systematic review and meta-analysis. J Affect Disord. 2023;335:204–15.

    Article  PubMed  Google Scholar 

  20. Tonhajzerova I, Visnovcova Z, Ondrejka I, Funakova D, Hrtanek I, Ferencova N. Major depressive disorder at adolescent age is associated with impaired cardiovascular autonomic regulation and vasculature functioning. Int J Psychophysiol. 2022;181:14–22.

    Article  PubMed  Google Scholar 

  21. De Maria B, Bari V, Cairo B, Vaini E, de Abreu RM, Perseguini NM, et al. Cardiac baroreflex hysteresis is one of the determinants of the heart period variability asymmetry. Am J Physiol Regul Integr Comp Physiol. 2019;317(4):R539–551.

    Article  PubMed  Google Scholar 

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

    Article  PubMed  Google Scholar 

  23. Guzik P, Piskorski J, Contreras P, Migliaro ER. Asymmetrical properties of heart rate variability in type 1 diabetes. Clin Auton Res. 2010;20(4):255–7.

    Article  PubMed  Google Scholar 

  24. Pawlowski R, Zalewski P, Newton J, Piatkowska A, Kozluk E, Opolski G, et al. An assessment of heart rate and blood pressure asymmetry in the diagnosis of vasovagal syncope in females. Front Physiol. 2023;13:1087837.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Shi B, Wang L, Yan C, Chen D, Liu M, Li P. Nonlinear heart rate variability biomarkers for gastric cancer severity: a pilot study. Sci Rep. 2019;9(1):13833.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Tonhajzerova I, Ondrejka I, Farsky I, Visnovcova Z, Mestanik M, Javorka M, et al. Attention deficit/hyperactivity disorder (ADHD) is associated with altered heart rate asymmetry. Physiol Res. 2019;63:S509–19.

    Google Scholar 

  27. Tonhajzerova I, Ondrejka I, Chladekova L, Farsky I, Visnovcova Z, Calkovska A, et al. Heart rate time irreversibility is impaired in adolescent major depression. Prog Neuropsychopharmacol Biol Psychiatry. 2012;39(1):212–7.

    Article  PubMed  Google Scholar 

  28. Piskorski J, Guzik P. Compensatory properties of heart rate asymmetry. J Electrocardiol. 2012;45(3):220–4.

    Article  PubMed  Google Scholar 

  29. Zalas D, Bobkowski W, Piskorski J, Guzik P. Heart rate asymmetry in healthy children. J Clin Med. 2023;12(3).

  30. Kemp AH, Quintana DS, Gray MA, Felmingham KL, Brown K, Gatt JM. Impact of depression and antidepressant treatment on heart rate variability: a review and meta-analysis. Biol Psychiatry. 2010;67(11):1067–74.

    Article  CAS  PubMed  Google Scholar 

  31. Koenig J, Schreiner MW, Klimes-Dougan B, Ubani B, Mueller B, Kaess M, Cullen KR. Brain structural thickness and resting state autonomic function in adolescents with major depression. Soc Cogn Affect Neurosci. 2018;13(7):741–53.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Moretta T, Kaess M, Koenig J. A comparative evaluation of resting state proxies of sympathetic and parasympathetic nervous system activity in adolescent major depression. J Neural Transm. 2023;130(2):135–44.

    Article  CAS  PubMed  Google Scholar 

  33. Riemann D, Krone LB, Wulff K, Nissen C. Sleep, insomnia, and depression. Neuropsychopharmacology. 2020;45(1):74–89.

    Article  PubMed  Google Scholar 

  34. American Psychiatric Association. Diagnostic and statistical manual of mental disorders. 5th ed. Arlington: American Psychiatric Publishing; 2013.

    Book  Google Scholar 

  35. Williams JB. A structured interview guide for the Hamilton depression rating scale. Arch Gen Psychiatry. 1988;45(8):742–7.

    Article  CAS  PubMed  Google Scholar 

  36. Hamilton M. A rating scale for depression. J Neurol Neurosurg Psychiatry. 1960;23(1):56–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Buysse DJ, Reynolds CF, Monk TH, Berman SR, Kupfer DJ. The Pittsburgh sleep quality index: a new instrument for psychiatry practice and research. Psychiatry Res. 1989;28(2):193–213.

    Article  CAS  PubMed  Google Scholar 

  38. Camm AJ, Malik M, Bigger JT, Breithardt G, Cerutti S, Cohen RJ, et al. Heart rate variability. Standards of measurement, physiological interpretation, and clinical use. Eur Heart J. 1996;17(3):354–81.

    Article  Google Scholar 

  39. Li K, Ruediger H, Ziemssen T. Spectral analysis of heart rate variability: time window matters. Front Neurol. 2019;10:545.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Piskorski J, Ellert J, Krauze T, Grabowski W, Wykretowicz A, Guzik P. Testing heart rate asymmetry in long, nonstationary 24 hour RR-interval time series. Physiol Meas. 2019;40(10):105001.

    Article  CAS  PubMed  Google Scholar 

  41. Cohen J. Statistical power analysis for the behavioral sciences. 2nd ed. Hillsdale, NJ: Lawrence Erlbaum Associates, 1988.

    Google Scholar 

  42. Otzenberger H, Gronfier C, Simon C, Charloux A, Ehrhart J, Piquard F, et al. Dynamic heart rate variability: a tool for exploring sympathovagal balance continuously during sleep in men. Am J Physiol. 1998;275(3):H946–50.

    CAS  PubMed  Google Scholar 

  43. Porges SW. The polyvagal perspective. Biol Psychol. 2007;74(2):116–43.

    Article  PubMed  Google Scholar 

  44. Thayer JF, Hansen AL, Saus-Rose E, Johnsen BH. Heart rate variability, prefrontal neural function, and cognitive performance: the neurovisceral integration perspective on self-regulation, adaptation, and health. Ann Behav Med. 2009;37(2):141–53.

    Article  PubMed  Google Scholar 

  45. Ozturk Y, Onat M, Ozyurt G, Mutlu C, Tufan AE, Akay AP. Relationships between family functioning, parenting and peer victimization in adolescent depression: a cross-sectional study. North Clin Istanb. 2021;8(3):212–21.

    PubMed  PubMed Central  Google Scholar 

  46. Maalouf FT, Brent D, Clark L, Tavitian L, McHugh RM, Sahakian BJ, et al. Neurocognitive impairment in adolescent major depressive disorder: state vs. trait illness markers. J Affect Disord. 2011;133(3):625–32.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Fröjd SA, Nissinen ES, Pelkonen MU, Marttunen MJ, Koivisto AM, Kaltiala-Heino R. Depression and school performance in middle adolescent boys and girls. J Adolesc. 2008;31(4):485–98.

    Article  PubMed  Google Scholar 

  48. Nook EC. Emotion differentiation and youth mental health: current understanding and open questions. Front Psychol. 2021;12:700298.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Vloet TD, Jans T, Frey A, Haeussler M, Vloet A, Geissler J, et al. Mean heart rate and parameters of heart rate variability in depressive children and the effects of antidepressant medication. Z Kinder Jugendpsychiatr Psychother. 2019;47(3):253–60.

    Article  PubMed  Google Scholar 

  50. Antelmi I, De Paula RS, Shinzato AR, Peres CA, Mansur AJ, Grupi CJ. Influence of age, gender, body mass index, and functional capacity on heart rate variability in a cohort of subjects without heart disease. Am J Cardiol. 2004;93(3):381–5.

    Article  PubMed  Google Scholar 

  51. Yeragani VK, Sobolewski E, Kay J, Jampala VC, Igel G. Effect of age on long-term heart rate variability. Cardiovasc Res. 1997;35(1):35–42.

    Article  CAS  PubMed  Google Scholar 

  52. Sibrecht G, Piskorski J, Krauze T, Guzik P. Heart rate asymmetry, its compensation, and heart rate variability in healthy adults during 48-h Holter ECG recordings. J Clin Med. 2023;12(3):1219.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Karmakar C, Khandoker A, Palaniswami M. Investigating the changes in heart rate asymmetry (HRA) with perturbation of parasympathetic nervous system. Australas Phys Eng Sci Med. 2012;35(4):465–74.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We are grateful to the adolescents, parents, and teachers who provided data for this study and Hangzhou Beneware Medical Equipment Co., Ltd for providing HolterSystem software free of charge.

Funding

This work was supported by Hangzhou Municipal General Medical and Health Plan (Grant No. A20240472).

Author information

Authors and Affiliations

Authors

Contributions

Hang Chen, Zhenghe Yu, and Shulin Chen conceived and designed the study; Wanlin Chen and Haisi Chen interviewed participants and collected the data; Wanlin Chen, Wenchen Jiang, Cheng Chen, Moya Xu, and Haoxuan Ruan analyzed the data; Wanlin Chen, Haisi Chen, and Wenchen Jiang wrote the paper. All authors contributed to and have approved the final manuscript.

Corresponding author

Correspondence to Shulin Chen.

Ethics declarations

Ethics approval and consent to participate

This study has been approved by Hangzhou Seventh People’s Hospital Ethics Committee (No. 2024033). We confirmed that written informed consent was obtained from all participants and their parents. All methods in this study were performed by the guidelines and regulations.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, W., Chen, H., Jiang, W. et al. Heart rate variability and heart rate asymmetry in adolescents with major depressive disorder during nocturnal sleep period. BMC Psychiatry 25, 497 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12888-025-06911-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12888-025-06911-3

Keywords