Acessibilidade / Reportar erro

Critical analysis of autoregressive and fast Fourier transform markers of cardiovascular variability in rats and humans

Abstract

The autonomic nervous system plays an important role in physiological and pathological conditions, and has been extensively evaluated by parametric and non-parametric spectral analysis. To compare the results obtained with fast Fourier transform (FFT) and the autoregressive (AR) method, we performed a comprehensive comparative study using data from humans and rats during pharmacological blockade (in rats), a postural test (in humans), and in the hypertensive state (in both humans and rats). Although postural hypotension in humans induced an increase in normalized low-frequency (LFnu) of systolic blood pressure, the increase in the ratio was detected only by AR. In rats, AR and FFT analysis did not agree for LFnu and high frequency (HFnu) under basal conditions and after vagal blockade. The increase in the LF/HF ratio of the pulse interval, induced by methylatropine, was detected only by FFT. In hypertensive patients, changes in LF and HF for systolic blood pressure were observed only by AR; FFT was able to detect the reduction in both blood pressure variance and total power. In hypertensive rats, AR presented different values of variance and total power for systolic blood pressure. Moreover, AR and FFT presented discordant results for LF, LFnu, HF, LF/HF ratio, and total power for pulse interval. We provide evidence for disagreement in 23% of the indices of blood pressure and heart rate variability in humans and 67% discordance in rats when these variables are evaluated by AR and FFT under physiological and pathological conditions. The overall disagreement between AR and FFT in this study was 43%.

Fast Fourier transform; Autoregressive model; Hypertension; SHR; Pharmacological blockade; Tilt test


Braz J Med Biol Res, April 2009, Volume 42(4) 386-396

Critical analysis of autoregressive and fast Fourier transform markers of cardiovascular variability in rats and humans

G.J.J. Silva1,3*, M.R. Ushizima2*, P.S. Lessa1, L. Cardoso3, L.F. Drager1, M.M. Atala1, F.M. Consolim-Colombo1, H.F. Lopes1, I.A. Cestari2, J.E. Krieger3 and Correspondence and Footnotes E.M. Krieger1

1Unidade de Hipertensão, 2Divisão de Bioengenharia, Instituto do Coração (InCor), Hospital das Clínicas, Faculdade de Medicina, Universidade de São Paulo, São Paulo, SP, Brasil

3Departamento de Medicina (LIM13), Faculdade de Medicina, Universidade de São Paulo, São Paulo, SP, Brasil

References Correspondence and Footnotes Correspondence and Footnotes Correspondence and Footnotes

Abstract

The autonomic nervous system plays an important role in physiological and pathological conditions, and has been extensively evaluated by parametric and non-parametric spectral analysis. To compare the results obtained with fast Fourier transform (FFT) and the autoregressive (AR) method, we performed a comprehensive comparative study using data from humans and rats during pharmacological blockade (in rats), a postural test (in humans), and in the hypertensive state (in both humans and rats). Although postural hypotension in humans induced an increase in normalized low-frequency (LFnu) of systolic blood pressure, the increase in the ratio was detected only by AR. In rats, AR and FFT analysis did not agree for LFnu and high frequency (HFnu) under basal conditions and after vagal blockade. The increase in the LF/HF ratio of the pulse interval, induced by methylatropine, was detected only by FFT. In hypertensive patients, changes in LF and HF for systolic blood pressure were observed only by AR; FFT was able to detect the reduction in both blood pressure variance and total power. In hypertensive rats, AR presented different values of variance and total power for systolic blood pressure. Moreover, AR and FFT presented discordant results for LF, LFnu, HF, LF/HF ratio, and total power for pulse interval. We provide evidence for disagreement in 23% of the indices of blood pressure and heart rate variability in humans and 67% discordance in rats when these variables are evaluated by AR and FFT under physiological and pathological conditions. The overall disagreement between AR and FFT in this study was 43%.

Key words: Fast Fourier transform; Autoregressive model; Hypertension; SHR; Pharmacological blockade; Tilt test

Introduction

The homeostasis of the cardiovascular system is efficiently maintained under control by feedback mechanisms that act to maintain perfusion pressure to the organs within a relatively narrow range of variation. The autonomic nervous system plays an important role in blood pressure regulation, producing adjustments in heart rate and myocardium contractility by altering sympathetic and parasympathetic activities to the heart, and peripheral vascular resistance and blood volume by constriction of arterial and venous vessels, respectively (1). These effects occur during each cardiac cycle and can be assessed by the variations in R-R intervals of the electrocardiogram.

