A high-temporal resolution algorithm for quantifying organization during atrial fibrillation

Atrial fibrillation (AF) has been described as a "random" or "chaotic" rhythm. Evidence suggests that AF may have transient episodes of temporal and spatial organization. We introduce a new algorithm that quantifies AF organization by the mean-squared error (MSE) in the linear prediction between two cardiac electrograms. This algorithm calculates organization at a finer temporal resolution. (approximately 300 ms) than previously published algorithms. Using canine atrial epicardial mapping data, we verified that the MSE algorithm showed nonfibrillatory rhythms to be significantly more organized than fibrillatory rhythms (p < .00001). Further, we compared the sensitivity of MSE to that of two previously published algorithms by analyzing AF with simulated noise and AF manipulated with vagal stimulation or by adenosine administration to alter the character of the AF. MSE performed favorably in the presence of noise. While all three algorithms distinguished between low and high vagal AF, MSE was the most sensitive in its discrimination. Only MSE could distinguish baseline AF from AF with adenosine. We conclude that our algorithm can distinguish different levels of organization during AF with a greater temporal resolution and sensitivity than previously described algorithms. This algorithm could lead to new ways of analyzing and understanding AF as well as improved techniques in AF therapy.     

Discrimination of atrial fibrillation from regular atrial rhythmsby spatial precision of local activation direction

This study tests the hypothesis that atrial fibrillation (AFib) can be discriminated from regular atrial rhythms by a measure of the variation in local activation direction. Human endocardial atrial recordings of AFib, sinus rhythm, atrial flutter, and supraventricular tachycardia were collected using a catheter with orthogonally placed electrodes, and the direction of each activation was calculated using methods previously described by our laboratory. Each recording was divided into segments containing 100 activations, and the spatial precision for each segment was calculated in three dimensions, as well as in each of the three two-dimensional (2-D) planes. The three-dimensional (3-D) spatial precision for 1161 segments of AFib in 11 recordings ranged from 0.09-0.85 (mean=0.45), whereas the spatial precision for 138 segments of regular rhythms in 28 recordings was ⩾0.91 in all but four instances. The 2-D spatial precision values overlapped for all rhythms. The results indicate that 3-D spatial precision of local activation direction is a useful discriminator of AFib     

A novel method for detection of the transition between atrial fibrillation and sinus rhythm

Automatic detection of atrial fibrillation (AF) for AF diagnosis, especially for AF monitoring, is necessarily desirable for clinical therapy. In this study, we proposed a novel method for detection of the transition between AF and sinus rhythm based on RR intervals. First, we obtained the delta RR interval distribution difference curve from the density histogram of delta RR intervals, and then detected its peaks, which represented the AF events. Once an AF event was detected, four successive steps were used to classify its type, and thus, determine the boundary of AF: 1) histogram analysis; 2) standard deviation analysis; 3) numbering aberrant rhythms recognition; and 4) Kolmogorov-Smirnov (K-S) test. A dataset of 24-h Holter ECG recordings (n = 433) and two MIT-BIH databases (MIT-BIH AF database and MIT-BIH normal sinus rhythm (NSR) database) were used for development and evaluation. Using the receiver operating characteristic curves for determining the threshold of the K-S test, we have achieved the highest performance of sensitivity and specificity (SP) (96.1% and 98.1%, respectively) for the MIT-BIH AF database, compared with other previously published algorithms. The SP was 97.9% for the MIT-BIH NSR database.     

Computer modeling of ventricular rhythm during atrial fibrillation and ventricular pacing

We propose a unified atrial fibrillation (AF)-ventricular pacing (VP) (AF-VP) model to demonstrate the effects of VP on the ventricular rhythm during atrial fibrillation AF. In this model, the AV junction (AVJ) is treated as a lumped structure characterized by refractoriness and automaticity. Bombarded by random AF impulses, the AVJ can also be invaded by the VP-induced retrograde wave. The model includes bidirectional conduction delays in the AVJ and ventricle. Both refractory period and conduction delay of the AVJ are dependent upon its recovery time. The electrotonic modulation by blocked impulses is also considered in the model. Our simulations show that, with proper parameter settings, the present model can account for most principal statistical properties of the RR intervals during AF. We further demonstrate that the AV conduction property and the ventricular rate in AF depend on both AF rate and the degree of electrotonic modulation in the AVJ. Finally, we show that multilevel interactions between AF and VP can generate various patterns of ventricular rhythm that are consistent with previous experimental observations.     

