ATRIAL FIBRILLATION
Analysis of surface electrocardiograms in atrial fibrillation: techniques, research, and clinical applications
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 |
|---|
|
|
|---|
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 |
|---|
|
|
|---|
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 guidelines1
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.
|
| Part I: Characterization of the fibrillatory process |
|---|
|
|
|---|
Background
Shortening of atrial refractoriness and maladaptation to rate are hallmarks of atrial electrical remodelling in AF.4
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.9
,10
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 baseline11
and greater cycle length prolongation12
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.13
15
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,16
,17
whereas two earlier studies proposed the use of this technique to subdivide AF and atrial flutter18
or automatically identify AF among different rhythms.19
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-lead16
,17
or (modified) orthogonal recordings.20
22
It has been suggested to analyse lead V1 when using a standard 12-lead ECG or, as proposed previously by Waktare et al.,23
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.24
It is essential to correct for baseline wander, preferably using linear phase filtering,25
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 TQ intervals.18
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 TQ 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,26
but has later been employed for the AF analysis.16
,17
,19
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).27
|
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)28
Since both approaches have been shown to produce similar performance,30
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.31
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.16
,17
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).16
,17
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 timefrequency 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.32
Timefrequency analysis
Power spectral analysis reflects the average signal behaviour during the analysed time intervalthe 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,16
,33
may be interpreted as spatial variation in fibrillatory rate. When these results are seen in the light of recent studies employing timefrequency analysis,34
,35
such peaks are most likely explained by temporal variation in fibrillatory rate rather than by spatial variation.
Timefrequency 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, timefrequency 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 timefrequency 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 timefrequency analysis was presented, in which the timefrequency 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).34
,35
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,36
dominant fibrillatory rates calculated from lead V1 substituted the right atrial free wall,16
,17
,37
and rates from an oesophageal lead reflected atrial septal and left atrial activity.16
Interestingly, rate differences and variability between ECG and electrograms increase with growing anatomical distance (right atriumcoronary sinuspulmonary 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).38
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).16
,17
,36
38
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.39
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.30
Although the gross atrial fibrillatory pattern, studied by multiple simultaneous epicardial recordings, is reproducible already in repeated measurements of only 8 s duration,40
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 15 min, yielding rate variation coefficients of 2.1%.16
In persistent AF, there is minor short-term rate variability within 30 min,16
,17
,39
whereas repeated daily frequency measurements on identical medication at the same time under similar conditions disclose an insignificant fibrillatory rate variability.41
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).42
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).43
The expression fibrillatory rate, with its unit fpm is therefore appropriate in clinical reports.21
,44
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.44
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.45
47
In the human atrium, beta stimulation has been found to predominate over vagal stimulation.46
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,48
,49
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 changeseither spontaneous (circadian) or by autonomic manoeuvres provokedon atrial electrophysiology. Indeed, the circadian variability of atrial fibrillatory rate has been explored in independent studies.50
,51
Fibrillatory rate obtained from Holter ECGs with persistent AF showed a significant decrease at night and an increase during the morning hours.50
,51
In six of 30 individuals studied by our group,50
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.52
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.53
In 12 patients, head-up tilt increased fibrillatory rate significantly, whereas there was no rate change in the remaining two (Figure 2).
|
With the availability of instantaneous fibrillatory rates obtained from timefrequency 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.54
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 episode56
and a rate decrease prior to termination.56
,57
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.58
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.59
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 timefrequency analysis60
,61
(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.62
|
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 ibutilde17
Previous investigations65
,66
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.67
Consequently, characterization of atrial electrophysiology has been suggested for identification of patients at risk for early AF recurrence.2
Previous studies reported higher rates in relapsed patients immediately prior to internal68
or external cardioversion22
,69
when compared with non-relapsed patients and also close relationships between fibrillatory rate and defibrillation thresholds.68
,70
So far, one study22
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.71
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.72
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.73
,74
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.75
|
| Part II: Characterization of ventricular response during AF |
|---|
|
|
|---|
Background
Given similar outcomes between rhythm and rate control strategies,73
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).
|
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 conduction77
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 authors80
85
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 Zimmerman84
suggested that an early increase in adrenergic tone was followed by late parasympathetic activation buffering the adrenergic predominance. More recently, Lombardi et al.85
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.66
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.86
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.85
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 regurgitation87
as well as in the setting of advanced heart failure,77
,88
but may also reflect a general risk during AF.89
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.89
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.90
93
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.94
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 20956490
,91
beats (Figure 3). Using this technique, a bimodal RR-interval distribution was found in 26 of 32 individuals (81%) with persistent AF.91
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.90
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.
|
The two different RR populations are suggested to correspond to conduction along two different atrio-nodal conduction routes.90
Some investigators,94
,98
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).95
Rarely, RR-interval histograms disclose more than two different and well-defined RR populations,93
which may be due to nodal escape rhythms91
and/or to any additional atrio-nodal inputs, for example from the left atrium.99
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 node100
and (ii) the degree of scatter as a measure of concealed conduction (Figure 4).
|
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).77
Lorenz plots of RR intervals from 24 h recordings sometimes exhibited two separate sectors of RR intervals.102
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.91
It has been suggested that Lorenz plots with two sectors hold information of the functional refractory periods of each of the two conduction routes.102
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.103
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,104
glimpses of underlying order have been recently documented in both long-term (hours)105
and short-term (few beats) dynamics of RR series.106
These findings support pioneering and sporadic observations of non-randomness in animals107
and humans108
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.109
Non-random features appear also when ventricular response dynamics are explored on a short-term basis. Stein et al.106
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.110
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.111
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.112
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 |
|---|
|
|
|---|
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 |
|---|
|
|
|---|
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 |
|---|
|
|
|---|
[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: 211850.
[2] Bollmann A. Quantification of electrical remodeling in human atrial fibrillation. Cardiovasc Res 2000; 47: 2079.
[3] Bollmann A. First comes diagnosis then comes treatment: an underappreciated paradigm in atrial fibrillation management. Eur Heart J 2005; 26: 24879.
[4] Wijffels MC, Kirchhof CJ, Dorland R, Allessie MA. Atrial fibrillation begets atrial fibrillation. A study in awake chronically instrumented goats. Circulation 1995; 92: 195468.
[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: 29617.
[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: 15315.[Abstract]
[7] Capucci A, Biffi M, Boriani G, et al. Dynamic electrophysiological behavior of human atria during paroxysmal atrial fibrillation. Circulation 1995; 92: 1193202.
[8] Konings KT, Kirchhof CJ, Smeets JR, Wellens HJ, Penn OC, Allessie MA. Highdensity mapping of electrically induced atrial fibrillation in humans. Circulation 1994; 89: 166580.
[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: 10338.[CrossRef][ISI][Medline]
[10] Boahene KA, Klein GJ, Yee R, Sharma AD, Fujimura O. Termination of acute atrial fibrillation in the WolffParkinsonWhite syndrome by procainamide and propafenone: importance of atrial fibrillatory cycle length. J Am Coll Cardiol 1990; 16: 140814.[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: 4298306.
[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: 104751.[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: 17997.[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: 24453.
[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: 4438.
[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: 6981.
[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: 143945.[CrossRef][ISI][Medline]
[18] Rosenbaum DS and Cohen RJ. Frequency based measures of atrial fibrillation in man. IEEE-EMBS 1990; 12: 5823.
[19] Slocum J, Sahakian A, Swiryn S. Diagnosis of atrial fibrillation from surface electrocardiograms based on computer-detected atrial activity. J Electrocardiol 1992; 25: 18.[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: 10114.[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: S15461.[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: S1625.[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: 82530.[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: 5406.[CrossRef][ISI][Medline]
[25] Sörnmo L and Laguna P. Bioelectrical Signal Processing in Cardiac and Neurological Applications. 2005; Amsterdam Elsevier.



