Skip Navigation

Europace 2006 8(11):911-926; doi:10.1093/europace/eul113
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (8)
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Bollmann, A.
Right arrow Articles by Sörnmo, L.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Bollmann, A.
Right arrow Articles by Sörnmo, L.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The European Society of Cardiology 2006. All rights reserved. For Permissions, please e-mail: journals.permissions@oxfordjournals.org


ATRIAL FIBRILLATION

Analysis of surface electrocardiograms in atrial fibrillation: techniques, research, and clinical applications

Andreas Bollmann1,*, Daniela Husser1, Luca Mainardi2, Federico Lombardi3, Philip Langley4, Alan Murray4, José Joaquín Rieta5, José Millet5, S. Bertil Olsson6, Martin Stridh6 and Leif Sörnmo6

1 Department of Cardiology, Otto-von-Guericke-University Magdeburg, Leipziger Str. 44, 39120 Magdeburg, Germany; 2 Polytechnic University of Milan, Milan, Italy; 3 University of Milan, Milan, Italy; 4 University of Newcastle upon Tyne, Newcastle, UK; 5 Polytechnic University of Valencia, Valencia, Spain; 6 Lund University, Lund, Sweden

Manuscript submitted 7 February 2006. Accepted after revision 11 June 2006.

* Corresponding author. Tel: +49 391 67 13229; fax: +49 391 67 13202. E-mail address: andreas.bollmann{at}medizin.uni-magdeburg.de


    Abstract
 Top
 Abstract
 Introduction
 Part I: Characterization of...
 Part II: Characterization of...
 Conclusions
 Acknowledgements
 References
 
Atrial fibrillation (AF) is the most common arrhythmia encountered in clinical practice. Neither the natural history of AF nor its response to therapy is sufficiently predictable by clinical and echocardiographic parameters. The purpose of this article is to describe technical aspects of novel electrocardiogram (ECG) analysis techniques and to present research and clinical applications of these methods for characterization of both the fibrillatory process and the ventricular response during AF. Atrial fibrillatory frequency (or rate) can reliably be assessed from the surface ECG using digital signal processing (extraction of atrial signals and spectral analysis). This measurement shows large inter-individual variability and correlates well with intra-atrial cycle length, a parameter which appears to have primary importance in AF maintenance and response to therapy. AF with a low fibrillatory rate is more likely to terminate spontaneously and responds better to antiarrhythmic drugs or cardioversion, whereas high-rate AF is more often persistent and refractory to therapy. Ventricular responses during AF can be characterized by a variety of methods, which include analysis of heart rate variability, RR-interval histograms, Lorenz plots, and non-linear dynamics. These methods have all shown a certain degree of usefulness, either in scientific explorations of atrioventricular (AV) nodal function or in selected clinical questions such as predicting response to drugs, cardioversion, or AV nodal modification. The role of the autonomic nervous system for AF sustenance and termination, as well as for ventricular rate responses, can be explored by different ECG analysis methods. In conclusion, non-invasive characterization of atrial fibrillatory activity and ventricular response can be performed from the surface ECG in AF patients. Different signal processing techniques have been suggested for identification of underlying AF pathomechanisms and prediction of therapy efficacy.

Key Words: Atrial fibrillation, ECG signal processing, Electrical remodelling, Antiarrhythmic drug monitoring, Cardioversion, Ablation


    Introduction
 Top
 Abstract
 Introduction
 Part I: Characterization of...
 Part II: Characterization of...
 Conclusions
 Acknowledgements
 References
 
In recent years, mechanisms leading to atrial fibrillation (AF) induction and maintenance have begun to be explored and rapidly evolving therapies including new antiarrhythmic drugs and catheter ablation techniques have been introduced into clinical practice. Nevertheless, outcome of AF patients is still unsatisfactory, which mainly can be attributed to the progressive nature of this arrhythmia and to only partially effective therapies. Moreover, AF treatment (e.g. rate vs. rhythm control, choice of antiarrhythmic drugs or device therapy, curative ablation) may be viewed as being based on trial and error, because no test is able to predict the natural history of this arrhythmia or its response to treatment. Consequently, current AF management guidelines1Go provide no treatment recommendations that ‘take the various mechanisms and patterns of AF into account’. Thus, it seems desirable to develop tests that quantify AF disease state and guide AF management.2Go,3Go

Virtually in every patient with AF, standard 12-lead surface electrocardiograms (ECG) and/or Holter recordings are acquired, the main purposes being confirmation of the presence of arrhythmia and the determination of ventricular rate. It is a common observation that fibrillatory waves have various appearances ranging from fine to coarse and from disorganized to organized and that the ventricular rate response varies in a rather unpredictable fashion. Interestingly, the mechanisms behind the fibrillatory process on the ECG as well as the ventricular response and, moreover, the possible prognostic information contained within the ECG have just recently begun to be explored.

The purpose of this article is (i) to describe technical aspects of novel ECG analysis techniques and (ii) to present research and clinical applications of these methods for characterization of both the fibrillatory process (Part I, Table 1) and the ventricular response (Part II) during AF.


View this table:
[in this window]
[in a new window]

 
Table 1 Methods for analysis of atrial activity

 

    Part I: Characterization of the fibrillatory process
 Top
 Abstract
 Introduction
 Part I: Characterization of...
 Part II: Characterization of...
 Conclusions
 Acknowledgements
 References
 
Background
Shortening of atrial refractoriness and maladaptation to rate are hallmarks of atrial electrical remodelling in AF.4Go During AF, re-excitation occurs during the repolarization phase of the preceding electrical wave, implying that local excitation almost always occurs without any obvious latency beyond the refractory period. Consequently, the average atrial fibrillatory rate is likely to reflect averaged refractoriness at any part of the tissue involved. This assumption has been independently verified in both animal5Go and human AF6Go,7Go by correlating both measures. Therefore, the length of the averaged atrial fibrillatory cycle (which is inversely related to fibrillatory frequency or rate) can be used as an index of the averaged atrial myocardial refractoriness and AF organization.8Go