In the early 70's, spectral analysis was applied to investigate cardiac interval time series (2,3) and to describe the rhythmic fluctuations that occur in the band of the respiratory frequency or high frequencies (Hering waves) and in lower frequencies (Mayer waves). In 1981, pharmacological studies showed for the first time that these variations were associated with the autonomic nervous system and hormonal mechanisms (4), suggesting a relationship between the spectral components and control mechanisms of the cardiovascular system. Taken together, the spectral analysis components facilitated the physical interpretation of the mechanisms involved in the regulation of the cardiovascular system. Later, these findings lent support to the sympathovagal balance hypothesis (5), whereby the normalized power of the low-frequency (LFnu) component of heart rate variability was associated with the sympathetic modulation and the high-frequency (HFnu) component was associated with parasympathetic modulation. The ratio between these two components established an index of the autonomic influence to the heart. This hypothesis was validated during postural changes (or orthostatic test) that impose perturbation on the cardiovascular system, revealing potential applications to assess physiological and physiopathological conditions.

The evaluation of autonomic markers on the basis of spectral analysis can be obtained by parametric (autoregressive model, AR, or autoregressive moving average) or non-parametric (fast Fourier transform, FFT) methods. In clinical research, both spectral techniques have been applied for the evaluation of the autonomic indices in myocardial infarction (6), Chagas disease (7) and diabetes (8), and in the validation of noninvasive blood pressure monitoring (9). However, a consensus regarding which is the most appropriate method for power spectral density interpretation of the cardiovascular system time series has not been established (10-12). The task force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology established the standardization of measurement, physiological interpretation, and clinical use of power spectral analysis (10), and the advantages and disadvantages of both parametric and non-parametric methods have been described by others (12).

In view of these considerations, we investigated the agreement of the AR and FFT methods regarding the extraction of autonomic markers during 1) postural changes (tilt test) in healthy volunteers, 2) pharmacological blockade of the autonomic nervous system to the heart in normal control rats, 3) in patients with arterial hypertension, and 4) in an experimental model of hypertension (spontaneously hypertensive rats, SHR).

Material and Methods

Animal model and subjects

All experimental procedures followed institutional guidelines for the care and use of laboratory animals (rats) and protocols used here were approved by the Institutional Review Board of the Medical School, University of São Paulo, São Paulo, SP, Brazil (humans).

Animals. Ten male normotensive Wistar rats and ten 18-week-old male SHR (Medical School, University of São Paulo) were fed standard laboratory chow and water ad libitum while housed (2-3 per cage) in a temperature-controlled room (22°C) with a 12-12-h dark-light cycle.

Human subjects. Sixteen normotensive volunteers and 8 subjects with mild hypertension and 17 subjects with severe hypertension were included in the study. Exclusion criteria were: age of less than 20 or more than 50 years, diabetes mellitus, cerebrovascular, aortic, or cardiac disease, smoking habit, and body mass index greater than 35 kg/m2 (Table 1). Hypertensive patients had their medication regimens changed with the progressive withdrawal of β-blocker, angiotensin-converting enzyme inhibitors and central α-blocker over a 2-week period.

Data acquisition and analysis

Experimental model. Under anesthesia with 30 mg/kg pentobarbitone sodium administered ip, polyethylene catheters (PE-10, 0.28 mm ID, 0.61 mm OD, Biocorp Australia, Australia) filled with heparin solution were inserted into the abdominal aorta and inferior vena cava through the left femoral artery and vein, respectively. The free ends of the cannulas were tunneled subcutaneously and exteriorized in the nape. Blood pressure was measured in conscious freely moving rats at least 24 h after the surgical procedures by connecting the arterial cannula to a pressure transducer (Statham P23Dd, Puerto Rico). The blood pressure waveform signal from the transducer was fed to an amplifier (GPA-4 model 2, Stemtech, Inc., USA), and to a 10-bit analog-to-digital converter (DataQ Instruments, Inc., USA) for measurement on a beat-to-beat basis at a sampling rate of 2.0 kHz.

Human subjects. All patients were instructed not to consume alcohol or caffeine 24 h before the protocol. During all experimental conditions, blood pressure was measured with a digital photoplethysmograph, which provides accurate beat-to-beat systolic and diastolic values (Finapress, Omeda 2300, Monitoring Systems, USA). This device uses a photoplethysmographic finger cuff to assess blood pressure continuously using the vascular unloading principle. Heart rate was monitored by electrocardiography (ECG Amplifier, model 13-4615-64, Gould, USA). The blood pressure waveforms and ECG signals were recorded on a beat-to-beat basis (DataQ Instruments) at a sampling rate of 1.0 kHz.