An automatic system for the analysis and classification of human atrial fibrillation patterns from intracardiac electrograms

This paper presents an automatic system for the analysis and classification of atrial fibrillation (AF) patterns from bipolar intracardiac signals. The system is made up of: 1) a feature-extraction module that defines and extracts a set of measures potentially useful for characterizing AF types on the basis of their degree of organization; 2) a feature-selection module (based on the Jeffries-Matusita distance and a branch and bound search algorithm) identifying the best subset of features for discriminating different AF types; and 3) a support vector machine technique-based classification module that automatically discriminates the AF types according to the Wells' criteria. The automatic system was applied on 100 intracardiac AF signal strips and on a selection of 11 representative features, demonstrating: a) the possibility to properly identify the most significant features for the discrimination of AF types; b) higher accuracy (97.7% using the seven most informative features) than the traditional maximum likelihood classifier; and c) effectiveness in AF classification also with few training samples (accuracy = 88.3% with only five training signals). Finally, the system identifies a combination of indices characterizing changes of morphology of atrial activation waves and perturbation of the isoelectric line as the most effective in separating the AF types.     

Convolutive blind source separation algorithms applied to the electrocardiogram of atrial fibrillation: study of performance

The analysis of the surface electrocardiogram (ECG) is the most extended noninvasive technique in medical diagnosis of atrial fibrillation (AF). In order to use the ECG as a tool for the analysis of AF, we need to separate the atrial activity (AA) from other cardioelectric signals. In this matter, statistical signal processing techniques, like blind source separation (BSS), are able to perform a multilead statistical analysis with the aim to obtain the AA. Linear BSS techniques can be divided in two groups depending on the mixing model: algorithms where instantaneous mixing of sources is assumed, and convolutive BSS (CBSS) algorithms. In this work, a comparison of performance between one relevant CBSS algorithm, namely Infomax, and one of the most effective independent component analysis (ICA) algorithms, namely FastICA, is developed. To carry out the study, pseudoreal AF ECGs have been synthesized by adding fibrillation activity to normal sinus rhythm. The algorithm performances are expressed by two indexes: the signal to interference ratio (SIRAA) and the cross-correlation (RAA) between the original and the estimated AA. Results empirically prove that the instantaneous mixing model is the one that obtains the best results in the AA extraction, given that the mean SIRAA obtained by the FastICA algorithm (37.6 +/- 17.0 dB) is higher than the main SIRAA obtained by Infomax (28.5 +/- 14.2 dB). Also the RAA obtained by FastICA (0.92 +/- 0.13) is higher than the one obtained by Infomax (0.78 +/- 0.16).     

Detection of atrial activity from high-voltage leads of implantableventricular defibrillators using a cancellation technique

The inability to detect atrial activity limits implantable ventricular cardioverter defibrillators (ICD) in discriminating tachycardias and can result in inappropriate therapy. This study attempted to detect atrial activity on the wide-spaced bipole signals formed by the high-voltage (HV) leads of the ICD during device implantation and to develop an algorithm for the detection of atrial fibrillation (AFib) from these signals. The authors used a method that cancelled ventricular and correlated atrial activity from the HV lead signals and measured frequency and amplitude distribution information to discriminate sinus rhythm (SR) and AFib segments. The authors analyzed 186 data segments from 21 patients (six AFib, 14 SR, one AFib and SR). For individual segments in this data set, the sensitivity of the algorithm was 78%, specificity 92.65%, positive and negative predictive values 79.59 and 91.97%, respectively. These results demonstrate that atrial activity is present in the HV lead signals, and AFib detection can be achieved in many, but not all cases, using information currently available to ICDs. Prior work from surface electrocardiograms suggests that this algorithm can function during ventricular tachycardias. However, specificity of the algorithm is not high enough for clinical use     

