Feasibility of multiple short, 40-s, intra-procedural ECG recordings to detect immediate changes in heart rate variability during catheter ablation for arrhythmias.

PURPOSE: This study aims to evaluate a method to detect heart rate variability (HRV) changes using short ECG segments during ablation for arrhythmias.
METHODS: HRV was averaged from sequentially shorter time windows from 5-min ECG recordings in 15 healthy volunteers. The 40-s window was identified as the shortest duration that yielded reproducible values in high frequency (HF) and low frequency (LF) HRV. This method was validated in patients undergoing tilt table testing to see if the expected modulation in HRV that occurs prior to syncope could be detected from multiple 40-s recordings. Lastly, this method was used to assess HRV changes in 75 patients undergoing ablation for atrial fibrillation (AF) and other arrhythmias, to see if autonomic modulation as a result of ablation could be detected. A further 14 patients had stepwise HRV measurements at different stages of the AF ablation procedure to determine whether intra-procedural HRV changes could be detected.
RESULTS: HRV, averaged from multiple 40-s recordings, demonstrated the expected increase immediately preceding syncope compared with baseline (LF: 341 ± 311-1,536 ± 1,368 ms(2), p< 0.05; HF: 342 ± 339-1,628 ± 1,755 ms(2), p < 0.05). AF ablation, particularly following right pulmonary vein circumferential ablation, produced immediately detectable reductions in LF (153 ± 251-50 + 116 ms(2), p < 0.001) and HF (86 ± 195-33 ± 83 ms(2), p < 0.001) without any change in RR interval (877 ± 191-843 ± 220 ms, p = 0.261). Ablation for atrial flutter did not change the mean RR interval, LF or HF HRV.
CONCLUSION: Averaging multiple 40-s windows give valid HF and LF HRV measurements that enable detection of intra-procedural changes. Left atrial ablation around the right-sided pulmonary veins is unique in producing reductions in HRV. This method has the potential for use as an endpoint marker for adjunctive autonomic ablation procedures.

Introduction

The intrinsic cardiac autonomic nervous system has been implicated in initiation and maintenance of atrial fibrillation (AF) by animal studies showing that autonomic stimulation leads to pulmonary vein ectopy [1] and increases susceptibility of atria to develop AF in the presence of ectopic activity [2].

In humans, modification of the autonomic nervous system has been suggested for prevention of AF [3, 4] by targeting presumed sites of dense atrial autonomic innervation [5]. Currently, there are no established techniques to monitor autonomic modification during ablation, so adjunctive autonomic ablation is performed either by empirical ablation at predetermined anatomical sites [6] or by targeting sites that cause atrioventricular block (AV) block [3, 4, 7, 8]. These are not ideal end-points because although they imply that autonomic ablation may have taken place, it is not possible to ascertain the extent of this autonomic modification in the context of the patient’s overall cardiac autonomic function.

Heart rate variability (HRV), as an index of cardiac autonomic modulation, might have potential for such monitoring. Oscillations in the interval between consecutive heart beats can be measured on fairly short recordings [9] by applying spectral analysis to the RR interval series, yielding HF (the “high frequency” oscillations of respiratory sinus arrhythmia) and LF (the lower frequency intrinsic cardiac oscillations with a cycle time of approximately 10 s) components of HRV. Autonomic ablation is known to attenuate LF and HF HRV. However, HRV has previously only been measured from recording times of many minutes or hours [8, 10, 11], preventing real-time use for intra-procedural monitoring during catheter ablation for cardiac arrhythmias. The catheter laboratory is a more complex recording environment than an ambulant outpatient setting, with frequent ectopy caused by intracardiac catheters, which necessitates a technique of HRV assessments from ECG segments short enough to be ectopy-free for adequate HRV analysis, without unduly prolonging procedure duration.

In this study, we hypothesised that the effects of ablation procedures on HRV can be detected using multiple ECG segments short enough to permit real-time monitoring of autonomic modulation.