Signal edition and artifact correction and rejection

Blood pressure waveforms were used to extract beat-by-beat time series of pulse (PI), and systolic (SBP) and diastolic (DBP) blood pressure. R-R interval time series were derived from ECG recordings. The signal pre-processing and spectral analysis were conducted by the same observer. Peak pressure is reported as SBP and the minimum pressure is reported as DBP. The tagged points were reviewed and edited if required after visual inspection with a special signal editor developed for Matlab (MATLAB 6.0, Mathworks, USA) (13). Pulse interval or R-R interval was calculated as the difference between the beginning and end points of the cycle (t1-t0).

Blood pressure and pulse/R-R interval variability

Time-domain analysis consisted of calculating mean PI or R-R intervals and SBP and DBP, followed by their variability as the variance from their respective time series. In the frequency-domain analysis, both AR and FFT methods were used to evaluate blood pressure and pulse interval and R-R interval variability. Both time- and frequency-domain parameters were determined by software. The spectral bands for humans (very low-frequency (VLF): 0.0-0.04 Hz; LF: 0.04-0.15 Hz; HF: 0.15-0.4 Hz) and rats (VLF: 0.0-0.2 Hz; LF: 0.2-0.8 Hz; HF: 0.8-2.8 Hz) were defined according to literature references (10,14). Spectral power for VLF, LF and HF bands was calculated by power spectrum density integration within each frequency band. The power density of each spectral component was calculated both as absolute values and as normalized units (nu) (10,15). The power at LF and HF for pulse interval was normalized and represents the relative value of each power component in proportion to the total power minus the VLF component. Data are reported as nu. The sympathovagal balance was defined by the LF/HF ratio. The LF component of the R-R and pulse intervals and of blood pressure variability was considered to be a marker of efferent sympathetic cardiac and vascular modulation, respectively, whereas the HF component of the R-R or pulse interval variability reflects respiratory-driven vagal modulation to the sinoatrial node (4,16).

For frequency-domain analysis, the whole 10-min time series of blood pressure and pulse or R-R intervals were evaluated under basal conditions using parametric and non-parametric methods, described in detail below.

Autoregressive method. Time series were analyzed using the Linear Analysis software (LA, version 8.5) (17), kindly provided by Alberto Porta (Department of Technologies for Health, Galeazzi Orthopaedic Institute, University of Milan, Italy). Time series were divided into segments of 300 beats that overlapped by 50%. A linear de-trending procedure was used, the spectra of each segment were calculated via the Levinson-Durbin recursion, and the order of the model was chosen according to Akaike's criterion.

Fast Fourier transform method. The linear trend was removed before the analysis of power spectral density obtained by FFT, using Welch's method. For the signal obtained in rats, we considered 10 min of data resampled at 13.8 Hz (first 213 points) and segments of 211 points with a Hanning window and 50% overlap. For the signal obtained in humans, we considered 10 min of data resampled at 5.0 Hz (first 211) and segments of 29 points with a Hanning window and 50% of overlap.

Experimental protocols

Tilt test. An established protocol currently used by Montano et al. (18) was adopted. During the entire experimental procedure, the normal volunteers (N = 8) were monitored by blood pressure and ECG measurements. Head-up tilt testing was performed with the subject in a non-sedated, post-absorptive state. After heart rate and blood pressure measurement under basal conditions (10 min), the subject was positioned upright at an angle of approximately 80 degrees for a maximum of 20 min. If syncope or pre-syncope with hypotension developed during the test, the table was lowered to the supine position.

Autonomic nervous system blockade. The experimental protocol was similar to that used in a previous study performed in our laboratory (19,20). Briefly, a 2-day protocol was used for the evaluation of the vagal and sympathetic effects on heart rate of 1 mg/kg methylatropine (Sigma Chemical Co., USA) and 1 mg/kg propranolol (Sigma) injected iv. On the first day of the study, heart rate was recorded for 10 min under basal conditions and after parasympathetic pharmacological blockade with methylatropine. On the second day of the protocol, the sympathetic influence on heart rate was evaluated by propranolol injection followed by methylatropine injection to assess the double sympathetic and parasympathetic blockade in order to determine the intrinsic heart rate. The effects of pharmacological blockade on heart rate were recorded for 10 min after both methylatropine and propranolol reached their maximum effects (~3 min, iv injection).

