Continuous Multi-Parameter Heart Rate Variability Analysis Heralds Onset of Sepsis in Adults

Abstract

Background

Early diagnosis of sepsis enables timely resuscitation and antibiotics and prevents subsequent morbidity and mortality. Clinical approaches relying on point-in-time analysis of vital signs or lab values are often insensitive, non-specific and late diagnostic markers of sepsis. Exploring otherwise hidden information within intervals-in-time, heart rate variability (HRV) has been documented to be both altered in the presence of sepsis, and correlated with its severity. We hypothesized that by continuously tracking individual patient HRV over time in patients as they develop sepsis, we would demonstrate reduced HRV in association with the onset of sepsis.

Methodology/Principal Findings

We monitored heart rate continuously in adult bone marrow transplant (BMT) patients (n = 21) beginning a day before their BMT and continuing until recovery or withdrawal (12±4 days). We characterized HRV continuously over time with a panel of time, frequency, complexity, and scale-invariant domain techniques. We defined baseline HRV as mean variability for the first 24 h of monitoring and studied individual and population average percentage change (from baseline) over time in diverse HRV metrics, in comparison with the time of clinical diagnosis and treatment of sepsis (defined as systemic inflammatory response syndrome along with clinically suspected infection requiring treatment). Of the 21 patients enrolled, 4 patients withdrew, leaving 17 patients who completed the study. Fourteen patients developed sepsis requiring antibiotic therapy, whereas 3 did not. On average, for 12 out of 14 infected patients, a significant (25%) reduction prior to the clinical diagnosis and treatment of sepsis was observed in standard deviation, root mean square successive difference, sample and multiscale entropy, fast Fourier transform, detrended fluctuation analysis, and wavelet variability metrics. For infected patients (n = 14), wavelet HRV demonstrated a 25% drop from baseline 35 h prior to sepsis on average. For 3 out of 3 non-infected patients, all measures, except root mean square successive difference and entropy, showed no significant reduction. Significant correlation was present amongst these HRV metrics for the entire population.

Conclusions/Significance

Continuous HRV monitoring is feasible in ambulatory patients, demonstrates significant HRV alteration in individual patients in association with, and prior to clinical diagnosis and treatment of sepsis, and merits further investigation as a means of providing early warning of sepsis.

About the Authors:

Saif Ahmad

Affiliation: Ottawa Hospital Research Institute, Ottawa, Ontario, Canada

Tim Ramsay

Affiliation: Ottawa Hospital Research Institute, Ottawa, Ontario, Canada

Lothar Huebsch

Affiliation: Department of Clinical Hematology, University of Ottawa, Ottawa, Ontario, Canada

Sarah Flanagan

Affiliation: Ottawa Hospital Research Institute, Ottawa, Ontario, Canada

Sheryl McDiarmid

Affiliation: Department of Clinical Hematology, University of Ottawa, Ottawa, Ontario, Canada

Izmail Batkin

Affiliation: Ottawa Hospital Research Institute, Ottawa, Ontario, Canada

Lauralyn McIntyre

Affiliations Ottawa Hospital Research Institute, Ottawa, Ontario, Canada, Department of Critical Care Medicine, University of Ottawa, Ottawa, Ontario, Canada

Sudhir R. Sundaresan

Affiliation: Division of Thoracic Surgery, University of Ottawa, Ottawa, Ontario, Canada

Donna E. Maziak

Affiliation: Division of Thoracic Surgery, University of Ottawa, Ottawa, Ontario, Canada

Farid M. Shamji

Affiliation: Division of Thoracic Surgery, University of Ottawa, Ottawa, Ontario, Canada

Paul Hebert

Affiliations Ottawa Hospital Research Institute, Ottawa, Ontario, Canada, Department of Critical Care Medicine, University of Ottawa, Ottawa, Ontario, Canada

Dean Fergusson

Affiliation: Ottawa Hospital Research Institute, Ottawa, Ontario, Canada

Alan Tinmouth

Affiliations Ottawa Hospital Research Institute, Ottawa, Ontario, Canada, Department of Clinical Hematology, University of Ottawa, Ottawa, Ontario, Canada

Andrew J. E. Seely

* E-mail: [email protected]

Affiliations Ottawa Hospital Research Institute, Ottawa, Ontario, Canada, Division of Thoracic Surgery, University of Ottawa, Ottawa, Ontario, Canada, Department of Critical Care Medicine, University of Ottawa, Ottawa, Ontario, Canada

Introduction

Severe sepsis and septic shock are major causes of mortality and costs in critically ill patients [1]–[2]. Acute neutropenia is a frequent and intended iatrogenic side effect of cytotoxic chemotherapy and radiotherapy, commonly employed in the management of malignant hematological diseases, most commonly leukemia and lymphoma, leading to increasing risk of opportunistic infections and sepsis [3]. This is particularly apparent in individuals undergoing induction chemotherapy and bone marrow transplantation (BMT) for acute leukemia due to severe prolonged neutropenia [4]. Thus, patients undergoing bone marrow transplantation comprise a group which is at a high risk of systemic infection (approximately 80%) and mortality (approximately 5%) [5]–[6].

