Heart rate variability as biomarker for bipolar disorder

Bipolar disorder (BD) is characterized by alterations in mood, energy levels and the ability to function. Accordingly, it is also associated with dysfunction of the autonomic nervous system (ANS), indexed by heart rate variability (HRV). Literature concerning differences in ANS functioning between mood states is still sparse. The main aim of the study was to investigate within-individual changes in HRV from manic to euthymic states in bipolar disorder (BD). This is the first study to do so using wrist-worn sensors. Seventeen patients with BD were equipped with photoplethysmography (PPG) sensor wristbands and provided 24-hour recordings both during a manic state and a euthymic state. We calculated mean heart rate and the commonly used HRV measures SDNN, RMSSD, HF, LF and Sample Entropy in 5-minute segments during rest at night. We compared HRV by mood state within individuals using paired t-tests and linear regression to control for age and sex. Recordings from 15 BD patients were analyzed. There were statistically significant increases in HRV measures SDNN, RMSSD, LF and Sample Entropy from mania to euthymia. Effect sizes were predominately large. Our findings reveal lower HRV in the manic state compared to the euthymic state. This indicates that HRV collected by wrist-worn PPG sensors is a possible biomarker for bipolar mood states. Movement artifacts were problematic and sampling during rest or in combination with actigraphy is recommended. Our findings can be further implemented to develop a monitoring device for bipolar patients.


42
Bipolar disorder (BD) is a severe mental disorder associated with intense mood fluctuations, 43 lifelong disability and increased mortality. It affects about 1-3% of the population, with a 44 lifetime suicide prevalence of 15-20% [1]. Patients typically need long term follow-up, mood 45 stabilizing medications and psychoeducation to reduce the frequency of new mood episodes 46 and achieve mood stability [2,3]. Clinical assessments of BD symptoms are limited to 47 subjective evaluations by clinicians and patient recall bias and poor illness insight [4]. 48 Consequently, objective methods for monitoring state changes are sought-after, and new tools 49 to improve early detection and intervention in BD are in demand [5]. 50 51 The autonomic nervous system (ANS) controls a multitude of bodily functions and links the 52 central nervous system to peripheral organs such as the heart [6][7][8]. Heart rate variability 53 (HRV) reflects fluctuations in time intervals between consecutive heartbeats and is a well-54 documented and studied measure of ANS activity. Studying HRV allows us to observe the 55 adaptive abilities of the heart and ANS to environmental changes [9]. HRV measures are 56 commonly calculated using time-domain, frequency-domain or non-linear approaches. 57 Through these methods, it is possible to observe the balance between the "fight-or-flight" 58 sympathetic and "rest-and-digest" parasympathetic branches of the ANS [10][11][12]. 59 60 The adaptiveness of the ANS is well-studied within psychology and psychopathology. High 61 HRV scores have been found correlated with positive traits such as sustained attention and 62 communication abilities. Reduced vagally mediated HRV has been found across several 63 different psychopathological conditions and is therefore suggested to be a possible 64 transdiagnostic biomarker for mental illness [13]. 65 66 . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted February 15, 2022. ;https://doi.org/10.1101https://doi.org/10. /2022 4 BD is characterized by changes in activation, and correspondingly, ANS activity and function 67 has been found to be altered and dysregulated in BD [10,[14][15][16][17][18][19]. A 2017 review article and 68 meta-analysis reported reduced HRV in BD compared to healthy controls [10]. Moreover,69 two recent studies report an inverse relationship between BD illness severity and HRV [19, 70 20]. This is of interest considering diagnostic tool development, general understanding of the 71 disorder, and as a candidate biomarker of cardiovascular risk [21]. 72 73 HRV differences between affective states have been sparsely studied, and review articles 74 reveal that most of the literature investigates group level differences [10,17]. Only one prior 75 study has investigated within-individual HRV changes related to mood states. This 2018 76 study compared 19 manic males to their euthymic selves, finding HRV changes suggestive of 77 ANS improvement [22]. A small observational study found increased HRV in mania 78 compared to euthymia and depression on a group level, opposing most of existing literature 79 [23]. Finally, a group has investigated HRV changes in the transition from pathological to 80 euthymic states in ten bipolar patients, finding a steady increase as symptoms subsided [24]. 81 The same group used innovative technology to accurately classify affective states in long-82 term HRV analysis of eight BD patients [25]. These findings indicate that HRV can be used 83 as an objective biomarker for distinguishing between affective states and monitoring the 84 course of bipolar illness. However, more studies on state dependent HRV changes within 85 individuals as are needed. Based on this existing literature, we hypothesized lower HRV in 86 mania compared to euthymia. 87 88 Aims of the study 89 To identify within-individual differences in heart rate variability measures between manic 90 and euthymic states in bipolar patients. We used high-quality wrist-worn pulse sensors in a 91 . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The participants' only instructions were not to shower or touch the device during the 137 recording. We used Empatica's online software, E4 Connect, to visualize gross loss of sensor 138 contact, as identified by flattened skin conductance and temperature measurements and peaks 139 . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted February 15, 2022. ; https://doi.org/10.1101/2022.02.14.22269413 doi: medRxiv preprint in heart rate, and to download the inter-beat-intervals (IBI) data as comma-separated value 140 files for post-processing. 141 142 IBI data was analyzed in Kubios HRV Premium software version 3.4.2 [28]. We set the 143 artifact correction threshold to automatic and applied the smoothing priors detrending method 144 (l=500) [28]. We located 5-minute samples of high quality manually using a systematic 145 approach. As artifact rates in PPG recordings are lowest during rest, we focused on night-146 time samples during sleep, cross-checking the accelerometer data from the device as 147 assurance of rest [27]. We set the artifact correction cut-off to the recommended 5% and 148 considered data with a higher correction rate as poor quality sections. All analyses of HRV 149 were performed on samples from as close to initiation of sleep as possible, ranging from 8.10 150 p.m. to 5.10 a.m. (mean 11.34 p.m., standard deviation 112 minutes.) 151