To test this, we first estimated the shortest ECG recording duration that still yielded reproducible average LF and HF HRV values in healthy volunteers. We then validated this technique by averaging HRV using multiple short-duration recordings in patients undergoing tilt table testing to see if the expected modulation in the autonomic nervous system that occurs prior to syncope could be detected using this method. Lastly, we used this method to assess HRV in patients undergoing ablation for AF and other arrhythmias or negative electrophysiological studies, to see if autonomic modulation as a result of ablation could be detected intra-procedurally.

Methods

Validation of short recordings

Fifteen healthy normal volunteers, who were not taking any drugs, were studied with ECG recordings (Hewlett-Packard 78351A) acquired from each subject for 20 min while lying still in a quiet room. The ECG was acquired from the limb lead with the largest R wave (usually lead II). Data was sampled at 1,000 Hz on a computer using an analogue-to-digital converter (National Instruments) and stored for off-line analysis of HRV using Matlab (Natick, MA, USA) [12]. RR intervals were determined and then converted into an evenly spaced temporal sequence at 5 Hz. Having identified ectopic-free segments, autoregressive spectral analysis was then carried out on the complete segment [12]. Low frequency (LF: 0.04–0.15 Hz) and high frequency (HF: 0.15–0.4 Hz) spectral power components were then calculated from the autoregression profiles by integration [12].

A 5-min (300 s) recording segment free of ectopy was selected for LF and HF HRV analysis. Further sampling of progressively shorter segments (2 × 150, 4 × 75, 6 × 50, 8 × 40, 10 × 30, and 15 × 20 s segments) of data within the 5-min recording was selected for HRV analysis to ascertain the shortest recording window that yielded reproducible average LF and HF values. The HRV indices for each of the recording durations were compared by repeated measures analysis of variance (ANOVA) with Newman–Keuls correction for comparison between the means. The HRV indices were further tested in random effects ANOVA and variance components models and intra-class correlation coefficients were calculated.

Further validation of the technique was conducted on ten patients undergoing tilt table testing for vasovagal syncope. LF and HF HRV were averaged from multiple 40-s ECG recordings within a 5-min window at the beginning of the tilt, and immediately preceding syncope to determine if the known changes in autonomic modulation that occur before syncope [13] were detected adequately using this technique.

Heart rate variability from short recordings during electrophysiological study and ablation

We studied 75 patients with symptomatic palpitations referred to our institution for electrophysiological (EP) study and ablation to see if the expected HRV changes could be identified following ablation. All patients gave full written informed consent for the procedure in-line with local institutional policy. There were 15 patients in each ablation arm, undergoing ablation for paroxysmal AF, atrial flutter, AV nodal re-entrant tachycardia (AVNRT) and AV re-entrant tachycardia (AVRT; Table 1). Fifteen patients had a negative EP study without ablation and served as controls. None of the patients received general anaesthesia for their procedures; however, patients were given sedation (15 of the AF patients and six of the non-AF patients) and analgesia (midazolam, morphine) as required if they were restless or in discomfort. Isoproterenol was given to provoke tachycardia if this was felt necessary. All patients were in sinus rhythm before and after their procedures.

Table 1

Patient characteristics and ablation procedure

	Age (years)	Sex (M/F)	Duration (min)	Procedure
AF	59 ± 11	6/9	149 ± 46	Segmental pulmonary vein isolation
AVNRT	48 ± 15	3/12	106 ± 37	Slow pathway modification
AVRT	34 ± 16	10/5	153 ± 72	Accessory pathway ablation
Atrial Flutter	64 ± 9	8/7	105 ± 37	Cavo-tricuspid ablation
Negative EP Study	43 ± 17	4/11	94 ± 34	No ablation

Radiofrequency energy was used for ablation in all patients with procedural details described in Table 1. Conventional end-points of ablation were used: wide area circumferential pulmonary vein isolation for atrial fibrillation (without any linear lesions or targeting of complex fractionated electrograms); bidirectional block across the cavo-tricuspid isthmus for atrial flutter; evidence of slow pathway modification for AVNRT and loss of pre-excitation or non-inducibility for AVRT. There was no deliberate attempt to modify or extend the ablation lesions to target areas of dense autonomic innervation [5].

