Person:

Goldberger, Ary

Major depression affects multiple physiologic systems. Therefore, analysis of signals that reflect integrated function may be useful in probing dynamical changes in this syndrome. Increasing evidence supports the conceptual framework that complex variability is a marker of healthy, adaptive control mechanisms and that dynamical complexity decreases with aging and disease. We tested the hypothesis that heart rate (HR) dynamics in non-medicated, young to middle-aged males during an acute major depressive episode would exhibit lower complexity compared with healthy counterparts. We analyzed HR time series, a neuroautonomically regulated signal, during sleep, using the multiscale entropy method. Our results show that the complexity of the HR dynamics is significantly lower for depressed than for non-depressed subjects for the entire night (P<0.02) and combined sleep stages 1 and 2 (P<0.02). These findings raise the possibility of using the complexity of physiologic signals as the basis of novel dynamical biomarkers of depression.

Background: The emergency department diagnosis of sinus versus nonsinus tachycardia is an important clinical challenge. The objective of this study was to identify subjects with a high prevalence of nonsinus tachycardia. Methods: Heart rate and cardiac rhythm were prospective reviewed in 500 consecutive patients with heart rate ≥ 100 beats/min in a busy emergency department. A predictive model based on age and heart rate was then developed to identify the probability of nonsinus tachycardia. Results: As age and heart rate increased, nonsinus tachycardias became more frequent. The probability of nonsinus tachycardia in a subject ≥ 71 years with heart rate ≥ 141 beats/minute was 93%, compared to only three percent in a subject ≤ 50 years with heart rate 100–120 beats/minute. A simple point score system based on age and heart rate helps predict the probability of sinus tachycardia versus nonsinus tachycardia. Conclusion: Nonsinus tachycardia is significantly more common than sinus tachycardia in elderly patients in the emergency department. The diagnosis of sinus tachycardia becomes much less likely as age and heart rate increase.

The loss of complexity in physiological systems may be a dynamical biomarker of aging and disease. In this study the effects of combined strength and endurance training compared with those of endurance training or strength training alone on heart rate (HR) complexity and traditional HR variability indices were examined in middle-aged women. 90 previously untrained female volunteers between the age of 40 and 65 years completed a 21 week progressive training period of either strength training, endurance training or their combination, or served as controls. Continuous HR time series were obtained during supine rest and submaximal steady state exercise. The complexity of HR dynamics was assessed using multiscale entropy analysis. In addition, standard time and frequency domain measures were also computed. Endurance training led to increases in HR complexity and selected time and frequency domain measures of HR variability (P<0.01) when measured during exercise. Combined strength and endurance training or strength training alone did not produce significant changes in HR dynamics. Inter-subject heterogeneity of responses was particularly noticeable in the combined training group. At supine rest, no training-induced changes in HR parameters were observed in any of the groups. The present findings emphasize the potential utility of endurance training in increasing the complex variability of HR in middle-aged women. Further studies are needed to explore the combined endurance and strength training adaptations and possible gender and age related factors, as well as other mechanisms, that may mediate the effects of different training regimens on HR dynamics.

Background: Physiologic signals, such as cardiac interbeat intervals, exhibit complex fluctuations. However, capturing important dynamical properties, including nonstationarities may not be feasible from conventional time series graphical representations. Methods: We introduce a simple-to-implement visualisation method, termed dynamical density delay mapping (“D3-Map” technique) that provides an animated representation of a system’s dynamics. The method is based on a generalization of conventional two-dimensional (2D) Poincaré plots, which are scatter plots where each data point, x(n), in a time series is plotted against the adjacent one, x(n + 1). First, we divide the original time series, x(n) (n = 1,…, N), into a sequence of segments (windows). Next, for each segment, a three-dimensional (3D) Poincaré surface plot of x(n), x(n + 1), h[x(n),x(n + 1)] is generated, in which the third dimension, h, represents the relative frequency of occurrence of each (x(n),x(n + 1)) point. This 3D Poincaré surface is then chromatised by mapping the relative frequency h values onto a colour scheme. We also generate a colourised 2D contour plot from each time series segment using the same colourmap scheme as for the 3D Poincaré surface. Finally, the original time series graph, the colourised 3D Poincaré surface plot, and its projection as a colourised 2D contour map for each segment, are animated to create the full “D3-Map.” Results: We first exemplify the D3-Map method using the cardiac interbeat interval time series from a healthy subject during sleeping hours. The animations uncover complex dynamical changes, such as transitions between states, and the relative amount of time the system spends in each state. We also illustrate the utility of the method in detecting hidden temporal patterns in the heart rate dynamics of a patient with atrial fibrillation. The videos, as well as the source code, are made publicly available. Conclusions: Animations based on density delay maps provide a new way of visualising dynamical properties of complex systems not apparent in time series graphs or standard Poincaré plot representations. Trainees in a variety of fields may find the animations useful as illustrations of fundamental but challenging concepts, such as nonstationarity and multistability. For investigators, the method may facilitate data exploration.