Results of a large single-centre randomized trial [7]–[8] indicate that, in early septic shock, mortality can be greatly reduced by employing early goal directed therapy (EGT). In addition, early effective antibiotic therapy is essential to minimize mortality secondary to septic shock [9]. Thus, early diagnosis of sepsis leading to aggressive resuscitation involving antibiotic administration, and maintaining adequate systemic oxygen delivery and tissue perfusion is vital to survival. Current clinical approaches for diagnosing sepsis are based on an increased absolute value of one or more vital signs in addition to laboratory tests such as blood cultures to look for evidence of a pathogen. Although this is the best clinical approach currently available, it represents a crude and potentially non-sensitive and non-specific means for early diagnosis of sepsis.

An alternative approach is to consider the host response to sepsis as a complex system, and utilize variability analysis as a means to better characterize the system [10]. Evidence from experimental and clinical studies indicates that heart rate variability (HRV) is altered in the presence of systemic infection [11]–[13] and correlates with its severity [14]–[15]. Investigators studying the prognostic value of HRV have shown it to be an early predictor of death [16] and multiple organ dysfunction syndrome (MODS) [17] in adult septic patients. Studies on infants and neonates have shown HRV to be predictive [18]–[19] and diagnostic [20] of sepsis.

For adult patients, prior studies have only assessed HRV for short intermittent epochs of 5–15 minutes. To our knowledge, there have been no published attempts to accomplish continuous HRV analysis in adult patients at high risk of sepsis, and no comprehensive analysis of HRV has been performed during the process of development of sepsis. The method that we propose examines a panel of HRV metrics computed continuously over time prior to and during the clinical diagnosis of sepsis (defined as systemic inflammatory response syndrome along with clinically suspected infection requiring treatment) in patients at high risk for systemic infection, namely neutropenic patients following bone marrow transplant. In this pilot study, we demonstrate that (a) it is feasible to enroll and perform continuous HRV analysis of prospectively collected non-stationary heart rate (HR) data in ambulatory patients, (b) continuous variability analysis demonstrates a reduction in individual patient HRV in association with the presence of sepsis, and (c) the alteration in HRV occurs within a clinically relevant period prior to when standard clinical measures lead to clinical diagnosis and treatment.

Methods

Ethics Statement

Written informed consent was obtained from all participants and the Ottawa Hospital Research Ethics Board authorized the study.

Design

This is a descriptive study in which prospective continuous HR recording, and retrospective analysis of both HRV and change in HRV were analyzed over time in ambulatory outpatients as they underwent BMT, and developed neutropenia associated with high risk for developing systemic infection.

Participants

Our subjects included patients (n = 21) undergoing BMT for hematological malignancy or other disorders (see Table 1) at the Ottawa Hospital (General Campus), Ontario, Canada. All patients were enrolled between 05/2007 to 05/2008. Inclusion criteria included treatment with myeloablative chemoradiotherapy followed by an allogeneic or autologous BMT, and informed consent. Exclusion criteria were pre-existing cardiopulmonary disease, taking beta-blockers or calcium-channel blockers, pre-existing arrhythmia (e.g. atrial fibrillation, atrial bigeminy), contraindication to electrocardiogram adhesives (e.g. allergy, severe psoriasis), and not being fluent in English or French (which precluded patients from understanding and signing the informed consent).

[Figure omitted. See PDF.]

Table 1. Study and patient characteristics, along with reason and type of BMT.

https://doi.org/10.1371/journal.pone.0006642.t001

Heart Rate (HR) Monitoring

We collected continuous Holter ECG data (average 12 [SD 4] days of HR monitoring) for all patients in the study, starting approximately 24 h before their BMT and continuing through neutropenia until its resolution or until withdrawal from the study. 3M Cavilon spray combined with daily replacement and relocation of gel ECG electrodes was utilized to reduce any skin irritation. We used a Zymed DigiTrak-Plus Holter system (Philips Healthcare, Markham, Ontario, Canada), which annotated all normal QRS peaks and arrhythmias, including premature atrial and ventricular beats. Only the beats that characterized normal sinus rhythm (NSR) were included, while all premature beats were excluded. RR intervals were derived from R wave annotations. Thus, for each patient, the input to our signal processing algorithms comprised time-stamped NSR-based RR intervals in seconds.

Feasibility

The feasibility measures evaluated included the number of patients dropping out of the study, the ability of patients to maintain a diary of clinical events (for example, temperature, diarrhea, vomiting, etc.), and compliance with Holter HR monitoring during enrollment. A patient was considered to be a dropout if he or she enrolled in the study, and subsequently discontinued Holter HR monitoring (i.e. due to discomfort, or other factors). A minimum of 72 h of monitoring was necessary to include patients in the study. We examined patient diaries for completeness and regularity in reporting clinical events as per given schedule and instructions. To quantify compliance with monitoring, we computed a ratio between the total time for which patients were monitored and the total time of lost data during the monitoring period. More specifically, percentage of data lost was computed as DL = 100*[TL/(TE−TS)], where DL is percentage of data lost, TS is monitoring start time, TE is monitoring end time, and TL is total time of lost data between TS and TE.