Procedures were defined as beginning at needle to skin and ending at catheter removal. All recordings were taken in the period between these events (i.e. intra-procedural). Sinus rhythm ECG recordings were taken for at least 5 min at the beginning and at the end of each procedure, with patients lying still on the operating table. The ECG recordings were stored on a Bard system (Bard EP, Lowell, MA, USA) for offline analysis. The ECG with the largest R wave amplitude was visually inspected to exclude ectopic beats before exporting segments for analysis.

Using the shortest duration that yielded valid HRV results from the study in normal volunteers, at least five 40-s ECG segments were exported at 1,000 Hz sampling frequency for HRV analysis to obtain average LF and HF values immediately before and after the procedure to ascertain if the expected changes in autonomic modulation was observed.

To prevent confounding by isoproterenol, we separately grouped patients who received isoproterenol during the procedure and those who did not. We performed multivariate analysis to establish whether patient age and procedure duration could contribute to the changes observed in patients who did not receive isoproterenol.

In a separate group of 14 patients undergoing ablation for AF, the LF and HF powers were calculated in a step-wise approach to see if ablation at any particular site was responsible for HRV changes observed. HRV indices were therefore measured at the following steps: baseline, after transseptal puncture, after right pulmonary vein circumferential ablation and after left pulmonary vein circumferential ablation. Seven patients underwent left-sided ablation before right-sided ablation and seven patients had right-sided ablation before left-sided ablation.

Statistical analysis

Data are expressed as mean ± standard deviation. The power of each frequency band was logarithmically transformed to normalise the distribution [8]. Pre- and post-ablation LF and HF HRV values were then compared using the two-tailed paired Student’s t test. A value of p < 0.05 indicates statistical significance. The potential confounding effects of age and procedure duration on changes in HRV were analysed in multivariate linear regression analyses. All statistical tests were performed using Stata v10.1 (StataCorp LP, TX, USA).

Results

Short duration recordings used for heart rate variability analysis in healthy volunteers

Fifteen healthy volunteers with a mean age of 43 ± 12 years were studied under control conditions with continuous ECG monitoring for 20 min. A 5-min ectopy-free recording segment was chosen and analysed in progressively shorter durations from 5 min down to 20 s.

The HF power, representing respiratory sinus arrhythmia, which has a cycle time typically of 3–4 s, was shown to have no consistent shift in its mean value as the analysis window was shortened, although as can be seen from the graph (Fig. 1a), there was increased scatter when the recording duration was 30 or 20 s.

Fig. 1

Validation of short recording duration for HRV calculation. Log HF (a) and log LF (b) values plotted against duration of analysis in 15 healthy volunteers. The HF power showed no shift in mean value in all recording durations, although scatter increased at 30 and 20 s recording durations. The LF power showed no shift in mean value down to 40 s. However, at window durations shorter than 40 s, there was a significant difference (p < 0.05), and increased scatter in the mean LF values. Mean HRV indices for each recording duration compared by one-way ANOVA with Newman–Keuls correction for comparison between the means. HF high frequency, LF low frequency, HRV heart rate variability

The LF power, representing oscillations whose typical cycle length is of the order of 10 s, showed no significant variability as the analysis window duration was shortened down to 40 s. However, at window durations shorter than 40 s, there was a significant difference and increased scatter in the mean LF values (Fig. 1b).

The intra-class correlations coefficients, calculated following variance components analyses were 0.80 for HF and 0.93 for LF, confirming that in these indices measured for different durations, within-subject variance was small compared with between-subject variance.

Short duration recordings used for heart rate variability analysis in patients undergoing tilt table testing

To validate this method, LF and HF HRV was calculated at baseline and immediately before onset of syncope using multiple 40 s windows in ten patients (mean age, 48 ± 22) undergoing tilt table testing for vasovagal syncope. The LF and HF HRV, averaged from 40-s windows, were found to increase significantly immediately preceding syncope (LF: 341 ± 311–1,536 ± 1,368 ms2, p < 0.01; HF: 342 ± 339–1,628 ± 1,755 ms2, p < 0.05). HRV assessed from conventional 5-min windows behaved in the same way (LF: 251 ± 242–735 ± 440 ms2, p < 0.01, HF: 297 ± 318–778 ± 739 ms2, p < 0.05).