Statistical analysis

Data are reported as means ± SEM. The agreement of all time- and frequency-domain data between AR and FFT was analyzed by two-way analysis of variance (ANOVA) for repeated measurements followed by the Bonferroni post hoc test. Demographic data were analyzed by one-way ANOVA followed by the Bonferroni post hoc test. Similarly, disagreement between AR and FFT was analyzed by the percentage of indices that presented statistical significance for the method factor (F1: AR vs FFT) in the two-way ANOVA. The level of significance was set at P < 0.05 in all analyses. All procedures were performed with the SigmaStat statistical software (Systat Software Inc., USA).

Results

Autonomic nervous system changes in response to postural and pharmacological blockade challenges

Normal healthy subjects (N = 8) submitted to postural challenge (orthostatic tilt test) exhibited hypotension (P = 0.039) accompanied by a significant decrease (P = 0.004) in R-R interval (Table 2 and Figure 1A). Neither method showed any significant reduction in R-R interval variance (P = 0.372). In addition, the FFT technique presented significantly higher values (P < 0.001) for R-R interval variance under both basal and tilt conditions compared to AR. Both the AR and FFT methods demonstrated a significant increase in LF for SBP (P = 0.022), and a decrease in HFnu for pulse interval (P = 0.003) induced by the tilt test; however, only the AR method was able to show changes in HF for SBP (P = 0.029) during the orthostatic tilt test. Moreover, a significant increase (P = 0.042) in LF/HF ratio after postural change was observed only with the AR method.

Pharmacological blockade (Table 3 and Figure 2) was used to assess the modulation of the autonomic nervous system produced by sympathetic and parasympathetic activities to the heart in an experimental model. Vagal blockade with methylatropine induced a significant reduction in pulse interval (P < 0.001). However, the blockade of this cardiac sympathetic activity using propranolol did not change pulse interval values (Figure 2A) compared to those of the basal condition. A significant difference in pulse interval variance was observed under basal conditions (P = 0.030) when AR and FFT were compared. The absolute values of LF (P = 0.363) and HF (P = 0.150), with power evaluated by both AR and FFT methods, were not affected by pharmacological blockade. However, the absolute value of the LF component, estimated by AR, was lower compared to FFT (P = 0.011) under basal conditions and after methylatropine blockade. Cardiac sympathetic blockade with propranolol (Figure 2B) tended to reduce the LFnu component by FFT, without any change in the AR method results. Furthermore, HFnu oscillation was significantly increased (P < 0.001) by propranolol blockade in both methods. Parasympathetic blockade with methylatropine (Figure 2C) did not change the LFnu component of either the AR or FFT method (P = 0.468). The reduction in HFnu oscillation of PI after parasympathetic blockade with methylatropine was detected by both the AR and FFT methods (P < 0.001). In addition, AR and FFT analysis was not concordant in the values for the LFnu (P < 0.001) and HFnu (P < 0.001) components under basal and methylatropine conditions. The LF/HF ratio for pulse interval (Figure 2D) was unchanged after both sympathetic and parasympathetic blockade when evaluated by the AR technique, but was increased after methylatropine injection when evaluated by the FFT method (P < 0.001). Furthermore, the AR and FFT techniques showed a disagreement in the LF/HF ratio of the pulse interval under basal and methylatropine conditions (P < 0.001).

Autonomic nervous system evaluation in hypertensive patients and in experimental hypertension

In patients, the autonomic changes commonly associated with arterial hypertension were used to evaluate the concordance of the FFT and AR methods (Table 4). As expected, SBP values increased significantly according to the severity of arterial hypertension (P < 0.001) compared to normotensive subjects. R-R interval values under basal conditions were significantly reduced in the hypertensive groups (P = 0.038). Both spectral estimation techniques demonstrated a significant increase in the values of SBP variance (P < 0.001) according to the severity of hypertension, and R-R interval variance decreased in both hypertensive (mild and severe) groups (P = 0.004). Neither the AR nor the FFT method showed any change in the LF (P = 0.688) or HF (P = 0.817) components of SBP in any of the hypertensive groups. The total power of SBP showed a nonsignificant increase according to the severity of hypertension only in the FFT method (P = 0.459). However, both total power (P < 0.001) and variance (P < 0.001) of SBP were significantly higher in FFT than in the AR method in all groups evaluated. The LFnu component for the R-R interval showed increasing values according to the severity of hypertension (P = 0.050) by both methods, whereas the HFnu component was significantly reduced (P = 0.012). The increase in the LFnu component was accompanied by a nonsignificant increase in LF/HF ratio (P = 0.065) in both methods, and presented a statistical difference between the AR and FFT methods in normotensive subjects (P = 0.031).