152
We selected commonly used HRV measures from all three domains. From the time domain 153 measures, we selected root mean square of the successive differences (RMSSD), standard 154 deviation of average R-R intervals (SDNN) and mean heart rate (HR). We viewed RMSSD as 155 our main HRV measure as it mirrors the variance in time between heartbeats, is widely used 156 for monitoring vagally mediated (i.e. parasympathetic) HRV changes, and is well-suited for 157 use on our 5-minute data segments [29]. In the frequency domain, we chose to examine high 158 frequency (HF) and low frequency (LF) oscillations (ms 2 /Hz). HF is perceived as reflective 159 of predominately parasympathetic activity, while LF interpretation remains debated and 160 involved in both parasympathetic and sympathetic activity [30]. The frequency measures 161 were presented in logarithmic form of results from an autoregressive approach. We selected 162 Sample Entropy (SampEn) as the primary measure from the non-linear domain. It reflects the 163 complexity and degree of chaos in the heart rate series [29]. 164 . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint We used paired two-tailed t-tests to explore the significance of mood state on mean HR and 175 HRV features of interest: RMSSD, SDNN, LF, HF and SampEn. Significance levels were set 176 to p = 0.008, by the Bonferroni method (a = 0.05), due to the small number of subjects (N = 177 15) relative to the number of tests (six). This was done to control for false positive findings 178 due to multiple testing. Prior to analysis, RMSSD values were log transformed due to 179 distribution skewness. Effect sizes of change by mood state were calculated using Hedges' g. 180 Sex and age are customary confounders of HRV [29]. We ran separate linear regression 181 models for each of our five HRV measures to examine confounding effects. The dependent 182 variables of the linear regression models were the manic-euthymic differences of the HRV 183 measures (e.g. RMSSD) by turn, and the independent variables were sex and age. Due to the 184 small number of subjects, we did not include mood state as an independent variable in the 185 linear regression models, using the models purely to examine confounding effects. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint

200
Analysis results

201
Considerable changes in heart rate and several HRV measures were observed from mania to 202 euthymia, with overall large effect sizes (Fig 1). We found a statistically significant decrease 203 in HR with a large effect size (g = 0.92), although the range was large in the manic state (49-204

beats per minute). Time domain HRV measures SDNN and RMSSD increased 205
significantly (g = 1.01; g = 1.13) from manic to euthymic state, as did the frequency measure 206 LF (g = 1.24) and the non-linear measure Sample Entropy (g = 0.80). HF also increased from 207 mania to euthymia, but the change was not statistically significant. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted February 15, 2022. ; https://doi.org/10.1101/2022.02.14.22269413 doi: medRxiv preprint measures (see S1 Table). 217 218