Diagnosis of Sepsis

In this study, sepsis was defined as systemic inflammatory response syndrome along with clinically suspected infection requiring treatment. Patients who were administered broad spectrum antibiotics by the physician (attending hematologist) were considered to have been clinically diagnosed with sepsis. Since exact time of diagnosis was uncertain, we assumed the time (hour and minute) of first administration of antibiotics as the time of clinical diagnosis of sepsis for consistency in procedure. It was not possible to identify precisely the actual time of physician diagnosis. Review of cases and procedures indicated that the diagnosis preceded the first administration of antibiotics by 30 to 90 minutes. The clinical indication to treat sepsis and the clinical diagnosis of sepsis are included in Table 2. Over 50% of the patients had sepsis diagnosed based on the presence of fever, defined a priori as one recording greater than 38.5 degrees centigrade or two recordings greater than 38.0 degrees centigrade within 12 h.

[Figure omitted. See PDF.]

Table 2. Indication for antibiotics and bacteriological diagnosis.

https://doi.org/10.1371/journal.pone.0006642.t002

Variability Analysis

We have developed a novel system for continuous individualized multiorgan variability analysis (CIMVA), with single-organ (the heart) application in this study. The CIMVA system (developed in Windows® Matlab®) comprises algorithms for computing and visualizing diverse measures of HRV. Details of diverse HRV measures computed by CIMVA are provided elsewhere [21], but in Table 3 we present a brief summary.

[Figure omitted. See PDF.]

Table 3. Summary of signal analysis capabilities of CIMVA.

https://doi.org/10.1371/journal.pone.0006642.t003

To accomplish continuous variability analysis over time, CIMVA employs a moving window approach, whereby a window (interval-in-time) of user specified width and step marches through the input signal, computing and time-stamping different variability metrics at each step, thus making it possible to monitor a change in HRV over time. A standard RR cleaning algorithm [22] is employed inside each window to detect gross artifact or noise, and HRV analysis is performed on the cleaned data. The cleaning algorithm excludes RR intervals less than 0.25 s and greater than 2.5 s, as well as those that differ by more than 15% from the previous one. The CIMVA system stores the number of samples lost due to RR cleaning in each window instance, thus keeping track of signal quality.

We employed a window width of 1200 samples (∼10 minutes) and step of 200 samples (∼2 minutes) to compute HRV over time. In Table 4, we present an example of the output from the CIMVA system.

[Figure omitted. See PDF.]

Table 4. Example of first five window instances of a CIMVA generated change in HRV over time matrix for diverse measures for a moving window width of 1200 samples and step of 200 samples.

https://doi.org/10.1371/journal.pone.0006642.t004

Subsequently, the variability time series (HRV over time) underwent a second data integration step or “smoothing” process. That is, the change in HRV over time was smoothed utilizing a moving average and hourly data interpolation technique. This smoothing approach stabilized the change in variability over time to characterize long-term behavior of various variability metrics and facilitated the study of lead times (ΔTP) for percentage changes (ΔHRV) from baseline HRV.

Continuous Analysis

We defined baseline variability as the mean variability for the first 24 h of recording, prior to transplantation. Given expected variation in baseline variability from patient to patient, for all measures of HRV, we computed an individualized percentage change in variability with respect to the baseline variability (ΔHRV = 100 * [baseline-current]/baseline). We studied the extent and timing of ΔHRV achieved by various metrics for infected and non-infected patients. To this end, we evaluated individual and between-patient average ΔHRV data. We defined a priori a significant ΔHRV as a 25% loss of variability. We also computed mean population correlation amongst all diverse HRV measures for the entire monitoring period.

Results

Monitoring was initiated in 21 patients. Four patients dropped out within 24 h of initiation of monitoring due to discomfort or other reasons, leaving 17 (81%) datasets for analysis. We observed that the requirement for patients to keep a diary of clinical events on a regular basis was difficult. We found general errors in the diaries. These included the inability to record the correct time of the occurrence of an event, and reporting inaccurate frequency ranges of occurrence of events like vomiting and loose motions. For n = 17 patients who completed the study, for a total monitoring time of 188 days, 12 days (6%) of data was lost due to non-compliance to monitoring. The reasons for non-compliance to monitoring included taking off the Holter due to discomfort or in order to take a bath, or arriving late at the hospital which led to a delay in attaching a new Holter after the old one had reached its recording limit. In terms of data quality, for n = 17, for 176 days (monitoring days - lost days of data = 188−12 = 176) of available RR-interval data, the data lost to RR cleaning was 1.06 days (0.6%).