Table 5 shows the time- and frequency-domain data obtained in an experimental model of hypertension (SHR). Similar to humans, basal SBP was significantly increased (P < 0.001) in SHR (N = 8) compared to normotensive Wistar rats (N = 8). Pulse interval values under basal conditions did not differ between groups (P = 0.597) and were not changed by the method (P = 0.714). High blood pressure levels were accompanied by significant increases in SBP (P = 0.004) and pulse interval (P = 0.005) variances, shown only by the FFT method. Both the AR and FFT methods demonstrated an increase in the power of LF oscillation for SBP (P = 0.012) in SHR compared to normotensive Wistar rats. Despite the relative concordance of the two methods, AR presented different variance values (P < 0.001) and total power estimation (P < 0.001) for SBP compared to FFT in the SHR group. Similarly, for pulse interval analysis, the AR and FFT methods presented a significant difference in variance (P < 0.001), absolute (P = 0.030) and normalized (P = 0.004) values of the LF component, absolute values of HF (P = 0.008), LF/HF ratio (P = 0.010), and total estimation power (P < 0.001).

Figure 1.
Mean pulse interval (PI, Panel A) and LFnu (Panel B) and HFnu (Panel C) oscillations, and LF/HF ratio (Panel D) for PI at basal condition and after tilt test in normal volunteers. LF = low-frequency; HF = high-frequency; LFnu = normalized power of LF; HFnu = normalized power of HF; AR = autoregressive method; FFT = fast Fourier transform method. *P ≤ 0.05, compared to respective basal group. +P ≤ 0.05, FFT compared to AR.

Figure 2.
Mean pulse interval (PI, Panel A) and LFnu (Panel B) and HFnu (Panel C) oscillations, and LF/HF ratio (Panel D) for PI at basal condition and after pharmacological blockade of sympathetic (propranolol) and parasympathetic (methylatropine) systems in rats. LF = low-frequency; HF = high-frequency; LFnu = normalized power of LF; HFnu = normalized power of HF; AR = autoregressive method; FFT = fast Fourier transform method. *P ≤ 0.05, compared to respective basal condition. +P ≤ 0.05, FFT compared to AR.

Table 1.
Demographic data of the normotensive and hypertensive subjects studied.

Table 2.
Data obtained for normal volunteers (N = 8) before (basal) and after (tilt) submission to the head-up tilt testing.

Table 3.
Time- and frequency-domain analysis of pulse interval of normotensive rats (Wistar, N = 8) submitted to pharmacological blockade.

Table 4.
Time- and frequency-domain analysis of systolic blood pressure (SBP) and R-R interval of normotensive subjects (N = 16), patients with mild hypertension (N = 8) and patients with severe hypertension (N = 17).

Table 5.
Time- and frequency-domain analysis of systolic blood pressure (SBP) and pulse interval (PI) of normotensive (Wistar, N = 8) and spontaneously hypertensive (SHR, N = 8) rats.

Discussion

The present study compared the agreement between the AR and FFT methods under several conditions: a physiological stimulus, postural changes induced in normal healthy volunteers or pharmacological blockade in an experimental model, and under pathological conditions, i.e., patients with mild and severe hypertension or a genetic experimental model of hypertension (SHR). Despite specific agreement between parametric and non-parametric methods in time- and frequency-domain indices, several disagreements were detected in the present study.