The importance of shortened refractoriness and its transformation into short atrial cycles or high rates, as assessed by direct intra-atrial measurements for AF progression and response to therapy, has been well established. First, spontaneous arrhythmia behaviour seems to be related to baseline fibrillatory cycle length, with sustained AF having shorter mean cycle lengths compared with non-sustained AF.9Go,10Go Secondly, response to antiarrhythmic drug therapy has been shown to be associated with both baseline AF cycle length before drug administration and drug-induced cycle length changes. Both longer right atrial cycle length at baseline11Go and greater cycle length prolongation12Go seem to be associated with antiarrhythmic drug-induced AF termination. Finally, there exists a close relation between the degree of electrical abnormalities and AF recurrence following cardioversion, with shorter atrial effective refractoriness and abnormal rate adaptation being predictors of AF relapses.13Go–15Go

Prompted by the likely usefulness of a non-invasive test that assesses the average fibrillatory rate from the surface ECG for exploration of AF pathophysiology and AF management, frequency analysis techniques were independently introduced by our groups,16Go,17Go whereas two earlier studies proposed the use of this technique to subdivide AF and atrial flutter18Go or automatically identify AF among different rhythms.19Go

Techniques for characterization of the fibrillatory process
ECG acquisition
Using digital recording techniques, traditional ECG signals are recorded. In most studies, atrial fibrillatory rate has been obtained by spectral analysis techniques of resting ECG recordings such as standard 12-lead16Go,17Go or (modified) orthogonal recordings.20Go–22Go It has been suggested to analyse lead V1 when using a standard 12-lead ECG or, as proposed previously by Waktare et al.,23Go a ‘bipolar modification of V1’ covering the atria (low C1) when applying different ECG recording systems, as the largest fibrillatory waves are usually found in these leads.

The sampling rate of the ECG is commonly determined by the highest frequency of the activity of interest: in our case, the atrial activity whose highest frequency components are well below 50 Hz. Hence, it may be tempting to conclude that the ECG signal can be acquired at a sampling rate of 100 Hz. However, extraction of the atrial signal usually involves signal processing steps which, if to be performed successfully, require a considerably higher sampling rate, i.e. 250 Hz or higher. This observation applies, in particular, to the extraction of atrial activity by QRST subtraction (given subsequently), where a too low sampling rate introduces large residual errors in the atrial signal. Consequently, it may be necessary to increase the original sampling rate to 1 kHz using digital interpolation techniques for QRST subtraction methods.24Go

It is essential to correct for baseline wander, preferably using linear phase filtering,25Go since such noise has serious repercussions on the performance of the atrial signal extraction method. The presence of baseline wander may also influence the lower frequency components of the atrial signal spectrum, though to a lesser extent. Following correction of baseline wander and QRS detection, it is necessary to identify ectopic beats. Although recordings free of ectopy are preferred since they facilitate atrial signal extraction, it is usually possible to handle the presence of ectopic beats within the framework of the selected approach.

Extraction of atrial signals during AF
Since the characterization of the fibrillatory process primarily revolves around the determination of rate, a straightforward approach would be to perform spectral analysis on samples in intervals without ventricular activity, i.e. the T–Q intervals.18Go However, such an approach does not fully exploit the on-going nature of the fibrillatory process; samples during the QRST interval can be equally used for determining rate, provided that the ventricular activity has been first removed. The availability of additional samples reflecting atrial activity contributes not only to improved accuracy of the spectral estimate but also circumvents the problem of vanishing T–Q intervals at high heart rates.

The extraction of an atrial signal during AF requires non-linear signal processing since atrial and ventricular activities overlap spectrally and therefore cannot be separated by linear filtering. Average beat subtraction is the most widespread technique for atrial signal extraction and relies on the fact that AF is uncoupled to the ventricular activity. Hence, subtraction of the average QRST complex produces a residual signal, which is the atrial fibrillatory signal. This technique was initially developed for identifying P-waves during ventricular tachycardia,26Go but has later been employed for the AF analysis.16Go,17Go,19Go Since average beat subtraction is performed in individual leads, it becomes sensitive to alterations in electrical axis, which are manifest as large QRS-related residuals. However, the effect of such alterations can be suppressed by using spatiotemporal QRST cancellation in which the average beats of adjacent leads are mathematically combined with the average beat of the analysed lead in order to produce optimal cancellation (Figure 1).27Go


