Atrial fibrillation (AF) is regarded as a complex arrhythmia, with one or more co-existing mechanisms, resulting in an intricate structure of atrial activations. Fractionated atrial electrograms (AEGs) were thought to represent arrhythmogenic tissue and hence have been suggested as targets for radiofrequency ablation. However, current methods for ablation target identification have resulted in suboptimal outcomes for persistent AF (persAF) treatment, possibly due to the complex spatiotemporal dynamics of these mechanisms. In the present work, we sought to characterize the dynamics of atrial tissue activations from AEGs collected during persAF using recurrence plots (RPs) and recurrence quantification analysis (RQA). 797 bipolar AEGs were collected from 18 persAF patients undergoing pulmonary vein isolation (PVI). Automated AEG classification (normal vs. fractionated) was performed using the CARTO criteria (Biosense Webster). For each AEG, RPs were evaluated in a phase space estimated following Takens' theorem. Seven RQA variables were obtained from the RPs: recurrence rate; determinism; average diagonal line length; Shannon entropy of diagonal length distribution; laminarity; trapping time; and Shannon entropy of vertical length distribution. The results show that the RQA variables were significantly affected by PVI, and that the variables were effective in discriminating normal vs. fractionated AEGs. Additionally, diagonal structures associated with deterministic behavior were still present in the RPs from fractionated AEGs, leading to a high residual determinism, which could be related to unstable periodic orbits and suggesting a possible chaotic behavior. Therefore, these results contribute to a nonlinear perspective of the spatiotemporal dynamics of persAF.

Objective: Atrial flutter (AFl) is a common arrhythmia that can be categorized according to different self-sustained electrophysiological mechanisms. The non-invasive discrimination of such mechanisms would greatly benefit ablative methods for AFl therapy as the driving mechanisms would be described prior to the invasive procedure, helping to guide ablation. In the present work, we sought to implement recurrence quantification analysis (RQA) on 12-lead ECG signals from a computational framework to discriminate different electrophysiological mechanisms sustaining AFl. Methods: 20 different AFl mechanisms were generated in 8 atrial models and were propagated into 8 torso models via forward solution, resulting in 1,256 sets of 12-lead ECG signals. Principal component analysis was applied on the 12-lead ECGs, and six RQA-based features were extracted from the most significant principal component scores in two different approaches: individual component RQA and spatial reduced RQA. Results: In both approaches, RQA-based features were significantly sensitive to the dynamic structures underlying different AFl mechanisms. Hit rate as high as 67.7% was achieved when discriminating the 20 AFl mechanisms. RQA-based features estimated for a clinical sample suggested high agreement with the results found in the computational framework. Conclusion: RQA has been shown an effective method to distinguish different AFl electrophysiological mechanisms in a non-invasive computational framework. A clinical 12-lead ECG used as proof of concept showed the value of both the simulations and the methods. Significance: The non-invasive discrimination of AFl mechanisms helps to delineate the ablation strategy, reducing time and resources required to conduct invasive cardiac mapping and ablation procedures.