In the time-domain analysis, the AR and FFT methods seemed to be concordant, with the exception of differences detected in the calculation of the time variance of both PI/R-R interval and SBP under all conditions evaluated. Our data partially agree with the study of Pitzalis et al. (21) who demonstrated that the time-domain measurements may be considered to be reproducible. Most of the differences observed between the AR and FFT methods are related to the frequency-domain analysis. Considering all the frequency-domain indices evaluated, the present study found that 43% of the data did not agree. Interestingly, most of the non-concordant time- and frequency-domain data were related to protocols using experimental animals (67%); in humans, 23% of non-concordant values were found in the present study. Moreover, we observed that the LF component was not concordant between methods in any of the experimental models, i.e., pharmacological blockade and genetic hypertension. The AR method presented 2- and 5-fold lower LFnu values for pulse interval compared to FFT in SHR and normotensive Wistar rats, respectively. Additionally, we observed significant differences in the HFnu component of the R-R interval for both normotensive and hypertensive patients. Indeed, several studies (22-24) have shown an overestimation of the HF component by the FFT method. Furthermore, a nonsignificant increase was observed in the LF/HF ratio for the R-R interval using the FFT method during the tilt test in normal healthy volunteers, whereas the AR method showed a significant increase in sympathovagal balance during the protocol involving postural changes. In the animal model of hypertension, we observed higher values of the total power for the R-R interval when analyzed by the FFT method. Pichon et al. (24) observed significant differences between the AR and FFT methods in the total power estimation for both sitting and standing positions. Moreover, these investigators showed that the power of the HF component was greater in the FFT method, suggesting a poor interchangeability in the heart rate variability indices in normal volunteers when parametric and non-parametric methods were compared. The disagreement between parametric and non-parametric methods was not limited to the analysis of heart rate variability. We were able to detect lower levels of the LF component for SBP obtained with the FFT method, compared to AR in normal healthy volunteers. Increased values of total power estimation for SBP in hypertensive patients were also observed using the FFT method. The difference in the total power estimation was consistent with the data obtained in the experimental model of hypertension (SHR) for SBP.

With regard to pathological conditions (arterial hypertension), we achieved a total of 33% non-concordant data for time- and frequency-domain indices. Indeed, Chemla et al. (25) compared both parametric and non-parametric methods and showed that FFT overestimated the LF (ms2 and nu) component and LF/HF ratio for the R-R interval in diabetic patients, while the AR method overestimated the total power, VLF and HFnu components. These data suggested that parametric and non-parametric methods of spectral analysis were not comparable in diabetes. In hypertensive patients, Badilini et al. (26) demonstrated that, under basal conditions, the FFT method provided a lower estimation of VLF and total power for the R-R interval, but presented higher values for the LF and HF components compared to the AR method. The increase in LFnu and HFnu values for the R-R interval obtained by the FFT method compared to AR agrees with our data. We showed that the FFT method overestimated the normalized values of the HF component for the R-R interval in both normotensive and severe hypertensive subjects, and overestimated the total power for SBP only in the normotensive subjects.

The Linear Analysis software used in the present investigation for the autoregressive technique showed several undetected or missing VLF, LF and HF values in the spectral analysis of both blood pressure and PI/R-R interval variability in most of the patients and animals studied. On the other hand, all spectral components obtained with the FFT analysis were detected in animal models and patients. Similarly, Chemla et al. (25) recently reported that the use of the AR method resulted in 50% missing values of both LF and HF components. These investigators showed that changing the model order of the analysis did not modify the null or missing data observed with the AR method.

Although a consensus regarding the most appropriate method for power spectral density estimation of cardiovascular system time series has not been established (10-12), several advantages and disadvantages regarding parametric and non-parametric methods have been described. The advantages of the non-parametric (FFT) methods are: the simplicity of the algorithm employed, good reproducibility, and high processing speed. On the other hand, the advantages of parametric (AR) methods are: good performance in time series with reduced number of points, smoother spectral components, which can be distinguished independently of predefined frequency bands, easy post-processing of the spectrum with an automatic calculation of low- and high-frequency components, easy identification of the central frequency of each component, and an accurate estimation of power spectrum density even for a small number of samples in which the signal is supposed to remain stationary. The basic disadvantage of parametric methods is the need to verify the suitability of the chosen order of model. On the other hand, methods based on the non-parametric (FFT) technique require some experience in dealing with the non-stationary segments of the time series and with the overlapping and windowing to filter the power spectral density. Considering the findings of the present investigation, in which spectral analysis was affected by different species, autonomic stimulation or blockade, blood pressure, and AR or FFT methodological approach, the present investigation strongly recommends that studies or groups that intend to use spectral analysis as a window to the autonomic nervous system use their signal database of well-known experimental protocols with autonomic nervous challenges to construct and share some pairs of data and spectral "fingerprints" to compare each new finding and conclusion within and between research groups.

Experimental investigations (27) have suggested that changes in autonomic balance might be helpful to identify high cardiovascular risk, i.e., life-threatening arrhythmias during acute myocardial ischemia. Based on this experimental evidence, in a prospective study, La Rovere et al. (28) demonstrated that time- and frequency-domain indices of heart rate variability have independent prognostic values after myocardial infarction. However, we are concerned with the use of autonomic markers derived from different spectral analysis methods to stratify cardiovascular risk. The dependence on the technique adopted and the absence of agreement between AR and FFT spectral analysis limit the direct application of these autonomic markers in clinical practice. In fact, the present study and others (21,22,24,25) have demonstrated that the AR and FFT methods are not exactly interchangeable.

