Time-Domain Analysis of Heart Rate Variability in Diabetic Patients with and without Autonomic Neuropathy
Introduction
The normal heart rate is determined by dynamic interaction between the spontaneous cardiac impulse generated by the sinoatrial (SA) node and conflicting influences of the sympathetic and parasympathetic nervous systems on the conducting tissue of the heart. The rate of spontaneous depolarization of the SA node is itself affected by its metabolic milieu and in the longer term, by hormonal influences. Normal resting heart rate is maintained by the tonic influence of the parasympathetic vagus nerve, and acceleration of the heart rate is affected both by the inhibition of vagal influence and the stimulation of the sympathetic nervous system. The activity of the autonomic is also governed by moment-to-moment changes in blood pressure and respiration, which alter heart rate continuously. The resultant heart rate is the summation of all these influences and thus inherently unpredictable on a beat-to-beat basis.
Recent advances in technology have enabled accurate recordings and the automated analysis of 24-hour ECG to detect beat-to-beat variability, providing not only more detailed, but much more accurate and precise information than the earlier tests.1 Heart rate variability (HRV) decreases with age,2 and shows a circadian variation, being maximum during sleep.3 It is also rate dependent, the heart rate showing more variability at lower heart rates. The loss of this beat-to-beat variability is a sign of disease.
It has long been known that cardiovascular autonomic diabetic neuropathy (CADN) is associated with a loss of heart rate variability.4,5 These patients have a poor cardiovascular prognosis,6,7 with a 5-year mortality greater than 50%.8 Some of this may be attributed to micro- or macrovascular disease,9 however, a recent study has shown the relatively poor prognosis of patients with CADN in the absence of clinically detectable micro- or macrovascular conditions.10 Clinically detectable autonomic failure is usually evident many years or decades after the onset of diabetes. It is likely that these patients develop subtle deficits in HRV much earlier, and these may include diminution in time-domain analysis. Detection of such changes may be used as markers of pathology, particularly to study the benefits of therapeutic interventions. Thus, the aim of this work was to study HRV in diabetic patients with clinical and sub-clinical autonomic neuropathy (AN), and to determine whether HRV in patients with sub-clinical AN is abnormal in comparison to normal controls.
SUBJECTS AND METHODS
Following Local Ethical Committee approval, we studied 10 consecutive, newly diagnosed non-insulin dependent diabetics mellitus (NIDDM), 10 with NIDDM with autonomic diabetic neuropathy (ADN), and 10 age-matched normal controls (Table 1). Informed consent was obtained from each subject. The patients were selected from King Abdulaziz University Hospital in Jeddah, Saudi Arabia, under supervision of a medical consultant. The analysis of HRV was done in the Department of Physiology. We found no gender differences in any measure of HRV (data not shown), which is consistent with previous published work.10 All subjects had no evidence of pre-existing cardiovascular disease, were normotensive, had no previous history of ill health, and were taking no regular medication except standard diabetic therapy.
![]() Table 1. Age and sex of D without AN, D with AN and normal subjects (n=10). Values are mean±1 SD and range. |
24-Hour ECG Recording
The 24-hour electrocardiogram (ECG) recordings were carried out using a two-channel tracker (DEL MAR Holter System). One channel records the ECG signal and the second channel records a time signal generated by the recorder. The time signal is used to make correction for variations in tape speed. Individuals went about normal daily activities. A new cassette tape and battery were used for each recording.
The electrodes were positioned in standard positions over the precordium. Reply, QRS detention and measurement of theRR interval was performed using DEL MAR ECG Analyzer. After analogue to digital conversion, a digital computer stored the time intervals (ms) between successive R waves as the interval tachogram (time series).
All intervals <200 ms, and >2000 ms were rejected as being artifacts. From the time series, using HRV analysis programs, time analysis was performed.
Time-Domain 24 Hours ECG Analysis
We measured several time domain variables (Table 2), which can be divided into three groups:
Measures of mean value and range, which includes mean RR and mean HR for the recording period;
Measures of variance, which include the standard deviation of all RR intervals around mean RR for the entire recording (SDRR), and St. George's index. The St. George's index was calculated from the frequency histogram of 24-h recordings. A triangle is fitted to the frequency histogram, the perpendicular height approximating to the model RR interval, and the area approximating to the total number of detected intervals. The index is calculated by dividing twice the area by height.11 The 24-hour SDRR and St. George's index are dependent not only on short-term variability but also on long-term variability and diurnal changes in mean RR interval. They, therefore depend not only on vagal modulation but also on cardiac sympathetic tone, neurohumoral activity and other factors. Hence, SDRR and St. George's index may be considered as a broadband measure of autonomic balance.
Parameters that analyze the differences between adjacent intervals including the standard deviation of the differences between adjacent RR intervals (root mean square of successive RR interval differences or RMSSD) and the Edinburgh index which is the number of times that the difference between adjacent normal RR intervals greater than 50 ms (NN50) is computed over the entire 24-hr ECG recording.12 NN50 and RMSSD quantitate short-term HRV and they reflect alterations that are almost wholly mediated by vagal tone. These parameters are considered as vagal indexes.13
![]() Table 2. Different time domain measures of heart rate variability. |
Statistical Analysis
Comparisons between data were made using Student's t-test for unpaired normally distributed values. However, because distributions of power spectrum values were skewed towards high values, these data were log transformed and then analyzed using Analysis Of Variance (ANOVA) test (general factorial model). P<0.05 was considered significant. Data were expressed as the mean ± SD and range.
RESULTS
The first step in any time series analysis is to plot the time series, so a graph of xt versus t was done. This was done to save us from rediscovering trivial aspects of the time series as well as suggest useful or interesting hypotheses to test. It has been reported that time series of RR intervals from patients with autonomic damage were noticeably smoother than those from normal subjects with large abrupt changes in heart rate.1 The beat-to-beat variation is represented as a time series (Figure 1A, B, C) in which the y axis is the time between beats (ms) and the x axis is the number of beats. The RR interval time series illustrates the normal pattern of fluctuations for the healthy adult (A) compared to the smooth pattern of beat-to-beat intervals for a diabetic patient with a symptomatic diabetic neuropathy (B) and a very smooth pattern of RR intervals for diabetic patient with clinical diabetic neuropathy (C). It was possible to observe a smaller RR variation in all diabetic subjects when compared to normal controls, and calculation confirmed this impression.

