Next Article in Journal
Patterns of Attention and Anxiety in Predicting Arithmetic Fluency among School-Aged Children
Next Article in Special Issue
Ultrasonic Deep Brain Neuromodulation in Acute Disorders of Consciousness: A Proof-of-Concept
Previous Article in Journal
Teaching Elaborative Reminiscing to Support Autobiographical Memory and Relationships in Residential and Community Aged Care Services
Previous Article in Special Issue
Pain in Persons with Disorders of Consciousness
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Brain Sciences
Volume 12
Issue 3
10.3390/brainsci12030375
brainsci-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Does the Heart Fall Asleep?—Diurnal Variations in Heart Rate Variability in Patients with Disorders of Consciousness

by
Monika Angerer
1,2,
Frank H. Wilhelm
3,
Michael Liedlgruber
3,
Gerald Pichler
4,
Birgit Angerer
5,
Monika Scarpatetti
4,
Christine Blume
6,7,† and
Manuel Schabus
1,2,*,†
1
Laboratory for Sleep, Cognition and Consciousness Research, Department of Psychology, University of Salzburg, 5020 Salzburg, Austria
2
Centre for Cognitive Neuroscience Salzburg (CCNS), University of Salzburg, 5020 Salzburg, Austria
3
Division of Clinical Psychology and Psychopathology, Department of Psychology, University of Salzburg, 5020 Salzburg, Austria
4
Apallic Care Unit, Albert Schweitzer Hospital, Geriatric Health Care Centres of the City of Graz, 8020 Graz, Austria
5
Private Practice for General Medicine and Neurology, 8430 Leibnitz, Austria
6
Centre for Chronobiology, Psychiatric Hospital of the University of Basel, 4002 Basel, Switzerland
7
Transfaculty Research Platform Molecular and Cognitive Neurosciences, University of Basel, 4055 Basel, Switzerland
*
Author to whom correspondence should be addressed.
†
Christine Blume and Manuel Schabus contributed equally to this work.
Brain Sci. 2022, 12(3), 375; https://doi.org/10.3390/brainsci12030375
Submission received: 30 January 2022 / Revised: 3 March 2022 / Accepted: 6 March 2022 / Published: 11 March 2022
(This article belongs to the Special Issue Emerging Treatments for Patients with Disorders of Consciousness)
Browse Figures
Review Reports Versions Notes

Abstract

:
The current study investigated heart rate (HR) and heart rate variability (HRV) across day and night in patients with disorders of consciousness (DOC). We recorded 24-h electrocardiography in 26 patients with DOC (i.e., unresponsive wakefulness syndrome (UWS; n = 16) and (exit) minimally conscious state ((E)MCS; n = 10)). To examine diurnal variations, HR and HRV indices in the time, frequency, and entropy domains were computed for periods of clear day- (forenoon: 8 a.m.–2 p.m.; afternoon: 2 p.m.–8 p.m.) and nighttime (11 p.m.–5 a.m.). The results indicate that patients’ interbeat intervals (IBIs) were larger during the night than during the day, indicating HR slowing. The patients in UWS showed larger IBIs compared to the patients in (E)MCS, and the patients with non-traumatic brain injury showed lower HRV entropy than the patients with traumatic brain injury. Additionally, higher HRV entropy was associated with higher EEG entropy during the night. Thus, cardiac activity varies with a diurnal pattern in patients with DOC and can differentiate between patients’ diagnoses and etiologies. Moreover, the interaction of heart and brain appears to follow a diurnal rhythm. Thus, HR and HRV seem to mirror the integrity of brain functioning and, consequently, might serve as supplementary measures for improving the validity of assessments in patients with DOC.

1. Introduction

Severe brain injury can cause coma, and, upon recovery, changes in consciousness often persist. These states are categorized under the term ‘disorders of consciousness’ (DOC). In a simplified approach, two major components are thought to be necessary for consciousness: wakefulness (i.e., the level of arousal) and awareness of the environment and the self (i.e., contents of consciousness) [1]. In patients living with DOC, wakefulness is preserved but awareness is only intermittently present or completely absent. More specifically, while patients with an unresponsive wakefulness syndrome (UWS) show some return of arousal (i.e., phases of sleep (closed eyes) and wakefulness (open eyes)) without signs of awareness during behavioral assessment, patients in a minimally conscious state (MCS) show inconsistent but reproducible signs of awareness that can be differentiated from reflexive behavior (e.g., response to commands, visual pursuit, intentional communication) [2,3]. If patients can communicate functionally and use objects adequately, their state is denoted as emergence from minimally conscious state (EMCS) [4]. Thus, while patients with UWS are assumed to be unconscious, patients in MCS and EMCS are assumed to be (at least minimally) conscious.
Clinical diagnoses are usually based on observations of the patients’ behavior using, for instance, the Glasgow Coma Scale [5] for acute situations or the Coma Recovery Scale-Revised (CRS-R) [6] for tracing their development during recovery. Unfortunately, behavioral assessments involve the risk of underestimating the level of consciousness [7]. This is because patients may be unable to respond behaviorally (and thus give evidence of their consciousness), for example, due to sensory or motor impairments. Finally, the fluctuating levels of consciousness carry the risk of examinations taking place during windows of unconsciousness [8]. Consequently, the absence of evidence for consciousness must not be mistaken for evidence of its absence. Hence, distinguishing between UWS and (E)MCS continues to be a challenge in clinical practice, and the rate of misdiagnoses is high (i.e., ~40%) [7] when comparing the medical consensus to the results obtained through clinical scales [9]. Approaches based on neuroimaging methods, such as functional magnetic resonance imaging (fMRI) or electroencephalography (EEG), have been used as additional tools to improve the validity of DOC diagnoses [10,11]. However, they require special expertise, come with high time and financial requirements, and often rely on tasks that are still challenging for patients [12]. Therefore, researchers and clinicians are looking for alternative or adjunct measures that (i) do not rely on patients’ behavioral responses, (ii) are time- and cost-efficient, (iii) and can be easily applied at beside. Lastly (iv), they should be useable during longer-term recordings, thereby taking the problem of fluctuating consciousness levels over the circadian day into account.
Among these measures, heart rate (HR) and heart rate variability (HRV) have been suggested to fulfill these criteria [13]. Specifically, HR indicates the average time interval between adjacent heartbeats (i.e., interbeat interval; IBI), while HRV quantifies the variability in these time intervals. These variations occur at different frequencies and reflect dynamics of the autonomic nervous system (ANS) regulation, being related to sympathetic and parasympathetic activity, breathing, thermo-, and blood pressure regulation, as well as changes in the vasomotor and renin-angiotensin system [14] (cf. Box 1). Furthermore, HR/HRV may mirror the interaction between the (injured) brain and the heart and thus represent a ‘peripheral’ window to ‘central’ functioning. The neural structure that enables the bidirectional communication of the heart and central nervous system has been described as the central autonomic network, which has been suggested to be involved in cognitive, emotional, and autonomic regulation, and linked to HR/HRV and cognitive performance [15].
In patients with severe brain injury, a decrease in HRV parameters in the time and frequency domain has been associated with a worse clinical outcome [16,17]. When looking at HRV entropy, it has been shown that patients with UWS show lower approximate entropy (ApEn) values than healthy controls [18]. Furthermore, higher complexity values, which were associated with higher CRS-R scores, were observed in patients in MCS as compared to patients with UWS [13]. Interestingly, no differences between patients with UWS and healthy controls have been found in any of the linear parameters (i.e., root mean square of successive differences between adjacent heartbeats (RMSSD), standard deviation of IBIs (SDRR), ratio between low and high frequencies (LF/HF ratio)) [18]. Furthermore, cognitive processes associated with ANS responses elicited by nociceptive stimuli have been identified in patients with MCS but not in patients with UWS. Specifically, following a repetitive nociceptive laser stimulation, patients in MCS showed a higher HR and SDRR as compared to before the stimulation, whereas patients with UWS did not show a significant change [19].
Importantly, like other physiological signals (such as body temperature and hormone secretion), HR/HRV in healthy populations shows strong alterations between sleep and wakefulness [20,21,22], so-called circadian variations (i.e., rhythms with a period length of approximately 24 h). In patients with severe brain injury, it has been shown that circadian temperature, melatonin, or motoric activity rhythms often deviate from the healthy norm and that a better integrity of the patients’ circadian rhythm is associated with a better clinical state [8,23,24]. Unfortunately, there is only limited research on how day vs. night variation affects the HR/HRV of patients with DOC. In patients with UWS, it has been shown that 6/11 patients had a higher HR during the night (i.e., 9 p.m.–5 a.m.) as compared to the day (5 a.m.–9 p.m.) [25]. Until now, it remains unclear to what degree the existence of diurnal variation in HR/HRV is linked to the patients’ clinical state.
Thus, our aim was to explore HR/HRV in patients with DOC in the time, frequency, and entropy domains. Specifically, as light is the most important zeitgeber for the internal biological clock [26], we also took the lighting conditions in the patients’ room into account and analyzed data from periods of definite day- and nighttime. Further, we investigated whether it is possible to differentiate between DOC diagnoses based on patients’ cardiac characteristics, and whether there is an association between the patients’ heart and brain activity.
Box 1. HRV parameters in the time, frequency, and entropy domains.
Methodologically, HRV can be evaluated in the frequency, time, and entropy domains. In the frequency domain, researchers usually compute a fast Fourier transform for three specific frequency bands that can be related to functions of the autonomic nervous system: high (HF; 0.15–0.4 Hz), low (LF; 0.04–0.15 Hz), and very low frequency (VLF; 0.003–0.04 Hz). While the HF band represents parasympathetic activity, the LF band reflects both parasympathetic and sympathetic activity and is particularly related to blood pressure regulation. The VLF band is mainly related to the thermoregulation, vasomotor, and renin–angiotensin system. In the time domain, variation in IBIs can be quantified by different parameters, such as the root mean square of successive differences between adjacent heartbeats (RMSSD), which is used to estimate variation mainly related to the HF band (pointing to parasympathetic activity; for other time-domain parameters, see Ref. [14]). By using non-linear methods (i.e., entropy domain), the complexity and irregularity of the IBI signal can be investigated. Approximate entropy (ApEn), detrended fluctuation analysis scaling exponent (DfaAlpha), Hurst exponent (Hurst), and sample entropy (SampEn) are amongst the most commonly used non-linear methods for HRV analysis [14,27,28]. Besides HRV, heart rate (HR) also provides information about autonomic nervous system regulation. Specifically, HR is a basic measure of overall cardiovascular arousal, which is influenced by both sympathetic activation and parasympathetic withdrawal.

