Skip to main content

ORIGINAL RESEARCH article

Front. Physiol., 14 November 2017
Sec. Computational Physiology and Medicine

Heart Rate Fragmentation: A Symbolic Dynamical Approach

  • 1Department of Medicine, Beth Israel Deaconess Medical Center, Margret and H. A. Rey Institute for Nonlinear Dynamics in Medicine, Harvard Medical School, Boston, MA, United States
  • 2Division of General Medicine and Primary Care, Department of Medicine, Beth Israel Deaconess Medical Center, Harvard Medical School, Boston, MA, United States

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.

Words, words, words. W. Shakespeare, Hamlet, Act 2, Scene 2

1. Introduction

Analysis of fluctuations in cardiac interbeat intervals, under the rubric of heart rate variability (HRV), continues to generate much interest as a uniquely accessible window into the complex network of regulatory mechanisms controlling the sino-atrial (SA) node (HRV, 1996; Billman, 2013). Particular emphasis has been placed on the analysis of short-term fluctuations, i.e., oscillatory patterns with cycle lengths ranging from approximately four to eight consecutive beats. Such fluctuations, termed respiratory sinus arrhythmia, are primarily ascribable to the coupling between heart rate and breathing, mediated by the parasympathetic (vagal) nervous system (Angelone and Coulter, 1964; Hirsch and Bishop, 1981).

However, short-term fluctuations in heart rate are not always a marker of healthy cardiopulmonary interactions (Makikallio et al., 1999; Domitrovich and Stein, 2002; Stein, 2002; Wiklund et al., 2008; Costa et al., 2017) (Figure 1, first three rows). They may also be associated with abnormalities in the function of the neuroautonomic system, the SA node and other electrophysiologic components (Geiger and Goerner, 1945; Binkley et al., 1995; Jalife, 2013). Anomalous short-term variability is important for two major reasons: (1) it may confound the assessment of vagal tone modulation using conventional time and frequency domain HRV measures, leading to inflated estimates of “healthy” autonomic function in the elderly and especially in those with clinical or pre-clinical organic heart disease; and (2) its presence, itself, may be a novel dynamical biomarker of pathology and increased risk of adverse cardiovascular outcomes.

FIGURE 1
www.frontiersin.org

Figure 1. Examples of respiratory sinus arrhythmia and anomalous (fragmented) sinus rhythm. Electrocardiograms (Holter lead) from a healthy subject (first row) and a patient with coronary artery disease (CAD) (second row), both from the present study. Normal-to-normal (NN) sinus interval time series from the healthy subject (third row, left) and the patient with CAD (third row, right). The fluctuation patterns of the former time series are characteristic of phasic (respiratory) sinus arrhythmia, while that of the latter are indicative of an abnormal, non-phasic sinus arrhythmia (Costa et al., 2017). Positive and negative changes in the value of the NN intervals, corresponding to heart rate decelerations and accelerations were mapped to symbols “−1” and “1,” respectively. Symbol “0” is used to represent intervals in which heart rate did not change. To assist in visual comparisons, pale gray backgrounds are used for data from the healthy subject and light red for data from the patient with CAD, respectively. The symbolic mapping of the differences between consecutive NN intervals for the ECG of the healthy subject (first 16 intervals) along with the first four words that were derived from this sequence are shown on the bottom left. The first word “−1−111” contains one hard inflection point. It belongs to the group W1 and, more specifically, to the subgroup W1H. The following three words, “−1110,” “110−1,” and “10−1−1” contain two inflection points. Therefore, they belong to group W2. However, the first word (“−1110”) belongs to the subgroup W2M since it contains one hard and one soft inflection point; the second (“110−1”) and the third (“10−1−1”) words belong to the subgroup W2S since they present two soft inflection points. The panels on the bottom right show the percentage of words in each group for the healthy subject (left) and patient with CAD (right). Note a substantially higher percentage of fragmented words for the patient with CAD than for the healthy subject. The abbreviation “a.u.” stands for arbitrary units.

To help further address these issues, we (Costa et al., 2017) recently introduced the concept of heart rate fragmentation, along with a set of metrics to quantify this property. The underlying framework is based on the observation that sustained physiologic changes in heart rate cannot persist at frequencies higher than those at which the intact parasympathetic nervous system operates. Although the maximal physiologic response frequency is difficult to pinpoint, anticorrelated, beat-to-beat changes in heart rate, characterized by frequent changes from acceleration to deceleration and vice versa, are clearly atypical or frankly abnormal. The fragmentation indices that we introduced (Costa et al., 2017) quantify the density of this type of pattern. The assumption was that the systems manifesting the highest degree of fragmentation (loss of “fluency”) were the most pathologic ones.

We showed (Costa et al., 2017) that: (i) the degree of fragmentation of the NN and RR time series, derived from 24-h Holter monitoring, varied directly as a function of cross-sectional age in cohorts of healthy young to elderly male and female subjects in sinus rhythm and of those with coronary artery disease (CAD); (ii) the degree of fragmentation was significantly higher in patients with CAD than in healthy subjects, both in unadjusted models and in those adjusted for age and sex; and (iii) fragmentation indices outperformed standard time and frequency domain measures, as well as, widely used non-linear measures, in separating healthy subjects from patients with CAD.

To gain additional insight into the temporal structure of heart rate fragmentation, we now introduce a symbolic dynamical approach to the quantification of this property. In general, symbolic mapping deliberately reduces the overall information content of a signal. At the same time, it provides a useful way of highlighting certain features deemed of interest, while deemphasizing others. In heart rate fragmentation studies, examples of features of interest are the changes in heart rate acceleration sign, whereas “details” that one may choose to ignore are the magnitudes of those changes. In this study, the general hypotheses were that the degree of fragmentation, quantified by a set of variables derived from the symbolic dynamical analysis described below, would be: (1) higher in older subjects than in their younger counterparts, and (2) higher in patients with CAD than in healthy subjects. In addition, we sought to explore whether different symbolic “phenotypes” could help in distinguishing physiologic aging from aging in the context of overt organic heart disease.