Example of RR interval series of 1000 RR intervals for: a) normal subject; b) diabetic patient with symptomatic autonomic neuropathy and; c) diabetic patient with clinical autonomic neuropathy.
Table 3. summarizes the result of time-domain analyses for normal controls, diabetic patients without autonomic neuropathy and patients with autonomic neuropathy. There was no significant difference in mean heart rate in diabetic patients compared to normal subjects [normal controls 79± 3 beats/min, D with AN 83±5(P=0.2) and D without AN 80±2 (P=0.6)]. All time-domain parameters (SDRR, SDRR index, SNN50, pNN50%, RMSSD and George index) were significantly decreased in patients with and without diabetic neuropathy compared to normal controls.
![]() Table 3. Time domain parameters of D without AN, D with AN, and normal subjects (n=10). Values are mean±SD (and range). |
There was no significant correlation between fasting blood sugar level and different HRV parameters in patients with and without diabetic autonomic neuropathy (results not shown).
DISCUSSION
Diabetes is known to be complicated by a dreadful disease, autonomic neuropathy, which affects about 40% of all diabetic patients and, in its severe form, offers a very poor prognosis. As HRV reflects the degree of autonomic control of the heart, it is widely used for the diagnosis of autonomic dysfunctions in noncardiological diseases on the assumption that if such a dysfunction in the heart is identified, it is a sign of a more widespread autonomic neuropathy affecting all organs. In the case of a systemic disease such as diabetes, this assumption is usually justified. Another reason why we used tests based on HRV as marker for cardiac autonomic diabetic neuropathy is because they are very accurate, noninvasive and generally reproducible.
In Saudi Arabia, the clinical importance of HRV tests is easy to appreciate if one considers the high incidence of diabetes in the Saudi population and the high incidence of diabetic neuropathy among diabetics. However, despite this high incidence of diabetes, a review of the literature revealed that our work is the first study that used tests based on HRV as markers for diabetic neuropathy in Saudi patients. Using time domain analysis of HRV, our study showed that despite the normal mean resting heart rate (as measured by RR interval) in diabetic subjects, there was a striking reduction in HRV when compared with normal controls. This reduction in HRV can be present even in patients without any clinical features of autonomic neuropathy. It has been reported1 that the measurement of risk RR interval changes, such as SDRR index, pNN50% and RMSSD aim to evaluate parasympathetic activity, which, of course, is directly responsible for rapid HR fluctuation. From the results that have been obtained in this study, it appears that vagal indices are decreased in diabetic patients compared to normal. However, the reduction is more pronounced in patients with autonomic neuropathy. Decreased parasympathetic activity is, of course, responsible for high heart rate usually reported in patients with diabetic autonomic neuropathy.
Although this method is very interesting, and has met with success in cardiological applications,1 it is not yet widely used by diabetologists, perhaps because it is not easy to carry out due to the large number of patients affected by diabetes mellitus.
ARTICLE REFERENCES:
1. . "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 Electrophysiology" . Circulation. 1996; 93: 1043–65.
2. . "Heart rate variability in healthy subjects: effect of age and the derivation of normal ranges for test of autonomic function" . Br Heart J. 1986; 55:348–54.
3. . "Circadian variation of heart rate variability" . Cardiovasc Res. 1990; 24:210–3.
4. . "R-R variation, a test of autonomic dysfunction" . Acta Neuro Scand. 1983; 67:285–93.
5. . "Heart rate variability in diabetes: relationship to age and duration of the disease" . Diabetes Care. 1985; 8:64–8.
6. . "Mortality in diabetic autonomic neuropathy" . Lancet. 1976; 1:601–3.
7. . "Cardiorespiratory arrest and diabetic autonomic neuropathy" . Lancet. 1978; 1:14–6.
8. . "The natural history of diabetic autonomic neuropathy" . Q J Med. 1980; 49:95–108.
9. . "Diabetic autonomic neuropathy and cardiovascular risk: Pittsburgh Epidemiology of Diabetes Complications Study" . Arch Intern Med. 1990; 150:1218–22.
10. . "Time and frequency-domain estimation of early diabetic cardiovascular autonomic neuropathy" . Clin Auto Res. 2001; 11:369–76.
11. . "Comparison of baroreflex sensitivity and heart period variability after myocardial infaction" . Journal of the Americal College of Cardiology. 1989; 14:1511–1518.
12. . "Evaluation of the Reynolds Pathfinder II system for 24 h heart rate variability analysis" . European Heart Journal. 1991; 12:1153–1162.
13. . "Twenty four hour heart rate variability: effects of posture, sleep and time of day in healthy controls and comparison with bedside tests of autonomic function in diabetic patient" . British Heart Journal. 1991; 65:239–244.