Background: This study evaluated the effects of stroke on regulation of cerebral blood flow in response to fluctuations in systemic blood pressure (BP). The autoregulatory dynamics are difficult to assess because of the nonstationarity and nonlinearity of the component signals. Methods: We studied 15 normotensive, 20 hypertensive and 15 minor stroke subjects (48.0 ± 1.3 years). BP and blood flow velocities (BFV) from middle cerebral arteries (MCA) were measured during the Valsalva maneuver (VM) using transcranial Doppler ultrasound. Results: A new technique, multimodal pressure-flow analysis (MMPF), was implemented to analyze these short, nonstationary signals. MMPF analysis decomposes complex BP and BFV signals into multiple empirical modes, representing their instantaneous frequency-amplitude modulation. The empirical mode corresponding to the VM BP profile was used to construct the continuous phase diagram and to identify the minimum and maximum values from the residual BP (BPR) and BFV (BFVR) signals. The BP-BFV phase shift was calculated as the difference between the phase corresponding to the BPR and BFVR minimum (maximum) values. BP-BFV phase shifts were significantly different between groups. In the normotensive group, the BFVR minimum and maximum preceded the BPR minimum and maximum, respectively, leading to large positive values of BP-BFV shifts. Conclusion: In the stroke and hypertensive groups, the resulting BP-BFV phase shift was significantly smaller compared to the normotensive group. A standard autoregulation index did not differentiate the groups. The MMPF method enables evaluation of autoregulatory dynamics based on instantaneous BP-BFV phase analysis. Regulation of BP-BFV dynamics is altered with hypertension and after stroke, rendering blood flow dependent on blood pressure.

Background: In older adults, depression has been associated with increased fall risk, but the reasons for this link are not fully clear. Given parallels between major depression and Parkinson's disease, we hypothesized that major depression and related affective disorders would be associated with impairment in the ability to regulate the stride-to-stride fluctuations in gait cycle timing. Methods: We measured stride-to-stride fluctuations of patients with two forms of mood disorders, unipolar major depressive disorder (MDD) and bipolar disorder, and compared their gait to that of a healthy control group. The primary outcomes were two measures of gait unsteadiness that have been associated with fall risk: stride time variability and swing time variability. Results: Compared to the control group, the two patient groups tended to walk more slowly and with decreased swing time and increased stride time. However, none of these differences was statistically significant. Compared to the control group, swing time variability was significantly larger in the subjects with bipolar disorder (p < 0.0001) and in the subjects with MDD (p < 0.0004). Conclusions: Patients with MDD and patients with bipolar disorder display gait unsteadiness. This perturbation in gait may provide a mechanistic link connecting depression and falls. The present findings also suggest the possibility that measurement of variability of gait may provide a readily quantifiable objective approach to monitoring depression and related affective disorders.

Background: Short-term heart rate variability (HRV) is most commonly attributed to physiologic vagal tone modulation. However, with aging and cardiovascular disease, the emergence of high short-term HRV, consistent with the breakdown of the neuroautonomic-electrophysiologic control system, may confound traditional HRV analysis. An apparent dynamical signature of such anomalous short-term HRV is frequent changes in heart rate acceleration sign, defined here as heart rate fragmentation. Objective: The aims were to: (1) introduce a set of metrics designed to probe the degree of sinus rhythm fragmentation; (2) test the hypothesis that the degree of fragmentation of heartbeat time series increases with the participants' age in a group of healthy subjects; (3) test the hypothesis that the heartbeat time series from patients with advanced coronary artery disease (CAD) are more fragmented than those from healthy subjects; and (4) compare the performance of the new fragmentation metrics with standard time and frequency domain measures of short-term HRV. Methods: We analyzed annotated, open-access Holter recordings (University of Rochester Holter Warehouse) from healthy subjects and patients with CAD using these newly introduced metrics of heart rate fragmentation, as well as standard time and frequency domain indices of short-term HRV, detrended fluctuation analysis and sample entropy. Results: The degree of fragmentation of cardiac interbeat interval time series increased significantly as a function of age in the healthy population as well as in patients with CAD. Fragmentation was higher for the patients with CAD than the healthy subjects. Heart rate fragmentation metrics outperformed traditional short-term HRV indices, as well as two widely used nonlinear measures, sample entropy and detrended fluctuation analysis short-term exponent, in distinguishing healthy subjects and patients with CAD. The same level of discrimination was obtained from the analysis of normal-to-normal sinus (NN) and cardiac interbeat interval (RR) time series. Conclusion: The fragmentation framework and accompanying metrics introduced here constitute a new way of assessing short-term HRV under free-running conditions, one which appears to overcome salient limitations of traditional HRV analysis. Fragmentation of sinus rhythm cadence may provide new dynamical biomarkers for probing the integrity of the neuroautonomic-electrophysiologic network controlling the heartbeat in health and disease.