Taken together, the data obtained in the present study demonstrate that the frequency-domain analysis, using parametric (AR) or non-parametric (FFT) methods do not agree in several aspects. The absence of agreement between methods extended to physiological conditions, i.e., postural changes in humans or pharmacological blockade in experimental animals. The typical changes in the autonomic markers observed during physiological conditions were not equally detected by AR and FFT. Considering the species studied, several differences were observed under basal conditions in both human and rat models when AR and FFT routines were applied for time- and frequency-domain analysis. Under pathological conditions in patients, such as arterial hypertension, both methods seem to be adequate for power spectral estimation and autonomic index estimation. However, in a genetic model of spontaneous hypertension (i.e., SHR) the parametric and non-parametric methods evaluated were not concordant. Furthermore, the present study showed that the disagreement between parametric and non-parametric methods was not restricted to the analysis of heart rate variability, but affected equally the spectral analysis of systolic blood pressure time series, and the overall disagreement between AR and FFT in this study was 43%.

Address for correspondence: E.M. Krieger, Unidade de Hipertensão, InCor, Av. Dr. Enéas C. Aguiar, 44, 2º andar, Bloco II, 05403-000 São Paulo, SP, Brasil. Fax: +55-11-3069-5048. E-mail: edkrieger@incor.usp.br