Frequency domain algorithm for quantifying atrial fibrillation organization to increase defibrillation efficacy

We hypothesized that frequency domain analysis of an interatrial atrial fibrillation (AF) electrogram would show a correlation of the variance of the signal and the amplitude of harmonic peaks with the periodicity and morphology (organization) of the AF signal and defibrillation efficacy. We sought to develop an algorithm that would provide a high-resolution measurement of the changes in the spatiotemporal organization of AF. AF was initiated with burst atrial pacing in ten dogs. The atrial defibrillation threshold (ADFT50) was determined, and defibrillation was repeated at the ADFT50. Bipolar electrograms from the shocking electrodes were acquired immediately preshock, digitally filtered, and a FFT was performed. The organization index (OI) was calculated as the ratio of the area under the first four harmonic peaks to the total area of the spectrum. For a 4-s window, the mean OI was 0.505 +/- 0.087 for successful shocks, versus 0.352 +/- 0.068 for unsuccessful shocks (p < 0.001). Receiver operator characteristic (ROC) curve analysis was used to determine the optimal sampling window for predicting successful shocks. The area of the ROC curve was 0.8 for a 1-s window, and improved to 0.9 for a 4-s window. We conclude that the spectrum of an AF signal contains information relating to its organization, and can be used in predicting a successful defibrillation.     

A detector for a chronic implantable atrial tachyarrhythmia monitor

Continuous long-term monitoring of atrial fibrillation (AF) and tachycardia (AT) is an unmet clinical need, which could be met with a chronically-implanted monitor. Improved therapeutic decisions based on accurate monitoring of parameters, such as daily AF/AT burden (hours/ day) may lead to improvements in clinical outcomes such as reduction in hospitalizations, symptoms, and strokes. This paper describes an AF/AT detector that detects AF as well as AT with an irregular ventricular response, and a supplementary AT detector for AT with more regular ventricular response. Seven databases with significant durations of AF, AT, and sinus rhythm were used to evaluate the performance of the detectors. All patient records with AF (N = 124) were detected by the AF/AT detector to have AF/AT burden with a mean, median, and 75 percentile of absolute error in burden detection of 8.8, 0, and 4 min, respectively. In patients having AF burden (= or > 10 min), the AF/AT detector was found to have burden accuracy within 20% of true burden in 96% of patients. The specificity was 94%, defined as follows: in patient records without AF/AT (N = 174), the percentage with AF/AT burden = or < 10 min in the 24-h recordings. The AF/AT detector underestimatesAT burden, thus degrading performance, in patients with significant amounts of AT with more regular ventricular response. The supplementary AT detector reduces the underestimation of AT while overestimating burden in patients without a significant amount of AT. The detectors described here could be implemented in an implantable monitor for accurate long-term AF/AT monitoring.     

Temporal and spatial phase analyses of the electrocardiogram stratify intra-atrial and intra-ventricular organization

We hypothesized that electrocardiogram (ECG) spatial phase analysis would define a spectrum of intracardiac organization from atrial fibrillation (AF), nonisthmus-dependent and isthmus-dependent atrial flutter (AFL) to supraventricular tachycardias (SVT), and similarly for ventricular arrhythmias. We analyzed arrhythmia ECGs of 33 patients with isthmus (n = 9) and nonisthmus (n = 5) dependent AFL and SVT: atrial (n = 3), atrioventricular nodal (n = 3), and orthodromic reciprocating (n = 3) tachycardias, as well as AF (n = 5), ventricular tachycardia (monomorphic, VT-MM; n = 7), and fibrillation (VF; n = 3). ECG spatial phase was considered coherent when the correlation coefficient of an atrial (or ventricular) template to its ECG over time maintained a constant relationship in XY, XZ, and YZ planes. Regularity was quantified spectrally from ECG and correlation series. Spatial coherence occurred in 9/9 cases of isthmus--but only 1/5 of cases of nonisthmus-dependent AFL (p < 0.01; chi2). All showed one dominant spectral peak (temporal coherence). In AF, spatial phase was inconsistent in all planes and spectra were broad band. Temporal and spatial coherence occurred in other SVT. VT-MM maintained spatial phase and a single spectral peak, while VF displayed neither. Our conclusions are that temporal and spatial phase analysis from the ECG stratifies intra-atrial and intra-ventricular organization and reveals subtle variability lost on visual inspection.     