Figure 1131
View larger version (28K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 1 Time–frequency analysis of AF (left), validation by comparison of surface ECG lead V1 with intracardiac electrograms (top right), and short-term reproducibility (bottom right). Left: 2 s (out of a 60 s recording) of an ECG signal from a patient with AF (upper panel), and the same interval after QRST cancellation (middle panel, amplitude scale is magnified five times). This fibrillatory signal is then subjected to spectral analysis techniques (lower panel). Time–frequency distribution (left box), power frequency spectrum in which the dominant fibrillatory rate is determined (middle box), and frequency trend (right box).27Go,34Go,35Go Top right: agreement of fibrillatory rates from the right atrium (left), coronary sinus (middle), superior pulmonary veins (right) with simultaneously recorded fibrillatory rates from surface ECG lead V1 in patients with persistent AF (Bland–Altman method).38Go Bottom right: repeated determination of fibrillatory rate over 24 h in 20 patients with persistent AF.42Go

 
Another approach to atrial signal extraction exploits the property that atrial and ventricular activities arise from different bioelectrical sources. These two sources, exhibiting different statistical properties, can be separated when multi-lead ECG recordings are available. Principal component analysis (PCA)28Go and independent component analysis (ICA)29Go are the two main representatives of the separation approach, which both have found their way into the AF analysis.

Since both approaches have been shown to produce similar performance,30Go the choice of extraction technique may instead be guided by their decidedly different modes of processing the ECG signal. Spatiotemporal QRST cancellation extracts the atrial signal in a specific lead and strives at not modifying the original f-wave morphology. On the other hand, PCA and ICA derive a global atrial signal with contributions from all leads, although dominated by the lead with the largest atrial activity. A recording length of at least 10 s is required for adequate computation of the average beat, whereas the recording length can be shorter when PCA or ICA is used. In addition to these techniques, PCA was recently suggested for QRST cancellation in single leads to account for beat-to-beat changes in morphology.31Go This technique can be viewed as a generalization of average beat subtraction, where more than one component can be used to create a template that best fits each QRST complex.

Fourier-based spectral analysis
Following extraction of the atrial signal, a power spectrum is obtained, which typically exhibits a distinct peak whose location determines the most common fibrillatory rate of nearby endocardial sites.16Go,17Go Most studies make use of non-parametric, Fourier-based spectral analysis in which the ECG signal is divided into shorter, overlapping segments, each segment subjected to windowing (e.g. Welch's method).16Go,17Go The desired power spectrum is obtained by averaging the power spectra of the respective segments. It is customary to pad each segment with zeros so that the location of the spectral peak can be determined more precisely (although zero padding per se does not improve spectral resolution).

Segment length is the most critical design parameter in non-parametric spectral analysis since it determines the estimation accuracy of fibrillatory rate and restricts spectral resolution. It is advisable that segment length is chosen to be at least a few seconds to produce an acceptable variance of the rate estimate. Concerning spectral resolution, it may seem warranted to use longer segments in order to assure that concurrent, slightly different rates are appropriately resolved; for example, 10 s segments of the atrial signal yield a spectral resolution of ~0.1 Hz. However, no evidence has been put forward, which supports the existence of concurrent fibrillatory rates, but time–frequency analysis of the atrial signal (given subsequently) rather suggests that rate often varies substantially over time. Hence, spectral resolution has secondary importance when choosing the segment length.

Since the frequencies characterizing the atrial signal are mostly confined to the interval well below 50 Hz, a sampling rate much lower than the one suggested earlier for the extraction of the atrial signal may be used to reduce the amount of computations. Although sampling rate decimation is dispensable when non-parametric spectral analysis is performed, it is crucial to perform when parametric spectral estimation techniques are employed in order to avoid the presence of spurious spectral peaks.32Go

Time–frequency analysis
Power spectral analysis reflects the average signal behaviour during the analysed time interval—the location of the dominant spectral peak representing the main information carrier with clinical value. However, power spectral analysis suffers from an inability to characterize temporal variations in fibrillatory rate. The presence of bi- or multimodal spectral peaks, reported in some early studies,16Go,33Go may be interpreted as spatial variation in fibrillatory rate. When these results are seen in the light of recent studies employing time–frequency analysis,34Go,35Go such peaks are most likely explained by temporal variation in fibrillatory rate rather than by spatial variation.

Time–frequency analysis is a powerful tool for unveiling the temporal variation of the atrial signal, whether such variation is spontaneous or due to an intervention. In particular, time–frequency analysis can be used to quantify changes in fibrillatory rate, which may be associated with an intervention. Similar to power spectral analysis, a plethora of methods has been developed for time–frequency analysis, and it remains an open issue to select the method that provides the most adequate characterization of the atrial signal.

Recently, a novel approach to time–frequency analysis was presented, in which the time–frequency distribution, defined by local power spectra from successive 1 s segments, is decomposed into a spectral profile and a set of functions that describes the second-to-second variations in fibrillatory rate and f-wave morphology (Figure 1).34Go,35Go The spectral profile differs from the conventional power spectrum in that the local spectra have been frequency-aligned before spectral averaging. As a result, the peaks of the spectral profile become more distinct, implying that determination of the harmonic structure of the atrial signal, i.e. the f-wave morphology, is facilitated.

Method validation and reproducibility
A direct comparison between endocardially recorded electrograms and body surface recordings clearly shows the validity of fibrillatory rate obtained from surface ECG as an index of the average atrial fibrillatory cycle. Although in one preliminary report median right atrial and surface ECG rates exhibited only moderate correlations,36Go dominant fibrillatory rates calculated from lead V1 substituted the right atrial free wall,16Go,17Go,37Go and rates from an oesophageal lead reflected atrial septal and left atrial activity.16Go Interestingly, rate differences and variability between ECG and electrograms increase with growing anatomical distance (right atrium—coronary sinus—pulmonary veins) to V1. In a first preliminary report, which assessed the relation between pulmonary venous and ECG fibrillatory rates in 25 patients with persistent AF, correlations were moderate (R=0.457), whereas they were higher for right atrial (R=0.885) and coronary sinus (R=0.558) rates (P<0.05 for all correlations).38Go However, non-standard surface ECG lead positions for characterization of left atrial/pulmonary vein activity or identifying spatial variation have not been studied.

So far, direct comparisons of endocardial and surface signals have only been reported when the surface ECG had been subjected to QRST cancellation (Figure 1).16Go,17Go,36Go–38Go This technique has, in addition, been validated by comparing ‘pure’ AF during brief ventricular asystole with atrial signals obtained with QRST cancellation also showing good agreement of atrial rates.39Go In a comparative study, in which spatiotemporal QRST cancellation, PCA, and ICA algorithms have been applied to the same surface ECG recordings, the dominant rate of the atrial signals extracted by each of the algorithms was within 1 Hz (60 fibrillations per minute), with larger differences due to poorly defined dominant peaks and residual ventricular activity corrupting the frequency response.30Go

Although the gross atrial fibrillatory pattern, studied by multiple simultaneous epicardial recordings, is reproducible already in repeated measurements of only 8 s duration,40Go calculation of fibrillatory rates during steady-state conditions from even longer time intervals discloses a true variability. Thus, increasingly longer recording times of up to 30 min result in a successive decrease in this variability. Prolonging the recording time during steady-state conditions, allowing integration of atrial fibrillatory activity over longer periods, thus enhances the reproducibility of the method. For practical reasons, steady-state recordings may be restricted to 1–5 min, yielding rate variation coefficients of 2.1%.16Go In persistent AF, there is minor short-term rate variability within 30 min,16Go,17Go,39Go whereas repeated daily frequency measurements on identical medication at the same time under similar conditions disclose an insignificant fibrillatory rate variability.41Go Very recently, surface ECG atrial fibrillatory wave characteristics have been shown to be highly reproducible over 24 h in clinically stable AF patients even in 10 s ECG segments, as acquired by standard 12-lead ECGs (Figure 1).42Go

To date, the dominant spectral peak has been presented as dominating atrial cycle length (in ms), fibrillatory rate [in fibrillations per minute (fpm)], or fibrillatory frequency (in Hz). Since it has become evident that the results of this method are of practical clinical importance, we believe that the results should be expressed in a way which is closest to the nomenclature of other surface ECG rate variables (e.g. sinus, atrial tachycardia, or atrial flutter rate).43Go The expression ‘fibrillatory rate’, with its unit fpm is therefore appropriate in clinical reports.21Go,44Go

Moreover, the calculation of cycle length from the frequency power spectrum with its standard unit Hz (cycle length in ms=1000/frequency in Hz) seems somewhat troublesome.44Go Although this calculation is appropriate for single measurements, its application to multiple simultaneous, or repeated measurements such as pre- and post-drug states may introduce errors. The reason is that the same frequency difference results in larger cycle length differences for low frequencies compared with high frequencies.

Research and clinical applications
Exploring autonomic modulations
Vagal and sympathetic stimulations have been shown to reduce atrial refractory periods and increase their heterogeneity.45Go–47Go In the human atrium, beta stimulation has been found to predominate over vagal stimulation.46Go Moreover, invasive studies performed in subjects with sinus rhythm have suggested a circadian pattern in atrial refractoriness, with longer refractory periods during night time and refractory period shortening during daytime,48Go,49Go which supports the role of the autonomic nervous system in modulating atrial electrophysiological properties.

Consequently, atrial fibrillatory rate seems to be ideal for monitoring the effect of autonomic tone changes—either spontaneous (circadian) or by autonomic manoeuvres provoked—on atrial electrophysiology. Indeed, the circadian variability of atrial fibrillatory rate has been explored in independent studies.50Go,51Go Fibrillatory rate obtained from Holter ECGs with persistent AF showed a significant decrease at night and an increase during the morning hours.50Go,51Go In six of 30 individuals studied by our group,50Go dominant nocturnal fibrillatory rate increased, however, concomitantly with a decrease in ventricular rate, whereas the opposite change occurred during the morning hours.

The second area of investigation concerns effects on atrial fibrillatory rate following vagal or sympathetic stimulation during experimental conditions. Carotid sinus massage, mainly supposed to induce vagal stimulation, resulted in a variable effect on fibrillatory rate in 19 patients.52Go A reproducible decrease was noted in nine individuals, whereas a rate increase occurred in eight, and no change was observed in two patients. Interestingly, calcium-channel blocker treatment was the only variable affecting the rate response to carotid sinus massage. Calcium-channel blockers were more frequently used in patients with a decrease in fibrillatory rate compared with patients with a rate increase. The effect of adrenergic stimulation via head-up tilting on fibrillatory rate was studied in 14 patients with long-lasting AF.53Go In 12 patients, head-up tilt increased fibrillatory rate significantly, whereas there was no rate change in the remaining two (Figure 2).


Figure 1132
View larger version (18K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 2 Monitoring of interventions during AF, such as response to antiarrhythmic drugs (top) or to autonomic manoeuvers (bottom). Top: the solid line indicates the actual fibrillatory rate during oral drug loading with flecainide (300 mg bolus, followed by 100 mg bid) over 72 h. Bottom: typical example illustrating changes in heart rate (HR) and dominant atrial cycle length (DACL, which is inversely related to fibrillatory rate) during head-up tilt testing. An increased sympathetic discharge and withdrawal of vagal activity are known to be the physiological responses to the tilt table test, which is associated with increases in both HR and fibrillatory rate (decrease in DACL).53Go

 
With the availability of instantaneous fibrillatory rates obtained from time–frequency analysis, their second-to-second variation can be explored by spectral analysis techniques, similar to those used for analysing heart rate variability (HRV). By applying this approach, controlled respiration caused cyclical fluctuations in fibrillatory rate in patients with long-duration AF, which was related to parasympathetic modulations of the atrial refractory period.54Go,55Go

Predicting spontaneous AF behaviour and therapeutic effects
In contrast to persistent AF, rate variability in paroxysmal AF seems to be related to its natural course, with a rate increase within the first 5 min of an AF episode56Go and a rate decrease prior to termination.56Go,57Go In fact, fibrillatory rate has been shown to be highly predictive of spontaneous termination of paroxysmal AF. In a study of 30 ECG recordings, a low fibrillatory rate was able to identify patients with spontaneously terminating AF with high accuracy.58Go

Class I and III antiarrhythmic drugs have been shown to increase atrial cycle length (decrease fibrillatory rate), which coincides with increased refractoriness and decreased conduction velocity.59Go A substantial reduction in atrial fibrillatory rates following several different intravenously or orally administered class I and III antiarrhythmic drugs as well as following verapamil or magnesium has been observed when applying Fourier-based spectral analysis and more recently time–frequency analysis60Go,61Go (Figure 2, Table 2). Of special note is the capability of these methods to monitor drug effects serially or continuously, which appear superior to monitoring plasma drug levels.62Go


View this table:
[in this window]
[in a new window]

 
Table 2 Frequency analysis of AF for monitoring and predicting class I and III antiarrhythmic drug action

 
Interestingly, it is also possible to identify suitable candidates for pharmacological cardioversion. A baseline fibrillatory rate of 360 fpm or lower was highly sensitive and specific for prediction of AF termination following intravenous ibutilde17Go or oral flecainide.20Go In contrast, no baseline fibrillatory rate difference between patients who converted to sinus rhythm and those who did not following oral bepridil administration has been noted. Instead, larger rate increases were associated with AF termination,63Go which has also been found with ibutilide treatment.64Go Although this concept is interesting and well established in experimental AF,59Go it needs to be emphasized that the authors of the former study63Go calculated atrial cycle length from the frequency spectrum, which may have introduced the aforementioned statistical errors.44Go

Previous investigations65Go,66Go have shown that most AF relapses occur within the first weeks after cardioversion, with decreased but constant recurrence rates thereafter. Early vulnerability to AF re-initiation within this time period is related to electrophysiological abnormalities, whereas structural abnormalities seem to be primarily responsible for later AF recurrences.67Go Consequently, characterization of atrial electrophysiology has been suggested for identification of patients at risk for early AF recurrence.2Go

Previous studies reported higher rates in relapsed patients immediately prior to internal68Go or external cardioversion22Go,69Go when compared with non-relapsed patients and also close relationships between fibrillatory rate and defibrillation thresholds.68Go,70Go

So far, one study22Go has investigated the combined predictive value of fibrillatory rate and echocardiographic left atrial parameters for predicting AF recidivism. It could be demonstrated that the combination of fibrillatory rate and systolic left atrial area predicted early AF recurrence after successful cardioversion with high accuracy and was able to provide individual risk estimates. Very recently, fibrillatory rate changes in response to linear left atrial ablation have been monitored non-invasively and the effect on fibrillatory rate of roof and mitral isthmus lines quantified.71Go Even though there was a trend to lower baseline rates with successful ablation, this study was not powered to predict outcome, although an invasive study supports this conclusion.72Go

Taken together, fibrillatory rate seems well suited to describe the individual atrial remodelling and monitor interventions (Table 3). The concept of ECG-guided AF therapy may gain even greater importance in the light of two studies showing the same mortality risk for rate- and rhythm-control strategies if adequate anticoagulation is performed.73Go,74Go This highlights the need to select candidates for cardioversion and other interventional procedures not only based on clinical judgement but also on measures that are able to determine the likelihood of maintaining sinus rhythm, such as rate parameters obtained from the surface ECG.75Go


View this table:
[in this window]
[in a new window]

 
Table 3 Frequency analysis of AF for monitoring spontaneous and intervention-induced fibrillatory rate changes

 

    Part II: Characterization of ventricular response during AF
 Top
 Abstract
 Introduction
 Part I: Characterization of...
 Part II: Characterization of...
 Conclusions
 Acknowledgements
 References
 
Background
Given similar outcomes between rhythm and rate control strategies,73Go,74Go ventricular rate control seems an attractive endpoint for the treatment of AF.76Go AF is accompanied by an irregular ventricular response, mostly with a higher average rate than during sinus rhythm, provided atrioventricular (AV) conduction is patent and unaffected by any drug. It is the result of fast irregular atrial excitations reaching the AV node and its inherent conduction properties.

On the one hand, ease of ventricular rate measurement using ECG recordings with variable length (e.g. standard resting ECG vs. Holter recordings) and under variable circumstances (e.g. exercise ECG) is an obvious attribute. On the other hand, more advanced ECG analysis methods for exploration of autonomic modulation and AV nodal physiology as well as determining outcome of different treatments have been proposed. In particular, this article provides an overview of (i) HRV during AF and before its onset and (ii) analyses of RR-interval histograms, Lorenz plots, and non-linear dynamics (Table 4).


View this table:
[in this window]
[in a new window]

 
Table 4 Methods for analysis of ventricular response

 
Techniques and applications for characterization of ventricular response
HRV during AF
The exploration of autonomic modulation of sinus node impulse activity via HRV analysis is a widely used method. Analysis of ambulatory ECG recordings during sinus rhythm has been proved effective to identify high-risk patients in different cardiovascular diseases. However, during AF, ventricular variability does not purely reflect sinus node modulation but is also dependent on AV nodal refractoriness as well as degree of concealed conduction77Go and left atrial size.78Go The typical RR-interval irregularity of ventricular response during AF has been considered a factor that might prevent an appropriate application of HRV analysis in these patients. Indeed, although time domain parameters can be easily computed also in AF patients, frequency domain analysis has been rarely performed as a consequence of the difficulty in interpreting the huge number of spectral peaks of the variability signal. Spectral analysis of tachograms of AF patients is commonly used to represent a white noise signal. However, in one study published a few years ago, a respiration-related component in addition to a predominant very low-frequency component has been recognized by subjecting patients with chronic AF to sequential autonomic blockade.79Go

These considerations of complex determinants of ventricular rate during AF may explain why HRV analysis has mostly been applied to evaluate autonomic modulation of the sinus node during sinus rhythm before AF onset and after cardioversion rather than on ventricular response during the arrhythmia.

Signs of enhanced parasympathetic activity or of sympathetic activation have been reported by different authors80Go–85Go who analysed short-term HRV before paroxysmal AF onset. One common finding was the detection of transient alterations of neural modulatory activities before the arrhythmic event. Timing of arrhythmias onset, age, and clinical characteristics of subjects enrolled in the studies are some of the factors that have been advocated to explain the difference in autonomic parameters before AF onset. In more detail, Bettoni and Zimmerman84Go suggested that an early increase in adrenergic tone was followed by late parasympathetic activation buffering the adrenergic predominance. More recently, Lombardi et al.85Go reported that signs of a predominant sympathetic and reduced vagal modulation of the sinus node were detectable in most paroxysmal AF episodes that started during the daytime and were characterized by the presence of atrial ectopies, whereas a more physiological autonomic modulation of the sinus node was associated with spontaneous recovery of sinus rhythm. Vagal predominance was detectable in ~30% of cases mainly occurring during night-time.

The prognostic value of HRV for predicting arrhythmia recurrence after successful cardioversion has also been analysed. Signs of increased sympathetic and reduced vagal modulation (LF/HF ratio >2) at the time of discharge after successful cardioversion were associated with early but not late AF recurrence.66Go In contrast, persistently higher values of 24 h time domain parameters reflecting vagal modulation have been reported in patients with AF recurrences after internal cardioversion performed under general anaesthesia.86Go Differences in patient population, pharmacological treatment, recording conditions, and/or type of analysis may explain why either sympathetic activation or enhanced vagal modulation may be associated with AF recurrences. It seems, therefore, that similar to what has been described to occur before paroxysmal AF onset, transient alteration of autonomic modulation rather than a specific autonomic change seems to favour both arrhythmia onset and recurrence.

As to the changes in RR intervals during AF, it must be recalled that in the electrophysiology laboratory, slowing or pseudo-regularization of atrial activity often precedes spontaneous termination of induced AF episodes. To evaluate whether such a phenomenon could also occur in the clinical setting, Lombardi et al.85Go analysed RR intervals during the initial, central, and final part of paroxysmal AF episodes, with the assumption that a less disorganized atrial activity could be reflected by a pseudo-regularization of ventricular response. No changes in mean RR intervals, variance, and coefficient of variation of RR intervals were observed in the 5 min preceding spontaneous recovery of sinus rhythm in comparison with the initial and central parts of the arrhythmic episode.

However, using different techniques for expressing HRV during AF, this method has been found useful for providing information about risk in patients with AF and severe mitral regurgitation87Go as well as in the setting of advanced heart failure,77Go,88Go but may also reflect a general risk during AF.89Go As with reduced HRV in sinus rhythm, reduced ventricular response irregularity during AF has been found to be associated with increased cardiac mortality in chronic AF patients.89Go

RR-interval histogram analysis
When constructing RR-interval histograms obtained from Holter recordings with persistent AF, uni-, bi-, or multimodal RR distribution patterns can be found.90Go–93Go In up to 55% of patients, the ventricular response during AF shows two separate populations of RR intervals already at naked eye inspection of the ECG.94Go A special technique for displaying and analysing RR-interval histograms is the so-called heart rate stratified histogram analysis (HRSH). RR-interval histograms are reconstructed on the basis of mean heart rate of the analysed ECG segment that contains 2095–6490Go,91Go beats (Figure 3). Using this technique, a bimodal RR-interval distribution was found in 26 of 32 individuals (81%) with persistent AF.91Go Prerequisites for this demonstration are the clustering of RR-interval distributions by average ventricular rate and enough observations over a span of ventricular rates reaching at least 120 b.p.m.90Go The advantage of this technique over simple RR-interval histogram analysis is first of all that the RR distribution is equally well visualized at all heart rate levels and comparable between them. Secondly, the identification of multiple RR populations is markedly improved by this technique.


Figure 1133
View larger version (37K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 3 Illustration of the HRSH method (top) and its application for monitoring drug effects (bottom).95Go Top: each QRS complex is identified, and all RR intervals in the recording are calculated. The entire series of RR is then divided into sequences of 20 individual RR intervals, and one-step movements (i.e. 1–21, 2–22) create the sequences in which the average heart rate is calculated. All individual RR intervals within sequences are then pooled into histograms at different heart-rate levels, calculated from the average of 20 consecutive RR intervals. Thus, RR intervals with the same average heart rate are pooled into the same histogram. At lower heart-rate levels, the peak corresponding to the longer RR-interval population is dominant. As the heart rate increases, the peak of the shorter RR-interval population will increase and become dominant. The heart-rate level at which a change in peak dominance (PDC) occurs is used for further analysis (left). By applying the formula shown (left), the original histogram smoothes and the position of the two RR-interval populations [PV(s) and PV(f)] and derived variables (PG, PVR) can be estimated (right). Bottom: the two superimposed histograms obtained at PDC exhibit changes in both RR populations induced by MgSO4 infusion. However, conduction delay in the longer RR population was more pronounced than in the shorter RR population, indicating that the two RR populations were differently affected by the drug. PDC, peak dominance change; PV(s), peak value of the longer RR-interval population; PV(f), peak value of the shorter RR-interval population; PG, peak gap (difference between the two RR populations); PVR, peak value ratio (quotient of the longer and shorter RR-interval populations.

 
The two different RR populations are suggested to correspond to conduction along two different atrio-nodal conduction routes.90Go Interestingly, in individuals with AF, in whom the HRSH method was performed before and following surgical incision in the right atrium, this intervention constantly changed the pattern of RR distribution during AF, including the transformation from a unimodal to a bimodal pattern.96Go A tempting interpretation of this finding is that the timing between the fibrillatory signals reaching the different atrio-nodal inputs has been changed by the incision. In fact, later pacing and conduction studies in atrio-nodal tissue preparations speak in favour of this interpretation.97Go RR-interval distribution had a sensitivity of 88% and a specificity of 80% for identifying dual AV node physiology, provided dual pathways could be independently confirmed.92Go

Some investigators,94Go,98Go moreover, analysed the correlation between RR-interval distribution pattern and efficacy of radiofrequency modification of the AV node in patients with chronic AF. These authors observed that this procedure was more effective in controlling ventricular response in patients with bimodal distribution of RR-interval pattern. Interference with posterior atrio-nodal input is considered to be the prevailing mechanism of rate control in these patients.

Finally, the HRSH analysis has been applied to analyse the effects of magnesium on AV nodal conduction. In this study, magnesium at high concentrations caused a delay in both the shorter and longer RR intervals, with a more pronounced conduction delay in the longer RR population, indicating a different effect on the two corresponding AV nodal pathways (Figure 3).95Go

Rarely, RR-interval histograms disclose more than two different and well-defined RR populations,93Go which may be due to nodal escape rhythms91Go and/or to any additional atrio-nodal inputs, for example from the left atrium.99Go

Lorenz plots
In Lorenz plots, individual RR intervals are plotted against the corresponding preceding RR interval. In this way, it is possible to identify (i) the lower envelope of the Lorenz plot as a measure of the cycle-length-dependent function of the functional refractory period of the AV node100Go and (ii) the degree of scatter as a measure of concealed conduction (Figure 4).


Figure 1134
View larger version (40K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 4 Illustration of Lorenz plot analysis (top) and its application for monitoring circadian rhythms (bottom).77Go Top: Lorenz plots of 512 successive RR intervals in AF and linear regression of the lower envelope of the plots. Dotted vertical lines indicate the separations between bins, each of which includes 64 points of the plots. Circled points indicate the minimal values for individual bins. The solid line indicates the regression line on the minimal values. The following variables can be obtained: (i) slope of the regression line on the lower envelope (LE slope); (ii) 1.0 s intercept of the regression line on the lower envelope (LE1.0); and (iii) the degree of scatter above the lower envelope calculated as the root mean square difference of RR intervals from the regression line on the lower envelope (scattering index). Bottom: circadian variations in variables obtained from the Lorenz plot analysis, with a time resolution of 5 min in a patient with chronic AF without heart failure (left) and with heart failure (right). Thick solid lines indicate the least square cosine curves fitted to the circadian variation of individual variables.

 
By applying Lorenz plot analysis to 24 h Holter recordings of 48 patients with chronic AF, it was suggested that both AV node refractoriness and the degree of concealed AV conduction during AF may show a circadian rhythm, but that these circadian rhythms may be attenuated in patients with heart failure (Figure 4).77Go In a similar vein, sleep apnoea has very recently been shown to entrain the ventricular response to AF by inducing rhythmic oscillations in AV node refractoriness and the degree of concealed conduction.101Go These findings point to the possibility of obtaining information concerning altered autonomic control of the RR intervals in patients with AF (and heart failure or other disease) with this simple technique.

Lorenz plots of RR intervals from 24 h recordings sometimes exhibited two separate sectors of RR intervals.102Go When this occurred, the RR-interval histogram disclosed a bimodal distribution in ~40% of patients. It should be noted, however, that these RR-interval histograms were not stratified for different average heart rates, a method which, as pointed out before, markedly increases the identification of more than one RR population.91Go It has been suggested that Lorenz plots with two sectors hold information of the functional refractory periods of each of the two conduction routes.102Go Interestingly, the circadian variability of the fast pathway functional refractoriness was more pronounced than that of the slow pathway.

By adding the number of occurrences of RR-interval pairs, a three-dimensional Lorenz plot can be constructed, in which clusters of RR intervals can be identified. Moreover, in AF patients whose RR intervals were clustered (47% of 66 patients), electrical cardioversion was more effective in restoring sinus rhythm, and also, those patients had a higher likelihood of remaining in sinus rhythm than their counterparts without well-defined clusters.103Go It is speculated that RR-interval clustering represents a relatively high degree of organization of atrial fibrillatory activity.

Non-linear dynamics
New insights into the inner dynamics of ventricular response during AF have recently been obtained using non-linear analysis methods. Even if ventricular responses during AF have been widely described as a non-harmonic, totally random process,104Go glimpses of underlying order have been recently documented in both long-term (hours)105Go and short-term (few beats) dynamics of RR series.106Go These findings support pioneering and sporadic observations of non-randomness in animals107Go and humans108Go during AF.

The spectral characteristics of 24 h RR-interval series in patients with chronic AF revealed the presence of two correlation ranges105: a flat spectrum was observed at higher frequencies (>0.005 Hz), whereas the long-term spectral components (long-range correlation) were characterized by a 1/f power law behaviour. The long-range correlation was similar to that observed in normal sinus rhythm, suggesting that during AF, long-term regulatory mechanisms are still effective in modulating ventricular responses. These findings are in agreement with the observation of circadian rhythms in patients with chronic AF.109Go

Non-random features appear also when ventricular response dynamics are explored on a short-term basis. Stein et al.106Go developed a non-linear forecasting algorithm to search for short-term predictability in the ventricular responses during AF. The authors observed a weak predictability (lasting for a few beats), which was explained as the evidence of autonomic modulation of the AV node buried by a random sequence of impulses due to multiple circulating wavefronts. Patterns of predictability were also assessed using approximate entropy (ApEn) and regularity measurements in non-terminating and terminating AF episodes, with a decreased regularity index and increased ApEn in self-terminating episodes.110Go

Altered complexity and correlation properties of ventricular activation were investigated before the onset of AF episodes. A reduction of ApEn was observed prior the onset of AF, and the values were significantly lower than in healthy controls.111Go More recently, analysis of heart rate turbulence (i.e. the short-term oscillation of heart rate after a premature beat) revealed different patterns before AF onset. In particular, the turbulence onset was significantly less negative during the 1 h preceding AF than during non-AF hours.112Go

Taken together, the different methods used for exploring and quantifying the ventricular response during AF have all shown a certain degree of usefulness, either in scientific explorations of AV nodal function or in selected clinical questions. Several additional issues remain to be further clarified. Most importantly, the RR distribution pattern is not only a likely sign of conduction via different AV nodal routes, but primarily affected by the atrial fibrillatory input to the AV node. Hence, further exploration of how these factors may be individually mirrored within different methods is needed.


    Conclusions
 Top
 Abstract
 Introduction
 Part I: Characterization of...
 Part II: Characterization of...
 Conclusions
 Acknowledgements
 References
 
Non-invasive characterization of atrial fibrillatory activity and ventricular response can be performed from the surface ECG in AF patients using different signal processing techniques. ECG analysis of AF has been suggested for identification of underlying AF pathomechanisms and prediction of therapy efficacy (drug-induced conversion, maintenance of sinus rhythm, selecting antiarrhythmic drugs, identifying candidates for non-pharmacological AF therapy). However, further larger studies are necessary to determine the role of these techniques in different AF management strategies in order to select and time the appropriate therapy for the individual patient.


    Acknowledgements
 Top
 Abstract
 Introduction
 Part I: Characterization of...
 Part II: Characterization of...
 Conclusions
 Acknowledgements
 References
 
This study has been supported by the NordForsk network ‘Electrocardiology in Atrial Fibrillation’ through travel grants to Bollmann, Husser, and Sörnmo. Husser has been supported by the German Cardiac Society. Bollmann and Husser have received travel grants from the Polytechnic University of Valencia (VIDI-UPV). Bollmann and Millet have been supported by travel grants from the Generalitat Valenciana (AEST05 and AEST06/171). Rieta has received grant support from the Generalitat Valenciana (GV06/299).

Conflict of interest: A.B., M.S. and L.S. have received consulting fees from Sanofi-Aventis (Frankfurt, Germany) for evaluation of antiarrhythmic drug effects using time-frequency analysis.


    References
 Top
 Abstract
 Introduction
 Part I: Characterization of...
 Part II: Characterization of...
 Conclusions
 Acknowledgements
 References
 
[1] Fuster V, Rydén LE, Asinger RW, et al. ACC/AHA/ESC guidelines for the management of patients with atrial fibrillation: executive summary. A Report of the American College of Cardiology/American Heart Association Task Force on Practice Guidelines and the European Society of Cardiology Committee for Practice Guidelines and Policy Conferences (Committee to Develop Guidelines for the Management of Patients With Atrial Fibrillation) Developed in Collaboration With the North American Society of Pacing and Electrophysiology. Circulation 2001; 104: 2118–50.[Free Full Text]

[2] Bollmann A. Quantification of electrical remodeling in human atrial fibrillation. Cardiovasc Res 2000; 47: 207–9.[Free Full Text]

[3] Bollmann A. First comes diagnosis then comes treatment: an underappreciated paradigm in atrial fibrillation management. Eur Heart J 2005; 26: 2487–9.[Free Full Text]

[4] Wijffels MC, Kirchhof CJ, Dorland R, Allessie MA. Atrial fibrillation begets atrial fibrillation. A study in awake chronically instrumented goats. Circulation 1995; 92: 1954–68.[Abstract/Free Full Text]

[5] Kim KB, Rodefeld MD, Schuessler RB, Cox JL, Boineau JP. Relationship between local atrial fibrillation interval and refractory period in the isolated canine atrium. Circulation 1996; 94: 2961–7.[Abstract/Free Full Text]

[6] Ramdat Misier AR, Opthof T, van Hemel NM, et al. Increased dispersion of refractoriness in patients with idiopathic paroxysmal atrial fibrillation. J Am Coll Cardiol 1992; 19: 1531–5.[Abstract]

[7] Capucci A, Biffi M, Boriani G, et al. Dynamic electrophysiological behavior of human atria during paroxysmal atrial fibrillation. Circulation 1995; 92: 1193–202.[Abstract/Free Full Text]

[8] Konings KT, Kirchhof CJ, Smeets JR, Wellens HJ, Penn OC, Allessie MA. High–density mapping of electrically induced atrial fibrillation in humans. Circulation 1994; 89: 1665–80.[Abstract/Free Full Text]

[9] Asano Y, Saito J, Matsumoto K, Kaneko K, Yamamoto T, Uchida M. On the mechanism of termination and perpetuation of atrial fibrillation. Am J Cardiol 1992; 69: 1033–8.[CrossRef][ISI][Medline]

[10] Boahene KA, Klein GJ, Yee R, Sharma AD, Fujimura O. Termination of acute atrial fibrillation in the Wolff–Parkinson–White syndrome by procainamide and propafenone: importance of atrial fibrillatory cycle length. J Am Coll Cardiol 1990; 16: 1408–14.[Abstract]

[11] Stambler BS, Wood MA, Ellenbogen KA. Antiarrhythmic actions of intravenous ibutilide compared with procainamide during human atrial flutter and fibrillation: electrophysiological determinants of enhanced conversion efficacy. Circulation 1997; 96: 4298–306.[Abstract/Free Full Text]

[12] Fujiki A, Nagasawa H, Sakabe M, et al. Spectral characteristics of human atrial fibrillation waves of the right atrial free wall with respect to the duration of atrial fibrillation and effect of class I antiarrhythmic drugs. Jpn Circ J 2001; 65: 1047–51.[CrossRef][Medline]

[13] Attuel P, Childers R, Cauchemez B, Poveda J, Mugica J, Coumel P. Failure in the rate adaptation of the atrial refractory period: its relationship to vulnerability. Int J Cardiol 1982; 2: 179–97.[CrossRef][ISI][Medline]

[14] Manios EG, Kanoupakis EM, Chlouverakis GI, Kaleboubas MD, Mavrakis HE, Vardas PE. Changes in atrial electrical properties following cardioversion of chronic atrial fibrillation: relation with recurrence. Cardiovasc Res 2000; 47: 244–53.[Abstract/Free Full Text]

[15] Biffi M, Boriani G, Bartolotti M, Bacchi Reggiani L, Zannoli R, Branzi A. Atrial fibrillation recurrence after internal cardioversion: prognostic importance of electrophysiological parameters. Heart 2002; 87: 443–8.[Abstract/Free Full Text]

[16] Holm M, Pehrson S, Ingemansson M, et al. Non-invasive assessment of the atrial cycle length during atrial fibrillation in man: introducing, validating and illustrating a new ECG method. Cardiovasc Res 1998; 38: 69–81.[Abstract/Free Full Text]

[17] Bollmann A, Kanuru NK, McTeague KK, Walter PF, DeLurgio DB, Langberg J.J. Frequency analysis of human atrial fibrillation using the surface electrocardiogram and its response to ibutilide. Am J Cardiol 1998; 81: 1439–45.[CrossRef][ISI][Medline]

[18] Rosenbaum DS and Cohen RJ. Frequency based measures of atrial fibrillation in man. IEEE-EMBS 1990; 12: 582–3.

[19] Slocum J, Sahakian A, Swiryn S. Diagnosis of atrial fibrillation from surface electrocardiograms based on computer-detected atrial activity. J Electrocardiol 1992; 25: 1–8.[ISI][Medline]

[20] Bollmann A, Binias KH, Toepffer I, Molling J, Geller C, Klein HU. Importance of left atrial diameter and atrial fibrillatory frequency for conversion of persistent atrial fibrillation with oral flecainide. Am J Cardiol 2002; 90: 1011–4.[CrossRef][ISI][Medline]

[21] Bollmann A, Husser D, Stridh M, et al. Frequency measures obtained from the surface electrocardiogram in atrial fibrillation research and clinical decision-making. J Cardiovasc Electrophysiol 2003; 14: S154–61.[CrossRef][ISI][Medline]

[22] Bollmann A, Husser D, Steinert R, et al. Echo- and electrocardiographic predictors for atrial fibrillation recurrence following cardioversion. J Cardiovasc Electrophysiol 2003; 14: S162–5.[CrossRef][ISI][Medline]

[23] Waktare JE, Gallagher MM, Murtagh AM, Camm AJ, Malik M. Optimum lead positioning for recording bipolar atrial electrocardiograms during sinus rhythm and atrial fibrillation. Clin Cardiol 1998; 21: 825–30.[ISI][Medline]

[24] Laguna P and Sörnmo L. Sampling rate and the estimation of ensemble variability for repetitive signals. Med Biol Eng Comput 2000; 38: 540–6.[CrossRef][ISI][Medline]

[25] Sörnmo L and Laguna P. Bioelectrical Signal Processing in Cardiac and Neurological Applications. 2005; Amsterdam Elsevier.