2. Materials and Methods

2.1. Patients

From a total of 38 non-sedated and spontaneously breathing patients with DOC, nine patients had to be excluded from the analyses due to severe cardiac arrhythmias. Three further patients were excluded because they were exposed to constant 24-h light. Thus, 26 patients (8 women) between 16 and 80 years old (mean = 50.29, sd = 19.96) with different etiologies (traumatic brain injury (TBI): n = 11, non-traumatic brain injury (NTBI): n = 15) from clinics in Austria (n = 12) and Belgium (n = 14) were included in the analyses. The patients’ behavioral state was assessed with the CRS-R [6]. While 16 patients were diagnosed with UWS, nine patients were in an MCS, and one was in an EMCS. Patients with MCS and EMCS were combined for statistical analyses, which allows us to analyze differences between unconscious patients with UWS and (at least minimally) conscious patients in (E)MCS. Note that the data of this study were collected in a broader research context with the aim to quantify sleep and circadian parameters. Specifically, data of both patient samples (i.e., Austria and Belgium) have been used in two previous publications, where we studied diurnal variations in EEG parameters [29], as well as circadian variations in skin temperature [30] of patients with DOC. However, these studies did not focus on diurnal variations in electrocardiogram (ECG) data. The studies have been approved by the local ethics committees and informed consent was obtained from the patients’ legal representatives. For details on the study sample, see Table 1.

2.2. Study Protocol

2.2.1. Austria

In Austria, the study protocol comprised two within-subject conditions, namely the (i) habitual light (HL) and (ii) dynamic daylight (DDL) condition. Each condition lasted one week, and orders were randomized. While patients were in a room with standard clinic room lighting in the HL condition, patients were in a room with a ‘biologically effective’ room light (i.e., light intensity and spectrum that mimics natural daylight) in the DDL condition. During both weeks, skin temperature and actimetry were assessed continuously. ECG data were acquired at the beginning of each study condition (i.e., with a minimum recording duration of 1.5 days and a maximum recording duration of 4.9 days). The patients’ behavioral repertoire was assessed twice with the CRS-R at the end of each week, once in the morning and once in the afternoon, in all patients (n = 12). The data from the DDL condition are beyond the scope of this study. For the current analyses, the CRS-R assessments from the HL condition as well as 24-h continuous ECG data recorded at the beginning of the HL condition (i.e., starting on the first or second day) were analyzed (for more details on the study protocol, see Ref. [30]).

2.2.2. Belgium

In Belgium, a 24-h polysomnography (PSG; including EEG, ECG, electromyogram (EMG), electrooculogram (EOG), and respiration) was performed in the patients’ usual clinical environment with standard clinic room lighting (i.e., HL). The CRS-R was performed twice (i.e., before and after the recording) in nine and once in five of the 14 patients.
For the current analyses, the data of the Belgian patient sample were combined with the HL data of the Austrian patient sample as lighting conditions were comparable in these two datasets. Thus, we are looking at HR/HRV under standard clinic room lighting conditions.

2.3. Behavioral Assessment and Data Analysis

2.3.1. Coma Recovery Scale-Revised

The behavioral state of the patients was assessed with the CRS-R [6], which is composed of six subscales reflecting auditory, visual, motor, oromotor, communication and arousal functions with a total of 23 items. While the lowest threshold item on each subscale represents reflexive behavior, the highest threshold item indicates cognitively mediated behavior. The CRS-R is performed in a hierarchical manner, which means that the examiner starts with the highest item of each subscale and moves down the scale until the patient’s response meets the criteria for one item. When two CRS-R assessments were available, we used the CRS-R assessment where the patients showed the highest behavioral reactivity (i.e., the best diagnosis or highest sum score) as this is thought to best represent the ‘true’ state of the patient. The highest CRS-R score and diagnosis of each patient is shown in Table 1 and the subscale scores in Table S1 of the Supplementary Materials.

2.3.2. Electrocardiography

For recording ECG in the patient sample from Austria, an ambulatory three-channel ECG device (eMotion Faros 180°, Mega Electronics Ltd., Kuopio, Finland) with self-adhesive Ambu® BlueSensor SP electrodes (AMBU A/S, Ballerup, Denmark) was used. Two electrodes were placed in an infraclavicular position on the right and left body side, and one on a rib on the lower left thoracic wall. The sampling rate was 1000 Hz.
In the patient sample from Belgium, the ECG was recorded as part of the PSG. For PSG recordings, BrainProducts amplifiers (BrainProducts, Gilching, Germany) were used. All PSG signals were recorded with a sampling rate of 500 Hz. For ECG recordings, two goldcup electrodes were used. One electrode was affixed in an infraclavicular position on the right body side and one on a rib on the lower left thoracic wall using EC2 electrode gel (Astro-Med® Inc., West Warwick, RI, USA).

2.3.3. Heart Rate and Heart Rate Variability