A new method for analysis of atrial activation during chronicatrial fibrillation in man

To further clarify the mechanisms maintaining chronic atrial fibrillation (CAF), a method identifying preferable activation patterns of the atria during fibrillation, by time averaging of multiple discrete excitation vectors, was developed. Repeated recordings, each of 56 atrial bipolar electrograms simultaneously acquired during 8 s, were made at multiple sites in the right atrial free wall and the left atrial appendage in 16 patients with CAF using a 2.17×3.54 cm electrode array. The local activation times (LAT's) in each recording were estimated as the median activation time at the respective measurement point. By calculating the time difference between the LAT's at adjacent measurement points in two spatial dimensions, a direction vector was created for each activation wave passing each set of measurement points, a total of 42 sets. By time averaging of the individual direction vectors (typically n=55) at each set of measurement points, preferable activation patterns were determined. Three types of activation patterns were found: 1) inconsistent activation (n=5), 2) consistent activation with preferential propagation directions (n=7) and 3) consistent activation with impulses originating from a localizable site within the recording area (n=4). All activation patterns were reproducible and the two latter patterns were proven significant using statistical tests. It is concluded that this new method is useful in further clarification of the mechanisms involved in the maintenance of atrial fibrillation     

Conduction velocity restitution of the human atrium--an efficient measurement protocol for clinical electrophysiological studies

Conduction velocity (CV) and CV restitution are important substrate parameters for understanding atrial arrhythmias. The aim of this work is to (i) present a simple but feasible method to measure CV restitution in-vivo using standard circular catheters, and (ii) validate its feasibility with data measured during incremental pacing. From five patients undergoing catheter ablation, we analyzed eight datasets from sinus rhythm and incremental pacing sequences. Every wavefront was measured with a circular catheter and the electrograms were analyzed with a cosine-fit method that calculated the local CV. For each pacing cycle length, the mean local CV was determined. Furthermore, changes in global CV were estimated from the time delay between pacing stimulus and wavefront arrival. Comparing local and global CV between pacing at 500 and 300 ms, we found significant changes in seven of eight pacing sequences. On average, local CV decreased by 20 ± 15% and global CV by 17 ± 13%. The method allows for in-vivo measurements of absolute CV and CV restitution during standard clinical procedures. Such data may provide valuable insights into mechanisms of atrial arrhythmias. This is important both for improving cardiac models and also for clinical applications, such as characterizing arrhythmogenic substrates during sinus rhythm.     

Spatiotemporal blind source separation approach to atrial activity estimation in atrial tachyarrhythmias

The analysis and characterization of atrial tachyarrhythmias requires, in a previous step, the extraction of the atrial activity (AA) free from ventricular activity and other artefacts. This contribution adopts the blind source separation (BSS) approach to AA estimation from multilead electrocardiograms (ECGs). Previously proposed BSS methods for AA extraction-e.g., independent component analysis (ICA)-exploit only the spatial diversity introduced by the multiple spatially-separated electrodes. However, AA typically shows certain degree of temporal correlation, with a narrowband spectrum featuring a main frequency peak around 3.5-9 Hz. Taking advantage of this observation, we put forward a novel two-step BSS-based technique which exploits both spatial and temporal information contained in the recorded ECG signals. The spatiotemporal BSS algorithm is validated on simulated and real ECGs from a significant number of atrial fibrillation (AF) and atrial flutter (AFL) episodes, and proves consistently superior to a spatial-only ICA method. In simulated ECGs, a new methodology for the synthetic generation of realistic AF episodes is proposed, which includes a judicious comparison between the known AA content and the estimated AA sources. Using this methodology, the ICA technique obtains correlation indexes of 0.751, whereas the proposed approach obtains a correlation of 0.830 and an error in the estimated signal reduced by a factor of 40%. In real ECG recordings, we propose to measure performance by the spectral concentration (SC) around the main frequency peak. The spatiotemporal algorithm outperforms the ICA method, obtaining a SC of 58.8% and 44.7%, respectively.     