219
There is a lack of disease monitoring options in bipolar disorder, and HRV is a promising 220 biomarker candidate. This study found a decrease in several HRV measures and increased 221 heart rate during mania when compared within individuals to their euthymic selves. It is one 222 of few studies examining state related changes longitudinally within bipolar subjects, and the 223 first to employ a user-friendly wrist worn device. These features are significant when 224 considering innovation possibilities. Our results are suggestive of decreased vagally mediated 225 parasympathetic activity during bipolar mania. 226 227 HRV differences were observed in measures from the time, frequency and non-linear 228 domains. From manic to euthymic state SDNN, RMSSD, LF and SampEn all increased 229 significantly with predominately large effect sizes. This is largely in agreement with a recent 230 study with a design and sample size resembling ours, although the only shared measurements 231 of statistical significance between our studies are RMSSD and SampEn [23]. Their study 232 approach was more controlled than our natural design, and there were several methodological 233 differences between our studies. First, they used an ECG chest strap for heart rate 234 measurement, instructing the participants to lie still for 20 minutes. Our study used a 235 lightweight wrist-worn PPG device and a natural approach, examining recordings made 236 . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted February 15, 2022. ;https://doi.org/10.1101https://doi.org/10. /2022 during sleep. Second, they excluded patients with high motor activity and all participants 237 were men. This could result in slightly different clinical profiles and a bias toward higher 238 HRV scores due to the male sex [37]. They also restricted medication intake to aripiprazole 239 and lithium, while we did not interfere with treatment. They did not control for possible 240 confounders. The similarity of our results, despite some methodological differences, 241 emphasize the robustness and biological underpinning of our shared findings. Lastly, we found a decline in heart rate in euthymia compared to mania (g = 0.82). This has 250 been described in previous literature [16,23,38]. Heart rate is regulated by several 251 mechanisms which are altered during a manic episode. Emotional and physiological stress 252 promote sympathetic upregulation in the limbic system, and secretion of stress hormones 253 [39]. Higher HR during mania is therefore not a surprising discovery. One could argue that 254 the higher HR during mania could be due to increased motor activity, and consequentially 255 higher HRV scores. However, we did not discover a significant difference in the amount of 256 movement, as measured by actigraphy in this sample [40]. 257 258 During the course of our study, methodological weaknesses became apparent. First, the 259 representability of our participants should be addressed. The design requires an active role 260 from treating physicians to refer manic inpatients capable of cooperating. This excludes 261 . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted February 15, 2022. ; https://doi.org/10.1101/2022.02.14.22269413 doi: medRxiv preprint subtypes of BD with paranoid traits or rapid mood changes, and those overly aggressive or 262 cognitively affected. Substance abuse is also a common comorbidity in BD, here excluded. 263 Due to the nature of scientific work, with ethical requirements and attempts at biological 264 classification, we believe that such a narrowing of scope is necessary at this stage. This is 265 also the case for comparable studies, as a certain degree of cooperation is required, and 266 comorbidities are commonly excluded to provide a clearer image of effects attributed solely 267 to BD. 268 269 Still, our sample contained older participants (43 ±13 years) compared to similar studies and 270 they were more educated than expected [41].We believe that this is partly due to our 271 recruitment method, appealing more to their inherent interest in their illness and contribution 272 to society. We experienced a higher degree of drop-out in younger participants as their manic 273 symptoms subsided. Such patients could for instance be lost due to a sudden repeal of 274 coerced treatment and subsequent ceased contact. A possible explanation is that both more 275 mature and educated BD patients are more likely to be interested in contributing to increased 276 knowledge and grounded enough to complete the study. 277

278
The effect of psychiatric drugs on HRV is controversial and study findings vary largely due 279 to heterogeneous methodology. Overall, the effect of psychiatric drugs on HRV has been 280 considered miniscule compared to the underlying HRV disturbances in psychiatric disorders. 281 The exceptions are tricyclic antidepressants and the atypical antipsychotic clozapine [38, 42, 282 43]. The patients in our sample were roughly on the same medication classes during manic 283 and euthymic recordings, see Table 1. Although dosages are not considered, this likely 284 implies smaller effects of medication differences on HRV by mood state. In addition to the 285 similarities between our findings and those of the previously mentioned study restricting 286 . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted February 15, 2022. ;https://doi.org/10.1101https://doi.org/10. /2022 medication use, similar results have been described during unmedicated mania, further 287 solidifying our results [16,23] instructions, we found that a selection of the 24-hour recordings were not of satisfactory 305 quality for HRV analyses. This was especially true for manic recordings. However, we 306 discovered a pattern of overall higher quality during rest at night. In order to keep as many 307 participants as possible in the analyses, we therefore focused on night recordings. This choice 308 provided a new caveat -the effect of circadian rhythms on HRV. In healthy individuals HRV 309 varies throughout the day and reaches a high-point around 03.00 a.m. [46]. In theory, due to 310 our participants falling asleep at different times (mean 11.34 p.m., standard deviation 112 311 . CC-BY 4.0 International license It is made available under a perpetuity. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted February 15, 2022. ; https://doi.org/10.1101/2022.02.14.22269413 doi: medRxiv preprint