To explore diurnal variations in HR/HRV, we used continuous 24-h ECG data from each patient and divided the ECG recording into periods of clear daytime (i.e., forenoon: 8 a.m.–2 p.m., afternoon: 2 p.m.–8 p.m.) and nighttime (i.e., 11 p.m.–5 a.m.) using BrainVision Analyzer 2.0 Software (Brain Products GmbH, Gilching, Germany). Importantly, all three periods were of equal length (i.e., 6 h). Periods with conditions of twilight (i.e., dawn: 5–8 a.m., dusk: 8–11 p.m.) were excluded from the analyses.
HRV analyses were conducted in ANSLAB 2.6 [31]. Accurate automatic R peak detections (obtained with the software’s algorithms) were carefully visually checked for the entire data of each patient and corrected whenever necessary. In a second and important preprocessing step, the IBI data were visually inspected for arrhythmias (e.g., ectopic beats) as single spikes in the IBI signal can seriously distort the spectral estimates of HRV [32]. While single spikes were interpolated with the software’s algorithm (i.e., mean of the adjacent IBIs), longer phases of arrhythmia were marked and excluded from the analyses (i.e., set to missing). In total, we removed per patient and daytime on average 0.3% (min = 0%, max = 7.2%) of the data points (for the amount of individual patients’ missing data, see Supplementary Materials: Table S2). Importantly, we paid special attention to adequate artefact correction because of the susceptibility of ApEn to noisy data, outliers, and missed heart beat detection [27].
Patients’ HRV data were analyzed in the frequency, time, and entropy domains. Due to obvious non-stationarities in the IBI signals (i.e., sudden shifts in the mean and/or standard deviation) we used complex demodulation (CDM) for quantifying the HRV frequency parameters. CDM produces equivalent results to fast Fourier transform but is less affected by non-stationarity (a prerequisite for accurate fast Fourier transform) [33,34].
We computed the mean IBI oscillation amplitude (in ms) for the very low (VLF; 0.003–0.04 Hz), low (LF; 0.04–0.15 Hz), and high (HF; 0.15–0.4 Hz) frequency bands for each minute. During the export of the mean IBI oscillation amplitude, ANSLAB automatically interpolates segments that were set to missing. Thus, we post hoc excluded all segments that were automatically interpolated using R version 3.6.1. [35]. Mean IBI (in ms; i.e., the time interval between two successive R-peaks), mean HR (equivalent to 60,000/IBI; in beats/min), and mean RMSSD (in ms) were also computed for each minute. For statistical analyses of HR/HRV parameters, all one-minute segments within each 6-h period (i.e., 8 a.m.–2 p.m., 2 p.m.–8 p.m., and 11 p.m.–5 a.m.) were averaged using R. Thus, we arrived at one value per patient and HR/HRV parameter for forenoon, afternoon, and night. For quantifying the complexity of the patients’ HR variations, ApEn, detrended fluctuation analysis scaling exponent (DfaAlpha), sample entropy (SampEn), and the Hurst exponent (Hurst) were computed for the whole 6-h segments (i.e., forenoon, afternoon, and night). ApEn describes the likelihood that patterns remain similar for subsequent comparisons, with higher ApEn values indicating higher irregularity and complexity in time-series data. DfaAlpha quantifies the presence or absence of fractal correlation properties in non-stationary time-series data, with higher DfaAlpha values suggesting fractal-like HR dynamics. SampEn describes the probability that two sequences of 1–5 consecutive data points (i.e., SampEn1-5) that are similar to each other will remain similar when one consecutive point is included. Higher SampEn indicates frequent incidence of dissimilarities in the time-series data. The Hurst exponent is a quantification of long-range dependence, with values of 0.5–1 indicating a long-term positive correlation, meaning that a high value in the series will probably be followed by another high value, values of 0–0.5 indicating time series with long-term negative correlation or switching between high and low values, and a value of 0.5 suggesting that the series are uncorrelated. For more details on non-linear HRV analyses, see Refs. [27,28].

2.3.4. Respiration

Since DOC may be related to respiratory alterations [36] and as low respiration rates (i.e., <9 cycles per minute (cpm) ≤ 0.15 Hz) can affect the HF estimates of HRV via shifts in the respiratory sinus arrhythmia to the LF band [37], we analyzed the respiratory signal in a subsample of patients (n = 12) for whom respiration was simultaneously recorded via a respiration belt around the thorax. Mean respiratory rate for each minute was computed in ANSLAB and averaged for 5-min segments using R. We found that only one 5-min segment (i.e., 1.39% of 6-h recordings) of two patients (i.e., P17: forenoon; P26: night) was <9 cpm. Thus, we decided that there was no need to control the HF parameter for respiration in our patient sample. Although respiration data were only available from a subset of patients, it is unlikely that breathing patterns in the other patients were appreciably different.

2.3.5. EEG Permutation Entropy

For analyzing heart–brain interaction, we correlated HRV entropy with EEG entropy in the Belgian patient sample (n = 14) where PSG was recorded. As the EEG data of the Belgian sample were already used in a previous study, where EEG permutation entropy (PE) had been computed, we used the EEG entropy values of the respective patients from that publication [29]. More specifically, PE of the entire EEG signal was computed for day (i.e., 8 a.m.–8 p.m.) and night (i.e., 11 p.m.–5 a.m.). PE quantifies the level of irregularity or unpredictability of an EEG signal. Higher PE values indicate more complex and/or random signals. For further information on the preprocessing of the EEG signal and the entropy analyses, see Ref. [29].

2.4. Statistical Analyses

Statistical analyses were done in R. In a first step, we evaluated whether the distribution of the data followed a normal distribution. As this was not the case for several variables (i.e., Shapiro–Wilk test for normality: p < 0.001; cf. Supplementary Materials: Table S3), we opted for advanced semi-parametrical and non-parametrical statistical tests. Specifically, for the analyses of differences in HR/HRV parameters (i.e., IBI, HR, RMSSD, VLF, LF, HF, ApEn, DfaAlpha, Hurst, and SampEn) and EEG entropy (i.e., PE) between different times of the day (i.e., within-subjects factor; forenoon, afternoon, night), diagnoses (i.e., between-subjects factor; (E)MCS vs. UWS), etiology (i.e., between-subjects factor; TBI vs. NTBI), or sex (i.e., between-subjects factor; female vs. male), we used advanced semi-parametric analyses for repeated measures designs as implemented in the ‘MANOVA.RM’ package available for R [38]. We report the resampled Wald-type statistic (WTS). As a resampling method, the function ‘perm’ was used, which randomly permutes all observations. The number of iterations used to calculate the resampled statistic was 10,000. To correct for multiple tests, p-values of post hoc comparisons were adjusted using the method of Benjamini and Hochberg (BH) [39] as implemented in the ‘p.adjust’ function in R. For correlation analyses of HR/HRV parameters, CRS-R sum score, and EEG entropy, we report Kendall’s Tau.
The significance level was set to α = 0.05 (two-sided) for all analyses. As suggested by Wasserstein et al. [40], we interpreted the overall pattern rather than focusing on individual p-values. Therefore, we also interpreted p-values 0.05 < p ≤ 0.10 if they were in line with the overall pattern of results and also ‘credible’ from a Bayesian point of view. More specifically, we additionally tested all trends of our post hoc comparisons and correlation analyses with Bayesian multilevel regression models via the Stan-based ‘brms’ package available for R [41,42]. We report regression coefficients and the 95% credible intervals (CIs; i.e., Bayesian confidence intervals). The CI describes the interval in which a parameter value falls with a 95% probability given the data observed, prior and model assumptions. Thus, an effect is considered to significantly differ from zero if zero is not included in the CI. We used weak- or non-informative default priors whose influence on results is negligible. For all computed regression models, no divergent transitions occurred, Rhat (i.e., potential scale reduction factor on split chains) was < 1.01, and effective sample size (ESS) was > 400.

3. Results

3.1. Interbeat Interval and Heart Rate

Analyses of the IBIs of 26 patients revealed a trend towards a main effect for time (FWTS(2) = 6.52, p = 0.068) and a significant effect for diagnosis (FWTS(1) = 5.8, p = 0.028) but no significant time   ×   diagnosis interaction (FWTS(2) = 1.36, p = 0.523). Specifically, the patients showed larger IBIs during the night as compared to forenoon (FWTS(1) = 7.76, p = 0.027) and afternoon (FWTS(1) = 4.74, p = 0.059; b = −37.44, 95%CI = [−67.46, −7.42]). No differences could be observed in the patients’ IBIs between fore- and afternoon (FWTS(1) = 0.21, p = 0.648; cf. Figure 1a). Further, patients in UWS showed larger IBIs as compared to patients in (E)MCS (FWTS(1) = 5.8, p = 0.028; cf. Figure 1b).
This effect was supported by correlation analyses of IBIs and CRS-R sum scores of 25 patients showing that lower CRS-R sum scores were associated with larger IBIs during forenoon (r τ (23) = −0.34, p = 0.02; cf. Figure 2a) and afternoon (r τ (23) = −0.38, p = 0.009; cf. Figure 2b) and by trend during the night (r τ (23) = −0.27, p = 0.07; cf. Figure 2c). However, we refrain from interpreting the effect during the night as it does not appear robust from a Bayesian point of view (b = −12.34, 95%CI = [−25.79, 1.16]).
Analyses of HR yielded similar results (cf. Supplementary Materials: Figures S1 and S2), which is expected as HR is inversely proportional to the IBI signal. In both IBI and HR, no differences for etiology (i.e., TBI vs. NTBI) or sex (i.e., male vs. female) were observed (cf. Supplementary Materials: Tables S4 and S5).

3.2. HRV Time Domain

Analyses of the RMSSD of 26 patients did not reveal significant main effects for time (FWTS(2) = 0.4, p = 0.828), diagnosis (FWTS(1) = 0.37, p = 0.554), and time   ×   diagnosis interaction (FWTS(2) = 3.91, p = 0.174; cf. Figure 3). No differences for etiology (i.e., TBI vs. NTBI) or sex (i.e., male vs. female) were observed (cf. Supplementary Materials: Tables S4 and S5).

3.3. HRV Frequency Domain

