Non-linear time series analysis: methods and applications to atrial fibrillation

We apply methods from non-linear statistical time series analysis to characterize electrograms of atrial fibrillation. These are based on concepts originating from the theory of non-linear dynamical systems and use the empirical reconstruction density in reconstructed phase space. Application of these methods is not restricted to deterministic chaos but is valid in a general time series context. We illustrate this by applying three recently proposed non-linear time series methods to fibrillation electrograms: 1) a test for time reversibility in atrial electrograms during paroxysmal atrial fibrillation in patients; 2) a test to detect differences in the dynamical behaviour during the pharmacological conversion of sustained atrial fibrillation in instrumented conscious goats; 3) a test for general Granger causality to identify couplings and information transport in the atria during fibrillation. We conclude that a characterization of the dynamics via the reconstruction density offers a useful framework for the non-linear analysis of electrograms of atrial fibrillation.

Spatial correlation analysis of the pharmacological conversion of sustained atrial fibrillation in conscious goats by cibenzoline

The nonlinear spatial redundancy and the linear spatial correlation function were used to investigate to what extent non-linearity was involved in the coupling of atrial regions and how organization in activation patterns of sustained atrial fibrillation (AF) had been modified by administration of the class IC agent cibenzoline in the experimental model of sustained AF in instrumented conscious goats. Electrograms were measured in five goats during sustained AF and when the fibrillation interval had been prolonged to about 25%, 50% and 85% (CIB25, CIB50, CIB85) with respect to control. The nonlinear association length and linear correlation length were estimated along the principal axes of two-dimensional correlation maps estimated from the spatial redundancy and the spatial correlation function, respectively. The estimated short axis association length in the right atrium increased already shortly after the start of infusion (CIB25, +61%), and remained significantly different from control during the experiment, including the effects of non-simultaneous interaction. At CIB85 the association length had almost become twice as long with respect to control (increase from 16 to 29 mm, 89%), while in the left atrium changes were less pronounced (increase from 9 to 12 mm, +32%). The linearized association length which was estimated using multivariate surrogate data increased more gradually and was less sensitive to changes in spatial organization. The results of the spatial correlation analysis suggest that the drug-induced nonlinearity in the spatio-temporal dynamics of sustained AF is related to activation patterns which are characterized by extended uniformly propagating fibrillation wavefronts (AF type I). We conclude that cibenzoline enhanced the spatial organization of sustained AF associated with a transition from type II to type I AF activation patterns. This may destabilize the perpetuation of AF since an increase in association length is equivalent to a reduction of atrial tissue mass available to support reentrant circuits. The results are consistent with the hypothesis that larger association lengths result from fewer and larger reentrant circuits. It is argued that effects of diminished curvature of fibrillation wavefronts are anti-arrhythmic under conditions of suppressed excitability imposed by cibenzoline. Termination of AF may be mediated by a mechanism resembling a bifurcation of the dynamics which sets in when the ends of fractionated wavefronts cannot sufficiently curve anymore to maintain a positive balance of newly generated wavelets needed to sustain AF.

Spatial correlation analysis of atrial activation patterns during sustained atrial fibrillation in conscious goats

In this study we applied both linear and nonlinear spatial correlation measures to characterize epicardial activation patterns of sustained atrial fibrillation in instrumented conscious goats. It was investigated if nonlinearity was involved in the spatial coupling of atrial regions and to what extent fibrillation was organized in the experimental model of sustained atrial fibrillation (AF) in instrumented goats. Data were collected in five goats during experiments to convert AF by continuous infusion of cibenzoline. Spatial organization during AF was quantified with the linear spatial cross correlation function and the nonlinear spatial cross redundancy which was calculated using the Grassberger-Procaccia correlation integral. Two different types of correlation were evaluated to distinguish simultaneous interaction from non-simultaneous interaction, for instance resulting from propagation of fibrillation waves. The nonlinear association length and the linear correlation length were estimated along the principal axes of iso-correlation contours in two-dimensional correlation maps of the nonlinear spatial redundancy and the linear spatial correlation function, respectively. To quantitatively assess the degree of nonlinearity, the association length was also estimated from the linearized spatial redundancy using multivariate surrogate data. The differences between the nonlinear and linearized association lengths indicated that a nonlinear component in the spatial organization of AF predominantly existed in the right atrium. The degree of organization characterized by association length along the short principal axis was higher in the right atrium (15 +/- 7 mm) than in the left atrium (8 +/- 4 mm). The spatial extension of coherent atrial patches was estimated from a surface of association equal to the area spanned by the principal axes of iso-correlation contours from the redundancy, including the effects from non-simultaneous interaction. Interpreting this area as the spatial domain of a fibrillation wavelet, the results suggest that the mapped region was activated on average by two wavelets in the left atrium and by one wavelet in the right atrium. Therefore, the activation pattern of sustained AF in goats was relatively organized, consistent with type II of AF. It is suggested that the surface of association is a measure of the number of independent wavelets present in the atria during sustained AF, and that larger association lengths result from fewer and larger reentrant circuits.

Nonlinear analysis of epicardial atrial electrograms of electrically induced atrial fibrillation in man

INTRODUCTION

We applied methods from the theory of nonlinear dynamics to characterize unipolar epicardial right atrial electrograms of electrically induced atrial fibrillation (AF) in man.

METHODS AND RESULTS

Electrograms were selected from a high-density mapping study, which confirmed the existence of at least 3 different types of induced AF (types I, II, and III) in patients undergoing open chest surgery. We analyzed sets of 5 electrograms (4 sec, sampling frequency 1 kHz, resolution 8 bits) in 9 patients (AF type I, n = 3; type II, n = 3; type III, n = 3). The Grassberger-Procaccia method was applied to estimate the correlation dimension and correlation entropy from the electrograms. In 2 patients (AF type I) some electrograms (2 of 5 and 3 of 5, respectively) showed scaling at normalized distances ranging from 0.2 to 0.5 in phase space. Correlation dimension D ranged from 1.8 to 3.2 and correlation entropy K from 2.2 to 3.8 nats/sec. The patients were ranked according to increasing coarse-grained correlation dimension Dcg (range 3.7 to 7.9) and coarse-grained correlation entropy Kcg (range 5.6 to 18.6 nats/sec). The method of surrogate data was applied to detect nonlinearity in the electrograms. Using the correlation integral as test statistic, it could be excluded that electrograms of type I AF have been generated by linear stochastic dynamics. Episodes of sinus rhythm (D ranging from 1.0 to 5.1 and K from 2.0 to 8.6 nats/sec) and induced atrial flutter (D ranging from 2.7 to 4.2 and K from 2.2 to 4.2 nats/sec) in 2 different patients showed features of low-dimensional chaos.

CONCLUSION

Nonlinear analysis discriminated between electrograms during electrically induced AF in humans. The results are consistent with a classification of AF into 3 types based on the spatiotemporal complexity of right atrial activation patterns.