We introduce a generalization of multiscale entropy (MSE) analysis. The method is termed MSEn, where the subscript denotes the moment used to coarse-grain a time series. MSEμ, described previously, uses the mean value (first moment). Here, we focus on MSEσ2, which uses the second moment, i.e., the variance. MSEσ2 quantifies the dynamics of the volatility (variance) of a signal over multiple time scales. We use the method to analyze the structure of heartbeat time series. We find that the dynamics of the volatility of heartbeat time series obtained from healthy young subjects is highly complex. Furthermore, we find that the multiscale complexity of the volatility, not only the multiscale complexity of the mean heart rate, degrades with aging and pathology. The “bursty” behavior of the dynamics may be related to intermittency in energy and information flows, as part of multiscale cycles of activation and recovery. Generalized MSE may also be useful in quantifying the dynamical properties of other physiologic and of non-physiologic time series.

Background: Poincaré delay maps are widely used in the analysis of cardiac interbeat interval (RR) dynamics. To facilitate visualization of the structure of these time series, we introduce multiscale Poincaré (MSP) plots. Methods: Starting with the original RR time series, the method employs a coarse-graining procedure to create a family of time series, each of which represents the system’s dynamics in a different time scale. Next, the Poincaré plots are constructed for the original and the coarse-grained time series. Finally, as an optional adjunct, color can be added to each point to represent its normalized frequency. Results: We illustrate the MSP method on simulated Gaussian white and 1/f noise time series. The MSP plots of 1/f noise time series reveal relative conservation of the phase space area over multiple time scales, while those of white noise show a marked reduction in area. We also show how MSP plots can be used to illustrate the loss of complexity when heartbeat time series from healthy subjects are compared with those from patients with chronic (congestive) heart failure syndrome or with atrial fibrillation. Conclusions: This generalized multiscale approach to Poincaré plots may be useful in visualizing other types of time series.

Background: We recently introduced the concept of heart rate fragmentation along with a set of metrics for its quantification. The term was coined to refer to an increase in the percentage of changes in heart rate acceleration sign, a dynamical marker of a type of anomalous variability. The effort was motivated by the observation that fragmentation, which is consistent with the breakdown of the neuroautonomic-electrophysiologic control system of the sino-atrial node, could confound traditional short-term analysis of heart rate variability. Objective: The objectives of this study were to: (1) introduce a symbolic dynamical approach to the problem of quantifying heart rate fragmentation; (2) evaluate how the distribution of the different dynamical patterns (“words”) varied with the participants' age in a group of healthy subjects and patients with coronary artery disease (CAD); and (3) quantify the differences in the fragmentation patterns between the two sample populations. Methods: The symbolic dynamical method employed here was based on a ternary map of the increment NN interval time series and on the analysis of the relative frequency of symbolic sequences (words) with a pre-defined set of features. We analyzed annotated, open-access Holter databases of healthy subjects and patients with CAD, provided by the University of Rochester Telemetric and Holter ECG Warehouse (THEW). Results: The degree of fragmentation was significantly higher in older individuals than in their younger counterparts. However, the fragmentation patterns were different in the two sample populations. In healthy subjects, older age was significantly associated with a higher percentage of transitions from acceleration/deceleration to zero acceleration and vice versa (termed “soft” inflection points). In patients with CAD, older age was also significantly associated with higher percentages of frank reversals in heart rate acceleration (transitions from acceleration to deceleration and vice versa, termed “hard” inflection points). Compared to healthy subjects, patients with CAD had significantly higher percentages of soft and hard inflection points, an increased percentage of words with a high degree of fragmentation and a decreased percentage of words with a lower degree of fragmentation. Conclusion: The symbolic dynamical method employed here was useful to probe the newly recognized property of heart rate fragmentation. The findings from these cross-sectional studies confirm that CAD and older age are associated with higher levels of heart rate fragmentation. Furthermore, fragmentation with healthy aging appears to be phenotypically different from fragmentation in the context of CAD.