Analyses of the VLF of 26 patients yielded a trend towards a main effect for time (FWTS(2) = 6.63, p = 0.064) but no significant effect for diagnosis (FWTS(1) = 1.74, p = 0.207) and the time   ×   diagnosis interaction (FWTS(2) = 1.08, p = 0.603). Post hoc comparisons of VLF between times of day did not yield significance anymore after correcting for multiple comparisons (i.e., fore- vs. afternoon: FWTS(1) = 1.15, p = 0.295; forenoon vs. night: FWTS(1) = 1.26, p = 0.295; afternoon vs. night: FWTS(1) = 4.73, p = 0.111 (before BH-correction: p = 0.03); cf. Figure 4a).
Analyses of the LF of 26 patients revealed a trend towards a main effect for time (FWTS(2) = 5.69, p = 0.088) but no significant effect for diagnosis (FWTS(1) = 0.46, p = 0.51) and the time   ×   diagnosis interaction (FWTS(2) = 0.36, p = 0.843). Post hoc comparisons of LF between times of day did not yield significance anymore after correcting for multiple comparisons (i.e., forenoon vs. afternoon: FWTS(1) = 1.65, p = 0.213; forenoon vs. night: FWTS(1) = 1.75, p = 0.213; afternoon vs. night: FWTS(1) = 4.86, p = 0.108 (before BH-correction: p = 0.028); cf. Figure 4b).
Analyses of the HF of 26 patients did not reveal significant main effects for time (FWTS(2) = 3.8, p = 0.191), diagnosis (FWTS(1) = 0.57, p = 0.460), and the time   ×   diagnosis interaction (FWTS(2) = 0.04, p = 0.980; cf. Figure 4c).
No etiology or sex differences were observed for any frequency band (cf. Supplementary Materials: Tables S4 and S5).

3.4. HRV Entropy Domain

Analyses of the ApEn of 26 patients revealed a significant main effect for time (FWTS(2) = 21.35, p = 0.001) but no significant effect for diagnosis (FWTS(1) = 1.79, p = 0.197) and the time   ×   diagnosis interaction (FWTS(2) = 1.99, p = 0.395). Specifically, the patients showed a higher ApEn during forenoon as compared to afternoon (FWTS(1) = 15.79, p < 0.001). No differences could be observed in the patients’ ApEn during the night as compared to forenoon (FWTS(1) = 2.19, p = 0.154) and afternoon (FWTS(1) = 2.77, p = 0.154; cf. Figure 5). Analyses of other entropy parameters did not yield significant main effects for time and diagnoses and the time   ×   diagnosis interaction (cf. Supplementary Materials: Table S6).
However, analyses of the DfaAlpha and SampEn1 yielded a significant main effect for etiology, with higher DfaAlpha (FWTS(1) = 4.59, p = 0.043; cf. Figure 6a) and SampEn1 (FWTS(1) = 4.36, p = 0.049; cf. Figure 6b) being observed in the patients with TBI as compared to patients with NTBI. No etiology and sex differences were observed for the other entropy parameters (cf. Supplementary Materials: Tables S4 and S5).

3.5. Correlation of EEG and HRV Entropy

Correlation analyses in 14 patients showed that a higher ApEn was associated with a higher PE during the night (i.e., 11 p.m.–5 a.m.; r τ (12) = 0.47, p = 0.019; cf. Figure 7b). No significant correlation was observed during the day (i.e., 8 a.m.–8 p.m., r τ (12) = 0.3, p = 0.157, cf. Figure 7a).

4. Discussion

In patients with DOC, variations in cardiac activity show a diurnal pattern. Specifically, we find preserved diurnal variations in the length of the patients’ IBIs. The IBIs were larger during the night than during the day, indicating that, as in healthy individuals [22], the heart slows down during the night due to parasympathetic dominance reflecting relaxation and sleep. Further, the complexity of patients’ HRV signal varies across wakefulness, with the signal being more irregular and complex (i.e., higher ApEn) during forenoon as compared to afternoon. This has also been found in healthy individuals [43] and might be due to an increased cardiac sympathovagal response in the morning after awakening. More specifically, it has been shown in healthy individuals that sleep-to-wake transitions in the morning are associated with higher sympathetic activation compared to those occurring during the rest of the day [22]. Thus, although patients with DOC often fluctuate between sleep and wake-like phases—hence experiencing several ‘awakenings’ throughout the day—the results indicate that the awakening in the morning is probably the one associated with the most prominent change in arousal. Further, it might also be the case that patients show less frequent or temporally more regular (i.e., systematic) alterations of arousal and/or awareness in the afternoon, which may be associated with lower entropy values. This would go in line with our findings from a previous study, where the patients tended to exhibit higher CRS-R sum scores at a later daytime (i.e., afternoon) [8], which also requires more stable arousal and awareness levels. Another explanation for the entropy drop from fore- to afternoon might be the recovery from stressful events that possibly took place more frequently during forenoon (e.g., therapies, medical rounds, nursing activities). It has been shown in healthy controls that a stressful task leads to a significant reduction in entropy in the succeeding relaxation period [44]. Thus, the variation in cardiac activity over the day might be necessary for an optimal interaction and adaptation to changing demands in the environment. Interestingly, while we found that the patients in UWS generally had larger IBIs (i.e., lower HRs) than the patients in (E)MCS, lower IBIs were associated with higher behavioral reactivity (i.e., higher CRS-R sum scores) in the patients with DOC only during the day but not during the night. Looking at the data, one can see that this effect is mainly driven by the nocturnal slowing of the heart in the patients in an (E)MCS, underlining that cardiac activity shows a diurnal pattern, particularly when patients have (partially) regained consciousness. Thus, our findings complement earlier research suggesting that better circadian rhythm integrity is associated with higher consciousness levels [8,23,24]. More specifically, we found in earlier studies that variations in peripheral biosignals, such as skin temperature, melatonin(-sulfate), and wrist actimetry, are better aligned to a healthy 24-h rhythm (i.e., circadian rhythm) and more pronounced in patients with a higher behavioral repertoire [8,23,24].
Concerning central biosignals, such as brain activity derived from EEG recordings, earlier studies have shown that more severely affected patients with UWS do not only show a stronger general slowing of the EEG but also no clear diurnal pattern [29,45]. This is in line with the findings we present here (i.e., lower HR/larger IBIs in the patients with UWS) and possibly reflects the interaction between peripheral and central biosignals. Thus, we additionally investigated whether there is an association between the heart and central measures of brain activity, that is, EEG entropy (i.e., a measure describing the level of irregularity or unpredictability of the brain signal). We found that heart and brain activity are coupled during the night but not during the day. Specifically, while during the night a higher EEG entropy (i.e., PE) was associated with a higher ECG entropy (i.e., ApEn), no such association was evident during the day. It might be the case that, during the night (i.e., habitual sleep), and due to less disturbance by external cues from the environment and a stronger focus on internal processes, brain and body rhythms are better connected. Conversely, during the day, the brain focuses more on the processing of sensory signals in the environment, leading to a decrease in the synchrony between brain and body rhythms. Interestingly, the effect in our data seems to be mainly driven by EEG entropy. Specifically, the patients showed lower EEG entropy values during the night as compared to the day (see Supplementary Materials: EEG Entropy), which might be an effect of the dominance of slow oscillations during the night (i.e., more synchronized brain activity, and thus less signal complexity) and probably indicates the existence of a ‘sleep-like state’ during the night. ECG entropy, however, does not differ between day and night. A reason for this could be that the patients spend most of the time in a lying position or a position where the upper body is raised in a 45° angle, which reduces events that usually influence HRV, such as changes in posture or physical movements. The brain signal, however, can still change independently—based on the patients’ clinical state. Specifically, while the patients in an (E)MCS showed higher EEG complexity during the day than during the night, no diurnal variation was evident in the patients with UWS (see Supplementary Materials: EEG Entropy).
Additionally, cardiac activity does not only inform regarding the integrity of diurnal variations and the degree of behavioral reactivity in patients with DOC but also differentiates between patients’ etiologies. Specifically, the patients with TBI showed a higher DfaAlpha and SampEn1 than the patients with NTBI, suggesting fractal-like (i.e., aperiodic) HR dynamics and higher dissimilarities in the HRV signal of the patients with TBI. In other words, cardiac activity is less complex and variable in patients with NTBI. Reduced HRV complexity has been shown to be a predictor for mortality [46,47]. This is in line with previous findings showing that patients with NTBI often have a less favorable prognosis than patients with TBI [48,49,50]. Importantly, while this study serves as a proof of principle study, additional clinical and multicenter studies are needed that investigate the extent of diagnostic quality (i.e., in terms of sensitivity and specificity) of the different HRV parameters.
When looking at the time (i.e., RMSSD) and frequency (i.e., VLF, LF, HF) domains of HRV, no differences between time of day, diagnoses, and etiologies were evident. One reason might be of a methodological nature. More specifically, the activity of the heart is not regular/periodic but rather fluctuates in complex/aperiodic patterns. Thus, it has frequently been argued that non-linear measures (i.e., measures of mathematical chaos/entropy) might be more appropriate for the analysis of HRV data [51,52]. This is in line with the findings from a study that showed a difference in HRV entropy between patients with UWS and healthy individuals (i.e., lower ApEn in patients with UWS) but no such differences in any of the linear parameters (i.e., IBI, SDRR, RMSSD, LF/HF ratio) [18].
Our results show that ApEn is of diagnostic relevance in our study sample. However, due to its special requirements and limitations, it has been discussed whether it is an ideal measure for quantifying entropy [27,53]. In our study, the often instead used ‘SampEn’ measure does not reveal the same results as ApEn (cf. Supplementary Materials: Table S6). A possible reason might be that ApEn correlates with RMSSD and HF power (see Table 1 in Ref. [27]), while SampEn correlates, rather, with LF power. ApEn thus combines aspects of HRV time and frequency domain measures with the nonlinear measure ‘entropy’ in a way that seems to be of particular diagnostic relevance for our research questions in severely brain-injured patients.