In the course of a median follow-up of 12 (IQR 9-14, SD 4) days for n = 17, we observed that 14 patients were diagnosed and treated for sepsis, whereas 3 were not. All patients recovered; no patient was admitted to the ICU or died. Out of the 14 infected patients, on an average, 12 (86%) showed a significant (25%) drop (ΔHRV) prior to sepsis for standard deviation (SD), root mean square successive difference (RMSSD), sample entropy (SampEn), multiscale entropy (MSE) area under curve (AUC), fast Fourier transform (FFT) low frequency (LF), FFT high frequency (HF), detrended fluctuation analysis (DFA) AUC, and wavelet spectral density AUC measures of variability. In addition, for 3 out of 3 non-infected patients, all measures, except RMSSD, SampEn and MSE, showed no significant reduction. Notably, these measures also showed a substantial correlation (>+30%) amongst each other (Table 5).

[Figure omitted. See PDF.]

Table 5. Mean population (n = 17) correlation amongst diverse HRV measures. Correlations greater than +30% are bolded and shown in parenthesis.

https://doi.org/10.1371/journal.pone.0006642.t005

Figure 1 shows the difference in individual patterns of variability for infected (n = 14, red plots) and non-infected (n = 3, green plots) patients for the wavelet measure (AUC) of ΔHRV. In addition, Figure 1 shows the timing of a 25% drop in wavelet ΔHRV (dashed vertical line and dot) in comparison to the time of sepsis (solid vertical line for infected patients) for all patients. Except two patients (Patient #s 4 & 16) out of the 14 infected (n = 14, red plots), all show a 25% drop in ΔHRV in the neighborhood of sepsis. For all non-infected patients (n = 3, green plots), no 25% drop was observed in wavelet ΔHRV. In Figure 1, a positive lead time ΔTP means the 25% drop occurred prior to sepsis, a negative ΔTP means it occurred after sepsis, and ΔTP = NaN (not a number) means no 25% drop occurred. Table 6 summarizes the lead time (ΔTP) for all patients (n = 17) studied for 25% drop in ΔHRV for SD, RMSSD, power law intercept, FFT HF, FFT LF, DFA, wavelet, SampEn, and MSE metrics of HRV. We note that except for power law analysis, the 25% drop in ΔHRV occurs prior to sepsis (positive ΔTP) for majority of the infected patients.

[Figure omitted. See PDF.]

Figure 1. Individual smoothed percentage change.

((HRV) in wavelet variability. Individual wavelet (HRV is studied every 1 h in infected (n = 14) and non-infected (n = 3) patients from 0 h (time after baseline variability) up to the end of the study. Green plots are non-infected patients whereas red plots are infected patients. Solid vertical line denotes the time of sepsis. Lead time (TP is studied for 25% drop (dashed vertical line and dot) from baseline (first 24 h) HRV.

https://doi.org/10.1371/journal.pone.0006642.g001

[Figure omitted. See PDF.]

Table 6. Lead time ΔTP (hours) at 25% drop from baseline (first 24 h) HRV.

https://doi.org/10.1371/journal.pone.0006642.t006

Figure 2 shows the difference in average patterns of variability every 12 h for the infected (n = 14) and non-infected (n = 3) population for five diverse HRV measures. For these measures, the infected population (plotted in red) shows a distinct and continuous drop in variability from 24–120 h whereas no such drop is observed for the non-infected population (plotted in green). Figure 3 shows the mean variation every 6 h in the above diverse HRV measures zoomed to (72 h around sepsis (time = 0 h) for infected patients (n = 14). Plots in Figure 2 and Figure 3 show a visually evident drop in diverse HRV measures before sepsis; led by entropy measures, followed by standard deviation, wavelet, and power law measures of HRV.

[Figure omitted. See PDF.]

Figure 2. Average percentage change (ΔHRV) in multi-parameter variability.

Average multi-parameter (HRV is studied every 12 h in infected (n = 14) and non-infected (n = 3) patients from 0 h (time after baseline variability) up to the end of the study. Green plots are non-infected patients whereas red plots are infected patients. The error bars represent the standard error mean (SEM). The grey horizontal bar at the bottom represents the range of hours where sepsis occurred for all infected patients. The blue horizontal line at the bottom represents the mean time ±standard deviation of sepsis for all infected patients.

https://doi.org/10.1371/journal.pone.0006642.g002

[Figure omitted. See PDF.]

Figure 3. Average percentage change (ΔHRV) in multi-parameter variability around sepsis.

Average multi-parameter ΔHRV is studied every 6 h in infected (n = 14) patients±72 h around sepsis (0 h). The error bars represent the standard error mean (SEM).

https://doi.org/10.1371/journal.pone.0006642.g003

Discussion

This study represents a pilot investigation of the methods and clinical relevance of tracking HRV continuously in high-risk patients prior to, during, and following clinical diagnosis and treatment for sepsis. We observed a significant loss (defined as a 25% reduction) in multi-parameter HRV in 12/14 or 86% of patients treated for sepsis. None of the patients (0/3) without sepsis demonstrated a significant loss in multi-parameter HRV. Evaluating individual patients as their own control, we observed that the occurrence of loss of HRV was temporally associated with clinical diagnosis and treatment of sepsis in individual patients. The onset of drop of HRV occurred prior to clinical diagnosis and treatment of sepsis, and importantly, recovery of HRV was observed as the patients improved (see Figures 1, 2 & 3). Although this study must be interpreted in light of the small sample size and observational nature, the temporal association demonstrated in this study utilizing continuous variability analysis is in keeping with past investigations utilizing intermittent variability analysis where reduced HRV was demonstrated in the presence of sepsis in adults [11]–[15] and infants [18]–[20].

For infected patients, we note that ΔTP values greater than 120 h are probably “false alarms” where the 25% drop does occur momentarily but is not followed by further sustained drops in ΔHRV. For example for Patient #16, the standard deviation metric shows a ΔTP of 237 h (Table 6). However, the variability pattern of standard deviation (not shown) is very similar to that of the wavelet metric for Patient #16 (Figure 1). This means that such drops in ΔHRV cannot be used for prediction of sepsis with appropriate confidence levels. That is, for infected patients, the drop in ΔHRV should be continuous and progressive for the prediction of sepsis to have more reliability and robustness. We found that amongst all metrics studied for prediction of sepsis, the wavelet analysis showed the highest reliability and stability whereby ΔTP was within a reasonable timeframe of sepsis, and “false alarms” were minimal with regard to both infected and non-infected patients (Figure 1 and Table 6). The main reason for the success of wavelet analysis seems to be the fact that it is a recursive filtering technique which is capable of analyzing a signal in both time and frequency domains simultaneously. These attributes make the wavelet analysis robust to both noise and non-stationarity in the data, thus allowing it to characterize signal variability in the most efficient and accurate manner.

While other studies have clearly documented population-based changes in HRV in association with sepsis, this study identified individual patient loss of HRV equal to 25% in association with, and prior to, clinical diagnosis and treatment of sepsis in all but two patients. Of note, the patient (Patient #16, Figure 1) who showed no drop in HRV but was diagnosed and treated for sepsis was a patient with insulin dependent diabetes. This finding requires further study, in particular since existing literature suggests that HRV may not be an effective tool for characterizing the autonomic nervous system control in patients with diabetes [23]–[24].

Although this investigation supports the feasibility of long term HR monitoring (up to 16 days) along with continuous HRV analysis in non-stationary HR data sets in ambulatory patients, the use of Holter monitors with gel electrodes proved demanding for patients recovering from BMT. A total of 4 patients out of 21 (19%) dropped out prematurely after initiating the study. For the 17 patients who completed the study, the Holter data lost due to non-compliance to monitoring was 6%. We note that some of the reasons for non-compliance were beyond the control of the patients. For example, if a patient had to bathe or was in discomfort, the Holter had to be taken off and sometimes a particular schedule for arrival at the hospital (resulting in delay in new Holter attachment and data loss) was based on bed availability, physician availability, and other similar factors. The quality of the Holter data recorded was clean and satisfactory; the data lost to RR cleaning was only 1.06 days (0.6%) in a total of 176 days of RR-interval data. We also visually inspected the recorded RR-interval data and found it to be satisfactory. We encountered difficulty in motivating patients to accurately and punctually maintain a diary of clinical events. This precluded us from rigorously correlating our HRV measures with other clinical data such as periodic temperature measurements, symptoms such as frequency ranges of diarrhea and vomiting, etc.

In 2 out of the 3 normal patients, relatively short Holter monitoring periods of 3 and 5 days, were observed. This was because these patients withdrew from the study after these timeframes. However, these patients did not develop sepsis during the monitoring period. Thus, they were considered as uninfected patients for the study.

In order to visually observe and computationally detect changes in HRV over time, we have pioneered a process of double temporal integration. The first step of this process involves computing a moving window based change in HRV over time. The second step involves smoothing of the change in HRV over time utilizing the moving average and data interpolation techniques. We observed that this process was instrumental in highlighting and studying stable and progressive changes in HRV hours prior to the onset of sepsis.

In addition to the general process of continuous variability analysis, we had to evolve the process of continuous individualized monitoring of individual variability metrics, where no guidelines exist. For example, in the performance of iteratively repeated sample and multi-scale entropy analysis, it remains an unanswered question as to whether the tolerance r, calculated as 15% of the standard deviation (SD), should be fixed, or continuously re-calculated based on the changing SD for every window instance. We empirically observed marked reduction in both entropy metrics by keeping the tolerance r fixed as 15% of the SD of the first 24 h (baseline) of HR data [r = 0.15*SD(HRbaseline)] as patients developed sepsis. However, re-calculating the entropy metrics based upon a variable value for r (a new r based on SD of each window instance), the reduction in entropy was not evident. This highlights the inter-dependence of the entropy metrics and the time-domain SD metric and requires further theoretical and empirical study.

We also observed that the DFA AUC and wavelet spectral density AUC metrics showed identical HRV patterns when computed continuously over time (99% correlation between DFA AUC and wavelet spectral density AUC – see Table 5). Our empirical observation and the correlation analysis highlight the interdependence between multiple measures of HRV, and the unresolved debate regarding the most efficient and effective means to characterize HRV.

The pathophysiology of altered HRV still remains uncertain, and was not the focus of this study. The loss of overall variation and complexity metrics are in keeping with existing hypotheses regarding augmentation of sympathetic tone, “decomplexification” of physiologic signals [25], as well as “uncoupling of biologic oscillators” [26] that have been proposed to occur in association with sepsis and critical illness. Regardless of the possible mechanisms, we interpret the alteration in HRV occurring in association with sepsis as reflective of a fundamental alteration to the underlying complex system producing the dynamics.

Although this study was not designed for a comprehensive lead-time analysis, we observed changes in HRV at least 24 h prior to the clinical diagnosis and treatment of sepsis. This analysis was limited by the fact that we were unable to identify the precise time of onset of sepsis, and rather identified the onset of antibiotic treatment as a temporal point for comparing the results of the variability analyses. Half of the patients were diagnosed with sepsis based upon detection of a fever 60–90 minutes prior to administration of antibiotics. Thus, intermittent sampling of temperature (about every 6 h) may have delayed diagnosis by several hours. Of note, we investigated the technology to monitor skin temperature continuously, and found it imprecise and unreliable. However, experimental models of sepsis have noted changes in HRV at least 6 h prior to the onset of fever [27]. Following these results, it will be worthwhile to correlate the diverse measures of HRV with circulating biomarkers of sepsis and inflammation. Our study group comprised ambulatory outpatients recovering from myeloablative therapy, who were already encumbered by continuous HR monitoring, and frequent blood sampling was not perceived as ethical in absence of these data.

Numerous unanswered questions remain. Although the clinical significance of this investigation lies in its logical extension, namely the potential of utilizing a system of prospective continuous HRV monitoring to provide an early warning system for sepsis, a larger validation cohort is clearly required, along with definitive investigations regarding generalizability, sensitivity and specificity. It is entirely unclear if hospitalized patients (e.g. surgical or critically ill patients) will demonstrate such profound loss of HRV in association with new onset of sepsis. Although a panel of HRV metrics has been recommended to optimally characterize variability, it remains to be seen which offer the optimal means to herald the presence of sepsis. Last, a more in depth evaluation of lead-time offered by HRV monitoring is warranted. As evident above from the discussion regarding double temporal integration, a process of prospective continuous HRV monitoring will necessarily lead to a delay inherent to the interval of time necessary to characterize variability and thus the underlying system. Akin to the uncertainty principle, there is a balance between the need for an adequate interval of time to characterize the system, and the need to characterize the system at any point in time. Optimization of this balance may depend on clinical applications. Nonetheless, the promise of a novel dimension of monitoring based on continuous variability analysis is supported by this observational trial, and merits the attention of future investigations.

Author Contributions

Conceived and designed the experiments: LH AJES. Performed the experiments: SA SF SM. Analyzed the data: SA TR LH SF IB LM SRS DEM FMS PH DF AJES. Contributed reagents/materials/analysis tools: LH SM. Wrote the paper: SA TR AJES. Contributed to the formulation of the research protocol, interpretation of the results and preparation of the manuscript: SA TR LH SF SM IB LM SRS DEM FMS PH DF AT. Principally involved in the analysis and drafting of the manuscript: SA AJES.

Citation: Ahmad S, Ramsay T, Huebsch L, Flanagan S, McDiarmid S, Batkin I, et al. (2009) Continuous Multi-Parameter Heart Rate Variability Analysis Heralds Onset of Sepsis in Adults. PLoS ONE4(8): e6642. https://doi.org/10.1371/journal.pone.0006642

References

1. Marik PE, Varon J (2001) Sepsis - State of the Art. Dis Mon 47: 462–532.PE MarikJ. Varon2001Sepsis - State of the Art.Dis Mon47462532

2. Angus DC, Linde-Zwirble WT, Lidicker J, Clermont G, Carcillo J, et al. (2001) Epidemiology of severe sepsis in the United States: analysis of incidence, outcome, and associated costs of care. Crit Care Med 29: 1303–1310.DC AngusWT Linde-ZwirbleJ. LidickerG. ClermontJ. Carcillo2001Epidemiology of severe sepsis in the United States: analysis of incidence, outcome, and associated costs of care.Crit Care Med2913031310

3. Neth OW, Bajaj-Elliott M, Turner MW, Klein NJ (2005) Susceptibility to infection in patients with neutropenia: the role of the innate immune system. Br J Haematol 129: 713–722.OW NethM. Bajaj-ElliottMW TurnerNJ Klein2005Susceptibility to infection in patients with neutropenia: the role of the innate immune system.Br J Haematol129713722

4. Lalami Y, Paesmans M, Muanza F, Barette M, Plehiers B, et al. (2005) Can we predict the duration of chemotherapy-induced neutropenia in febrile neutropenic patients, focusing on regimen-specific risk factors? A retrospective analysis. Ann Oncol. Y. LalamiM. PaesmansF. MuanzaM. BaretteB. Plehiers2005Can we predict the duration of chemotherapy-induced neutropenia in febrile neutropenic patients, focusing on regimen-specific risk factors? A retrospective analysis.Ann Oncol

5. Celebi H, Akan H, Akcaglayan E, Ustun C, Arat M (2000) Febrile neutropenia in allogeneic and autologous peripheral blood stem cell transplantation and conventional chemotherapy for malignancies. Bone Marrow Transplant 26: 211–214.H. CelebiH. AkanE. AkcaglayanC. UstunM. Arat2000Febrile neutropenia in allogeneic and autologous peripheral blood stem cell transplantation and conventional chemotherapy for malignancies.Bone Marrow Transplant26211214

6. Regazzoni CJ, Khoury M, Irrazabal C, Myburg C, Galvalisi NR, et al. (2003) Neutropenia and the development of the systemic inflammatory response syndrome. Intensive Care Med 29: 135–138.CJ RegazzoniM. KhouryC. IrrazabalC. MyburgNR Galvalisi2003Neutropenia and the development of the systemic inflammatory response syndrome.Intensive Care Med29135138

7. Rivers EP, Nguyen HB, Havstad S, Ressler J, Muzzin A, et al. (2001) Early Goal-Directed Therapy in the Treatment of Severe Sepsis and Septic Shock. N Engl J Med 345: 1368–1377.EP RiversHB NguyenS. HavstadJ. ResslerA. Muzzin2001Early Goal-Directed Therapy in the Treatment of Severe Sepsis and Septic Shock.N Engl J Med34513681377

8. Otero RM, Nguyen HB, Huang DT, Gaieski DF, Goyal M, et al. (2006) Early Goal-Directed Therapy in Severe Sepsis and Septic Shock Revisited: Concepts, Controversies, and Contemporary Findings. Chest 130: 1579–1595.RM OteroHB NguyenDT HuangDF GaieskiM. Goyal2006Early Goal-Directed Therapy in Severe Sepsis and Septic Shock Revisited: Concepts, Controversies, and Contemporary Findings.Chest13015791595

9. Kumar A, Roberts D, Wood KE, Light B, Parrillo JE, et al. (2006) Duration of hypotension before initiation of effective antimicrobial therapy is the critical determinant of survival in human septic shock. Crit Care Med 34: 1589–1596.A. KumarD. RobertsKE WoodB. LightJE Parrillo2006Duration of hypotension before initiation of effective antimicrobial therapy is the critical determinant of survival in human septic shock.Crit Care Med3415891596

10. Seely AJE, Christou NV (2000) Multiple organ dysfunction syndrome: exploring the paradigm of complex nonlinear systems. Crit Care Med 28: 2193–2200.AJE SeelyNV Christou2000Multiple organ dysfunction syndrome: exploring the paradigm of complex nonlinear systems.Crit Care Med2821932200

11. Korach M, Sharshar T, Jarrin I, Fouillot JP, Raphaël JC, et al. (2001) Cardiac variability in critically ill adults: influence of sepsis. Crit Care Med 29: 1380–1385.M. KorachT. SharsharI. JarrinJP FouillotJC Raphaël2001Cardiac variability in critically ill adults: influence of sepsis.Crit Care Med2913801385

12. Piepoli M, Garrard CS, Kontoyannis DA, Bernardi L (1995) Autonomic control of the heart and peripheral vessels in human septic shock. Intensive Care Med 21: 112–119.M. PiepoliCS GarrardDA KontoyannisL. Bernardi1995Autonomic control of the heart and peripheral vessels in human septic shock.Intensive Care Med21112119

13. Annane D, Trabold F, Sharshar T, Jarrin I, Blanc AS, et al. (1999) Inappropriate sympathetic activation at onset of septic shock: a spectral analysis approach. Am J Respir Crit Care Med 160: 458–465.D. AnnaneF. TraboldT. SharsharI. JarrinAS Blanc1999Inappropriate sympathetic activation at onset of septic shock: a spectral analysis approach.Am J Respir Crit Care Med160458465

14. Barnaby D, Ferrick K, Kaplan DT, Shah S, Bijur P, et al. (2002) Heart Rate Variability in Emergency Department Patients with Sepsis. Acad Emerg Med 9: 661–670.D. BarnabyK. FerrickDT KaplanS. ShahP. Bijur2002Heart Rate Variability in Emergency Department Patients with Sepsis.Acad Emerg Med9661670

15. Garrard CS, Kontoyannis DA, Piepoli M (1993) Spectral analysis of heart rate variability in the sepsis syndrome. Clin Auton Res 3: 5–13.CS GarrardDA KontoyannisM. Piepoli1993Spectral analysis of heart rate variability in the sepsis syndrome.Clin Auton Res3513

16. Chen WL, Chen JH, Huang CC, Kuo CD, Huang CI, et al. (2008) Heart rate variability measures as predictors of in-hospital mortality in ED patients with sepsis. Am J Emerg Med 26: 395–401.WL ChenJH ChenCC HuangCD KuoCI Huang2008Heart rate variability measures as predictors of in-hospital mortality in ED patients with sepsis.Am J Emerg Med26395401

17. Pontet J, Contreras P, Curbelo A, Medina J, Noveri S, et al. (2003) Heart rate variability as early marker of multiple organ dysfunction syndrome in septic patients. J Crit Care 18: 156–163.J. PontetP. ContrerasA. CurbeloJ. MedinaS. Noveri2003Heart rate variability as early marker of multiple organ dysfunction syndrome in septic patients.J Crit Care18156163

18. Griffin MP, Lake DE, Bissonette EA, Harrell FE, O'Shea TM, et al. (2005) Heart rate characteristics: novel physiomarkers to predict neonatal infection and death. Pediatrics 116: 1070–1074.MP GriffinDE LakeEA BissonetteFE HarrellTM O'Shea2005Heart rate characteristics: novel physiomarkers to predict neonatal infection and death.Pediatrics11610701074

19. Cao H, Lake DE, Griffin MP, Moorman JR (2004) Increased nonstationarity of neonatal heart rate before the clinical diagnosis of sepsis. Ann Biomed Eng 32: 233–244.H. CaoDE LakeMP GriffinJR Moorman2004Increased nonstationarity of neonatal heart rate before the clinical diagnosis of sepsis.Ann Biomed Eng32233244

20. Ellenby MS, McNames J, Lai S, McDonald BA, Krieger D, et al. (2001) Uncoupling and recoupling of autonomic regulation of the heart beat in pediatric septic shock. Shock 16: 274–277.MS EllenbyJ. McNamesS. LaiBA McDonaldD. Krieger2001Uncoupling and recoupling of autonomic regulation of the heart beat in pediatric septic shock.Shock16274277

21. Seely AJE, Macklem PT (2004) Complex systems and the technology of variability analysis. Crit Care 8: 367–384.AJE SeelyPT Macklem2004Complex systems and the technology of variability analysis.Crit Care8367384

22. Clifford GD, Francisco A, McSharry P (2006) Advanced Methods And Tools for ECG Data Analysis. London: Artech House Publishers. GD CliffordA. FranciscoP. McSharry2006Advanced Methods And Tools for ECG Data Analysis.LondonArtech House Publishers384

23. Töyry JP, Niskanen LK, Mäntysaari MJ, Länsimies EA, Haffner SM, et al. (1997) Do High Proinsulin and C-Peptide Levels Play a Role in Autonomic Nervous Dysfunction? Power Spectral Analysis in Patients With Non–Insulin-Dependent Diabetes and Nondiabetic Subjects. Circulation 96: 1185–1191.JP TöyryLK NiskanenMJ MäntysaariEA LänsimiesSM Haffner1997Do High Proinsulin and C-Peptide Levels Play a Role in Autonomic Nervous Dysfunction? Power Spectral Analysis in Patients With Non–Insulin-Dependent Diabetes and Nondiabetic Subjects.Circulation9611851191

24. Low PA, Benrud-Larson LM, Sletten DM, Opfer-Gehrking TL, Weigand SD, et al. (2004) Autonomic Symptoms and Diabetic Neuropathy: A population-based study. Diabetes Care 27: 2942–2947.PA LowLM Benrud-LarsonDM SlettenTL Opfer-GehrkingSD Weigand2004Autonomic Symptoms and Diabetic Neuropathy: A population-based study.Diabetes Care2729422947

25. Goldberger AL, Peng CK, Lipsitz LA (2002) What is physiologic complexity and how does it change with aging and disease? Neurobiol Aging 23: 23–26.AL GoldbergerCK PengLA Lipsitz2002What is physiologic complexity and how does it change with aging and disease?Neurobiol Aging232326

26. Godin PJ, Buchman TG (1996) Uncoupling of biological oscillators: a complementary hypothesis concerning the pathogenesis of multiple organ dysfunction syndrome. Crit Care Med 24: 1107–1116.PJ GodinTG Buchman1996Uncoupling of biological oscillators: a complementary hypothesis concerning the pathogenesis of multiple organ dysfunction syndrome.Crit Care Med2411071116

27. Lin MT, Wang JK, Lu FL, Wu ET, Yeh SJ, et al. (2006) Heart rate variability monitoring in the detection of central nervous system complications in children with enterovirus infection. J Crit Care 21: 280–286.MT LinJK WangFL LuET WuSJ Yeh2006Heart rate variability monitoring in the detection of central nervous system complications in children with enterovirus infection.J Crit Care21280286