Abstract
Objective
Suicide risk in bipolar disorder (BD) is estimated to be up to 20 times higher than in the general population. While there is a large body of evidence suggesting that increased sympathetic activation is associated with disease and death, there is a paucity of research on the role of autonomic nervous system (ANS) dysfunction in patients with BD who have attempted suicide.
Methods
Fifty-three participants with BD used a wearable device to assess the association between history of suicide attempt, current suicidal ideation, and ANS dysfunction, including measures of heart rate variability (HRV) and respiratory rate. Data were analyzed in a series of unadjusted and adjusted bivariate models of association controlling for relevant variables.
Results
A history of suicide attempts was significantly associated with an increase in respiratory rate (p < 0.01). These results remained significant after adjusting for age, BMI, and current mood state. There was no association between current suicidal ideation and heart rate or respiratory rate. In the frequency domain, HRV parameters suggest reduced parasympathetic (i.e., vagal) activity in participants with a history of suicide attempts and in those with current suicidality, suggesting changes in sympathicovagal balance in BD.
Conclusions
Our results suggest that changes in the ANS in patients with BD and a history of suicide attempt are not restricted to pure vagally mediated HRV parameters, but rather signal a general ANS dysregulation. This ANS imbalance may be contributing to illness burden and cardiovascular disease. Further research on the relationship between ANS and suicidality in BD is needed.
Introduction
While predictive models for suicide have been successfully developed at the level of populations, current demographic and clinical variables are neither sensitive nor specific enough for making individual clinical predictions. A recent meta-analysis concluded that the ability to predict suicide has not improved much over the past 50 years, 1 partially because analyses have only used static (e.g., demographic) variables. This is particularly important in bipolar disorder (BD), a mood disorder that has been associated with premature death due to both cardiovascular disease2–5 and suicide.6–8Strategies to improve the treatment of patients with BD have only partially decreased the incidence of suicide.9,10 This is in part because it is difficult to identify individual patients who are at the highest risk for suicide and to intervene to prevent it. New analytical approaches coupled with new technologies that passively capture densely sampled data may allow us to address this problem, which has been intractable so far.
There is a large body of evidence suggesting that autonomic nervous system (ANS) dysfunction (i.e., typically there is increased activation in the sympathetic system compared to the parasympathetic system) is associated with disease and death.11,12 Heart rate variability (HRV) has been used as an index of autonomic function, and its predictive power for morbidity and mortality in patients with or without cardiovascular disease has been confirmed. 13 Healthy people exhibit a high degree of HRV, reflecting the ability to quickly adapt to the physical or psychological demands of the environment 14 ; conversely, low HRV has been associated with arrhythmias, myocardial infarction, and sudden death. 13 In BD, previous studies have shown that illness burden is associated with lower HRV15–19 and that those BD patients with a family history of suicide present with HRV measures associated with lower parasympathetic activation. 20
HRV is regulated in part by respiratory rate.21,22 Different respiratory patterns have been associated with increased sympathetic (e.g., hyperventilation) or parasympathetic stimulation (e.g., slow, “deep” diaphragmatic breathing). 23 In turn, these different respiratory patterns can modify gas exchange, intracellular cascades,24,25 and brain activity, 26 which can further influence other components of the ANS, including heart rate and HRV, in a feedback loop. Moreover, recent studies have described an association between ANS dysfunction and suicidal ideation,27–34 reinforcing the notion that ANS dysfunction is a maladaptive mechanism that might be compromising the ability of the individual to adequately respond to the demands of the environment. Furthermore, baseline respiratory abnormalities have been described in patients with panic disorder35,36 and BD37,38 and more so in those with both diagnoses.39–41 In parallel, the prevalence of depressive symptoms in patients with chronic obstructive pulmonary disease, asthma, and bronchitis is between 44 and 80%,42–46 with suicide risk almost 4 times higher in patients with 3 or more markers of chronic hypoxia, 47 despite antidepressant treatment. 48 These findings have led some authors to conceptualize the association of pulmonary disease and suicide as a function of hypoxia.49,50
Based on these studies, we conducted a pilot study to assess whether integrating demographic and clinical data with densely sampled physiologic measures of sympathetic and parasympathetic activation (i.e., heart rate, respiratory rate, and HRV) may differentiate BD patients with or without a history of suicide attempt. We hypothesize that patients with BD and a history of a suicide attempt would have higher sympathetic activation (i.e., both higher heart rate and respiratory rate and lower HRV) than patients with BD without a history of a suicide attempt.
Materials and Methods
The study protocol has been previously described 19 and is briefly summarized below.
Participants
The protocol was approved by local research ethics board (REB # 151/2018). After providing informed consent in writing, each participant was interviewed to collect demographic and clinical information. Outpatients were recruited in 2 academic psychiatric hospitals (the Centre for Addiction and Mental Health, Toronto, Canada, and the Royal Ottawa Hospital, Ottawa, Canada) between April 2016 and December 2019. Participants were men or women, 18 years or older, with a primary diagnosis of BD I or II according to the Diagnostic and Statistical Manual of Mental Disorders (DSM) 5 criteria, confirmed with the Schedule for Affective Disorders and Schizophrenia – Lifetime version (SADS-L), 51 in any phase of the illness. Exclusion criteria included active substance abuse or dependence according to DSM 5, mood disorder secondary to a general medical condition, and conditions that may confound HRV measurements including atrial fibrillation or flutter, high-grade heart block, QTc > 500 ms, permanent pacemaker implantation, hemodynamically significant valvular or congenital heart disease, angina, history of acute myocardial infarction, bypass, angioplasty, diabetes mellitus type 1 or 2, or cardiac medications –including β-blockers, calcium channel blockers, and thiazides.
Participants received treatment as usual, according to standard local practices. We classified current psychotropic medications into 4 groups (lithium carbonate, antipsychotics, anticonvulsants, others); 1 participant was not taking any psychotropic medication at the time of measurements.
Psychiatric Measures
We obtained the following clinical information from the SADS-L to determine 12 factors characterizing burden of illness, as described in 19 : (1) lifetime number of depressive episodes; (2) lifetime number of manic, hypomanic, or mixed episodes; (3) lifetime number of all types of episodes; durations (in weeks) of the most severe (4) depressive and (5) hypomanic/manic episode; (6) if applicable, duration (in weeks) of current episode; (7) illness duration calculated by subtracting from the current age their age at onset (defined as the age at which the participant first experienced a depressive, manic, or hypomanic episode, according to DSM 5 criteria); history of (8) suicide attempts or (9) psychotic symptoms; (10) co-morbid anxiety disorders; family history of (11) mood disorders (MDD, BD); or (12) completed suicide. A history of suicide attempt was defined based on the SADS-L section on suicidal behavior, including questions on lifetime history of suicide attempts, intent, and their lethality. Durations of the most severe episodes and of the current episode were transformed logarithmically to better accommodate the assumption of linearity. An index of burden of illness (IBI) was calculated using a weighted geometric mean of the 12 factors characterizing illness burden, which is appropriate when comparing different items, and each item has multiple properties that have numeric ranges. 52 Upon enrolling in the study, clinicians completed the Montgomery–Asberg Depression Rating Scale (MADRS) 53 and the Young Mania Rating Scale (YMRS). 54 Euthymia was defined as both a MADRS score of <7 and YMRS score of <5 for at least 4 weeks.
Data Collection
Cardiovascular and Respiratory Measures
Each participant was provided with a compression shirt fitted with a battery-charged BioModule™ from Zephyr Technology 55 that they wore for 1 period of 24 h. Participants were instructed on how to turn on and off the device and wore it for 24 h, including while sleeping at night. We asked participants not to take a shower during these 24 h, as the device is not waterproof. The BioModule™ is worn against the skin, weighs 35 g, and is attached to the front of a shirt, acting as a data logger. It simultaneously measures 5 physiological and activity-related parameters in real time: respiratory rate, ECG, skin temperature, tri-axial accelerometry, and posture. Data were stored onto the BioModule™ and later transferred to a computer via USB port. The ECG was recorded at 250 Hz and stored onto the BioModule™ and later transferred to a computer via USB port. All ECG data were uploaded into and analyzed using Kubios HRV Premium software. 56
Respiratory rate was detected by a pressure sensor in the strap, which detects torso expansion and contraction due to breathing, and it is reported as the number of breaths per minute. As part of our data preprocessing, a trained rater identified time periods during which the data were very noisy due to movement artifacts. These data segments were not used in the HRV analysis to prevent noisy data from distorting time-varying analysis results. The program we used to analyze the data provided individual streams for each type of data (i.e., ECG data and breathing data), and it has a built-in QRS detector. Once the ECG was imported into the program, the R-wave time instants are automatically detected by applying the built-in QRS algorithm. This detection algorithm is based on the Pan–Tompkins detection algorithm, 57 which consists of applying a series of filters to highlight the frequency content of rapid heart depolarization and removes background noise of any kind. Subsequently, it squares the signal to amplify the QRS. Finally, it applies adaptive thresholds to detect the peaks of the filtered signal. The detector consisted of a preprocessing component, including bandpass filtering of the ECG (to reduce power line noise, baseline wander, and other noise components) and squaring the data samples (to highlight peaks and moving average filtering, thus smoothing close by peaks). The decision rules in this algorithm include amplitude threshold and comparison to expected value between adjacent R-waves. Both these results are adjusted every time a new R-wave is detected. Before R-wave time instant extraction, the R-wave is interpolated at 2,000 Hz to improve the time resolution of the detection. Finally, we used the automatic feature of the “artifact correction algorithm,” in which artifacts are detected from a time series consisting of differences between successive RR intervals. This time series provides a robust way to separate misplaced beats from the normal sinus rhythm. To ensure optimal correction of artifacts, the threshold for correction is estimated from the time-varying distribution of this time series. Detected ectopic beats were corrected by replacing these with interpolated RR values; extra beats were corrected by removing extra R-wave detection and recalculating RR interval series.
Ectopic beats were corrected by linear interpolation with the adjacent complexes; ECG tracings with more than 1% premature beats were eliminated from the analysis. The remaining ECG samples were analyzed with an R-peak detection algorithm and thus produced RR interval (milliseconds) data for each participant, with an accuracy of 98%. The RR data were detrended using the smoothness priors approach (lambda = 500), and an automatic artifact correction algorithm was applied to the RR data, resulting in 1.9% of the beats corrected in the sample.
Based on our prior work,19,20 we analyzed the following HRV measures: (a) 3 HRV measures in the time domain – (i) the root mean square of successive differences (RMSSD), which is associated with short-term, rapid changes in heart rate; (ii) the standard deviation of all the normal RR intervals (SDNN), which is considered the gold standard for medical stratification of cardiac risk when recorded over a 24-h period58,59; (iii) the standard deviation of the average RR intervals (SDANN), which is calculated across 5 min blocks of time, instead of 24-h series – (b) 3 HRV measures in the frequency domain: high frequency (HF), associated with parasympathetic activity; low frequency (LF), associated with sympathetic activity; and LF/HF ratio.
Statistical Analysis
Due to the small sample size and large number of variables, our analysis should be considered hypothesis generating. Thus, we focused on the interpretation of patterns of the bivariate association and their 95% confidence intervals. We calculated odds ratio (OR) between history of suicide attempt (exposure) and HRV and respiratory rate parameters (outcome) using logistic regression and nonparametric tests in SPSS. 60 We conducted unadjusted and adjusted analyses (accounting for age, current mood state, and body mass index (BMI)) in the model. 95% confidence intervals for correlation coefficients were calculated using the inverse hyperbolic tangent transformation of the correlation coefficient, which has approximate normal distribution with variance 1.06/ n – 3, where n is the sample size. 61 We also explored the interaction between HRV and respiratory rate, being mindful of the small number of events and possible instability of these models.
Considering the precision of the estimates, a sample of 50 subjects, the width of a 95% confidence interval would be approximately 0.57 for a Spearman rank correlation coefficient of size 0. This width decreases to around 0.52 for a correlation of 0.30 and to 0.38 for a correlation of 0.60. Thus, if statistical tests for single correlation coefficient are considered, a sample of 50 subjects provides 80% power to detect correlation coefficients of size 0.38 assuming 2 tailed tests, significance level of 0.05, and no adjustment for multiple tests.
Results
Fifty-three participants were enrolled in the study, and none dropped out. Table 1 presents their demographic and clinical characteristics: 32 (60%) were euthymic, 14 (26%) were in a depressive episode, 6 (11%) were in a hypomanic or manic episode, and 1 (2%) was in a mixed episode. Mean (SD) for respiratory rate in the sample was 17.4 (5.1), which is in keeping with normal parameters in adults (12–20 breaths per minute). Almost half of the sample (47%) smoked cigarettes, but none had been diagnosed with a primary pulmonary disease. See Table S1 in the Supplementary Material for details on intent and lethality.
Demographic and Clinical Characteristics of the Sample.
AAP, atypical antipsychotic; BD, bipolar disorder; ECT, electroconvulsive therapy; HRV, heart rate variability; MDD, major depressive disorder; SD, standard deviation
Adjusted and Unadjusted Models for the Association Between History of Suicidal Ideation and Respiratory Rate.
Adjusted and Unadjusted Models for the Association Between Current Suicidal Ideation and Respiratory Rate.
Adjusted and Unadjusted Models for the Interaction Between History of Suicidal Ideation, HRV, and Respiratory Rate.
Adjusted and Unadjusted Models for the Interaction Between Current Suicidal Ideation, HRV, and Respiratory Rate.
HR, heart rate; HRV, heart rate variability; LF, low frequency; HF, high frequency; RMSSD, root mean square of successive differences; SDANN, standard deviation of the average RR intervals; SDNN, standard deviation of all the normal RR intervals.
A history of suicide attempts was significantly associated with a 26% increase in respiratory rate (U = 14.09; p < 0.01). These results remained significant after adjusting for age, current mood state, and BMI. There was no association between current suicidal ideation and heart rate or respiratory rate (p > 0.05).
A history of suicide attempts was also associated with decreased global HRV in all time domain variables (RMSSD, SDANN, and SDNN); none of these variables were significant in the adjusted or unadjusted models (p > 0.05). See Table 2 and Figure 1 for further details. However, in the frequency domain, after adjusting for age, current mood state, and BMI, a history of suicide attempt was associated with an 82% increase in LF/HF ratio (Wald = 4.1; p = 0.04; [1.02, 3.24]; OR = 1.82). Similarly, after controlling for age, BMI and current mood state, participants with current suicidal ideation showed a 67% increase in LF/HF ratio (Wald X2 = 3.84; p <0.05; [1.00, 2.79]; OR = 1.67). See Table 3 for details.

Interaction between HRV, RR, and history of suicidal behavior in 53 patients with BD. Legend: Error bars are the estimates of the difference in association between HRV and RR between those participants with a history of suicide attempt compared to those without. Estimates higher than 0 indicate that the association is higher among those participants with a history of suicide. HR, heart rate; HF, high frequency; LF, low frequency; Min, minimum; Max, maximum; RR, respiratory rate; SDANN, standard deviation of the average RR intervals; SDNN, standard deviation of all the normal RR intervals; SE, standard error; STD: standard.
Furthermore, when we investigated the interaction between HRV and respiratory rate among those with and without history of suicide attempts, we found a significant interaction for LF (t = 4.52; df = 23, p < 0.001), indicating that the association was significantly higher among those with a history of suicide attempts (i.e., an increase of 1 unit in LF was associated with an increase of 8.28 (SE = 1.60; t = 5.21; p < 0.001) units in respiratory rate among those participants with a history of suicide attempt). These results remained significant after adjusting for age, current mood state, and BMI (t = 2.79; df = 20, p = 0.01; 1 unit increase in LF was associated with 8.52 unit increase in respiratory rate (SE = 3.24; t = 2.62; p = 0.01) see Table 4 for details. Conversely, there was no significant difference when comparing participants with current suicidal ideation and those without (Table 5).
While our study was not powered to analyze the potential additive effects of both family history of suicide and a history of suicide attempts on ANS dysfunction, participants with both conditions showed the lowest RMSSD, SD1, and HF values, all indicative of decreased parasympathetic activation (p > 0.05).
Discussion
The aim of this hypothesis-generating study was to examine differences in ANS in patients with BD with and without a history of suicide attempt and with and without current suicidal ideation. In accordance with our previous findings, 20 our results indicate that individuals with BD have lower HRV. Lower RMSSD and SDNN suggest reduced parasympathetic (i.e., vagal) activity and total variability, respectively, in individuals with BD. Higher LF/HF ratio suggests changes in sympathicovagal balance in BD, and increased HR and respiratory rate may be indicative of higher sympathetic activity. Our results suggest that changes in the ANS in individuals with a history of suicide attempt or with current suicidal ideation are not restricted to pure vagally mediated HRV parameters, but rather indicate a general ANS dysregulation.
While our sample size is not large enough to draw conclusions on the association between family history of suicide, history of suicide attempts, and ANS parameters, we consider that our results confirm prior hypotheses related to maladaptive physiological responses to stressors in suicide attempters.34,62 We hypothesize that there is a withdrawal of an already “weakened” parasympathetic branch (e.g., in those with family history of suicide) via decreased RSA, contributing to decreased inhibition of the sympathetic response. This autonomic imbalance is further stimulated by an increased respiratory rate. Our work adds new knowledge to the field by identifying the role of respiratory rate and the potential implications of modulating it to improve mood regulation.
Further research is needed to understand the mechanistic implications of respiratory rate on ANS regulation (other than the importance of RSA 63 ) and the role of hypocapnia (as a result of higher respiratory rate) in BD. Several studies have shown that changes in carbon dioxide levels change the intracellular and extracellular concentrations of H+ very quickly, a change that can only be normalized after 120 h of compensatory (renal) mechanisms. 24 Thus, changes in acidity affect the charge on enzymes, the rate of metabolic reactions, the binding of molecules, as well as pharmacokinetic reactions, pulmonary perfusion, and cardiac activity.64,65 Although a thorough review of acid-base disturbances is outside the scope of this paper, it is important to mention that alkalosis (as a result of hypocapnia) has been associated with sudden death in the context of arrhythmias. 66
Our main limitation is the modest sample size. Other limitations include that this is a cross-sectional analysis and the lack of a control group. While we did not control for changes in metabolic requirements or posture in this study, both of which can affect ANS activity, this is a hypothesis-generating study that sheds light on respiratory rate as an important aspect to be considered in BD. Moreover, although an increase in LF (associated with standing) is often cited as evidence for its association with sympathetic activity, only one-third of healthy subjects show this change, 67 further strengthening our findings. Similarly, while we excluded participants treated with medications that could potentially affect ANS markers, many psychotropics have targets at different receptor levels. Finally, while there is still some controversy between HRV indices and its relationship with the sympathetic and parasympathetic systems, 68 our findings add new knowledge to the field by describing the association between ANS dysfunction and a history of suicide attempt in patients with BD.
Conclusion
Our findings add to the body of knowledge showing that a history of suicide attempt and current suicidal ideation is associated with ANS dysfunction. Together with our previous results on sympathetic activation in BD patients with higher illness burden, our results provide a plausible mechanistic hypothesis to explain the association between increased respiratory rate and further increase in sympathetic stimulation in patients with a history of suicide attempt. Longitudinal studies are needed to understand the role of slow (diaphragmatic) breathing on parasympathetic stimulation and its repercussions on mood regulation in patients with BD.
Supplemental Material
sj-docx-1-cpa-10.1177_07067437231194334 - Supplemental material for A History of Suicide Attempt Is Associated with Increased Sympathetic Activation in Bipolar Disorder
Supplemental material, sj-docx-1-cpa-10.1177_07067437231194334 for A History of Suicide Attempt Is Associated with Increased Sympathetic Activation in Bipolar Disorder by Abigail Ortiz, Yunkyung Park, Stephane MacLean, M. Ishrat Husain, Marcos Sanches, Arun Ravindran and Benoit H. Mulsant in The Canadian Journal of Psychiatry
Footnotes
Acknowledgments
We would like to thank Ms. Pooja Moorti for her help with data preprocessing and harmonization.
Data Availability Statement
The data that support these findings are not publicly available due to ethical restrictions (i.e., containing information that could compromise the privacy of research participants).
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
This work was funded by a grant from the University of Ottawa Department of Psychiatry Research Funds (AO) and was supported in part by Academic Scholar Awards from the Department of Psychiatry, University of Toronto (AO, MIH) and by the Labatt Family Chair in Biology of Depression in Late-Life Adults at the University of Toronto (BHM).
Supplemental Material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