An atrioventricular node model for analysis of the ventricular response during atrial fibrillation

This paper introduces a model of the atrioventricular node function during atrial fibrillation (AF), and describes the related ECG-based estimation method. The proposed model is defined by parameters that characterize the arrival rate of atrial impulses, the probability of an impulse choosing either one of the two atrioventricular nodal pathways, the refractory periods of these pathways, and the prolongation of the refractory periods. These parameters are estimated from the RR intervals using maximum likelihood estimation, except for the shorter refractory period which is estimated from the RR interval Poincaré plot, and the mean arrival rate of atrial impulses by the AF frequency. Simulations indicated that 200-300 RR intervals are generally needed for the estimates to be accurate. The model was evaluated on 30-min ECG segments from 36 AF patients. The results showed that 88% of the segments can be accurately modeled when the estimated probability density function (PDF) and an empirical PDF were at least 80% in agreement. The model parameters were estimated during head-up tilt test to assess differences caused by sympathetic stimulation. Both refractory periods decreased as a result of stimulation, and the likelihood of an impulse choosing the pathway with the shorter refractory period increased.     

Evaluation of techniques for recognition of ventricular arrhythmiasby implanted devices

Implantable devices that provide antitachycardia and defibrillation capability currently have limited ability to distinguish among different cardiac rhythms. We have investigated three methods of electrogram analysis: rate, irregularity, and amplitude distribution. In 35 episodes in 19 patients, we applied these three algorithms to 15 s recorded passages of ventricular electrograms during supraventricular tachycardia (N = 11), ventricular tachycardia (N = 11), and ventricular fibrillation (N = 13). Each was individually paired with a recording of sinus rhythm from the same patient. All recordings were obtained during standard electrophysiologic testing. Each algorithm was successful at distinguishing the tachyarrhythmias from sinus rhythm at one or more levels of algorithm parameterization. Rate alone discriminated supraventricular tachycardia from ventricular fibrillation but did not distinguish between supraventricular and ventricular tachycardia. Rate combined with irregularity distinguished between ventricular tachycardia and ventricular fibrillation, but did not discriminate between ventricular and supraventricular tachycardia. Although the amplitude distribution algorithm was unable to separate perfectly any of the three tachyarrhythmias, it provided the best performance in separating supraventricular and ventricular tachycardia (82 percent sensitivity and specificity). We conclude that algorithms based on rate, irregularity, and amplitude distribution analysis of ventricular electrograms may distinguish sinus rhythm from tachyarrhythmias, but may not distinguish among tachyarrhythmias.     

Preventive ablation strategies in a biophysical model of atrial fibrillation based on realistic anatomical data.

Ablation strategies to prevent episodes of paroxysmal atrial fibrillation (AF) have been subject to many clinical studies. The issues mainly concern pattern and transmurality of the lesions. This paper investigates ten different ablation strategies on a multilayered 3-D anatomical model of the atria with respect to 23 different setups of AF initiation in a biophysical computer model. There were 495 simulations carried out showing that circumferential lesions around the pulmonary veins (PVs) yield the highest success rate if at least two additional linear lesions are carried out. The findings compare with clinical studies as well as with other computer simulations. The anatomy and the setup of ectopic beats play an important role in the initiation and maintenance of AF as well as the resulting therapy. The computer model presented in this paper is a suitable tool to investigate different ablation strategies. By including individual patient anatomy and electrophysiological measurement, the model could be parameterized to yield an effective tool for future investigation of tailored ablation strategies and their effects on atrial fibrillation.     

Development of an automatic implanted drug infusion system for themanagement of cardiac arrhythmias

Conventional management of cardiac arrythmias relies on oral drug therapy which minimizes recurrence of the arrhythmia, but risks unpleasant side effects and even long-term toxicity. The authors propose acute management instead, from an implanted drug pump which automatically senses the onset of arrhythmia, delivers a pharmacokinetically-based infusion to terminate the episode, and discontinues drug delivery until the next occurrence. A bedside system consisting of a personal computer and conventional intravenous pump has been developed and tested in five dogs and 24 patients during a catheter electrophysiologic study. After detection of the arrhythmia plasma levels of the antiarrhythmia drug rose immediately to the therapeutic range and were subsequently well-controlled for 30 to 60 minutes. In all five dogs and in seven of the eight patients in whom atrial fibrillation was induced during the study, conversion to normal rhythm occurred within fifteen minutes     

