Anti-Aging Research: Heart Rate Variability
----------------------------------------------------------------------------------------------------------------------------------------------------------
*** SPECIAL NOTE FROM DR. WICHMAN ***
For information regarding Heart Rate Variability (HRV) and Type 1 Diabetes, please read http://care.diabetesjournals
Also, the following excellent article is reproduced from ISRN Endocrinology at https://www.ncbi.nlm.nih.gov/p
Cardiac Autonomic Neuropathy Measured by Heart Rate Variability and Markers of Subclinical Atherosclerosis in Early Type 2 Diabetes
1. Introduction
Cardiac autonomic neuropathy (CAN) is a serious complication of type 2 diabetes mellitus (T2DM) which carries an approximately 5-fold increased risk of mortality in these patients [1, 2]. Damage to the autonomic innervations of the heart is a harbinger of dire outcomes such as lethal arrhythmias and sudden cardiac death [3].
Heart rate variability (HRV) analysis is a noninvasive tool to assess cardiac autonomic function. The value of HRV is that it can detect CAN before conventional tests of cardiovascular autonomic function like the Ewing battery [4–6].
In the Framingham Heart Study, both frequency and time domain measures of HRV were found to be inversely associated with the risk of mortality [7]. In the Atherosclerosis Risk In Communities (ARIC) Study, decreased HRV was independently associated with the risk of developing coronary heart disease [8, 9].
On the other hand, T2DM is associated with accelerated atherosclerosis, with markedly increased incidence of cardiovascular morbidity and mortality [10].
Thus, early detection of atherosclerosis is of utmost importance in T2DM to prevent negative outcomes.
Framingham Risk Score has been proved useful for detection of coronary heart disease risk before its serious negative end points become apparent [11, 12]. Notwithstanding, due to limitations of this scoring system surrogate markers of atherosclerosis have been devised to enhance prevention and earlier diagnosis of coronary heart disease in T2DM [13, 14].
We performed this study to evaluate whether CAN as measured by decreased HRV is associated with increased risk of subclinical atherosclerosis in T2DM.
2. Materials and Methods
2.1. Subjects
This case-control study was performed in Dr. Shariati hospital, Tehran University of Medical Sciences. A total of 57 diabetic and 54 nondiabetic control subjects were recruited in this study. Participants were drawn from the list of diabetics who attended diabetes outpatient clinics in Dr. Shariati Hospital, Tehran, Iran, for regular followup between September 2010 and January 2011. The control group was healthy relatives in law of diabetic participants. The study protocol was approved by the ethics committee of Tehran University of Medical Sciences, and the participants signed their informed consent at the time of recruitment.
The inclusion criteria were defined as both sexes, aged between 30 and 65 years. Current or previous smokers, subjects with diabetic foot, renal failure (GFR < 90), history of malignancy, and cirrhosis were excluded from the study. In addition, none of the participants had clinical coronary heart disease as evidenced by Rose questionnaire and electrocardiographic (ECG) criteria (Minnesota codes 1.1–1.3, 4.1–4.4, 5.1–5.3, and 7.1) at the time of the study.
Patients with proliferative retinopathy were excluded from the study according to their results of retinoscopy which was performed by two ophthalmologists. Peripheral neuropathy was diagnosed using the10-gram Semmes-Weinstein monofilament test on the plantar foot surface in diabetic patients. Failure to sense the filament was considered as diabetic neuropathy and resulted in exclusion from the study.
Height was measured with a Stadiometer, and weight was assessed by a calibrated beam balance. Body mass index (BMI) was calculated as weight (Kg) divided by height (M) squared. Blood pressure was measured twice (5 minutes apart) using a standard calibrated mercury sphygmomanometer on both right and left arms after the participants had been sitting calm for at least 10 minutes. The highest blood pressure of two sides was considered as participant's blood pressure. The two groups were matched based on age.
2.1.1. Definitions
Diabetes mellitus was defined as patients with FBS ≥ 126 mg/dL, or 2-h postload glucose ≥200 mg/dL or else were using oral hypoglycemic agents [14]. Hypertensive patients were those with systolic blood pressure of ≥ than 140 or diastolic blood pressure of ≥90 mmHg or were using antihypertensive drugs [15].
2.2. Laboratory Data
Venous blood samples were collected in the morning after 12 h fasting. The blood samples were centrifuged, and then serum was collected for measuring the biomedical parameters.
In all the participants, fasting plasma glucose (FBS) and 2-h postload glucose levels were measured.
Plasma levels of glucose, triglyceride (TG), total cholesterol, HDL cholesterol, LDL cholesterol, creatinine, blood urea nitrogen (BUN), and highly sensitive C-reactive Protein (hs-CRP) were measured by a colorimetric method using Pars Azmoon kit with an autoanalyzer (Hitachi 902, Boehringer Mannheim Germany). Serum Insulin concentration was assessed by immunoassay (ELISA) using a Bioscience kit (Monobind kit, USA). HbA1c was detected by High-performance liquid chromatography (HPLC) (Knauer, Germany), coupled with fluorescence detector. The method was validated over a linearity range of 1–100 μmol/L of the plasma. The intra- and interassay coefficients of variation (CVs) for all these measurements were <4%, which was less than allowed CVs. A morning clean catch midstream urine sample was collected for measuring creatinine and evaluation of microalbuminuria.
Homeostasis model assessment (HOMA) index was calculated as the product of the fasting plasma insulin level (μIU/mL) and the fasting plasma glucose level (mmol/L), divided by 22.5.
2.2.1. Surrogate Atherosclerosis Markers
We measured the following markers of subclinical atherosclerosis (SCA) in our participants.
Carotid Intima Media Thickness (CIMT) —
Ultrasonographic analysis of the carotid artery was performed with a high-resolution ultrasound scanner, equipped with a linear array 13 MHz transducer (MyLab 70 XVision, Biosound Esaote, USA). One physician who was blind to the nature of the group performed CIMT measurements. A rapid cross sectional scanning was made in the first step to pinpoint the possible plaques. The scan was started from the proximal part of the common carotid artery (CCA) toward the bifurcation, followed by scanning the internal and then the external carotid arteries. This process was followed by a longitudinal scanning of the CCA. In this step, the dynamic sequence images were stored for the following measurements of CIMT. A segment of the artery (usually where the vessel walls were most clearly seen throughout the recording) was magnified to identify a distinct lumen-intima and media-adventitia interface. CIMT was defined as the distance between the leading edge of the lumen-intima interface and the leading edge of the media-adventitia interface. For detection of CIMT, a special software (Vascular tools 5 (Medical Imaging Applications LLC, USA)) was employed. The regions of interest were defined as 1.0 cm distal to the bifurcation, the bifurcation and 1.0 cm proximal to the internal carotid artery in both near and far walls. Then for each subject CIMT was reported as the average of 12 measurements (6 measurements from the right and 6 from the left carotid artery).
Coronary Artery Calcium Score (CAC) —
Anterior-posterior and lateral chest scout views were first obtained for planning. Calcium score images were obtained with a Phillips 64 MDCT scanner using 64 × 2.5 mm × 400 ms with 120 KVP and 50–75 mAs to cover the entire heart and proximal ascending aorta. A positive calcium score was defined by 130 HU with an area of 1 mm2 or greater. The amount of calcium was quantified using the Agatston scoring method [16].
Flow-Mediated Dilation (FMD) —
FMD of the brachial artery was measured according to the American College of Cardiology guidelines [17]. The diameter of the right brachial artery was measured 3–5 cm above the antecubital fossa. Then a blood pressure cuff was inflated around the right forearm to at least 50 mmHg above the systemic blood pressure for 4-5 minutes. 60 seconds after cuff release, the diameter of the brachial artery was measured again. The brachial FMD was calculated as the percentage of change in the maximum postocclusion diameter of the brachial artery relative to the mean baseline diameter. All the measurements were performed in the end-diastolic phase coinciding with the R-wave on an electrocardiograph monitor. Every measurement was taken as the average of 3 consecutive cardiac cycles. FMD values greater than 5–10% were considered normal [17].
Heart Rate Variability —
The evaluation of HRV was performed in a quiet and temperature-controlled room according to the guidelines of the Task Force for Pacing and Electrophysiology [18]. Participants were advised to abstain from caffeinated food and beverages on the day of their assessments. Repeat assessments were performed at precisely the same time of day after 48 hours. After 15 minutes of supine rest with a regular and calm breathing pattern, a continuous 10-minute ECG recording was collected using an applanation tonometer interface with HRV software (SphygmoCor, AtCor Medical Pty, Sydney, Australia). The high-frequency (HF band: 0.15–0.45 Hz), low-frequency (LF band: 0.04–0.15 Hz), and very-low-frequency (VLF band: 0.01–0.04 Hz) components of HRV (measured in absolute units; i.e., ms2) were obtained. Total power (TP) of HRV was also calculated to be used in regression analysis as a global marker of cardiac autonomic function. Normalized HF and LF powers were determined by dividing their absolute powers by the total power minus the VLF component and multiplied by 100 [18–20]. From the electrocardiographic recording, the following statistical and geometric time domain indices were calculated from RR intervals: standard deviation of the NN intervals (SDNN), the square root of the mean squared difference of successive NNs (RMSSD), and the triangular index (TI). Frequency domain variables including total, HF, and LF powers and LF : HF ratio were derived from spectral analysis of successive R-R intervals [18].
Assessment of Cardiac Autonomic Function —
HRV measurement was performed after Valsalva and standing maneuvers in addition to supine state, using SphygmoCor software (SphygmoCor, AtCor Medical Pty, Sydney, Australia).
For valsalva maneuver, the participant was requested to blow into the mouthpiece of the device manometer to a pressure of 40 mmHg for 15 seconds. Then the valsalva ratio was calculated as the relationship between the longest and shortest R-R intervals after strain. For standing maneuver, the participant was requested to breathe at a normal pace for 5 minutes in the supine state. Then he/she was asked to go from supine to a full upright position and remain standing until the end of measurement. The standing ratio was calculated as longest R-R interval around the 30th beat after standing up to the shortest R-R interval around the 15th beat during standing.
We also performed deep breathing test by calculating the ratio of maximum and minimum heart rates during six cycles of paced deep breathing (E/I index) [21].
2.3. Data Analysis
For data analysis, the SPSS (version 18.0) was used. Normality of data distribution was evaluated by the Kolmogorov-Smirnov test. For data not normally distributed, the Mann-Whitney U test was applied. TP, LFnorm, HFnorm, LF : HF, RMSSD, and valsalva ratio were log-transformed before analysis due to their nonnormal distribution. For comparing data with normal distribution, unpaired t-test was used. Adjustment for confounding factors was performed using univariate analysis of variance (ANOVA). Correlation of variables was demonstrated using Pearson's and Spearman's correlation coefficients in normal distributed parametric and nonparametric variables, respectively. For assessment of association of variables, linear regression and logistic regression were used for parametric and binary variables, respectively.
3. Results
Table 1 summarizes the general characteristics and subclinical atherosclerosis markers of the two groups. Mean duration of diabetes was 8.6 years in T2DM participants. Systolic blood pressure was significantly higher and HDL-C levels were significantly lower in diabetics. Total- and LDL-cholesterols were significantly higher in diabetics. Fasting insulin and). HbA1c levels were significantly higher in diabetics than in controls. CIMT and CAC scores were significantly higher, while brachial FMD was significantly lower in diabetics.
Table 2 compares the autonomic and HRV indices between T2DM and control groups. SDNN and RMSSD were significantly lower in diabetics than in controls. Significant reduction of spectral power in HF band (expressed as normalized units) and in total power was also observed in our T2DM participants relative to controls.
Valsalva ratio, E/I index and heart rate response to standing were significantly lower in diabetics relative to controls.
HF power was highly correlated to total power. In addition, the observed associations for total power were similar to those of HF power. So we performed our analyses with total power for convenience.
Unadjusted correlations between total power and conventional CHD risk factors revealed that diastolic BP in all participants and BMI in diabetics were negatively associated with total power (Table 3).
Multiple logistic regression of the association between surrogate atherosclerosis markers and total power spectra revealed that CIMT was inversely and independently associated with total power both in diabetics and controls (Table 4). This relationship between total power and CIMT showed a dose-response pattern throughout the distribution of HRV. On the other hand, although CAC score also was inversely associated with HRV in a stepwise manner, this relation lost its significance after adjustment for diabetes (Table 2). The relationship between FMD% and total power was lost after multiple regression (Table 4).
Figure 1 shows the median values of the total power in normal participants (controls) and diabetics divided by tertiles according to HbA1c levels. We found an inverse association between total power and median HbA1c levels.
4. Discussion
We found that substantial autonomic dysfunction was present in our early T2DM patients even before overt clinical symptoms of CVD became apparent. It is noteworthy that involvement of the vagal parasympathetic component of autonomic nervous system was obvious in our T2DM patients. This is evidenced by increased resting heart rate and decreased Valsalva ratio; E/I index and standing ratio in diabetics relative to controls. These findings are in line with those of Freccero et al. who reported a high frequency of parasympathetic and sympathetic neuropathy in both type 1 and type 2 diabetic patients [22]. They suggested that severe damage to large myelinated nerve fibers in addition to the widespread neurological degeneration which usually affects the small nerve fibers of the autonomic nervous system was responsible for profound parasympathetic neuropathy in patients with T2DM. Other researchers also found appreciable degree of autonomic neuropathy in patients with T2DM [23, 24].
In fact significantly reduced HRV measures in T2DM patients compared to controls have been previously verified in large-population-based studies [25–27]. In the Hoorn Town of The Netherlands HF, and LF powers and SDNN were lower among T2DM participants compared to those with normal fasting glucose [25]. The same results had also been documented in the Framingham Heart Study and in the Atherosclerosis Risk in Communities (ARIC) cohort [26, 28].
The significance of our findings is that HRV reduction in diabetics was present since the early stages of diabetes even before clinical atherosclerotic cardiovascular disease became evident. Thus, it is imperative to screen for autonomic neuropathy as early as possible in T2DM to prevent or retard its serious consequences.
In addition, we found that surrogate atherosclerosis markers were associated with lower HRV. Especially, increased CIMT in our T2DM participants was significantly associated with decreased HRV. This association was independent of conventional risk factors such as hypertension, BMI, and dyslipidemia. Also, this association was of greater magnitude in diabetic than nondiabetic participants. A stepwise increase in the odds of CIMT with decreasing quartile of total power was observed in those with T2DM.
Gottsäter et al. in a longitudinal study of T2DM patients found a similar correlation between decreasing HRV in the form of total power and increasing CIMT [28]. In their study, this relation was observed both at baseline and with progression of diabetes during 3 years of followup. The association between decreased total power and E/I ratio with increased CIMT has also been observed by Eller et al. previously [29]. In addition, our findings are concordant with those of researchers in the ARIC cohort who had previously reported significant association between lower HRV and coronary heart disease in diabetic subjects [30].
Furthermore, increased CAC score was also associated with lower HRV in our T2DM participants; however, after adjustment for diabetes and other conventional risk factors this relationship lost its significance. Notwithstanding, Rodrigues et al. in their recent study of type 1 diabetes patients found that reduced HRV prospectively predicted progression of coronary artery calcium as a powerful cardiovascular disease risk marker [31]. Previously Colhoun et al. had reported clustering of HRV with coronary calcification and other cardiovascular risk factors in asymptomatic young type 1 diabetes patients [32].
We did not find an independent association between reduced HRV and endothelial dysfunction in this study after adjusting for conventional risk factors. However, Pinter et al. in a recent study of young healthy male volunteers found a significant positive correlation between vagal HRV indices (RMSSD, PNN50, and HF power) and normalized flow-mediated dilation [33]. They suggested that endothelial mediators especially nitric oxide released locally from the capillary endothelium could enhance the effects of vagal inputs during respiratory cycle. Vagal discharge increases during expiration and decreases during inspiration, producing respiratory sinus arrhythmia. As a result the major component of short-term HRV which is respiratory related is influenced by endothelial dysfunction.
The association of decreased HRV with atherosclerosis markers could be due to ischemic damage to cardiac nerves [34], which may have been occurred before these markers became apparent. In line with this hypothesis, Gautier et al. in their study of participants in the Pittsburg Study and Kuopio Ischemic Heart Disease Risk Factor Trial found a significant negative relation between mean CIMT and LF power. In other words, LF power, which reflects sympathovagal balance and baroreceptor buffering, was correlated with CIMT. They suggested that decreased LF power could be attributed to chronic downregulation of vascular sympathetic receptors over time. This may potentially cause vascular smooth muscle cell differentiation with subsequent intimal migration and extracellular matrix production [35]. Gottsäter et al. also concluded based on their findings that decreased LF power could predict the progression of atherosclerosis in T2DM [28]. These changes are characteristic of atherosclerotic progression. In our study although LF power was lower in diabetics than controls, this difference hardly became significant (P < 0.10). The current study cannot establish this array of events in the course of atherosclerosis development. Further experimental prospective investigations are needed to resolve this cause-and-effect relationship.
Advanced glycation end products may play a critical role in damage to cardiac nerves and subsequent autonomic dysfunction, which promotes progression of atherosclerosis. The inverse association between total power and median HbA1c levels in our study is concordant with this hypothesis. Colhoun and colleagues found that in young type 1 diabetic patients HbA1c was a very good proxy marker of diabetes which could explain differences in HRV indices between diabetics and nondiabetics [32]. On the other hand, Mäkimattila et al. found that total power was the factor most influenced by chronic hyperglycemia in type 1 diabetes patients. In addition, they showed that total power and retinopathy score were the most sensitive measures of diabetes complications [36]. In conclusion, we found that autonomic dysfunction, especially parasympathetic neuropathy, was present since early stages in T2DM. This was strongly related to subclinical atherosclerosis markers in these patients. Specifically decreased total power as a global measure of parasympathic neuropathy was independently associated with increased CIMT.
Limitations —
A limitation of this study is its cross-sectional nature. So it is not possible to establish a definite cause-and-effect relationship between findings as one is not able to differentiate whether decreased HRV precedes atherosclerosis markers or appears subsequently. In addition, the sample size was relatively limited, and we reported P values < 0.10.
Association of cardiac autonomic neuropathy with arterial stiffness in type 2 diabetes mellitus patients
- Ataollah Bagherzadeh,
- Afshin Nejati-Afkham,
- Yaser Tajallizade-Khoob,
- Akbar Shafiee,
- Farshad Sharifi,
- Morteza Abdar Esfahani,
- Zohre Badamchizade,
- Sudabeh Alatab and
- Hossein FakhrzadehThis email address is being protected from spambots. You need JavaScript enabled to view it.
DOI: 10.1186/2251-6581-12-55
© Bagherzadeh et al.; licensee BioMed Central Ltd. 2013
Received: 31 July 2013
Accepted: 30 September 2013
Published: 20 December 2013
Abstract
Background
Diabetic patients are at the risk of cardiac autonomic neuropathy (CAN) and arterial stiffness. This study aimed to investigate the association of heart rate variability (HRV) as an index for CAN and pulse wave velocity (PWV) as an index for arterial stiffness.
Methods
Uncomplicated diabetes type-2 patients who had no apparent history of cardiovascular condition underwent HRV and PWV measurements and the results were compared with the control group consisting of non-diabetic peers. Also, the findings were adjusted for the cardiovascular risk factors and other confounding factors.
Results
A total of 64 diabetic patients (age= 52.08±8.50 years; males=33 [51.6%]) were compared with 57 controls (age= 48.74±6.18 years; males=25 [43.9%]) in this study. Hypertension, dyslipidemia, and thereby systolic blood pressure and statin use were significantly more frequent in the diabetic group, while the serum levels of cholesterol, HDL-C and LDL-C were significantly higher in the controls. Pulse wave was significantly increased in the diabetic patients (p<0.001). Main HRV parameters were significantly lower in diabetics than in controls. After adjustment for the confounders, PWV and HRV remained significantly different between the groups (p=0.01 and p=0.004, respectively). Multiple logistic regression of the association between pulse wave velocity and HRV index was independently significant both in diabetics and controls.
Conclusions
There exists a significant relationship between heart rate variability and arterial stiffness as a measure for atherosclerosis in diabetic patients, although the role of the confounding factors is noteworthy.
Introduction
More than half of the mortality in diabetic patients occurs due to cardiovascular disease [1]. Although most of this risk stems from atherosclerosis along with vascular ageing, there are other conditions involved [2]. Cardiac autonomic neuropathy (CAN) is a serious complication of diabetes mellitus, shown to influence both mortality and cardiovascular events in these patients [3]. As the heart rate and vascular tone are regulated via the autonomic system, CAN can increase heart rate variability and decrease myocardial perfusion [4].
On the other hand, heart rate variability (HRV) can be assessed easily to measure CAN, much earlier than its clinical appearance [5]. It has been shown that low HRV is associated with increased mortality in patients with ischemic heart disease or diabetes mellitus [3, 6]. HRV is also a sensitive indicator of baroreflex control, specifically the vagal control [7, 8]. Therefore, arterial stiffness may affect baroreceptor function and thereby, HRV. Increased arterial stiffness evaluated by pulse wave velocity (PWV) and/or augmentation index has been associated with the presence of coronary atherosclerosis and worse cardiovascular prognosis both in general population [9, 10] and specific disease groups, including diabetes mellitus [11]. Therefore, understanding the exact relationship of heart rate variability, as an index for CAN, with arterial stiffness is crucial.
The aims of the present study were to investigate (i) the association of disturbed PWV, as an index for increased arterial stiffness, and HRV, as an index for CAN, in the presence or absence of uncomplicated diabetes mellitus and (ii) the relation of PWV with HRV in uncomplicated diabetic patients.
Methods
In this case–control study, uncomplicated diabetic type 2 (T2DM) patients who were referred to the Diabetes clinic of Dr. Shariati hospital, Tehran University of Medical Sciences, Tehran, Iran, between August 2011 and December 2012 were recruited. The control group consisted of their spouse or relatives who were proved to be non-diabetic. The inclusion criteria were defined as both sexes, aged between 30 and 65 years. Individuals with clinically proven coronary heart disease, diabetic foot, diabetic retinopathy or nephropathy, renal failure (GFR < 90), history of malignancy, and cirrhosis were excluded from the study. Height was measured with a Stadiometer, and weight was assessed by a calibrated beam balance. Body mass index (BMI) was calculated as weight (Kg) divided by height (M) squared. Blood pressure was measured twice (5 minutes apart) using a standard calibrated mercury sphygmomanometer on both right and left arms after the participants had been sitting calm for at least 10 minutes. The highest blood pressure of two sides was considered as participant’s blood pressure.
All the participants signed an informed consent at the time of recruitment and the study protocol was approved by the ethics committee of Tehran University of Medical Sciences, and the board of research at Dr. Sahriati hospoital.
Diabetes mellitus was defined as patients with fasting blood sugar (FBS) ≥ 126 mg/dL, or 2-h postprandial glucose ≥200 mg/dL or those who were using insulin or oral hypoglycemic agents. Hypertensive patients were those with systolic blood pressure of ≥than 140 or diastolic blood pressure of ≥90 mmHg or were using antihypertensive drugs.
Venous blood samples were drawn from the anticubital vein in the morning after 12 h fasting. The blood samples were centrifuged, and then serum was collected for measuring the biochemical parameters.
In all the participants, FBS and 2-h postprandial glucose levels were measured. Plasma levels of glucose, triglyceride (TG), total cholesterol, high density lipoprotein cholesterol (HDL-C), low density lipoprotein cholesterol (LDL-C), creatinine, and blood urea nitrogen (BUN) were measured by a colorimetric method using Pars Azmoon® kit with an autoanalyzer (Hitachi 902, Boehringer Mannheim Germany). HbA 1 C was detected by High-performance liquid chromatography (HPLC) (Knauer, Germany), coupled with fluorescence detector. The method was validated over a linearity range of 1–100 μmol/L of the plasma. The intra- and interassay coefficients of variation (CVs) for all of the measurements were <4%, which was less than allowed CVs.
Evaluation of HRV was performed in a quiet and temperature-controlled room according to the guidelines of the Task Force for Pacing and Electrophysiology [5]. Participants were advised to abstain from caffeinated food and beverages on the day of assessment. Repeat assessments were performed at precisely the same time of day after 48 hours. After 15 minutes of supine rest with a regular and calm breathing pattern, a continuous 10-minute ECG recording was collected using an applanation tonometer interface with HRV software (SphygmoCor, AtCor Medical Pty, Sydney, Australia). The high-frequency (HF band: 0.15–0.45 Hz), low-frequency (LF band: 0.04– 0.15 Hz), and very-low-frequency (VLF band: 0.01–0.04 Hz) components of HRV (measured in absolute units; i.e., ms 2 ) were obtained. Total power (TP) of HRV was also calculated to be used in regression analysis as a global marker of cardiac autonomic function.
Normalized HF and LF powers were determined by dividing their absolute powers by the total power minus the VLF component and multiplied by 100 [5, 12].
From the electrocardiographic recording, the following statistical and geometric time domain indices were calculated from RR intervals: standard deviation of the NN intervals (SDNN), and the square root of the mean squared difference of successive NNs (RMSSD). Frequency domain variables including total, HF, and LF powers and LF:HF ratio were derived from spectral analysis of successive R-R intervals [5].
HRV measurement was performed after Valsalva and standing maneuvers in addition to supine state, using SphygmoCor software (SphygmoCor, AtCor Medical Pty, Sydney, Australia). For Valsalva maneuver, the participant was requested to blow into the mouthpiece of the device manometer to a pressure of 40 mmHg for 15 seconds. Then the Valsalva ratio was calculated as the relationship between the longest and shortest R-R intervals after strain. For the standing maneuver, the participant was requested to breathe at a normal pace for 5 minutes in the supine state. Then, the participant was asked to change position from supine to a full upright and remain erected until the end of the test. The standing ratio was calculated as longest R-R interval around the 30th beat after standing up to the shortest R-R interval around the 15th beat during standing.
PWV was measured using the SphygmoCor System (AtCor Medical Pty Ltd Head Office, West Ryde, Australia), with the individual in the supine position. The pulse waves of the carotid and femoral arteries were analyzed, estimating the delay according to the ECG wave and calculating PWV. PWV was calculated as the ratio of the distance travelled (calculated as distance in mm of distal minus proximal, where measurements are performed from the supra-sternal notch to the sampling site) and the foot-to-foot time delay between the pulse waves and expressed in meters per second (m/sec).
Statistical analysis
Continuous data are presented as mean ± SD. For comparing data with normal distribution, unpaired t-test was used. Correlation of variables was demonstrated using Pearson’s and Spearman’s correlation coefficients in normal distributed parametric and nonparametric variables, respectively. For assessment of the association of variables, linear regression and logistic regression were used for parametric and binary variables, respectively. P-values were always 2-sided, and P < 0.05 was considered significant. The SPSS statistical software package (version 18.0 for Windows; SPSS Inc. Chicago, IL) was used for data analysis.
Results
General characteristics of the study population
|
Characteristics |
Normal (N=57) |
Diabetic (N=64) |
P-value |
|---|---|---|---|
|
Age (years) |
48.74 ± 6.18 |
52.08 ± 8.50 |
0.15 |
|
Sex (male) % |
25 (43.9) |
33 (51.6) |
0.39 |
|
BMI (kg/m2) |
29.12 ± 5.06 |
27.93 ± 4.46 |
0.12 |
|
Hypertension |
4 (7.0) |
23 (35.9) |
<0.001 |
|
Dyslipidemia |
9 (15.8) |
36 (57.1) |
<0.001 |
|
Smoking, n (%) |
1 (1.8) |
8 (12.5) |
0.02 |
|
Statin use |
2 (3.5) |
27 (42.9) |
<0.001 |
|
FBS (mg/dl) |
93.84 ± 13.22 |
162.60 ± 56.32 |
<0.001 |
|
Hb A1c (%) |
5.31 ± 0.68 |
7.95 ± 1.70 |
<0.001 |
|
Cholesterol (mg) (mg/dl) |
201.80 ± 31.04 |
174.12 ± 38.50 |
<0.001 |
|
Triglyceride (mg/dl) (mg/dl) |
173.52 ± 99.17 |
194.08 ± 120.45 |
0.42 |
|
HDL-C (mg/dl) |
46.25 ± 10.74 |
40.34 ± 8.23 |
0.001 |
|
LDL-C (mg/dl) |
114.52 ± 22.43 |
95.12 ± 23.13 |
<0.001 |
|
Creatinine (mg/dl) |
0.95 ± 0.13 |
0.94 ± 0.15 |
0.62 |
|
BUN |
12.88 ± 3.19 |
12.36 ± 4.15 |
0.45 |
|
Hemoglobulin |
14.87 ± 4.37 |
14.30 ± 1.44 |
0.34 |
|
SBP (mmHg) |
123.46 ± 14.09 |
131.00 ± 17.33 |
0.01 |
|
DBP (mmHg) |
77.47 ± 10.11 |
76.03 ± 8.73 |
0.4 |
Comparing pulse wave velocity and heart rate variability indices between the study groups
|
Variables |
Normal (N=57) |
Diabetic (N=64) |
P-value |
|---|---|---|---|
|
Pulse wave velocity (m/s) |
8.00 ± 1.61 |
10.11 ± 2.45 |
<0.001 |
|
RMSSD (ms) |
31.78 ± 21.27 |
20.13 ± 19.51 |
0.02 |
|
HRV index |
7.78 ± 2.94 |
5.57 ± 2.22 |
<0.001 |
|
Heart rate(bpm) |
68.95 ± 9.44 |
73.64 ± 8.41 |
0.005 |
|
PNN 50 |
9.58 ± 14.56 |
3.01 ± 7.10 |
0.002 |
|
LF norm (ms2) |
55.98 ± 18.79 |
61.65 ± 20.99 |
0.12 |
|
HF norm (ms2) |
44.02 ± 18.78 |
38.34 ± 20.99 |
0.12 |
|
LF:HF ratio |
1.83 ± 1.70 |
2.72 ± 2.73 |
0.03 |
|
Total power (ms2) |
1133.17 ± 1266.26 |
703.73 ± 1279.76 |
0.06 |
|
Valsalva ratio |
1.58 ± 0.28 |
1.44 ± 0.31 |
0.01 |
|
Standing ratio |
1.27 ± 0.13 |
1.23 ± 0.20 |
0.29 |
After adjustment for the confounding variables, including age, BMI, hypertension, dyslipidemia, smoking, statin use, FBS, total cholesterol and systolic blood pressure, PWV and HRV remained significantly different between the groups (p=0.01 and p=0.004, respectively).
Partial and total correlations of the total power with PWV and other cardiovascular parameters
|
Characteristics |
Nondiabetic subjects |
Diabetic subjects |
All subjects |
|---|---|---|---|
|
Age (years) |
−0.42 † |
−0.16 |
−0.14 |
|
Sex (male) % |
0.17 |
0.06 |
−0.07 |
|
BMI (kg/m2) |
0.07 |
−0.23 |
−0.06 |
|
Hypertension |
0.2 |
0.28 * |
0.15 |
|
Dyslipidemia |
0.14 |
0.19 |
0.2 * |
|
Current smoker (%) |
0.02 |
−0.17 |
−0.17 |
|
Statin use |
−0.11 |
−0.14 |
−0.12 |
|
FBS (mg/dl) |
−0.23 |
−0.003 |
−0.02 |
|
Hb A1c (%) |
0.004 |
−0.09 |
0.09 |
|
Cholesterol (mg) (mg/dl) |
−0.04 |
0.04 |
−0.03 |
|
Triglyceride (mg/dl) (mg/dl) |
−0.02 |
0.25 * |
−0.08 |
|
HDL-C (mg/dl) |
−0.13 |
−0.02 |
−0.10 |
|
LDL-C (mg/dl) |
−0.04 |
−0.10 |
−0.001 |
|
Creatinine (mg/dl) |
−0.10 |
0.18 |
0.06 |
|
BUN |
−0.06 |
−0.05 |
−0.12 |
|
Hemoglobulin |
−0.20 |
0.21 |
0.02 |
|
SBP (mmHg) |
−0.10 |
−0.31 * |
−0.12 |
|
DBP (mmHg) |
−0.35 † |
−0.36 † |
−0.32 † |
|
Pulse wave velocity (m/s) |
−0.28 * |
0.37 † |
−0.19 * |
Association between pulse wave velocity and quartiles of HRV index in multivariate logistic regression model
|
Unadjusted OR (CI) |
Diabetes adjusted OR (CI) |
Full adjusted OR (CI) | |
|---|---|---|---|
|
Fist quartile of HRV index† |
|||
|
Second quartile of total power |
0.73 (0.57, 0.92) |
0.81 (0.63, 1.04) |
0.93 (0.81, 1.06) |
|
Third quartile of total power |
0.51 (0.38, 0.71) |
0.57 (0.41, 0.80) |
0.55 (0.37, 0.84) |
|
Fourth quartile of total power |
0.55 (0.40, 0.74) |
0.64 (0.47, 0.89) |
0.91 (0.61, 1.37) |
|
P trends |
<0.001 |
0.001 |
0.01 |
Discussion
In this study, we observed that considerable cardiac autonomic dysfunction exists in the uncomplicated diabetic patients in comparison with the normal controls, as well as an increase in the arterial stiffness, measured by pulse wave velocity.
Current evidence confirms that HRV is a good measure of cardiac autonomic neuropathy in diabetic patients and its decrease is accompanied by increased mortality and morbidity [13]. In our study, increased resting heart rate and decreased Valsalva ratio and standing ratio in the diabetic patients illustrate the parasympathetic involvement of the autonomous system as compared with the normal controls. This has also been shown previously both in type 1 and type 2 diabetes patients [14, 15].
Decreased HRV in the uncomplicated diabetes patients of our study and previous studies highlights the obscure process of autonomic neuropathy in diabetic patients that begins even before clinical atherosclerotic cardiovascular disease becomes apparent [16].
It has also been shown that surrogate atherosclerosis markers were associated with lower HRV, and increased carotid intima media thickness (CIMT) in T2DM participants was significantly associated with decreased HRV, independent from conventional cardiovascular risk factors [16]. Therefore, the presence of cardiac autonomic neuropathy should be considered much earlier in the course of diabetes, rather than after the development of clinical cardiovascular disease.
PWV is known as a potentially applicable atherosclerotic risk marker irrespective of classical cardiovascular risk factors and ethnicity [17]. In previous studies, arterial stiffness assessed by pulse wave analysis had a prognostic value for cardiovascular morbidity and mortality, mostly in hypertensive patients [18, 19, 20]. One study demonstrated that increased aortic pulse wave velocity was associated with the presence of angiographic coronary artery disease in overweight and obese patients, although the arterial stiffness indices were not consistently associated with obesity [21]. Similarly, it has been shown that cardiac parasympathetic function is a strong predictor of large arterial stiffness, in young type 1 diabetes patients without macrovascular and renal complications [22]. In the Pittsburgh Epidemiology of Diabetes Complications study, cardiac autonomic neuropathy was associated with increased arterial stiffness indices, in patients with childhood-diagnosed type 1 diabetes [23]. Moreover, a novel relationship between arterial stiffness, hyperinsulinaemia and autonomic neuropathy in a Type 2 diabetic population has been shown in a study which signifies their pathogenic roles in the development of cardiovascular disease in diabetic patients [13].
Based on our findings and previous works, one could suggest that atherosclerosis, both as a result of diabetes and increased age is influenced by the cardiac autonomic neuropathy, which in turn results in increased risk of cardiovascular diseases and related mortality in the type 2 diabetes patients.
Study limitations
Among limitation to the study, we can mention that our measurements were cross-sectional and we could not assess how the changes in both HRV and PWV through time, as well as the fluctuations in the serum glucose levels and glycosated hemoglobin, affect each other. We also did not perform glucose tolerance test in the normal controls, so there is a probability that those with glucose intolerance may have some degrees of disturbed HRV and PWV. Also the effects of the type of diabetes treatment and other prescribed medications on HRV and PWV need to be investigated in future studies.
Conclusion
In this study, we observed increased arterial stiffness and decreased heart rate variability in the uncomplicated type 2 diabetes patients as compared with normal controls. The relationship between heart rate variability indices and pulse wave velocity was significant after adjustment for diabetes; however, this effect was lost after adjustment for confounders. Based on the findings of this study, it seems that there exists a relationship between heart rate variability and arterial stiffness as a measure for atherosclerosis in diabetic patients, although the role of the confounding factors should be taken into account.
Copyright
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
----------------------------------------------------------------------------------------------------------------------------------------------------------
Heart Rate Variability (HRV) is the measure of the variability between heartbeats and the health of the autonomic nervous system. Imbalances between the parasympathetic nervous system (maintains balance body functions) and sympathetic nervous system ( which controls emergency energy in stressful situations) is an indicator of health issues such as heart arrhythmia or diabetes.
Also denoted as the RR interval, an ECG can measure the patient’s HRV against the patient’s resting heart rate. Low HRV rates are a key indicators of heart disease, stress and the chance of heart attack survivability. Using this data, a doctor can diagnose existing issues and proscribe a treatment in response.
Sources and Research:
Heart Rate Variability: Standards of Measurement, Physiological Interpretation, and Clinical Use