5. Conclusions

To summarize, patients with severe brain injuries—particularly those who (partially) regained consciousness—still had preserved diurnal variations, as characterized by a heart rate slowing during the night. This suggests preserved integrity of circadian rhythms in the autonomic nervous system activity. Further, cardiac activity differentiated between patients’ etiologies and diagnoses. The patients with UWS had larger IBIs (i.e., lower heart rate) than the patients in an (E)MCS, and the patients with NTBI had a less complex HRV signal than the patients with TBI. Thus, cardiac activity and its variations might represent a peripheral window to central (brain) functioning. Indeed, we found an interaction of heart and brain signal complexity, which also followed a diurnal pattern. Specifically, while a more complex brain signal was associated with a more complex heart signal during the night, no such association was found during the day. In conclusion, HR and HRV seem to mirror the integrity of brain functioning and, consequently, might serve as supplementary measures that aid the differentiation between clinical states. Ultimately, this has the potential to improve the validity of assessments in patients with DOC.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/brainsci12030375/s1, Table S1. CRS-R scores; Table S2. Amount of missing data from each subject separately for time (i.e., forenoon [8 a.m.–2 p.m.], afternoon [2 p.m.–8 p.m.], night [11 p.m.–5 a.m.]); Table S3. Shapiro-Wilk tests (p-values) for normality for each variable separately for time (i.e., forenoon [8 a.m.–2 p.m.], afternoon [2 p.m.–8 p.m.], night [11 p.m.–5 a.m.]); Figure S1. Heart rate (HR) separately for time and diagnosis contrasts; Figure S2. Correlation between heart rate (HR) and CRS-R sum score separately for time; Table S4. Wald-type statistic (WTS) of the different HRV parameters separately for the main effects ‘time’ (i.e., forenoon [8 a.m.–2 p.m.], afternoon [2 p.m.–8 p.m.], night [11 p.m.–5 a.m.]) and ‘etiology’ (i.e., traumatic vs. non-traumatic brain injury), and the ‘time × etiology’ interaction (N = 26); Table S5. Wald-type statistic (WTS) of the different HRV parameters separately for the main effects ‘time’ (i.e., forenoon [8 a.m.–2 p.m.], afternoon [2 p.m.–8 p.m.], night [11 p.m.–5 a.m.]), ‘sex’ (i.e., male vs. female), and the ‘time × sex’ interaction (N = 26); Table S6. Wald-type statistic (WTS) of the different entropy parameters separately for the main effects ‘time’ (i.e., forenoon [8 a.m.–2 p.m.], afternoon [2 p.m.–8 p.m.], night [11 p.m.–5 a.m.]) and ‘diagnosis’ (i.e., E/MCS vs. UWS), and the ‘time × diagnosis’ interaction (N = 26).

Author Contributions

Conceptualization, M.S. (Manuel Schabus), C.B. and M.A.; methodology, M.S. (Manuel Schabus), M.A., C.B., M.L. and F.H.W.; software, M.A., M.L. and F.H.W.; formal analysis, M.A.; investigation, M.S. (Manuel Schabus), G.P., B.A. and M.S. (Monika Scarpatetti); resources, M.S. (Manuel Schabus), F.H.W. and M.L.; data curation, M.A. and M.S. (Manuel Schabus); writing—original draft preparation, M.A. and C.B.; writing—review and editing, M.A., C.B., F.H.W., M.S. (Manuel Schabus) and M.L.; visualization, M.A.; supervision, M.S. (Manuel Schabus) and C.B.; project administration, M.S. (Manuel Schabus) and M.A.; funding acquisition, M.S. (Manuel Schabus). All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Austrian Science Fund (FWF; Y-777), and the Belgian National Funds for Scientific Research (FNRS/FRIA). M.A. was additionally supported by the Doctoral College ‘Imaging the Mind’ (FWF; W1233-B), and another grant from the FWF (P-33630). C.B. was supported by grants from the University of Basel, the Austrian Science Fund (FWF; J-4243), the Psychiatric Hospital of the University of Basel, the Freiwillige Akademische Gesellschaft Basel, and the Novartis Foundation for biological–medical research.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committees of the Medical University of Liège and Graz.

Informed Consent Statement

Informed consent was obtained from the legal representatives of all patients involved in the study.

Data Availability Statement

The data that support the findings of this study are available from the corresponding authors upon reasonable request.

Acknowledgments

We thank Sarah Haberl, Kerstin Dokter, Marlen Dido Grand, and Julius Köppen, as well as the staff at the Albert Schweitzer Hospital in Graz for their constant support and help with the data collection. Further, we thank Marina Meier and Jasmin Preiß for their support with the data preparation. Additionally, we thank Malgorzata Wislowska for sharing patients’ EEG entropy data with us. Finally, we also thank the Coma Science Group (Liège) for contributing data, and the patients’ families for their trust and cooperation in this study. Open Access Funding by the Austrian Science Fund (FWF).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ANSAutonomic nervous system
ApEnApproximate entropy
BHBenjamini–Hochberg correction for multiple comparisons
CDMComplex demodulation
CICredibility interval
CpmCycles per minute
CRS-RComa Recovery Scale-Revised
DDLDynamic daylight
DfaAlphaDetrended fluctuation analysis scaling exponent
DOCDisorders of consciousness
ECGElectrocardiography
EEGElectroencephalography
EMCSEmergence from minimally conscious state
EMGElectromyography
EOGElectrooculography
fMRIFunctional magnetic resonance imaging
HFHigh frequencies (0.15–0.4 Hz)
HLHabitual light
HRHeart rate
HRVHeart rate variability
HurstHurst exponent
IBIInterbeat interval
LFLow frequencies (0.04–0.15 Hz)
LF/HF ratioRatio between low and high frequencies
MCSMinimally conscious state
NTBINon-traumatic brain injury
PEPermutation entropy
PSGPolysomnography
RMSSDRoot mean square of successive differences between adjacent heartbeats
SampEnSample entropy
SDRRStandard deviation of interbeat intervals
TBITraumatic brain injury
UWSUnresponsive wakefulness syndrome
VLFVery low frequencies (0.003–0.04 Hz)
WTSWald-type statistic