Ablation procedures and success

Radiofrequency ablation was deemed a primary success in all 60 patients undergoing ablation. None of the patients experienced any acute complication during the procedure.

Acute change in heart rate variability analysed using multiple 40-s windows

The AF ablation procedure (which in no case involved isoproterenol use) significantly reduced LF by 67% and HF by 61% (example shown in left panel of Fig. 3). There was no change in the mean RR interval (877 ± 191–843 ± 220 ms, p = 0.26) or the LF/HF ratio.

Fig. 3

Power spectral density and RR tachogram plots pre- and post-ablation power spectral density and RR tachogram plots for example patients in AF and atrial flutter

The other EP procedures, in which no isoproterenol was used (15 atrial flutter, five AVNRT, nine AVRT, and five negative EP studies) had no effect on the mean RR interval, LF, HF or LF/HF values (Table 2, Fig. 2). An example of LF and HF values following ablation for atrial flutter is shown in the right panel of (Fig. 3).

Table 2

Immediate effects of electrophysiological procedures on heart rate variability

	Number of Subjects	RR interval (ms)			LF (ms2)			HF (ms2)			LF/HF ratio
		Pre-ablation	Post-ablation	P value	Pre-ablalion	Post-ablation	P value	Pre-ablation	Post-ablation	P value	Pre-ablation	Post-ablation	P value
No isoprenaline given
PAF	15	877 ± 191	843 ± 220	NS	153 ± 251	50 ± 116	<0.001	86 ± 195	33 ± 83	<0.001	1.8 ± 2.3	15 ± 3.0	NS
Atrial flutter	15	802 ± 187	809 ± 145	NS	188 ± 450	86 ± 120	NS	96 ± 184	47 ± 100	NS	2.0 ± 2.4	1.8 ± 6.6	p < 0.05
AVNRT	5	841 ± 74	795 ± 103	NS	205 ± 108	257 ± 122	NS	71 ± 37	65 ± 39	NS	2.9 ± 1.3	4.0 ± 2.8	NS
AVRT	9	723 ± 121	719 ± 115	NS	272 ± 186	248 ± 205	NS	82 ± 65	102 ± 77	NS	3.3 ± 2.9	2.4 ± 2.7	NS
Negative EP study	5	827 ± 240	728 ± 193	NS	49 ± 39	76 ± 77	NS	38 ± 63	27 ± 32	NS	1.3 ± 4.7	2.8 ± 4.7	NS
Isoprenaline given
AVNRT	10	773 ± 171	604 ± 126	0.02	189 ± 147	106 ± 103	0.03	88 ± 47	30 ± 38	<0.01	2.2 ± 1.2	3.5 ± 5.6	NS
AVRT	6	840 ± 97	660 ± 75	<0.001	657 ± 567	344 ± 289	NS	250 ± 123	62 ± 47	0.03	2.6 ± 3.3	5.6 ± 6.4	<0.05
Negative EP Study	10	841 ± 202	700 ± 128	0.01	475 ± 377	204 ± 187	<0.05	258 ± 232	126 ± 184	<0.05	1.8 ± 1.5	1.6 ± 3.8	NS

NS indicates not significant

Fig. 2

Changes in RR interval and HRV indices after ablation. Non-AF includes patients having ablation for atrial flutter, atrioventricular nodal re-entry tachycardia, atrioventricular re-entrant tachycardia and negative electrophysiological studies. NS nonsignfiiciant

The remaining procedures, which were performed with isoproterenol administration (ten AVNRT, six AVRT and ten negative EP studies), showed a reduction in the mean RR interval, accompanied by the expected reductions in LF and HF HRV and an increase in the LF/HF ratio (Table 2, Fig. 2) detected using the multiple 40-s window technique.