2. Methods

2.1. Databases

We employed the same two long-term (~24-h) ECG ambulatory databases from the Intercity Digital Electrocardiogram Alliance (IDEAL) study previously analyzed (Costa et al., 2017). The deidentified recordings are made available via the University of Rochester Telemetric and Holter ECG Warehouse (THEW) archives (http://thew-project.org/databases.htm).

1. Healthy Subjects Database (THEW identification: E-HOL-03-0202-003) The database comprises 24-h Holter recordings from 202 ostensibly healthy subjects (102 males). The ECG signals were recorded at a sampling frequency of 200 Hz. Automated beat annotations were manually reviewed and adjudicated. We excluded 45 subjects with more than 2% non-sinus beats, 37 younger than 25 years old, 10 with BMI >30 Kg/m2 and one with <12 h of data. Overall, our analyses comprised 109 healthy adult subjects (60 male), age (median, 25th–75th) 40, 33–49 years.

2. Coronary Artery Disease Subjects Database (THEW identification E-HOL-03-0271-002).

This database comprises 24-h Holter recordings from 271 patients (223 males). Subjects had an abnormal coronary angiogram (at least one vessel with luminal narrowing >75%) and either exercise-induced ischemia or a documented previous myocardial infarction. For our analysis, we also excluded 11 subjects whose Holter recordings contained ≥20% non-sinus beats and four with less than 12 h of data. Overall, we analyzed 256 subjects (208 male), age (median, 25th–75th): 60; 51–67 years; left ventricular ejection fraction 56.5%, 50–66.

As previously described Costa et al. (2017), presumed waking and sleeping periods were estimated as the six consecutive hours of highest and lowest heart rates, respectively. These periods were calculated from the NN interval time series using a 6-h moving average window shifted 15 min at a time. From the continuous ECG of each subject, the time series of the RR and NN intervals were derived. The former is the sequence of intervals between consecutive QRS complexes. The latter, is the subset of intervals between consecutive normal sinus to normal sinus QRS complexes.

2.2. Symbolic Mapping and Dynamical Analysis

The original interbeat interval time series, {si}, 1 ≤ iN (N, time series length) was mapped to a ternary symbolic sequence as follows: “−1” if ΔNNi < 0, “0” if ΔNNi = 0, and “1” if ΔNNi > 0. Of note, since the ECG signals were sampled at 200 Hz, the resolution (τ) of both the NN interval and the increment time series was 5 ms (1/200 s). Taking the sampling frequency into consideration, the symbolic mapping rules can be alternatively written as: “1” if ΔNNi ≤ −5 ms, “0” if −5 < ΔNNi < 5 ms, and “−1” if ΔNNi ≥ 5 ms. Next, the percentages of different segments of l consecutive symbols, wi = (si, si+1, …, si+l−1), 1 ≤ iNl+1, termed “words,” were calculated. (With an alphabet of n symbols, the number of different words of length l is nl.) Words derived from the NN interval time series were termed NN words. Words derived from the RR interval time series were termed RR words.

Since we were interested in the analysis of short-term dynamical patterns occurring at the respiratory frequency, we chose a word length of four, which corresponds to time scales of approximately 3–5 s, depending on the heart rate. Subsequently, the words were grouped according to the number and type of transition between consecutive symbols (Figure 1). Reversals in heart rate acceleration (ΔNNi × ΔNNi+1 < 0), i.e., transitions from symbol “1” to “−1” or vice versa, were termed hard (H) inflection points. Transitions to or from zero acceleration (ΔNNi × ΔNNi+1 = 0, ΔNNi ≠ ΔNNi+1), i.e., transitions from symbols “1” or “−1” to “0”, or vice versa, were termed soft (S) inflection points. The higher the number of inflection points in a word the more fragmented it was. Words of length four can contain no more than three inflection points. Word groups with only hard, only soft and a combination of hard and soft inflection points were, respectively, labeled WjH, WjS, and WjM (where “M” stands for “mixed” and j indicates the number of inflection points). The word groups with more than one inflection point, Wj (2 ≤ j ≤ 3), for which the type of inflection point was not specified, comprised the words from subgroups WjH, WjS, and WjM. The word group W1 comprised the words from subgroups W1H and W1S. Figure 2 shows a schematic representation of all the different words (n = 81).

FIGURE 2
www.frontiersin.org

Figure 2. Schematic representation of all (#81) different words of length 4 with an alphabet of 3 symbols. The symbols “/”, “\”, and “−” represent heart rate acceleration, deceleration and no change, respectively. Words were grouped by the number and type of inflection points. The labels, 0–80, shown in parentheses, are the decimal value of the ternary representation of each pattern using the symbols “2” if ΔNNi < 0, “1” if ΔNNi > 0 and “0” if ΔNNi = 0. For example, the label for the word comprising 4 consecutive accelerations, i.e., the word 2222, is 80 (= 2 × 33+2 × 32+2 × 31+2 × 30). Abbreviations: W, word subgroup. The subscript and superscript of W indicate, respectively, the number and the type of inflection points, hard (H), soft (S) or a combination of hard and soft (M, mixed) that the words in that subgroup contain.

Of note, to calculate the percentage of each NN word, two different denominators can be used: the total number of NN words and the total number of RR words. The former is not affected by the presence of ectopic beats, while the latter takes them into consideration. Here, the percentages of W0, Wj, WjH, WjS, and WjM, 1 ≤ j < 3 were computed using the total number of NN words. In addition, we calculated the percentages of hard (soft) NN words that contained one, two and three hard (soft) inflection points. In these cases, the denominators were the total number of hard (soft) NN words with at least one inflection point. These word subgroups were labeled WjH* (WjS*). Thus, while WjH (WjS) represents the overall percentage of words with j hard (soft) inflection points, WjH* (WjS*) represents the percentage of hard (soft) words with j inflection points. They were calculated as follows: WjH* = WjH/i=13WiH (WjS* = WjS/i=13WiS).

The percentages of hard (PIPH) and soft (PIPS) inflection points were also computed. PIPH and PIPS are subcategories of the fragmentation index, PIP, previously introduced (Costa et al., 2017).

We analyzed how PIPH, PIPS and the different group of words changed with the participants' age and with disease in unadjusted and adjusted [for age and sex, and age, sex and average NN interval (AVNN)] logistic models. Taking into consideration that heart rate fragmentation has been shown (Costa et al., 2017) to increase with cross-sectional age and with CAD in these databases, we specifically hypothesized that the percentages of words in groups W0 and W0 (least fragmented), would decrease with the participants' age and with disease, while the percentages of words in groups W0 and W0 (most fragmented), would increase, regardless of the type of inflection points.

2.3. Statistical Analysis

Outcome variables were summarized by their median, 25th and 75th percentile values.

Linear regression models were used to quantify the dependence of each of the outcome variables (y: Wj, WjH, WjH*, WjS, WjS*, WjM, PIPH, and PIPS) with the participants' age. These models included the interaction term between age and sample population to assess whether the regression slopes were the same in the two populations. In addition, these models included AVNN to control for the effects of this variable on each of the outcome variables (y = c + β1 × age + β2 × population + β3 × age × population + β4 × AVNN, where c is a constant). (The values of β1 and of β1 + β3 along with their confidence intervals are provided in Table 1 for the groups of healthy subjects and those with CAD, respectively.) Statistical significance was set at a p < 0.05.

TABLE 1
www.frontiersin.org

Table 1. Slope and [95%CI] of the association between each outcome measure and the participants' age for the group of healthy subjects and those with CAD, for the 24-h and putative awake and sleep periods.

Logistic regression analysis was used to assess the relationships between presence of CAD and each of the outcome measures in unadjusted models, and models adjusted for age and sex, and age, sex, and AVNN. To facilitate comparisons among various outcome variables, we report normalized odds ratios (i.e., the odds ratio for a one standard deviation change in a given measure).

The area under the receiver operating characteristic (AUC) curve was used to assess the discrimination of each model. The likelihood ratio test was used to compare the fit of two nested models. All analyses were performed using raw measures.

3. Results

3.1. Relationship between Heart Rate Fragmentation Indices and Participants' Age

As previously reported in Costa et al. (2017), the overall percentage of inflection points (soft and hard combined) significantly increased with the participants' age in both healthy subjects and those with CAD. However, the percentages of only soft and only hard inflection points changed with the participants' age in different ways for each of the groups (Table 1). In healthy subjects, the percentage of soft inflection points significantly increased with the participants' age (slope of the regression line, [95% CI]: 0.26[0.18, 0.35] %/yr), while the percentage of hard inflection points did not (0.04[−0.06, 0.13] %/yr). In patients with CAD, both, the percentages of soft and hard inflection points increased with the participants' age at identical rates, 0.16[0.10, 0.22] %/yr and 0.16[0.10, 0.23] %/yr, respectively. These findings were consistent across all time periods with only one exception. During the putative sleep period, the increase in the percentage of hard inflection points did not reach statistical significance for those with CAD.

Overall, the percentage of words W0 and W0, i.e., the least fragmented (most “fluent”), decreased with the participants' age both in the groups of healthy subjects and those with CAD. All relationships were significant with the exception of the one with W0 during the putative sleep period for the healthy subjects. Complementarily, the percentage of words W0, the most fragmented (least “fluent”) significantly increased with the participants' age in both sample populations for all time periods (Figure 3 and Table 1). The percentage of words W0 that capture patterns of transitions between fluent and fragmented dynamics also tended to increase with the participants' age. However, only some of these relationships reached significance.

FIGURE 3
www.frontiersin.org

Figure 3. Relationship between the percentage of words with no inflection points (W0), one (W0), two (W0) and three (W0) inflection points and the participants' age for the healthy subjects (blue) and those with coronary artery disease (CAD, red) during the 24-h (A), putative awake (B) and putative sleep (C) periods. Symbols and lines represent, respectively, word percentages for each subject and the regression lines derived from linear regression analyses controlled for the average NN interval. In each plot, the rates of change of the outcome variables per year of age for the healthy subjects and the patients with CAD are indicated in blue and red, respectively.

In both healthy subjects and patients with CAD, the word groups with the strongest association with the participants' age were W1, among the most fluent, and W0, among the most fragmented. In both cases, the magnitude of the rate of change was above 0.2%/yr (negative rate for W1 and positive rate for W3).

In general, the number of inflection points in a given word subgroup, not the type of inflection points (H, S or M) determined the directionality of the changes in its density with the participants' age. Among a total of 84 relationships (14 word subgroups, three time periods and two sample populations) only 7 did not change with the participants' age in the expected direction (Table 1); all except one of these relationships (W2S for the group of patients with CAD) were for the putative sleep period.

The most notable difference between healthy subjects and patients with CAD in the symbolic analysis related to the words with three hard inflection points (W3H). While the percentage of these words did not significantly change with cross-sectional age in the group of healthy subjects for any time periods, it significantly increased in the group of patients with CAD, for all time periods, at rates varying between 0.13 and 0.20%/year.

3.2. Changes in Heart Rate Fragmentation Indices with Coronary Artery Disease

A 1-year increase in age was associated with an increase of 14% in the odds of having CAD (odds ratio [95%CI]: 1.14[1.11, 1.17] <0.0001). The AUC for the model with age as the only covariate was 0.853. Male sex carried a 3.54-fold increase in the odds of CAD (3.54[2.17, 5.78], p < 0.0001). The AUC for the null model with age and sex as the sole independent variables was 0.882. The AUC for the null model with age, sex and AVNN was 0.910.

3.2.1. Unadjusted Analyses

Summary statistics of the calculated indices for the groups of healthy subjects and patients with CAD are presented in Table 2. In unadjusted analyses (Table 3), the percentage of hard inflection points, PIPH, was significantly higher in those with CAD than in healthy subjects for the three time periods. The percentage of soft inflection points, PIPS, showed a similar behavior. However, statistical significance was only reached for the putative sleep time.

TABLE 2
www.frontiersin.org

Table 2. Measures of heart rate fragmentation/fluency in healthy subjects and those with coronary artery disease.

TABLE 3
www.frontiersin.org

Table 3. Logistic regression analysis and area under the ROC curve for unadjusted models of CAD.

In these models, lower percentages of words without (W0) and with only one (W0) inflection point, and higher percentages of words with two (W0) and three (W0) inflection points were significantly associated with the presence of CAD, for all time periods. Similar results were obtained for the subgroups of hard (WjH and WjH*) and mixed (WjM) words with any number of inflection points, for all time periods.

The percentages of words with one and two soft inflection points were significantly higher in healthy subjects than in patients with CAD for the 24-h and putative awake periods. For the sleep period, the differences were not significant. The percentage of soft words with three inflection points, the most fragmented of this class, tended to be higher in patients with CAD than in healthy subjects. However, significance was only reached for the sleep period.

3.2.2. Analyses Adjusted for Age and Sex

The percentage of hard inflection points, PIPH, remained positively associated with CAD in models adjusted for age and sex (Table 4) for the three time periods. The odds of CAD more than tripled, for each one-standard deviation increase in PIPH during the 24-h and putative awake periods. For the putative sleep period, the odds doubled. The percentage of soft inflection points, PIPS tended to be higher in healthy subjects than in those with CAD but the difference did not reach statistical significance for any of the time periods. Adding PIPH to a model with only age and sex significantly improved its performance, whereas PIPS did not.

TABLE 4
www.frontiersin.org

Table 4. Logistic regression analysis and AUC for models of CAD adjusted for age and sex.

After adjusting for age and sex, the proportion of fragmented words, W2 and W3, remained positively associated with CAD (Table 4) for all time periods, while the proportion of fluent words, W0 and W1, remained negatively associated with CAD (Table 4) for all periods. Specifically, for the 24-h period, a one-standard deviation increase in W2 and W3 was associated with an increase of 180 and 75% in the odds of CAD, respectively. In addition, a one-standard deviation increase in W0 and W1 was associated with a drop of 63 and 61% in the odds of CAD, respectively. All word groups Wj, 0 ≤ j ≤ 3, significantly improved the performance of a model with age and sex alone.

Hard word subgroups, WjH, and WjH* changed with disease in the same way as the groups Wj, (1 ≤ j ≤ 3). Specifically, W1H and W1H* were lower in those with CAD than in healthy subjects, for all time periods. All comparisons were statistically significant except the one with the variable W1H for the putative awake period. In addition, W2H, W3H, W2H* and W3H* were significantly higher in those with CAD than in healthy subjects, for all time periods.

Soft word subgroups with one inflection point, W1S and W1S*, changed with disease in the same way as the fluent, hard word subgroups. Specifically, they were more frequent in healthy subjects than in patients with CAD. In contrast, the percentages of soft words with two and three inflection points were lower in patients with CAD than healthy subjects. The comparison with W2S were significant for all time periods. For W3S, only the comparison for the putative awake period reached significance. Mixed words with three inflection points were more discriminatory than those with two.

For the majority of cases (44 out of 54), adding a word group or subgroup to a model with age and sex significantly improved its performance. The exceptions were: W2M, W3S, and W3S* for the 24-h period, W1H, W2M, W3S* and W3M for the putative awake time, and W2S*, W2M, and W3S for the putative sleep time.

3.2.3. Analyses Adjusted for Age, Sex and AVNN

The major difference between the results of the analyses adjusted for age and sex and those adjusted for age, sex and AVNN, concerned the variable PIPS. When AVNN was added to the models, the differences in PIPS between patients with CAD and healthy subjects became strongly significant: a one-standard deviation increase in PIPS was associated with an increase in the odds of CAD ranging from 100 to 200%, for the different time periods.

In these models, PIPH remained positively associated with CAD (Table 5) for the 24-h and putative awake periods, despite a decrease in the values of the odds ratio. However, for the sleep period, statistical significance was lost.

TABLE 5
www.frontiersin.org

Table 5. Logistic regression analysis and AUC for models of CAD adjusted for age, sex and the average value of the NN intervals.

For all of the time periods, fully adjusted models with word groups without inflection points, with one and three inflection points were more discriminatory than those with two inflection points (Table 5). Specifically, a one-standard deviation increase in the percentage of fluent words W0 and W1 was associated with a decrease in the odds of CAD ranging from 38 to 64% and from 50 to 67%, respectively. A one-standard deviation increase in the percentage of fragmented words W3 was associated with an increase in the odds of CAD ranging from 132 to 278%.

Of note, in these fully adjusted models, the number of inflection points, not their type (H, S or M), i.e., the degree of fragmentation of the words, determined the directionality of the effects in the odds of CAD. The majority of the word subgroups significantly improved the performance of a null model with age, sex and AVNN. Within group 1, the most discriminatory word subgroups were W1H and W1H*. They appeared in significantly higher densities in healthy subjects than in patients with CAD, for all time periods. Within group 2, W2H* and W2M were the most discriminatory variables. For all time periods, significantly higher percentages of these words were observed in patients with CAD than in healthy subjects. Within group 3, all word subgroups were highly discriminatory of the two sample population. The only exception was W3H during the putative sleep period. As expected, CAD was associated with a significant increase in the density of these words.

4. Discussion

Recently, we described a property of short-term HRV termed fragmentation and introduced a set of metrics to quantify this feature (Costa et al., 2017). The key marker of heart rate fragmentation is an overall increase in the frequency of changes in heart rate acceleration sign.

The purpose here was to further explore the property of heart rate fragmentation using symbolic dynamical analysis with the same databases previously studied. The findings were consistent with those reported in Costa et al. (2017), indicating an increase in heart rate fragmentation with the participants' age and with the presence of CAD. In addition, the symbolic analyses suggested a potentially important dynamical difference between aging in the healthy population and in those with coronary disease. This difference was not anticipated prior to this analysis.

Briefly, the notable findings of the cross-sectional study of heart rate fragmentation with the participants' age were that: (i) in the cohort of healthy subjects, the percentage of soft but not of hard inflection points significantly increased as a function of age; (ii) the percentages of both soft and hard inflection points significantly increased with the participants' age in the cohorts of subjects with CAD; (iii) overall, the density of fluent words tended to decrease and the density of the fragmented ones tended to increase with the participants' age in both populations, for all time periods (Figure 3); (iv) the percentage of words with three hard inflection points, the most fragmented, changed with age differently in the cohorts of healthy subjects and patients with CAD. For the latter, these words markedly increased for all time periods. For the former, the increases did not reach significance. The results from the symbolic analysis adjusted for age, sex and AVNN were consistent with the finding that the overall percentage of transitions from acceleration/deceleration to deceleration/acceleration did not change with the participants' age in the healthy group.

The key findings from the symbolic dynamical analysis comparing healthy subjects and those with CAD were that: (i) the percentages of fluent words, W0, W1, W1H* were significantly higher in healthy subjects than in patients with CAD, for all time periods, in both unadjusted and adjusted models; (ii) the percentages of fragmented words, W2, W2H*, W3, and W3H*, were significantly lower in healthy subjects than in those with CAD, for all time periods and models; and (iii) although not all word subgroups were statistically different in the two sample populations for all time periods and models, importantly none of the fluent word subgroups, W1H, W1S, or W1S* was significantly higher in patients with CAD than in healthy subjects for any time period or for any model. Similarly, none of the fragmented word subgroups with three inflection points, W3H, W3H*, W3M, W3S, and W3S* was significantly higher in healthy subjects than in patients with CAD for any of the time periods in any model.

Overall, word subgroups with two inflection points were less discriminatory than the other ones, particularly in the fully adjusted models. This finding is not entirely surprising in light of the fact that words with two inflection points encode patterns that represent a transition between dynamical short-term fluency and fragmentation.

Qualitatively similar results to those presented here were obtained from the analysis of RR intervals time series (instead of NN) and words of length 5 (not presented here). Taken together these results robustly support the notion that heart rate fragmentation increases as a function of the participants' age and in the presence of overt CAD.

In this study, words without inflection points included segments of four consecutive accelerative, decelerative and zero acceleration intervals. Excluding the latter, i.e., the segments with no heart rate variability (neither fragmented nor fluent) from the word group W0, and quantifying their density separately, could potentially allow for a better characterization of a given study population, for example, one with chronic heart failure. In the present study, the results for the word group W0, including or excluding the word “0000,” were very similar. Therefore, we reported solely the results for which that word was included.

We wish to emphasize that the interpretation of the results for the word group W0 (with or without the inclusion of the word “0000”) can be dependent on the physiologic context. A deficit of these words is likely a consequence of a high degree of heart rate fragmentation. However, an excess is likely a consequence of long-term (above the normal respiratory frequency) trends in the data. These trends can be pathologic, as seen, for example, with sleep apnea syndromes (Guilleminault et al., 1984; Lipsitz et al., 1995; Guzik et al., 2013; Jiang et al., 2017), or physiologic, e.g., when associated with even mild bouts of exercise and recovery. The former conjecture is supported by the work of Guzik et al. (2013), who in a study of heart rate variability in subjects with various degrees of obstructive sleep apnea, found that an increased number of long (>5 intervals) deceleration and acceleration runs were most common in patients with severe sleep apnea.

Symbolic mapping of both the NN interval time series (Ravelo-Garcia et al., 2014; Cysarz et al., 2015) and of its increments have been used in many other studies, both in our laboratory (Yang et al., 2003; Costa et al., 2005, 2008) and others (Cysarz et al., 2000, 2015; Kantelhardt et al., 2002; Piskorski and Guzik, 2011; Guzik et al., 2012, 2013; Jiang et al., 2017). For example, Ashkenazy et al. (2001) used a binary map of the increment time series to analyze the correlation properties of the sign and magnitude heart rate time series of healthy subjects and patients with heart failure. For the shortest time scale explored, 6–16 NN intervals, they found that the dynamics of the sign time series of healthy subjects were closer to brown noise than those of patients with heart failure. This finding supports the hypothesis that long (>5 intervals) deceleration and acceleration runs are more common in healthy subjects than in patients with heart failure.

Guzik et al. (2012, 2013) and Piskorski and Guzik (2011) specifically analyzed the percentages of acceleration and deceleration runs of various lengths in a population of postinfarction patients. Overall, they found that decelerations runs of 2–10 intervals were significantly less frequent in non-survivors and used the runs of lengths 2, 4, and 8 to stratify all-cause mortality risk. The frequency of occurrence of runs of different lengths can be related to the concept of fragmentation. A higher percentage of short (<3) and a lower percentage of longer runs are expected in more fragmented than less fragmented time series. However, there is no direct correspondence between runs of a given length and a specific word. For example, runs of length 3 are necessarily part of the word group W1, but this word group also includes runs of lengths 1 and 2. Runs of length 2 are part of word groups W1 and W2; and runs of length 1 are part of all word groups but W0.

Cysarz et al. (2000, 2015) and Porta et al. (2007) used a binary map of the increment time series (“1” if ΔRRi+1 > RRi; “0” if ΔRRi+1 ≤ RRi) to analyze putative sympathetic/parasympathetic changes in neuroautonomic control under different conditions. However, for the types of fragmentation analyses proposed here, binary maps of heart rate increments are not recommended. In fact, if positive (or negative) and zero increments are mapped to the same symbol, soft inflection points are either “ignored” (when an accelerative interval is preceded or followed by an interval in which heart rate does not change), or “transformed” into hard inflection points (when a decelerative interval is preceded or followed by an interval in which heart rate does not change). Consequently, the word groups will contain words with different numbers of inflection points. For example, with the ternary map we used, the word group W0 contained only three words, specifically those labeled 0, 80, and 40 in Figure 2. However, with the binary map, this word group would also include the words 2, 8, 26, 54, 72, and 78 from W1S (Figure 2), the words 6, 18, 28, 56, 62, and 74 from subgroup W2S, and the words 20 and 60, from subgroup W3S, with one, two and three soft inflection points, respectively. The same would be true for other word groups. Thus, the binary mapping of the NN interval time series does not preserve all the information necessary for assessing heart rate fragmentation.

Our analyses were based on the definition of acceleration/deceleration as a decrease (increase) in consecutive NN (RR) intervals of ≤ (≥) 5 ms. We could have chosen any multiple of 5 ms, but not a lower value, since, as mentioned in the Methods section, the ECG signals were recorded at 200 Hz. ECG signals recorded with a higher sampling frequency (SF) would permit other choices, specifically, any multiple of 1/SF. However, higher resolution/lower thresholds may not necessarily translate into an enhanced ability to discriminate different populations. In fact, the lower the threshold the more likely the results are to be affected by both biological and instrumental noise. On the other hand, the larger the threshold the higher the number of significant changes in acceleration/deceleration that will not be detected. Future studies will help determine an “optimal” range of thresholds for fragmentation analysis.

As noted, increased fragmentation under free-running conditions is not directly attributable to variations in sympathetic or parasympathetic activity (Costa et al., 2017). These autonomic effectors do not operate on fast enough time scales to account for sustained beat-to-beat changes in heart rate acceleration sign. However, such rapid heart rate acceleration changes have been noted with a variety of pathophysiologic alterations, including subtle premature supraventricular extrasystoles coming from the SA node itself (or from nearby areas), SA exit block variants, modulated sinus parasystole, and possibly mechanical atrial stretch effects (Friedman, 1956; Nazir and Lab, 1996). These conditions are most likely to occur with the breakdown of the sinus regulatory control, possibly due to inflammation or fibrosis at various anatomic sites (Ghiassian et al., 2016). In fact, the most common clinical settings of atrial fibrillation, which may represent an end stage of supraventricular fragmentation, are seen with aging and chronic heart disease, conditions in which vagal tone is usually diminished, SA node size is reduced and intercellular coupling may be impaired (Moghtadaei et al., 2016).

A notable but unanticipated finding of this study was the difference in the pattern of fragmentation seen in the cross-sectional analysis of older healthy subjects vs. those with organic heart disease (advanced atherosclerosis). Fragmentation in older healthy subjects was mostly due to the increase in the percentage of transitions from acceleration/deceleration to zero acceleration or vice versa (soft inflection points). Fragmentation in those with CAD, fragmentation was also due to the increase in the percentage of transitions from acceleration to deceleration or vice versa (hard inflection points).

Speculatively, the increase in hard inflection points with disease, i.e., the emergence of beat-to-beat reversals in heart rate acceleration, might also relate to higher degrees of fibrosis and inflammation, substrates for the development of conduction and/or pacemaker abnormalities. The increase in soft inflection points likely relate, in part, to the well-documented decrease in the variance of the NN interval high-frequency fluctuations with aging (Pikkujamsa et al., 1999). In fact, if the structure of the variability is sufficiently preserved (a sign of health), a decrease in the amplitude of the time series, would translate into an increase in the likelihood of having consecutive NN intervals with the same value, that is, of zero accelerations and thus of soft inflection points (assuming that the temporal resolution of the time series does not change).

Although a benign increase in heart rate fragmentation should be rare, it might arise with vagally induced prominent sinus bradycardia with SA Wenckebach, a condition sometimes seen in very healthy (athletic) young subjects. Future studies in well-characterized, larger databases, with outcome data related to incident atrial fibrillation and advanced sinus node disease, should also help ascertaining the translational value of the symbolic analysis of heart rate fragmentation proposed here and the utility of heart rate fragmentation as a quantifiable descriptor of HRV.

Finally, it may be of interest to explore the utility of the concept of dynamical fragmentation and adapt the symbolic dynamic analysis introduced here to the study of changes in repolarization parameters, such as those described under the heading of T wave alternans.

5. Conclusion

A symbolic dynamical approach to the analysis of heart rate increment time series in a cross-sectional study of healthy subjects and those with CAD, provides evidence supporting the conjecture that fragmentation increases with age and disease. In addition, our results suggest that fragmentation in ostensibly healthy aging is different from fragmentation in the context of overt disease. Future studies analyzing larger databases with outcome measures are needed to confirm these findings and to assess their translational value.

Author Contributions

MC and AG developed the fragmentation concept and symbolic approach to its quantification. RD directed the statistical analysis. All three authors contributed to the interpretation of the findings and worked collaboratively on the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We gratefully acknowledge support from G. Harold and Leila Y. Mathers Charitable Foundation; the James S. McDonnell Foundation; and the National Institutes of Health grants R01GM104987 and R24HL114473.

References

Angelone, A., and Coulter, N. A. Jr. (1964). Respiratory sinus arrhythmia: a frequency dependent phenomenon. J. Appl. Physiol. 19, 479–482.

Google Scholar

Ashkenazy, Y., Ivanov, P. C., Havlin, S., Peng, C. K., Goldberger, A. L., Stanley, H. E., et al. (2001). Magnitude and sign correlations in heartbeat fluctuations. Phys. Rev. Lett. 86, 1900–1903. doi: 10.1103/PhysRevLett.86.1900

PubMed Abstract | CrossRef Full Text | Google Scholar

Billman, G. E. (2013). The effect of heart rate on the heart rate variability response to autonomic interventions. Front. Physiol. 4:222. doi: 10.3389/fphys.2013.00222

PubMed Abstract | CrossRef Full Text | Google Scholar

Binkley, P. F., Eaton, G. M., Nunziata, E., Khot, U., and Cody, R. J. (1995). Heart rate alternans. Ann. Intern. Med. 122, 115–117. doi: 10.7326/0003-4819-122-2-199501150-00007

PubMed Abstract | CrossRef Full Text | Google Scholar

Costa, M., Goldberger, A. L., and Peng, C. K. (2005). Broken asymmetry of the human heartbeat: loss of time irreversibility in aging and disease. Phys. Rev. Lett. 95:198102. doi: 10.1103/PhysRevLett.95.198102

PubMed Abstract | CrossRef Full Text | Google Scholar

Costa, M. D., Davis, R. B., and Goldberger, A. L. (2017). Heart rate fragmentation: a new approach to the analysis of cardiac interbeat interval dynamics. Front. Physiol. 8:255. doi: 10.3389/fphys.2017.00255

PubMed Abstract | CrossRef Full Text | Google Scholar

Costa, M. D., Peng, C. K., and Goldberger, A. L. (2008). Multiscale analysis of heart rate dynamics: entropy and time irreversibility measures. Cardiovasc. Eng. 8, 88–93. doi: 10.1007/s10558-007-9049-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Cysarz, D., Bettermann, H., and van Leeuwen, P. (2000). Entropies of short binary sequences in heart period dynamics. Am. J. Physiol. Heart Circ. Physiol. 278, H2163–H2172.

PubMed Abstract | Google Scholar

Cysarz, D., Van Leeuwen, P., Edelhauser, F., Montano, N., Somers, V. K., and Porta, A. (2015). Symbolic transformations of heart rate variability preserve information about cardiac autonomic control. Physiol. Meas. 36, 643–657. doi: 10.1088/0967-3334/36/4/643

PubMed Abstract | CrossRef Full Text | Google Scholar

Domitrovich, P. P., and Stein, P. K. (2002). A new method to detect erratic sinus rhythm in RR-interval files generated from Holter recordings. Comput. Cardiol. 26, 665–668. doi: 10.1109/CIC.2002.1166860

CrossRef Full Text | Google Scholar

Friedman, B. (1956). Alternation of cycle length in pulsus alternans. Am. Heart J. 51, 701–712. doi: 10.1016/S0002-8703(56)80006-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Geiger, A., and Goerner, J. (1945). Premature beats of sinus origin: electrocardiographic demonstration of a clinical case. Am. Heart J. 30, 284–291. doi: 10.1016/S0002-8703(45)90007-X

CrossRef Full Text | Google Scholar

Ghiassian, S. D., Menche, J., Chasman, D. I., Giulianini, F., Wang, R., Ricchiuto, P., et al. (2016). Endophenotype network models: common core of complex diseases. Sci. Rep. 6:27414. doi: 10.1038/srep27414

PubMed Abstract | CrossRef Full Text | Google Scholar

Guilleminault, C., Winkle, R., Connolly, S., Melvin, K., and Tilkian, A. (1984). Cyclical variation of the heart rate in sleep apnoea syndrome: mechanisms, and usefulness of 24 h electrocardiography as a screening technique. Lancet 323, 126–131. doi: 10.1016/S0140-6736(84)90062-X

CrossRef Full Text | Google Scholar

Guzik, P., Piskorski, J., Awan, K., Krauze, T., Fitzpatrick, M., and Baranchuk, A. (2013). Obstructive sleep apnea and heart rate asymmetry microstructure during sleep. Clin. Auton. Res. 23, 91–100. doi: 10.1007/s10286-013-0188-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Guzik, P., Piskorski, J., Barthel, P., Bauer, A., Muller, A., Junk, N., et al. (2012). Heart rate deceleration runs for postinfarction risk prediction. J. Electrocardiol. 45, 70–76. doi: 10.1016/j.jelectrocard.2011.08.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Hirsch, J. A., and Bishop, B. (1981). Respiratory sinus arrhythmia in humans: how breathing pattern modulates heart rate. Am. J. Physiol. 241, H620–H629.

PubMed Abstract | Google Scholar

HRV (1996). Heart rate variability: standards of measurement, physiological interpretation and clinical use. Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. Circulation 93, 1043–1065.

Jalife, J. (2013). Modulated parasystole: still relevant after all these years. Heart Rhythm 10, 1441–1443. doi: 10.1016/j.hrthm.2013.06.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, J., Chen, X., Zhang, C., Wang, G., Fang, J., Ma, J., et al. (2017). Heart rate acceleration runs and deceleration runs in patients with obstructive sleep apnea syndrome. Sleep Breath 21, 443–451. doi: 10.1007/s11325-016-1437-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Kantelhardt, J. W., Ashkenazy, Y., Ivanov, P. C., Bunde, A., Havlin, S., Penzel, T., et al. (2002). Characterization of sleep stages by correlations in the magnitude and sign of heartbeat increments. Phys. Rev. E 65:051908. doi: 10.1103/PhysRevE.65.051908

PubMed Abstract | CrossRef Full Text | Google Scholar

Lipsitz, L. A., Hashimoto, F., Lubowsky, L. P., Mietus, J., Moody, G. B., Appenzeller, O., et al. (1995). Heart rate and respiratory rhythm dynamics on ascent to high altitude. Heart 74, 390–396. doi: 10.1136/hrt.74.4.390

PubMed Abstract | CrossRef Full Text | Google Scholar

Makikallio, T. H., Koistinen, J., Jordaens, L., Tulppo, M. P., Wood, N., Golosarsky, B., et al. (1999). Heart rate dynamics before spontaneous onset of ventricular fibrillation in patients with healed myocardial infarcts. Am. J. Cardiol. 83, 880–884. doi: 10.1016/S0002-9149(98)01068-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Moghtadaei, M., Jansen, H. J., Mackasey, M., Rafferty, S. A., Bogachev, O., Sapp, J. L., et al. (2016). The impacts of age and frailty on heart rate and sinoatrial node function. J. Physiol. (Lond.) 594, 7105–7126. doi: 10.1113/JP272979

PubMed Abstract | CrossRef Full Text | Google Scholar

Nazir, S. A., and Lab, M. J. (1996). Mechanoelectric feedback in the atrium of the isolated guinea-pig heart. Cardiovasc. Res. 32, 112–119. doi: 10.1016/S0008-6363(96)00077-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Pikkujamsa, S. M., Makikallio, T. H., Sourander, L. B., Raiha, I. J., Puukka, P., Skytta, J., et al. (1999). Cardiac interbeat interval dynamics from childhood to senescence. Circulation 100, 393–399. doi: 10.1161/01.CIR.100.4.393

PubMed Abstract | CrossRef Full Text | Google Scholar

Piskorski, J., and Guzik, P. (2011). The structure of heart rate asymmetry: deceleration and acceleration runs. Physiol. Meas. 32, 1011–1023. doi: 10.1088/0967-3334/32/8/002

PubMed Abstract | CrossRef Full Text | Google Scholar

Porta, A., Tobaldini, E., Guzzetti, S., Furlan, R., Montano, N., and Gnecchi-Ruscone, T. (2007). Assessment of cardiac autonomic modulation during graded head-up tilt by symbolic analysis of heart rate variability. Am. J. Physiol. Heart Circ. Physiol. 293, H702–H708. doi: 10.1152/ajpheart.00006.2007

PubMed Abstract | CrossRef Full Text | Google Scholar

Ravelo-Garcia, A. G., Saavedra-Santana, P., Julia-Serda, G., Navarro-Mesa, J. L., Navarro-Esteva, J., Alvarez-Lopez, X., et al. (2014). Symbolic dynamics marker of heart rate variability combined with clinical variables enhance obstructive sleep apnea screening. Chaos 24:024404. doi: 10.1063/1.4869825

PubMed Abstract | CrossRef Full Text | Google Scholar

Stein, P. K. (2002). Heart rate variability is confounded by the presence of erratic sinus rhythm. Comput. Cardiol. 26, 669–672. doi: 10.1109/CIC.2002.1166861

CrossRef Full Text | Google Scholar

Wiklund, U., Hornsten, R., Karlsson, M., Suhr, O. B., and Jensen, S. M. (2008). Abnormal heart rate variability and subtle atrial arrhythmia in patients with familial amyloidotic polyneuropathy. Ann. Noninvasive Electrocardiol. 13, 249–256. doi: 10.1111/j.1542-474X.2008.00228.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, A. C., Hseu, S. S., Yien, H. W., Goldberger, A. L., and Peng, C. K. (2003). Linguistic analysis of the human heartbeat using frequency and rank order statistics. Phys. Rev. Lett. 90:108103. doi: 10.1103/PhysRevLett.90.108103

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: aging, coronary artery disease, fragmentation, heart rate variability, symbolic dynamics, vagal tone

Citation: Costa MD, Davis RB and Goldberger AL (2017) Heart Rate Fragmentation: A Symbolic Dynamical Approach. Front. Physiol. 8:827. doi: 10.3389/fphys.2017.00827

Received: 04 August 2017; Accepted: 06 October 2017;
Published: 14 November 2017.

Edited by:

Zbigniew R. Struzik, The University of Tokyo, Japan

Reviewed by:

Dirk Cysarz, Witten/Herdecke University, Germany
Heikki Veli Huikuri, University of Oulu, Finland

Copyright © 2017 Costa, Davis and Goldberger. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Madalena D. Costa, mcosta3@bidmc.harvard.edu

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.