References

  1. Laureys, S. The neural correlate of (un)awareness: Lessons from the vegetative state. Trends Cogn. Sci. 2005, 9, 556–559. [Google Scholar] [CrossRef] [PubMed]
  2. Laureys, S.; Celesia, G.G.; Cohadon, F.; Lavrijsen, J.; Leon-Carrion, J.; Sannita, W.G.; Sazbon, L.; Schmutzhard, E.; von Wild, K.R.; Zeman, A.; et al. Unresponsive wakefulness syndrome: A new name for the vegetative state or apallic syndrome. BMC Med. 2010, 8, 1–4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Giacino, J.T.; Malone, R. The vegetative and minimally conscious states. Handb. Clin. Neurol. 2008, 90, 99–111. [Google Scholar] [CrossRef] [PubMed]
  4. Bruno, M.A.; Vanhaudenhuyse, A.; Thibaut, A.; Moonen, G.; Laureys, S. From unresponsive wakefulness to minimally conscious PLUS and functional locked-in syndromes: Recent advances in our understanding of disorders of consciousness. J. Neurol. 2011, 258, 1373–1384. [Google Scholar] [CrossRef] [PubMed]
  5. Teasdale, G.; Jennett, B. Assessment of coma and impaired consciousness. A practical scale. Lancet 1974, 2, 81–84. [Google Scholar] [CrossRef]
  6. Kalmar, K.; Giacino, J.T. The JFK Coma Recovery Scale—revised. Neuropsychol. Rehabil. 2005, 15, 454–460. [Google Scholar] [CrossRef] [PubMed]
  7. Andrews, K.; Murphy, L.; Munday, R.; Littlewood, C. Misdiagnosis of the vegetative state: Retrospective study in a rehabilitation unit. Brit. Med. J. 1996, 313, 13–16. [Google Scholar] [CrossRef] [Green Version]
  8. Blume, C.; Angerer, M.; Raml, M.; del Giudice, R.; Santhi, N.; Pichler, G.; Kunz, A.B.; Scarpatetti, M.; Trinka, E.; Schabus, M. Healthier rhythm, healthier brain? Integrity of circadian melatonin and temperature rhythms relates to the clinical state of brain-injured patients. Eur. J. Neurol. 2019, 26, 1051–1059. [Google Scholar] [CrossRef] [Green Version]
  9. Schnakers, C.; Vanhaudenhuyse, A.; Giacino, J.; Ventura, M.; Boly, M.; Majerus, S.; Moonen, G.; Laureys, S. Diagnostic accuracy of the vegetative and minimally conscious state: Clinical consensus versus standardized neurobehavioral assessment. BMC Neurol. 2009, 9, 35. [Google Scholar] [CrossRef] [Green Version]
  10. Demertzi, A.; Antonopoulos, G.; Heine, L.; Voss, H.U.; Crone, J.S.; de Los Angeles, C.; Bahri, M.A.; Di Perri, C.; Vanhaudenhuyse, A.; Charland-Verville, V.; et al. Intrinsic functional connectivity differentiates minimally conscious from unresponsive patients. Brain 2015, 138, 2619–2631. [Google Scholar] [CrossRef] [Green Version]
  11. Stender, J.; Gosseries, O.; Bruno, M.A.; Charland-Verville, V.; Vanhaudenhuyse, A.; Demertzi, A.; Chatelle, C.; Thonnard, M.; Thibaut, A.; Heine, L.; et al. Diagnostic precision of PET imaging and functional MRI in disorders of consciousness: A clinical validation study. Lancet 2014, 384, 514–522. [Google Scholar] [CrossRef]
  12. Monti, M.M.; Vanhaudenhuyse, A.; Coleman, M.R.; Boly, M.; Pickard, J.D.; Tshibanda, L.; Owen, A.M.; Laureys, S. Willful Modulation of Brain Activity in Disorders of Consciousness. N. Engl. J. Med. 2010, 362, 579–589. [Google Scholar] [CrossRef] [Green Version]
  13. Riganello, F.; Larroque, S.K.; Bahri, M.A.; Heine, L.; Martial, C.; Carriere, M.; Charland-Verville, V.; Aubinet, C.; Vanhaudenhuyse, A.; Chatelle, C.; et al. A Heartbeat Away From Consciousness: Heart Rate Variability Entropy Can Discriminate Disorders of Consciousness and Is Correlated With Resting-State fMRI Brain Connectivity of the Central Autonomic Network. Front. Neurol. 2018, 9, 769. [Google Scholar] [CrossRef]
  14. Shaffer, F.; Ginsberg, J.P. An Overview of Heart Rate Variability Metrics and Norms. Front. Public Health 2017, 5, 258. [Google Scholar] [CrossRef] [Green Version]
  15. Thayer, J.F.; Hansen, A.L.; Saus-Rose, E.; Johnsen, B.H. Heart Rate Variability, Prefrontal Neural Function, and Cognitive Performance: The Neurovisceral Integration Perspective on Self-regulation, Adaptation, and Health. Ann. Behav. Med. 2009, 37, 141–153. [Google Scholar] [CrossRef]
  16. Gujjar, A.R.; Sathyaprabha, T.N.; Nagaraja, D.; Thennarasu, K.; Pradhan, N. Heart rate variability and outcome in acute severe stroke: Role of power spectral analysis. Neurocrit. Care 2004, 1, 347–353. [Google Scholar] [CrossRef]
  17. Rapenne, T.; Moreau, D.; Lenfant, F.; Vernet, M.; Boggio, V.; Cottin, Y.; Freysz, M. Could heart rate variability predict outcome in patients with severe head injury? A pilot study. J. Neurosurg. Anesthesiol. 2001, 13, 260–268. [Google Scholar] [CrossRef]
  18. Sara, M.; Sebastiano, F.; Sacco, S.; Pistoia, F.; Onorati, P.; Albertini, G.; Carolei, A. Heart rate non linear dynamics in patients with persistent vegetative state: A preliminary report. Brain Inj. 2008, 22, 33–37. [Google Scholar] [CrossRef]
  19. Leo, A.; Naro, A.; Cannavo, A.; Pisani, L.R.; Bruno, R.; Salviera, C.; Bramanti, P.; Calabro, R.S. Could autonomic system assessment be helpful in disorders of consciousness diagnosis? A neurophysiological study. Exp. Brain Res. 2016, 234, 2189–2199. [Google Scholar] [CrossRef]
  20. Garbarino, S.; Lanteri, P.; Feeling, N.R.; Jarczok, M.N.; Quintana, D.S.; Koenig, J.; Sannita, W.G. Circadian Rhythms, Sleep, and the Autonomic Nervous System. J. Psychophysiol. 2020, 34, 1–9. [Google Scholar] [CrossRef]
  21. Bilan, A.; Witczak, A.; Palusinski, R.; Myslinski, W.; Hanzlik, J. Circadian rhythm of spectral indices of heart rate variability in healthy subjects. J. Electrocardiol. 2005, 38, 239–243. [Google Scholar] [CrossRef] [PubMed]
  22. Boudreau, P.; Yeh, W.H.; Dumont, G.A.; Boivin, D.B. A circadian rhythm in heart rate variability contributes to the increased cardiac sympathovagal response to awakening in the morning. Chronobiol. Int. 2012, 29, 757–768. [Google Scholar] [CrossRef] [PubMed]
  23. Blume, C.; Lechinger, J.; Santhi, N.; del Giudice, R.; Gnjezda, M.T.; Pichler, G.; Scarpatetti, M.; Donis, J.; Michitsch, G.; Schabus, M. Significance of circadian rhythms in severely brain-injured patients A clue to consciousness? Neurology 2017, 88, 1933–1941. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Angerer, M.; Schabus, M.; Raml, M.; Pichler, G.; Kunz, A.B.; Scarpatetti, M.; Trinka, E.; Blume, C. Actigraphy in brain-injured patients—A valid measurement for assessing circadian rhythms? BMC Med. 2020, 18, 1–10. [Google Scholar] [CrossRef]
  25. Riganello, F.; Napoletano, G.; Cortese, M.D.; Arcuri, F.; Solano, A.; Lucca, L.F.; Tonin, P.; Soddu, A. What impact can hospitalization environment produce on the ANS functioning in patients with Unresponsive Wakefulness Syndrome? —24-hour monitoring. Brain Inj. 2019, 33, 1347–1353. [Google Scholar] [CrossRef]
  26. Czeisler, C.A.; Gooley, J.J. Sleep and circadian rhythms in humans. Cold Spring Harb. Symp. Quant. Biol. 2007, 72, 579–597. [Google Scholar] [CrossRef] [Green Version]
  27. Voss, A.; Schulz, S.; Schroeder, R.; Baumert, M.; Caminal, P. Methods derived from nonlinear dynamics for analysing heart rate variability. Philos. Trans. Math. Phys. Eng. Sci. 2009, 367, 277–296. [Google Scholar] [CrossRef]
  28. Sheng, H.; Chen, Y.Q.; Qiu, T. On the robustness of Hurst estimators. IET Signal Process. 2011, 5, 209–225. [Google Scholar] [CrossRef]
  29. Wislowska, M.; Del Giudice, R.; Lechinger, J.; Wielek, T.; Heib, D.P.J.; Pitiot, A.; Pichler, G.; Michitsch, G.; Donis, J.; Schabus, M. Night and day variations of sleep in patients with disorders of consciousness. Sci. Rep. 2017, 7, 266. [Google Scholar] [CrossRef]
  30. Angerer, M.; Pichler, G.; Angerer, B.; Scarpatetti, M.; Schabus, M.; Blume, C. From Dawn to Dusk—Mimicking Natural Daylight Exposure Improves Circadian Rhythm Entrainment in Patients with Severe Brain Injury, Sleep 2022, In press. Available online: https://osf.io/vk2d6/ (accessed on 29 January 2022).
  31. Blechert, J.; Peyk, P.; Liedlgruber, M.; Wilhelm, F.H. ANSLAB: Integrated multichannel peripheral biosignal processing in psychophysiological science. Behav. Res. Methods 2016, 48, 1528–1545. [Google Scholar] [CrossRef]
  32. Berntson, G.G.; Stowell, J.R. ECG artifacts and heart period variability: Don’t miss a beat! Psychophysiology 1998, 35, 127–132. [Google Scholar] [CrossRef]
  33. Wilhelm, F.H.; Grossman, P.; Roth, W.T. Assessment of heart rate variability during alterations in stress: Complex demodulation vs. spectral analysis. Biomed. Sci. Instrum. 2005, 41, 346–351. [Google Scholar]
  34. Hayano, J.; Taylor, J.A.; Yamada, A.; Mukai, S.; Hori, R.; Asakawa, T.; Yokoyama, K.; Watanabe, Y.; Takata, K.; Fujinami, T. Continuous Assessment of Hemodynamic Control by Complex Demodulation of Cardiovascular Variability. Am. J. Physiol. 1993, 264, H1229–H1238. [Google Scholar] [CrossRef]
  35. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2015. [Google Scholar]
  36. Asehnoune, K.; Roquilly, A.; Cinotti, R. Respiratory Management in Patients with Severe Brain. Inj. Crit. Care 2018, 22, 1–6. [Google Scholar] [CrossRef] [Green Version]
  37. Grossman, P.; Karemaker, J.; Wieling, W. Prediction of Tonic Parasympathetic Cardiac Control Using Respiratory Sinus Arrhythmia—The Need for Respiratory Control. Psychophysiology 1991, 28, 201–216. [Google Scholar] [CrossRef]
  38. Friedrich, S.; Konietschke, F.; Pauly, M. Resampling-Based Analysis of Multivariate Data and Repeated Measures Designs with the R Package MANOVA. RM. R J. 2019, 11, 380–400. [Google Scholar]
  39. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate—A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. B 1995, 57, 289–300. [Google Scholar] [CrossRef]
  40. Wasserstein, R.L.; Schirm, A.L.; Lazar, N.A. Moving to a World Beyond “p < 0.05”. Am. Stat. 2019, 73, 1–19. [Google Scholar] [CrossRef] [Green Version]
  41. Bürkner, P.C. brms: An R Package for Bayesian Multilevel Models Using Stan. J. Stat. Softw. 2017, 80, 1–28. [Google Scholar] [CrossRef] [Green Version]
  42. Carpenter, B.; Gelman, A.; Hoffman, M.D.; Lee, D.; Goodrich, B.; Betancourt, M.; Riddell, A.; Guo, J.Q.; Li, P.; Riddell, A. Stan: A Probabilistic Programming Language. J. Stat. Softw. 2017, 76, 1–29. [Google Scholar] [CrossRef] [Green Version]
  43. Cugini, P.; Curione, M.; Cammarota, C.; Bernardini, F.; Cipriani, D.; Rosa, R.D.; Francia, P.; Laurentis, T.D.; Marco, E.D.; Napoli, A.; et al. Is a Reduced Entropy in Heart Rate Variability an Early Finding of Silent Cardiac Neurovegetative Dysautonomia in Type 2 Diabetes Mellitus. J. Clin. Basic Cardiol. 2001, 4, 289–294. [Google Scholar]
  44. Byun, S.; Kim, A.Y.; Jang, E.H.; Kim, S.; Choi, K.W.; Yu, H.Y.; Jeon, H.J. Entropy analysis of heart rate variability and its application to recognize major depressive disorder: A pilot study. Technol. Health Care 2019, 27, S407–S424. [Google Scholar] [CrossRef] [Green Version]
  45. Lechinger, J.; Bothe, K.; Pichler, G.; Michitsch, G.; Donis, J.; Klimesch, W.; Schabus, M. CRS-R score in disorders of consciousness is strongly related to spectral EEG at rest. J. Neurol. 2013, 260, 2348–2356. [Google Scholar] [CrossRef]
  46. Batchinsky, A.L.; Cancio, L.C.; Salinas, J.; Kuusela, T.; Cooke, W.H.; Wang, J.J.; Boehme, M.; Convertino, V.A.; Holcomb, J.B. Prehospital loss of R-to-R interval complexity is associated with mortality in trauma patients. J.Trauma 2007, 63, 512–518. [Google Scholar] [CrossRef] [Green Version]
  47. Gao, L.; Smielewski, P.; Czosnyka, M.; Ercole, A. Cerebrovascular Signal Complexity Six Hours after Intensive Care Unit Admission Correlates with Outcome after Severe Traumatic Brain Injury. J. Neurotrauma 2016, 33, 2011–2018. [Google Scholar] [CrossRef] [Green Version]
  48. Bagnato, S.; Boccagni, C.; Sant’Angelo, A.; Fingelkurts, A.A.; Fingelkurts, A.A.; Galardi, G. Longitudinal Assessment of Clinical Signs of Recovery in Patients with Unresponsive Wakefulness Syndrome after Traumatic or Nontraumatic Brain Injury. J. Neurotrauma 2017, 34, 535–539. [Google Scholar] [CrossRef]
  49. Luaute, J.; Maucort-Boulch, D.; Tell, L.; Quelard, F.; Sarraf, T.; Iwaz, J.; Boisson, D.; Fischer, C. Long-term outcomes of chronic minimally conscious and vegetative states. Neurology 2010, 75, 246–252. [Google Scholar] [CrossRef]
  50. Katz, D.I.; Polyak, M.; Coughlan, D.; Nichols, M.; Roche, A. Natural history of recovery from brain injury after prolonged disorders of consciousness: Outcome of patients admitted to inpatient rehabilitation with 1–4 year follow-up. Prog. Brain Res. 2009, 177, 73–88. [Google Scholar] [CrossRef]
  51. Echeverria, J.C.; Woolfson, M.S.; Crowe, J.A.; Hayes-Gill, B.R.; Croaker, G.D.; Vyas, H. Interpretation of heart rate variability via detrended fluctuation analysis and alphabeta filter. Chaos 2003, 13, 467–475. [Google Scholar] [CrossRef]
  52. Goldberger, A.L. Is the Normal Heartbeat Chaotic or Homeostatic. News Physiol. Sci. 1991, 6, 87–91. [Google Scholar] [CrossRef] [Green Version]
  53. Richman, J.S.; Moorman, J.R. Physiological time-series analysis using approximate entropy and sample entropy. Am. J. Physiol. Heart C 2000, 278, H2039–H2049. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Brainsci 12 00375 g001 550