Waveform Characterization of Atrial Fibrillation Using Phase Information

A novel method for characterization of f-wave morphology in atrial fibrillation (AF) is presented. The method decomposes atrial activity into fundamental and harmonic components, dividing each component into short blocks for which the amplitudes, frequencies, and phases are estimated. The phase delays between the fundamental and each of the harmonics, here referred to as harmonic phase relationships, are used as features of f-wave morphology. The estimated waves are clustered into typical morphologic patterns. The performance of the method is illustrated by simulated signals, ECG signals recorded from 36 patients with organized AF, and an ECG signal recorded during drug loading with flecainide. The results show that the method can distinguish a wide variety of f-wave morphologies, and that typical morphologies can be established for further analysis of AF.      

Propagation pattern analysis during atrial fibrillation based on sparse modeling

In this study, sparse modeling is introduced for the estimation of propagation patterns in intracardiac atrial fibrillation (AF) signals. The estimation is based on the partial directed coherence function, derived from fitting a multivariate autoregressive model to the observed signal using least-squares (LS) estimation. The propagation pattern analysis incorporates prior information on sparse coupling as well as the distance between the recording sites. Two optimization methods are employed for estimation of the model parameters, namely, the adaptive group least absolute selection and shrinkage operator (aLASSO), and a novel method named the distance-adaptive group LASSO (dLASSO). Using simulated data, both optimization methods were superior to LS estimation with respect to detection and estimation performance. The normalized error between the true and estimated model parameters dropped from 0.20 ± 0.04 for LS estimation to 0.03 ± 0.01 for both aLASSO and dLASSO when the number of available data samples exceeded the number of model parameters by a factor of 5. For shorter data segments, the error reduction was more pronounced and information on the distance gained in importance. Propagation pattern analysis was also studied on intracardiac AF data, the results showing that the identification of propagation patterns is substantially simplified by the sparsity assumption.     

Comparison of atrial signal extraction algorithms in 12-lead ECGs with atrial fibrillation

Analysis of atrial rhythm is important in the treatment and management of patients with atrial fibrillation. Several algorithms exist for extracting the atrial signal from the electrocardiogram (ECG) in atrial fibrillation, but there are few reports on how well these techniques are able to recover the atrial signal. We assessed and compared three algorithms for extracting the atrial signal from the 12-lead ECG. The 12-lead ECGs of 30 patients in atrial fibrillation were analyzed. Atrial activity was extracted by three algorithms, Spatiotemporal QRST cancellation (STC), principal component analysis (PCA), and independent component analysis (ICA). The amplitude and frequency characteristics of the extracted atrial signals were compared between algorithms and against reference data. Mean (standard deviation) amplitude of QRST segments of V1 was 0.99 (0.54) mV, compared to 0.18 (0.11) mV (STC), 0.19 (0.13) mV (PCA), and 0.29 (0.22) mV (ICA). Hence, for all algorithms there were significant reductions in the amplitude of the ventricular activity compared with that in V1. Reference atrial signal amplitude in V1 was 0.18 (0.11) mV, compared to 0.17 (0.10) mV (STC), 0.12 (0.09) mV (PCA), and 0.18 (0.13) mV (ICA) in the extracted atrial signals. PCA tended to attenuate the atrial signal in these segments. There were no significant differences for any of the algorithms when comparing the amplitude of the reference atrial signal with that of the extracted atrial signals in segments in which ventricular activity had been removed. There were no significant differences between algorithms in the frequency characteristics of the extracted atrial signals. There were discrepancies in amplitude and frequency characteristics of the atrial signal in only a few cases resulting from notable residual ventricular activity for PCA and ICA algorithms. In conclusion, the extracted atrial signals from these algorithms exhibit very similar amplitude and frequency characteristics. Users of these algorithms should be observant of residual ventricular activities which can affect the analysis of the fibrillatory waveform in clinical practice.