Intra-procedural changes in heart rate variability assessed in a step-wise approach

Fourteen patients with a mean age of 59 ± 8 years and a mean AF duration of 56 ± 30 months were studied for step-wise changes in HRV. In both groups, a significant reduction in both LF and HF HRV were only observed following right pulmonary vein, but not left pulmonary vein wide area circumferential ablation (Fig. 4).

Fig. 4

Changes in HRV indices measured in a step-wise fashion Step-wise HRV changes in 14 patients. Seven patients had left before right pulmonary vein ablation and seven patients had right before left pulmonary vein wide area circumferential ablation. LF and HF HRV indices were measured in a step-wise fashion at each stage of the procedure, including: baseline, after transseptal puncture, after left-sided circumferential PV isolation and after right-sided circumferential PV isolation. HF high frequency, LF low frequency, HRV heart rate variability, PV pulmonary veins

The effect of patient age and procedure duration on heart rate variability changes

In order to determine whether procedure duration had an effect on the changes in HRV, we carried out a multivariate analysis to compare changes recorded in the AF ablation (n = 15) and the non-AF procedures in which no isoproterenol was used (n = 34) (Table 3). There was only a non-significant trend towards falling LF and HF after longer procedures. The effect of AF ablation is much larger, is clearly statistically significant, is equally evident in both LF and HF HRV, and persists even after adjusting for any effect of age or procedure duration. The unadjusted regression coefficients for the between group differences in log (LF) HRV and log (HF) HRV were −0.61 (95%CI: −0.99, −0.22, p = 0.003) and −0.40 (−0.73, −0.06, p = 0.023), and following adjustment for age and procedure duration: −0.51 (−0.88, −0.14, p = 0.008) and −0.33 (−0.61, −0.05, p = −0.021), respectively.

Table 3

Univariate and multivariate regression coefficients

	Regression coefficient	SE	p
Dependent variable = delta log LF
Univariate
Age	−0.003	0.005	0.62
Procedure duration	−0.005	0.002	0.061
AF ablation vs. non-AF ablation	−0.608	0.191	0.003
Mutivariate
AF ablation vs. non-AF ablation adjusted for age and procedure	−0.512	0.183	0.008
Dependent variable = delta log HF
Univariate
Age	0.000	0.005	0.93
Procedure duration	−0.003	0.002	0.137
AF ablation vs. non-AF ablation	−0.395	0.169	0.023
Multivariate
AF ablation vs. non-AF ablation adjusted for age and procedure	−0.333	0.139	0.021

Dependent variables: delta log LF and delta log HF in participants who did not receive isoproterenol (n = 49)

Discussion

This study presents and validates a technique for measuring HRV using shorter windows of time than are typically used. Averaging windows of 40 s (a duration that is easily visualised on a single screen during a procedure to establish that it is ectopic-free) gives data comparable to 5-min values which are conventionally recommended for ambulant outpatients undergoing routine HRV assessment, and readily detects intra-procedural autonomic changes within individuals.

Although HRV changes have been previously documented following AF ablation, the use of this technique intra-procedurally has, for the first time in humans, demonstrated that significant HRV reduction occurs following right, but not left-sided pulmonary vein wide area circumferential ablation.

Real-time intra-procedural recordings: validity and potential utility

We showed that average LF and HF indices of HRV remained stable with recording durations as brief as 40 s compared with 5 min recordings in healthy volunteers, despite the increase scatter seen with the 40-s window. The five 40-s recordings used encompassed a total of 200 s, which is approximately 20 cycles worth of the typical 10-s cycles of LF oscillations and many more cycles of HF oscillations: from basic mathematical principles, this is adequate to quantify oscillatory amplitude.

This method of monitoring HRV reliably detected a significant increase in LF and HF HRV immediately preceding syncope in patients undergoing tilt table testing [13]. This short recording window identified HRV changes in patients undergoing catheter ablation, which were consistent with previous studies using longer term recordings [11, 14].

In addition, short-segment HRV measurements taken at sequential stages of an AF ablation procedure identified significant intraprocedural changes in HF and LF HRV following right-sided circumferential pulmonary vein ablation, but not left-sided circumferential ablation or transseptal puncture.

Our study shows that the principle of HRV analysis using multiple short-term windows is valid. Ectopy-free segments of this duration are easy to collect, and clearly the effect size is large enough to be detectable in such recordings. Thus short-term recordings may have potential value during or immediately after ablation procedures.

Specific immediate effect on heart rate variability of AF ablation

An important consideration when measuring changes in HRV is to ensure a constant RR interval, since any cause of reduced RR interval (i.e. sinus tachycardia) is necessarily accompanied by a smaller absolute amplitude of variability of the RR interval. We detected a significant reduction in HF and LF HRV following AF ablation, without a reduction in the RR interval, suggesting that the reduction is due to modifications of autonomic nervous modulation.

Ablation of other arrhythmias did not affect HF or LF. However, if isoproterenol was given during the procedure, mean RR interval fell as expected, as did the HF and LF parameters of HRV.

Autonomic modulation of the sinus node

Autonomic modulation of the sinus node causes physiological oscillations in heart rate leading to LF and HF HRV [9]. Hou et al. [15] showed in animal studies that a significant amount of sympathetic and parasympathetic input to the sinus node traverses the right anterior ganglionated plexus located in the antral region of the right superior pulmonary vein. We postulate that in humans, a similarly sited ganglionated plexus is responsible for most of the nervous traffic to the sinus node and in turn can influence LF and HF variations in the heart rate. Our study showed HRV reduction occurred only following ablation around the right-sided pulmonary veins, which could conceivably disrupt autonomic nervous inputs into the sinus node.

Implications for autonomic ablation

Advances in understanding of functional cardiac neuroanatomy [16, 17] have led to development of autonomic ablation procedures that have been performed in humans in an attempt to treat cardiac arrhythmias including AF [3, 4, 6, 8] and neurocardiogenic syncope [18]. Autonomic ablation has been performed based on fixed anatomical landmarks [6], identification of distinct spectral patterns by fast Fourier transforms of endocardial signals [18], and identification and ablation of sites guided by high frequency stimulation [3, 4, 7] resulting in vagal responses. However, these methods do not provide a quantifiable measure of autonomic modification. Real-time measurement of HRV has the potential to be developed into a tool to assess intraprocedural autonomic modification.

Limitations

Healthy volunteers, typical of the general population in which conventional HRV measurements were developed, were used in the validation arm of this study to ascertain the minimum recording period that still yielded valid HRV results. In such subjects, ectopy-free ECG segments are easy to obtain. This is not always the case in patients undergoing catheter ablation, making a reason to study the usability of shorter (40 s) recordings for which ectopy-free segments are more easily recorded. Detectable changes in HRV in such short recordings increase the potential practical clinical value of the findings during ablation procedures.

The changes in HRV during AF ablation appear to be related only to right-sided circumferential pulmonary vein ablation. There may, however, by important effects of ablation of autonomic ganglia in other regions of the heart (including adjacent to the left pulmonary veins), which are not reflected by HRV changes, but which might still affect ablation outcomes.

These findings may stimulate further development into acute real-time longer recordings, which would give greater confidence in the values obtained. We are currently developing a suitable processing system, contained within a “black box” that will permit on-line processing of the data to obtain intra-procedural HRV measurements during adjunctive autonomic ablation procedures. We envisage that the ability to perform such intra-procedural HRV measurements will facilitate further studies to improve our understanding of autonomic mechanisms of arrhythmias, which may ultimately lead to development of clinically relevant autonomic ablation procedures.

Conclusions

Multiple windows as short as 40 s give valid measurements of mean HF and LF HRV. Such measurements are easily obtainable during routine EP procedures.

Left atrial ablation around the right-sided pulmonary veins is unique in producing reductions in LF and HF HRV, which is immediately detectable from multiple 40-s recordings. The ability to detect LF and HF HRV modification may, in the future, permit valuable intra-procedural measurements, ultimately able to establish end-point markers for adjunctive autonomic ablation procedures.