*These authors contributed equally to this study. Research supported by FINEP, FAPESP and the Zerbini Foundation. Received June 18, 2008. Accepted February 26, 2009.

  • 1. Cowley AW Jr. Long-term control of arterial blood pressure. Physiol Rev 1992; 72: 231-300.
  • 2. Sayers BM. Analysis of heart rate variability. Ergonomics 1973; 16: 17-32.
  • 3. Hyndman BW, Kitney RI, Sayers BM. Spontaneous rhythms in physiological control systems. Nature 1971; 233: 339-341.
  • 4. Akselrod S, Gordon D, Ubel FA, Shannon DC, Berger AC, Cohen RJ. Power spectrum analysis of heart rate fluctuation: a quantitative probe of beat-to-beat cardiovascular control. Science 1981; 213: 220-222.
  • 5. Pagani M, Lombardi F, Guzzetti S, Rimoldi O, Furlan R, Pizzinelli P, et al. Power spectral analysis of heart rate and arterial pressure variabilities as a marker of sympatho-vagal interaction in man and conscious dog. Circ Res 1986; 59: 178-193.
  • 6. Kleiger RE, Miller JP, Bigger JT Jr, Moss AJ. Decreased heart rate variability and its association with increased mortality after acute myocardial infarction. Am J Cardiol 1987; 59: 256-262.
  • 7. Guzzetti S, Iosa D, Pecis M, Bonura L, Prosdocimi M, Malliani A. Impaired heart rate variability in patients with chronic Chagas' disease. Am Heart J 1991; 121: 1727-1734.
  • 8. Bianchi A, Bontempi B, Cerutti S, Gianoglio P, Comi G, Natali Sora MG. Spectral analysis of heart rate variability signal and respiration in diabetic subjects. Med Biol Eng Comput 1990; 28: 205-211.
  • 9. Parati G, Casadei R, Groppelli A, Di Rienzo M, Mancia G. Comparison of finger and intra-arterial blood pressure monitoring at rest and during laboratory testing. Hypertension 1989; 13: 647-655.
  • 10. Heart rate variability. Standards of measurement, physiological interpretation, and clinical use. Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. Eur Heart J 1996; 17: 354-381.
  • 11. Malik M. Sympathovagal balance: a critical appraisal. Circulation 1998; 98: 2643-2644.
  • 12. Parati G, Saul JP, Di Rienzo M, Mancia G. Spectral analysis of blood pressure and heart rate variability in evaluating cardiovascular regulation. A critical appraisal. Hypertension 1995; 25: 1276-1286.
  • 13. Ushizima MR. Desenvolvimento de ferramentas para análise de sinais biológicos nos domínios tempo, freqüência e tempo-freqüência: aplicação ao estudo da regulação da pressão arterial. [Doctoral thesis]. Campinas: Universidade Estadual de Campinas; 2000.
  • 14. Ushizima MR, Moreira ED, Costa ET, Castiglioni P, Krieger EM, Di Rienzo M, et al. Effects of sinoaortic denervation on blood pressure, pulse interval, sympathetic nerve activity and diaphragmatic electromyogram variability in conscious rats. IEEE Proceedings of Computers in Cardiology 2001; 485-488.
  • 15. Pagani M, Montano N, Porta A, Malliani A, Abboud FM, Birkett C, et al. Relationship between spectral components of cardiovascular variabilities and direct measures of muscle sympathetic nerve activity in humans. Circulation 1997; 95: 1441-1448.
  • 16. Japundzic N, Grichois ML, Zitoun P, Laude D, Elghozi JL. Spectral analysis of blood pressure and heart rate in conscious rats: effects of autonomic blockers. J Auton Nerv Syst 1990; 30: 91-100.
  • 17. Rubini R, Porta A, Baselli G, Cerutti S, Paro M. Power spectrum analysis of cardiovascular variability monitored by telemetry in conscious unrestrained rats. J Auton Nerv Syst 1993; 45: 181-190.
  • 18. Montano N, Ruscone TG, Porta A, Lombardi F, Pagani M, Malliani A. Power spectrum analysis of heart rate variability to assess the changes in sympathovagal balance during graded orthostatic tilt. Circulation 1994; 90: 1826-1831.
  • 19. De Angelis K, Wichi RB, Jesus WR, Moreira ED, Morris M, Krieger EM, et al. Exercise training changes autonomic cardiovascular balance in mice. J Appl Physiol 2004; 96: 2174-2178.
  • 20. Gava NS, Veras-Silva AS, Negrao CE, Krieger EM. Low-intensity exercise training attenuates cardiac beta-adrenergic tone during exercise in spontaneously hypertensive rats. Hypertension 1995; 26: 1129-1133.
  • 21. Pitzalis MV, Mastropasqua F, Massari F, Forleo C, Di Maggio M, Passantino A, et al. Short- and long-term reproducibility of time and frequency domain heart rate variability measurements in normal subjects. Cardiovasc Res 1996; 32: 226-233.
  • 22. Badilini F, Maison-Blanche P, Coumel P. Heart rate variability in passive tilt test: comparative evaluation of autoregressive and FFT spectral analyses. Pacing Clin Electrophysiol 1998; 21: 1122-1132.
  • 23. Fagard RH, Pardaens K, Staessen JA, Thijs L. Power spectral analysis of heart rate variability by autoregressive modelling and fast Fourier transform: a comparative study. Acta Cardiol 1998; 53: 211-218.
  • 24. Pichon A, Roulaud M, Antoine-Jonville S, de Bisschop C, Denjean A. Spectral analysis of heart rate variability: interchangeability between autoregressive analysis and fast Fourier transform. J Electrocardiol 2006; 39: 31-37.
  • 25. Chemla D, Young J, Badilini F, Maison-Blanche P, Affres H, Lecarpentier Y, et al. Comparison of fast Fourier transform and autoregressive spectral analysis for the study of heart rate variability in diabetic patients. Int J Cardiol 2005; 104: 307-313.
  • 26. Badilini F, Maison-Blanche P, Champomier P, Provost JC, Coumel P, Milon H. Frequency-domain heart rate variability in 24-hour Holter recordings: role of spectral method to assess circadian patterns and pharmacological autonomic modulation. J Electrocardiol 2000; 33: 147-157.
  • 27. Schwartz PJ, Vanoli E, Stramba-Badiale M, De Ferrari GM, Billman GE, Foreman RD. Autonomic mechanisms and sudden death. New insights from analysis of baroreceptor reflexes in conscious dogs with and without a myocardial infarction. Circulation 1988; 78: 969-979.
  • 28. La Rovere MT, Bigger JT Jr, Marcus FI, Mortara A, Schwartz PJ. Baroreflex sensitivity and heart-rate variability in prediction of total cardiac mortality after myocardial infarction. ATRAMI (autonomic tone and reflexes after myocardial infarction) investigators. Lancet 1998; 351: 478-484.
  • Correspondence and Footnotes

  • Publication Dates

    • Publication in this collection
      20 Mar 2009
    • Date of issue
      Apr 2009

    History

    • Accepted
      26 Feb 2009
    • Received
      18 June 2008
    Associação Brasileira de Divulgação Científica Av. Bandeirantes, 3900, 14049-900 Ribeirão Preto SP Brazil, Tel. / Fax: +55 16 3315-9120 - Ribeirão Preto - SP - Brazil
    E-mail: bjournal@terra.com.br