Figure 1. Interbeat interval (IBI) separately for time and diagnosis contrasts. (a) While patients’ IBIs were larger during the night as compared to the day (i.e., forenoon, afternoon), they did not differ between fore- and afternoon. (b) Patients in UWS showed larger IBIs than patients in (E)MCS. Error bars represent the mean and 95% confidence interval. * p < 0.05; + p ≤ 0.1; ns = not significant. Abbreviations: (E)MCS = (emergence from) minimally conscious state; UWS = unresponsive wakefulness syndrome; ms = milliseconds.
Figure 1. Interbeat interval (IBI) separately for time and diagnosis contrasts. (a) While patients’ IBIs were larger during the night as compared to the day (i.e., forenoon, afternoon), they did not differ between fore- and afternoon. (b) Patients in UWS showed larger IBIs than patients in (E)MCS. Error bars represent the mean and 95% confidence interval. * p < 0.05; + p ≤ 0.1; ns = not significant. Abbreviations: (E)MCS = (emergence from) minimally conscious state; UWS = unresponsive wakefulness syndrome; ms = milliseconds.
Brainsci 12 00375 g001
Brainsci 12 00375 g002 550
Figure 2. Correlation between interbeat interval (IBI) and CRS-R sum score separately for time. Larger IBIs were associated with lower CRS-R sum scores throughout the (a,b) day (i.e., forenoon, afternoon) and (c) night. Please note that the effect during the night is no longer ‘credible’ from a Bayesian point of view, and will not be interpreted. ** p < 0.01; * p < 0.05. Abbreviations: (E)MCS = (emergence from) minimally conscious state; UWS = unresponsive wakefulness syndrome; ms = milliseconds.
Figure 2. Correlation between interbeat interval (IBI) and CRS-R sum score separately for time. Larger IBIs were associated with lower CRS-R sum scores throughout the (a,b) day (i.e., forenoon, afternoon) and (c) night. Please note that the effect during the night is no longer ‘credible’ from a Bayesian point of view, and will not be interpreted. ** p < 0.01; * p < 0.05. Abbreviations: (E)MCS = (emergence from) minimally conscious state; UWS = unresponsive wakefulness syndrome; ms = milliseconds.
Brainsci 12 00375 g002
Brainsci 12 00375 g003 550
Figure 3. Root mean square of successive differences between adjacent heartbeats (RMSSD). Patients’ RMSSD did not differ between time ranges and diagnoses. Error bars represent the mean and 95% confidence interval. Abbreviations: (E)MCS = (emergence from) minimally conscious state; UWS = unresponsive wakefulness syndrome; ms = milliseconds.
Figure 3. Root mean square of successive differences between adjacent heartbeats (RMSSD). Patients’ RMSSD did not differ between time ranges and diagnoses. Error bars represent the mean and 95% confidence interval. Abbreviations: (E)MCS = (emergence from) minimally conscious state; UWS = unresponsive wakefulness syndrome; ms = milliseconds.
Brainsci 12 00375 g003
Brainsci 12 00375 g004 550
Figure 4. Very low (VLF), low (LF), and high frequency (HF) oscillation amplitude of the IBI signal computed by complex demodulation (CDM). Analyses of the (a) VLF and (b) LF band revealed significant main effects for time (i.e., forenoon, afternoon, night). However, post hoc comparisons between times of day did not yield significance anymore after correction for multiple comparisons. (c) Patients’ HF also did not differ between times and diagnoses. Error bars represent the mean and 95% confidence interval. # = significant before correction for multiple comparisons; ns = not significant. Abbreviations: (E)MCS = (emergence from) minimally conscious state; UWS = unresponsive wakefulness syndrome; ms = milliseconds.
Figure 4. Very low (VLF), low (LF), and high frequency (HF) oscillation amplitude of the IBI signal computed by complex demodulation (CDM). Analyses of the (a) VLF and (b) LF band revealed significant main effects for time (i.e., forenoon, afternoon, night). However, post hoc comparisons between times of day did not yield significance anymore after correction for multiple comparisons. (c) Patients’ HF also did not differ between times and diagnoses. Error bars represent the mean and 95% confidence interval. # = significant before correction for multiple comparisons; ns = not significant. Abbreviations: (E)MCS = (emergence from) minimally conscious state; UWS = unresponsive wakefulness syndrome; ms = milliseconds.
Brainsci 12 00375 g004
Brainsci 12 00375 g005 550
Figure 5. Approximate entropy (ApEn). Patients showed a higher ApEn during forenoon as compared to afternoon. No significant differences were evident between night and day (i.e., night vs. forenoon or afternoon). Error bars represent the mean and 95% confidence interval. *** p < 0.001; ns = not significant. Abbreviations: (E)MCS = (emergence from) minimally conscious state; UWS = unresponsive wakefulness syndrome; arb. units = arbitrary units.
Figure 5. Approximate entropy (ApEn). Patients showed a higher ApEn during forenoon as compared to afternoon. No significant differences were evident between night and day (i.e., night vs. forenoon or afternoon). Error bars represent the mean and 95% confidence interval. *** p < 0.001; ns = not significant. Abbreviations: (E)MCS = (emergence from) minimally conscious state; UWS = unresponsive wakefulness syndrome; arb. units = arbitrary units.
Brainsci 12 00375 g005
Brainsci 12 00375 g006 550
Figure 6. Detrended fluctuation analysis scaling exponent (DfaAlpha) and sample entropy 1 (SampEn1). Patients with TBI showed a higher (a) DfaAlpha and (b) SampEn1 as compared to patients with NTBI. Error bars represent the mean and 95% confidence interval. * p < 0.05. Abbreviations: NTBI = non-traumatic brain injury; TBI = traumatic brain injury; arb. units = arbitrary units.
Figure 6. Detrended fluctuation analysis scaling exponent (DfaAlpha) and sample entropy 1 (SampEn1). Patients with TBI showed a higher (a) DfaAlpha and (b) SampEn1 as compared to patients with NTBI. Error bars represent the mean and 95% confidence interval. * p < 0.05. Abbreviations: NTBI = non-traumatic brain injury; TBI = traumatic brain injury; arb. units = arbitrary units.
Brainsci 12 00375 g006
Brainsci 12 00375 g007 550
Figure 7. Correlation between EEG permutation entropy (PE) and HRV approximate entropy (ApEn) separately for day and night. (a) While EEG PE and HRV ApEn did not correlate during the day (i.e., 8 a.m.–8 p.m.), (b) a higher EEG PE was associated with a higher HRV ApEn during the night (i.e., 11 p.m.–5 a.m.). * p < 0.05; ns = not significant. Abbreviations: (E)MCS = (emergence from) minimally conscious state; UWS = unresponsive wakefulness syndrome; arb. units = arbitrary units.
Figure 7. Correlation between EEG permutation entropy (PE) and HRV approximate entropy (ApEn) separately for day and night. (a) While EEG PE and HRV ApEn did not correlate during the day (i.e., 8 a.m.–8 p.m.), (b) a higher EEG PE was associated with a higher HRV ApEn during the night (i.e., 11 p.m.–5 a.m.). * p < 0.05; ns = not significant. Abbreviations: (E)MCS = (emergence from) minimally conscious state; UWS = unresponsive wakefulness syndrome; arb. units = arbitrary units.
Brainsci 12 00375 g007
Table
Table 1. Demographic information and highest CRS-R sum score/diagnosis.
Table 1. Demographic information and highest CRS-R sum score/diagnosis.
Patient IDAgeSexEtiologyTime Since
Injury (Months)		DiagnosisCRS-R Sum Score
P164MNTBI8.8UWS3
P277MNTBI15.7UWS3
P335MNTBI14.6UWS4
P471MNTBI20.7EMCS21
P555FNTBI23.8UWS4
P659MNTBI25.8UWS4
P780MTBI21.9MCS9
P822FNTBI48.3UWS-
P948MTBI7.7UWS5
P1076MTBI3.5UWS4
P1170MNTBI3.4UWS2
P1271FTBI4.7UWS4
P1316FTBI0.8MCS16
P1421MTBI7.0UWS6
P1548MNTBI1.4MCS18
P1666MNTBI3.2MCS10
P1761MNTBI2.0MCS10
P1836MTBI6.0MCS6
P1974FTBI0.5UWS3
P2031FNTBI1.4MCS11
P2143FTBI6.0MCS21
P2261FNTBI0.9UWS6
P2337MNTBI9.4UWS5
P2434MNTBI240.0UWS8
P2532MTBI6.1UWS5
P2620MTBI36.0MCS 13
M = male; F = female; NTBI = non-traumatic brain injury; TBI = traumatic brain injury; UWS = unresponsive wakefulness syndrome; MCS = minimally conscious state; EMCS = emergence from MCS; CRS-R = Coma Recovery Scale-Revised (cf. Behavioral Assessment and Data Analysis). Please note that we could not obtain valid CRS-R assessments in one patient (P8) because not all subscales could be evaluated (i.e., due to eyes being closed and it being impossible to induce eye-opening even when physically stimulating the patient). In this case, we used the diagnosis that was obtained during a CRS-R assessment ten days earlier. While the diagnosis usually stays quite stable over time, CRS-R scores can slightly change. Thus, we did not use the CRS-R sum score of this earlier assessment for analyses.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Angerer, M.; Wilhelm, F.H.; Liedlgruber, M.; Pichler, G.; Angerer, B.; Scarpatetti, M.; Blume, C.; Schabus, M. Does the Heart Fall Asleep?—Diurnal Variations in Heart Rate Variability in Patients with Disorders of Consciousness. Brain Sci. 2022, 12, 375. https://doi.org/10.3390/brainsci12030375

AMA Style

Angerer M, Wilhelm FH, Liedlgruber M, Pichler G, Angerer B, Scarpatetti M, Blume C, Schabus M. Does the Heart Fall Asleep?—Diurnal Variations in Heart Rate Variability in Patients with Disorders of Consciousness. Brain Sciences. 2022; 12(3):375. https://doi.org/10.3390/brainsci12030375

Chicago/Turabian Style

Angerer, Monika, Frank H. Wilhelm, Michael Liedlgruber, Gerald Pichler, Birgit Angerer, Monika Scarpatetti, Christine Blume, and Manuel Schabus. 2022. "Does the Heart Fall Asleep?—Diurnal Variations in Heart Rate Variability in Patients with Disorders of Consciousness" Brain Sciences 12, no. 3: 375. https://doi.org/10.3390/brainsci12030375

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop