Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Association of Diabetes Related Complications with Heart Rate Variability among a Diabetic Population in the UAE

  • Ahsan H. Khandoker,

    Affiliation Department of Biomedical Engineering, Khalifa University, Abu Dhabi, United Arab Emirates

  • Haitham M. Al-Angari,

    Affiliation Khalifa University Center of Biotechnology, Abu Dhabi, United Arab Emirates

  • Kinda Khalaf,

    Affiliation Department of Biomedical Engineering, Khalifa University, Abu Dhabi, United Arab Emirates

  • Sungmun Lee,

    Affiliation Department of Biomedical Engineering, Khalifa University, Abu Dhabi, United Arab Emirates

  • Wael Almahmeed,

    Affiliations Institute of Cardiac Science, Sheikh Khalifa Medical City, Abu Dhabi, United Arab Emirates, Heart and Vascular Institute, Cleveland Clinic, Abu Dhabi, United Arab Emirates

  • Habiba S. Al Safar,

    Affiliations Department of Biomedical Engineering, Khalifa University, Abu Dhabi, United Arab Emirates, Khalifa University Center of Biotechnology, Abu Dhabi, United Arab Emirates

  • Herbert F. Jelinek

    hjelinek@csu.edu.au

    Affiliations School of Community Health, Charles Sturt University, Albury, New South Wales, Australia, Australian School of Advanced Medicine, Macquarie University, Sydney, New South Wales, Australia

Abstract

Microvascular, macrovascular and neurological complications are the key causes of morbidity and mortality among type II diabetes mellitus (T2DM) patients. The aim of this study was to investigate the alterations of cardiac autonomic function of diabetic patients in relation to three types of diabetes-related complications. ECG recordings were collected and analyzed from 169 T2DM patients in supine position who were diagnosed with nephropathy (n = 55), peripheral neuropathy (n = 64) and retinopathy (n = 106) at two hospitals in the UAE. Comparison between combinations of patients with complications and a control diabetic group (CONT) with no complication (n = 34) was performed using time, frequency and multi-lag entropy measures of heart rate variability (HRV). The results show that these measures decreased significantly (p<0.05) depending on the presence and type of diabetic complications. Entropy, (median, 1st- 3rd interquartile range) for the group combining all complications (1.74,1.37–2.09) was significantly lower than the corresponding values for the CONT group (1.77, 1.39–2.24) with lag-1 for sequential beat-to-beat changes. Odds ratios (OR) from the entropy analysis further demonstrated a significantly higher association with the combination of retinopathy and peripheral neuropathy versus CONT (OR: 1.42 at lag 8) and an even OR for the combination of retinopathy and nephropathy (OR: 2.46 at lag 8) compared to the other groups with complications. Also, the OR of low frequency power to high frequency power ratio (LF/HF) showed a higher association with these diabetic-related complications compared to CONT, especially for the patient group combining all complications (OR: 4.92). This study confirms that the type of microvascular or peripheral neuropathy complication present in T2DM patients have different effects on heart rate entropy, implying disorders of multi-organ connectivity are directly associated with autonomic nervous system dysfunction. Clinical practice may benefit from including multi-lag entropy for cardiac rhythm analysis in conjunction with traditional screening methods in patients with diabetic complications to ensure better preventive and treatment outcomes in the Emirati Arab population.

Introduction

Type II diabetes mellitus (T2DM) is a multifactorial disease characterized by altered glucose metabolism that can affect organ function either directly or indirectly through oxidative stress and inflammatory mechanisms linked to hyperglycemia [1]. Oxidative stress and inflammation are considered as the main causative factors of diabetes complications including renal, retinal vascular and neurological impairments [2, 3]. However, the mechanisms for the development of diabetic related complications are most likely multifactorial and involve in addition to oxidative stress and inflammation, environmental and lifestyle factors, as well as genetic predisposition [4].

T2DM complications do not always occur in isolation but are often found as a group in patients, especially in those who do not have good glucose control, and most often involve autonomic dysfunction. All major organs of the body are modulated by the autonomic nervous system (ANS) and ANS dysfunction is a major cause of increased morbidity and mortality in T2DM. Approximately half of the nephropathy patients attending secondary and tertiary care experience autonomic impairment and cardiopathy [5]. Similarly, diabetic retinopathy has been associated with cardiac autonomic dysfunction in both type 1 [6] and type 2 diabetes mellitus patients [7]. Diabetic peripheral neuropathy (DPN) and peripheral vascular disease have also been associated with symptoms of distal sympathetic autonomic neuropathy [8].

The sympathetic and parasympathetic parts of the ANS are both affected by increased blood glucose levels, oxidative stress, and inflammation processes, potentially leading to multiple organ dysfunction and cardiac autonomic neuropathy (CAN) [9]. CAN is characterized by an altered cardiac rhythm due to the initial changes in the parasympathetic followed by sympathetic modulation of the cardiac rhythm [1012]. Meta-analyses of published data have demonstrated that reduced cardiovascular autonomic function as measured by heart rate variability (HRV) is strongly (i.e., relative risk is doubled) associated with an increased risk of silent myocardial ischemia and mortality in diabetes patients [13]. However, up to date there are no reports on how single or multiple diabetic complications (DPN, nephropathy and retinopathy) affect cardiac autonomic dysfunction.

HRV analysis is a noninvasive method to assess variations in autonomic nervous system modulation of the heart and has been used previously to evaluate impairments in the sympathetic and vagal tone associated with diabetes, cardiovascular disease, depression, schizophrenia and Parkinson’s disease [1418].

Our hypothesis is that multi-lag HRV parameters are associated with different combinations of T2DM related complications. Multi-lag heart rate analysis has the advantage of identifying possible contributions of the parasympathetic and sympathetic components of the ANS and therefore may be more sensitive to changes in the ANS associated with different complications. For example, the parasympathetic nervous system influences the heart rate on a beat-by-beat basis up to 10 beats in sequence. Sympathetically driven changes in the heart rate are slower and occur most likely at longer time lags. The aim of this study was twofold: 1) to investigate whether the various diabetic complications and combinations affect HRV and entropy at multiple lags (1–8); and 2) to study if the observed HRV results differ depending on the type of diabetic complication and the combinations of these complications present in the current group of patients. Our findings have confirmed that the types of microvascular complications present in T2DM patients lead to differences in heart rate variability and multi-lag entropy results at lower to higher lags, which suggests that this method may be a useful adjunct to clinical decision making for treatment options related to diabetes complications.

Materials and Methods

A total of 169 (119 females and 50 males), unrelated diabetic patients who were enrolled during a routine visit to the endocrinology and cardiology clinics at Shaikh Khalifa Medical Center (SKMC) and Mafraq Hospital in Abu Dhabi, United Arab Emirates participated in the study during the period between July 2014 and May 2015. Each volunteer agreed to take part in this study after a briefing session and upon, signing an informed consent form that had been approved by the Institutional Ethics Committee of both hospitals (REC-04062014 and R292 respectively). Inclusion criteria were UAE-born national, diagnosed diabetic patient at one of the hospitals with and without complications, able to give consent, does not have coronary artery disease and above 18 years of age. Exclusion criteria were, current pregnancy, and other pathophysiology such as cancer.

The patients included in the study were all diagnosed with T2DM either from the medical records or reported medication used at the time of the hospital visit. The presence of diabetic associated complications was confirmed by a qualified physician, based on the criteria outlined by the World Health Organization (WHO) consultation group report [19].

The patients in this study were categorized into 8 groups: 1) ALL-C group which included patients with retinopathy, nephropathy and DPN, 2) RNp: included patients with retinopathy and nephropathy, 3) NNp: patients with DPN and nephropathy, 4) RN: patients with retinopathy and DPN, DPNn: patients with DPN, 6) R: patient with retinopathy, 7) Np: patients with nephropathy, and 8) a control group (CONT), which included patients with none of these three complications. A patient was diagnosed with DPN if they presented with (1) foot ulcers, (2) loss of sensation/numbness/burning/tingling in the feet, (3) loss of toe, foot or leg due to diabetes, (4) pain in calf muscles while walking, or (5) peripheral vascular disease in the legs [20]. The presence of nephropathy was determined by urine albumin level higher than 20 μg/min for microalbuminuria and higher than 200 μg/min for macroalbuminuria or if the estimated glomerular filtration rate (eGFR) was less than 60 ml/min/1.73m2. Retinopathy was defined as either white or red lesions or both present in the retina according to WHO criteria [21, 22].

2.1 Data Collection

Demographic data collection and clinical assessment for each participant were completed at the clinics. ECG signals were recorded for 10 minutes under supine rest and edited using the MLS310 HRV module (version 1.0, ADInstruments, Australia) included in the Chart software package. High frequency noise was removed with a 45Hz low-pass filter and a 0.5 Hz high pass filter adjusted for wandering baseline. Ectopic beats were selected visually and deleted manually. Linear interpolation was used to replace ectopic beats that occur immediately before and after the ectopic interval. Intervals between successive R waves of the QRS complex (i.e., RR intervals in seconds) were calculated using the algorithm developed by Pan and Tompkins [23]. The HRV analysis described in the following sections was performed on 500 RR intervals selected from the middle part of the 10 minute recording [24].

2.2 Time and Frequency Domain Analysis

The mean RR (mRR), the standard deviation of normal RR intervals (SDNN), the standard descriptors of the Poincaré plot (SD1 and SD2) [25] were computed for time domain analysis. Spectral analysis was performed using the fast Fourier transform on linearly resampled (1 Hz) time series using Welch’s method with a 50% overlap between adjacent segments [26]. A Hanning window was applied to avoid spectral leakage. Subsequently, spectral powers were binned in the low frequency (LF) band (0.04–0.15 Hz), the high frequency (HF) band (0.15–0.40 Hz). From these bands HF, LF and the LF/HF were determined [27].

2.3 Multi-lag Entropy Analysis

In entropy analysis, the RR interval lengths associated with the continuous heartbeat are recorded as a percentile change. The entropy can then be determined from the probability distribution of the percentile change of the RR intervals using information theory such as the information content suggested by Shannon [28].

In conventional entropy analysis, the percentile change of the successive RR intervals with respect to the first RR interval is expressed as the percentage index (PI) defined as: (1)

Entropy is defined from the probability distribution of PI using Shannon’s formula: (2) where p(i) is a probability of PI(n) having values in the range i ≤ PI(n)≤ i+1, where i is an integer. Multi-lag entropy differs as it analyses the probability changes of the RR intervals within a set of predefined heart beats, which in the current analysis is set between 1 and 10 beats. For multi-lag entropy analysis, a lag of m points was introduced in Eq (1). Hence, in the multi-lag entropy analysis, PI is expressed as the percentile change of the i-th and i+m-th RR intervals with respect to the i-th RR interval and is defined as: (3) Where m is an integer and m = 1 represents the conventional entropy analysis. More details on the conventional entropy method can be found in previous work [29, 30].

An example of the derived probability of RR intervals of certain lengths occurring in the heart rate tachogram and the derived histogram (for the multi-lag entropy computation) from the CONT and DPNn group is shown in Fig 1.

thumbnail
Fig 1. Deriving probability and histogram from RR intervals.

RR interval data from a CONT patient (a) and a DPNn patient (b), their PI (c,d) and PI histogram (e,f).

https://doi.org/10.1371/journal.pone.0168584.g001

2.4 Statistical Analysis

The Lilliefors test of normality was performed on all features to test for the normal distribution of data. A Kruskal-Wallis nonparametric test of variance was applied for comparing the CONT group and each of the other seven groups mentioned above. We used a parametric analysis of covariance (ANCOVA) to correct for the age and years of diabetes and test feature significance between the CONT group and each of the other groups with complications. A nonparametric ANCOVA (Quade test) was used when the assumptions for parametric ANCOVA were violated [31, 32] Spearman’s Rank-Order Correlation was performed between the HRV features and the demographic data. A binary logistic regression model was used to check association with the HRV features for each of the complication groups compared to CONT and odds ratio (OR) for each of the HRV features was determined. MatLab (R2014b) statistical toolbox was used to perform the statistical analysis.

Results

The data set used in this study included 135 T2DM patients who were diagnosed with retinopathy, DPN, nephropathy, or combinations of these three complications as well as 34 T2DM patients with no complications. Table 1 summarizes the clinical variables of the patient groups in this study.

Significant differences in the demographic variables (using Kruskal-Wallis test) were observed for age, diabetes mellitus years (DM years), total cholesterol level, high and low density lipoprotein between the groups with different diabetic complications. Age was significantly lower in the CONT group compared to the remaining groups except for N and NNp, which also had a lower mean age but were not significant. DM years was significantly lower in the CONT group as compared to all groups-except for the Np group. Cholesterol levels differed significantly between the control and the diabetic complications and combinations (Table 1). Medication information is categorized into 4 groups: a) Combination drug use (Biguanides, dipeptidyl peptidase-4 inhibitors, sulfonureas and sodium glucose transport inhibitors), b) Incretin mimetics, c) Insulin analogs, d) Antianxiety, antidepressive and antipsychotic medications (Table 1). The values in the table indicate the number of patients who were using that type of medication in each group.

Without correcting for confounding effects, only the RN group entropy results with lags 5–8 were significantly lower than CONT (p<0.05). No other comparisons were significant. To correct for any possible confounding effects we applied ANCOVA and included age and DM years as covariates in the model and show the results in Table 2. Comparisons between the HRV time, frequency and entropy features for the CONT group and multiple combinations of diabetes complications are shown in Table 2. SD1 was significantly lower in the groups that included nephropathy patients as compared to CONT. In contrast, SD2 was significantly lower for the groups that presented with retinopathy and combined with either nephropathy or DPN. DPN as a single complication was not significantly different from the CONT group for any HRV features. However in combination with retinopathy or nephropathy, significantly lower HRV results were observed.

thumbnail
Table 2. Results for 8 T2DM groups using time, frequency and multi-lag entropy features.

https://doi.org/10.1371/journal.pone.0168584.t002

Groups that included retinopathy patients showed a significant difference in LF and LF/HF as compared to CONT. Overall, spectral indices were lower than those found for the ALL-C group. Nephropathy as a single complication had a significantly lower LF, while no HRV feature differed significantly between DPN and the control group. LF was lower for all groups investigated as compared to the control and was significant for five of the seven groups investigated.

For the multi-lag entropy analysis, significant differences were only obtained when the groups with multiple diabetic complications (ALL-C, RNp, NNp, and RN) were compared to the CONT group. However, the RNp group showed significant difference from the CONT at lags higher than 2. Fig 2 shows the multi-lag entropy for the 8 T2DM groups.

thumbnail
Fig 2. Multi-lag entropy of the 8 T2DM groups for lags 1–8.

https://doi.org/10.1371/journal.pone.0168584.g002

In addition, we investigated whether any of the HRV measures were correlated with demographic variables using Spearman’s Rank-Order correlation (Tables 3, 4 and 5). Our results indicate that HRV features were negatively correlated with age apart from the ALL-C HRV features, which were positively correlated with age (Table 3). The highest correlations were found for the Np group (-0.65 to -0.85, p<0.05) with entropy lags above 2. Significant correlations between age and the multi-lag entropy (at all lag) were also found for the CONT, and RNp groups. An example for the Spearman’s Rank-Order correlation for multi-lag entropy (lag 1) with age for the eight patient groups is shown in Fig 3. A significant correlation with age was observed for patients with all three complications present (Fig 3a), retinopathy and nephropathy (Fig 3c) and nephropathy (Fig 3h).

thumbnail
Table 3. HRV feature correlations with age for the 8 T2DM groups.

https://doi.org/10.1371/journal.pone.0168584.t003

thumbnail
Table 4. HRV feature correlations with DM years for the 8 T2DM groups.

https://doi.org/10.1371/journal.pone.0168584.t004

thumbnail
Table 5. HRV feature correlations with systolic blood pressure, high density lipoprotein level and total cholesterol level for CONT group.

https://doi.org/10.1371/journal.pone.0168584.t005

thumbnail
Fig 3. Correlation between age and entropy at lag 1 for the 8 groups.

a) CONT, b) ALL-C, c) RNp, d) NNp, e) RN, f) DPNn, g) R, n) Np. Only CONT, RN and Np groups showed significant correlation (p<0.05).

https://doi.org/10.1371/journal.pone.0168584.g003

Years of diabetes did not significantly correlate with all HRV features (Table 4). Only the entropy features from lag 2 to lag 7 and presence of retinopathy and DPN had any degree of correlation with years of diabetes. The highest significant correlation was found between the presence of DPN and entropy at lag 4 (-0.62) followed by a correlation of -0.43 for retinopathy and entropy at lag 6. Some of the clinical variables and HRV features were also significantly correlated but only in the CONT group (Table 5).

High density lipoprotein (HDL) showed the highest correlation with LF power. All correlations were negative and in agreement with clinical observations, where abnormally high systolic blood pressure or cholesterol levels affects heart rhythm negatively.

The OR in Table 6 indicate that the time and frequency features had an OR of 1, thereby suggesting no association or increased risk of diabetes complications, except for the LF/HF feature, which had OR values between 1.33 and 4.92. The highest OR was found for the ALL-C group. In contrast the OR for the multi-lag entropy had values less than 1 except for RN and entropy lags between 4 and 8. OR between 0.4 and 0.7 were obtained for the R group and between 0.2 and 0.4 for the DPNn group. The combination of retinopathy and nephropathy (RNp group) lowered the OR (0.15–0.4), whereas the combination of all three complications again lowered the OR (0.5–0.7). The OR between the RN group and the other complication groups (Table 7) shows that the differences in ANS dysfunction between the retinopathy-DPN and retinopathy-nephropathy complications are manifested more on the multi-lag entropy measure than the LF/HF ratio. To study the influence of gender on the results we stratified data from the CONT and the RN group (the group that showed most significant difference from CONT) by sex. The male groups (CONT: 12, RN: 6 subjects) had matched age and DM years while the female CONT group had significantly lower DM years, p<0.01 (CONT: 22, RN: 21 subjects). However, no significant differences were found in any of the HRV features between the male and female groups.

thumbnail
Table 6. Odds ratio of the HRV features between the CONT and DM with complications groups.

https://doi.org/10.1371/journal.pone.0168584.t006

thumbnail
Table 7. Odds ratio of the HRV features between the RN group and the other DM with complications groups.

https://doi.org/10.1371/journal.pone.0168584.t007

Discussion

Clinical studies have shown that development of T2DM complications is definitely linked to hyperglycemia, but the pathophysiology of retinopathy, nephropathy and DPN are different [4, 33]. DPN which occurs in up to 60% of long-term T2DM, can lead to cardiac arrhythmias and sudden cardiac death [34]. Similarly diabetic retinopathy and nephropathy have been linked with cardiovascular events [35, 36]. Our current study investigated the relationship between DPN, diabetic retinopathy and nephropathy specific T2DM complications and combinations of these complications with the cardiac autonomic nervous system by using traditional and multi-lag HRV analysis in a previously not studied Arab Emirati T2DM cohort. The prevalence of DPN in this study was 38%, which is comparable with data from other T2DM populations from Saudi Arabia (38%) [37], Mexico (40%) [38], UK (36%) [39] and Egypt (22%) [40]. Nephropathy population in this study (33%) was also comparable to other studies from Jordan (33%) [41], Egypt (42%) [40] and Saudi Arabia (40%) [42]. Retinopathy prevalence (63%) was higher than reported results from Saudi Arabia (31%) [43], Oman (14.4%) [44] and Yemen (55%) [40].

The pathophysiology of diabetic complications also includes dysfunction of the ANS, which also affects cardiac function and is the main cause of cardiovascular morbidity and mortality in T2DM [34]. Diabetic associated complications are typically considered to be a function of disease duration, which is reflected in ANS dysfunction, where the parasympathetic division is usually affected earlier compared to sympathetic dysfunction. However the role of age, years of diabetes and the interaction between different types of diabetic complications with cardiac autonomic dysfunction has not been investigated previously. Traditional clinical features including age and the number of years of T2DM differed significantly between the diabetic complication groups suggesting a multifactorial model for progression in diabetes related complications. Lifestyle factors, such as poor glycemic control, diet, obesity and hypertension have been shown to lead to earlier onset of diabetic complications [45]. Our results confirm these findings and further bring attention to the differences in HRV and some anthropometric and clinical data parameters with respect to different diabetic complications. The DPNn group was the youngest and the nephropathy group the oldest in our cohort. The older age of the nephropathy group may reflect age related changes in addition to diabetes, which play a role in development of nephropathy, especially as the mean years of diabetes in this group was the lowest and similar to the control group with no complications. These age-related changes could explain the significantly lower LF power in the CONT as compared to the other groups having only one complication. Blood pressure, BMI and HbA1c were all above normal in the nephropathy group, which also play an important role in the pathogenesis of diabetic nephropathy. Patients with all three complications had on average over 15 years of disease duration. Our findings show that retinopathy is more likely to occur as a single lesion later during disease progression, which differs to previous studies. In the Early Treatment Diabetic retinopathy Study approximately 20% of T2DM patients already had retinopathy at the onset of diabetes [46]. When combined with DPN, the average age of the cohort was younger bringing the age closer to the DPN cohort (Table 1). This suggests a strong link between DPN and retinopathy, with possibly large vessels being more susceptible to the effects of hyperglycemia followed by small vessel disease such as retinopathy and diabetic nephropathy in the current cohort. Single diabetes complications were significantly correlated with age (nephropathy) and DM years (retinopathy) in the current study. This agrees with previous clinical research on T1DM [47]. However associations between retinopathy and nephropathy in T2DM have not been fully explored [48]. Combinations of diabetic complications resulted in a different set of significant correlations with age. Of these only the ALL-C group was positively correlated with age. The positive correlation with age and the most severe affected diabetes complication groups explains why higher HRV indices do not always reflect better cardiac autonomic function and are similarly to lower HRV values signs of autonomic pathology [49]. All groups had increased HbA1c above recommended levels, indicating that glucose control is an important factor in developing diabetic complications [50]. Significant correlations between HRV features and clinical variables as found in this study and previously in a Caucasian cohort suggest that HDL cholesterol may also be contributing to HRV, followed by systolic blood pressure [51, 52].

The mean age of the DPNn group was the lowest despite a relative long T2DM duration similar to the remaining diabetic complications groups. Of interest is that this group had the highest mean HbA1c suggesting that peripheral nerve function is more susceptible to uncontrolled diabetes but may be also exacerbated by the type of footwear most often worn in the Middle East, where wearing open toes sandals can lead to foot deformity and loss of microvascular flow to the foot [53]. DPN is also associated with autonomic nervous system dysfunction. Symptoms include postural hypotension, abnormal sweating, gastrointestinal and urinary dysfunction, which often is not noticed by the patients [54]. Including tests for autonomic neuropathy such as HRV analysis of the cardiac rhythm are therefore an important option for diagnosis of diabetic complications and identification of ANS associated pathology [55, 56].

With respect to the type of HRV analysis most useful in clinical practice, we have recently reported that multi-lag entropy analysis addresses limitations of single lag (beat-to-beat) entropy analysis in HRV studies [57]. Multi-lag analysis considers that a heartbeat influences not only the beat immediately following it and identified by single-lag analysis, but also up to 6–10 beats downstream, which requires multi-lag analysis [58]. Our HRV results show that impairment of the ANS, characterized by decreasing time and frequency domain HRV results as well as the entropy HRV features, clearly increases as patients have more than one complication and provides a better understanding of whether the changed heart rate variability is primarily based on vagal (lags 1–6) or a combination of sympathetic and parasympathetic dysfunction (lags 7–10).

The majority of the correlation results of diabetic complications with years of DM were negative suggesting that with increasing years of diabetes, the risk of ANS dysfunction increases, but was only significant in the retinopathy group. In terms of single complications, it is interesting to observe that DPN patients showed the least difference in HRV compared to the CONT group, despite recent results that DPN is strongly associated with ANS dysfunction and suggests large fiber involvement with vascular pathology including cardiovascular disease (CVD) [54]. This lack of association found in the current study may be due to small fibre disease having a different aetiology to large fibre disease. Retinopathy patients had significantly lower LF, and LF/HF results compared to the CONT group indicating a possible sympathovagal dysfunction. The effect of the increase in the number of diabetic complications was also manifested in the multi-lag entropy, which decreased (except for the case of nephropathy). The entropy evaluates total acceleration-inhibition activities of the heart rate, or total heart period variation. The multi-lag entropy analysis gives additional information on the short-term (beat-to-beat) and long-term variations in HRV associated with parasympathetic and sympathetic function related to diabetes complications.

The relationship between autonomic impairment and diabetic complications was also seen when odds ratios were determined. The OR for all complications being present determined and LF/HF was nearly 4 times greater compared to the OR for presence of any single or multiple complications. When retinopathy was combined with DPN, the OR at higher lags of the entropy measure was significant suggesting that long-term changes in heart rate was more associated with sympathetic dysfunction. The accurate detection of sympathetic dysfunction is crucial but conventional HRV features are not sensitive to the sympathetic component [11]. Multi-lag entropy however is able to resolve sympathetic function by including higher lags in the analysis that are unlikely to contain a parasympathetic component. This was supported by the higher OR in the multi-lag measure (compared to the conventional LF/HF ratio) between the RN and the RNp groups.

Conclusion

The results of this study offer a novel approach to detecting both parasympathetic and sympathetic changes that characterize differences in patients with single or multiple diabetes complications. Our findings further characterize a local group of UAE citizens with diabetes attending a hospital diabetes outpatient clinic and hence provide a database to evaluate the influence of diabetic complications in this ethnic group.

Author Contributions

  1. Conceptualization: AHK HMA KK SL WA HSA HFJ.
  2. Data curation: AHK HMA KK SL WA HSA HFJ.
  3. Formal analysis: HMA.
  4. Investigation: AHK HMA KK SL WA HSA HFJ.
  5. Methodology: AHK HMA KK SL WA HSA HFJ.
  6. Project administration: AHK HMA KK SL WA HSA HFJ.
  7. Supervision: AHK HMA KK SL WA HSA HFJ.
  8. Validation: AHK HMA KK SL WA HSA HFJ.
  9. Visualization: AHK HMA KK SL WA HSA HFJ.
  10. Writing – original draft: AHK HMA KK SL WA HSA HFJ.
  11. Writing – review & editing: AHK HMA KK SL WA HSA HFJ.

References

  1. 1. Brownlee M. Biochemistry and molecular cell biology of diabetic complications. Nature. 2001;414(6865):813–20. pmid:11742414
  2. 2. Brennan E, McEvoy C, Sadlier D, Godson C, Martin F. The genetics of diabetic nephropathy. Genes. 2013;4(4):596–619. pmid:24705265
  3. 3. Wilson N, Wright D. Inflammatory mediators in diabetic neuropathy. Journal of Diabetes. 2011:S5:004.
  4. 4. Fowler M. Microvascular and macrovascular complications of diabetes. Clinical Diabetes. 2008;26(2):77–82.
  5. 5. Krishnan AV, Kiernan MC. Uremic neuropathy: clinical features and new pathophysiological insights. Muscle Nerve. 2007;35(3):273–90. pmid:17195171
  6. 6. Krolewski AS, Barzilay J, Warram JH, Martin BC, Pfeifer M, Rand LI. Risk of early-onset proliferative retinopathy in IDDM is closely related to cardiovascular autonomic neuropathy. Diabetes. 1992;41(4):430–7. pmid:1607070
  7. 7. Smith SE, Smith SA, Brown PM. Cardiac autonomic dysfunction in patients with diabetic retinopathy. Diabetologia. 1981;21(6):525–8. pmid:7338291
  8. 8. Vinik AI, Mehrabyan A. Diabetic neuropathies. The Medical Clinics of North America. 2004;88(4):947–99, xi. pmid:15308387
  9. 9. Vinik A, Maser R, Nakave A. Diabetic cardiovascular autonomic nerve dysfunction. US Endocrine Disease. 2007;2:66–74.
  10. 10. Ewing DJ, Borsey DQ, Bellavere F, Clarke BF. Cardiac autonomic neuropathy in diabetes: comparison of measures of R-R interval variation. Diabetologia. 1981;21(1):18–24. pmid:7274613
  11. 11. Bellavere F, Balzani I, De Masi G, Carraro M, Carenza P, Cobelli C, et al. Power spectral analysis of heart-rate variations improves assessment of diabetic cardiac autonomic neuropathy. Diabetes. 1992;41(5):633–40. pmid:1568534
  12. 12. Jelinek HF, Cree MJ. Automated Image Detection of Retinal Pathology. Boca Raton: CRC Press; 2009.
  13. 13. Vinik AI, Maser RE, Mitchell BD, Freeman R. Diabetic autonomic neuropathy. Diabetes Care. 2003;26(5):1553–79. pmid:12716821
  14. 14. Khandoker AH, Jelinek HF, Palaniswami M. Identifying diabetic patients with cardiac autonomic neuropathy by heart rate complexity analysis. Biomedical Engineering Online. 2009;8:3. pmid:19178728
  15. 15. Khandoker AH, Jelinek HF, Moritani T, Palaniswami M. Association of cardiac autonomic neuropathy with alteration of sympatho-vagal balance through heart rate variability analysis. Medical Engineering & Physics. 2010;32(2):161–7.
  16. 16. Khandoker AH, Fujibayashi M, Moritani T, Palaniswami M, editors. Assessing sympatho-vagal balance in schizophrenia through tone-entropy analysis. Computing in Cardiology; 2010; Belfast, UK.
  17. 17. Khandoker AH, Luthra V, Abouallaban Y, Saha S, Ahmed KI, Mostafa R, et al., editors. Reduced variability in pulse wave velocity and heart rate in depressed patients with suicidal ideation. Computing in Cardiology; 2015.
  18. 18. Marzbanrad F, Khandoker AH, Hambly BD, Ng E, Tamayo M, Lu Y, et al. Methodological comparisons of heart rate variability analysis in patients with type 2 diabetes and angiotensin converting enzyme polymorphism. IEEE Journal of Biomedical and Health Informatics. 2016;20(1):55–63. pmid:26529791
  19. 19. James PA, Oparil S, Carter BL, Cushman WC, Dennison-Himmelfarb C, Handler J, et al. Evidence-based guidelines for the management of high blood pressure in adults: report from the panel members appointed to the Eighth Joint National Committee (JNC 8). JAMA. 2014;311:507–20. pmid:24352797
  20. 20. Aring AM, Jones DE, Falko JM. Evaluation and prevention of diabetic neuropathy. American Family Physician. 2005;71(11):2123–8. pmid:15952441
  21. 21. McCarty CA, Taylor KI, McKay R, Keeffe JE, Working Group on Evaluation of NDRG. Diabetic retinopathy: effects of national guidelines on the referral, examination and treatment practices of ophthalmologists and optometrists. Clinical & Experimental Ophthalmology. 2001;29(2):52–8.
  22. 22. Jelinek HF, Khandoker A, Palaniswami M, McDonald S, editors. Tone-Entropy analysis as a cardiac risk stratification tool. Computing in Cardiology; 2010; Belfast, Ireland: IEEE Press.
  23. 23. Pan J, Tompkins WJ. A real-time QRS detection algorithm. IEEE Transactions on Biomedical Engineering. 1985;32(3):230–6. pmid:3997178
  24. 24. Jelinek HF, Alothman T, Cornforth DJ, Khalaf K, Khandoker A, editors. Effect of biosignal preprocessing and recording length on clinical decision making for cardiac autonomic neuropathy. ESGCO 2014; 2014; Trento, Italy.
  25. 25. Brennan M, Palaniswami M, Kamen P. Do existing measures of Poincare plot geometry reflect nonlinear features of heart rate variability? IEEE Transactions on Biomedical Engineering. 2001;48(11):1342–7. pmid:11686633
  26. 26. Welch PD. The use of fast fourier transform for the estimation of power spectra: A method based on time averaging over short, modified periodograms. IEEE Transactions on Audio and Electroacoustics. 1967;15(2):70–3.
  27. 27. Malik M. 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. Circulation. 1996;93(5):1043–65. pmid:8598068
  28. 28. Shannon CE. A mathematical theory of communication. Bell System Technical Journal. 1948;27:379–423.
  29. 29. Oida E, Moritani T, Yamori Y. Tone-entropy analysis on cardiac recovery after dynamic exercise. Journal of Applied Physiology. 1997;82(6):1794–801. pmid:9173943
  30. 30. Amano M, Oida E, Moritani T. Age-associated alteration of sympatho-vagal balance in a female population assessed through the tone—entropy analysis. European Journal of Applied Physiology. 2005;94:602–10. pmid:15942770
  31. 31. Quade D. Rank analysis of covariance. Journal of the American Statistical Association. 1967;62(320):1187–200.
  32. 32. Maxwell SE, O’Callaghan MF, Delaney HD. Analysis of covariance In: Edwards LK, editor. Applied Analysis of Variance in Behavioral Science. New York: Marcel Dekker. Inc.; 1993. p. 63–104.
  33. 33. Cade WT. Diabetes-related microvascular and macrovascular diseases in the physical therapy setting. Physical Therapy. 2008;88(11):1322–35. pmid:18801863
  34. 34. Vinik AI, Ziegler D. Diabetic cardiovascular autonomic neuropathy. Circulation. 2007;115(3):387–97. pmid:17242296
  35. 35. Sasso FC, Chiodini P, Carbonara O, De Nicola L, Conte G, Salvatore T, et al. High cardiovascular risk in patients with Type 2 diabetic nephropathy: the predictive role of albuminuria and glomerular filtration rate. The NID-2 Prospective Cohort Study. Nephrol Dial Transplant. 2012;27(6):2269–74. pmid:22090446
  36. 36. Targher G, Bertolini L, Zenari L, Lippi G, Pichiri I, Zoppini G, et al. Diabetic retinopathy is associated with an increased incidence of cardiovascular events in Type 2 diabetic patients. Diabetic medicine: a journal of the British Diabetic Association. 2008;25(1):45–50.
  37. 37. Nielsen JV. Peripheral neuropathy, hypertension, foot ulcers and amputations among Saudi Arabian patients with type 2 diabetes. Diabetes Res Clin Pract. 1998;41(1):63–9. pmid:9768374
  38. 38. Paisey RB, Arredondo G, Villalobos A, Lozano O, Guevara L, Kelly S. Association of differing dietary, metabolic, and clinical risk factors with microvascular complications of diabetes: a prevalence study of 503 Mexican type II diabetic subjects. II. Diabetes Care. 1984;7(5):428–33. pmid:6499636
  39. 39. Young MJ, Boulton AJ, MacLeod AF, Williams DR, Sonksen PH. A multicentre study of the prevalence of diabetic peripheral neuropathy in the United Kingdom hospital clinic population. Diabetologia. 1993;36(2):150–4. pmid:8458529
  40. 40. Herman WH, Aubert RE, Engelgau MM, Thompson TJ, Ali MA, Sous ES, et al. Diabetes mellitus in Egypt: glycaemic control and microvascular and neuropathic complications. Diabetic medicine: a journal of the British Diabetic Association. 1998;15(12):1045–51.
  41. 41. Jbour AS, Jarrah NS, Radaideh AM, Shegem NS, Bader IM, Batieha AM, et al. Prevalence and predictors of diabetic foot syndrome in type 2 diabetes mellitus in Jordan. Saudi medical journal. 2003;24(7):761–4. pmid:12883610
  42. 42. Hussein MM, Mooij JM, Roujouleh H, el-Sayed H. Observations in a Saudi-Arabian dialysis population over a 13-year period. Nephrol Dial Transplant. 1994;9(8):1072–6. pmid:7800203
  43. 43. El-Asrar AM, Al-Rubeaan KA, Al-Amro SA, Kangave D, Moharram OA. Risk factors for diabetic retinopathy among Saudi diabetics. International ophthalmology. 1998;22(3):155–61. pmid:10548460
  44. 44. Khandekar R, Al Lawatii J, Mohammed AJ, Al Raisi A. Diabetic retinopathy in Oman: a hospital based study. The British journal of ophthalmology. 2003;87(9):1061–4. pmid:12928265
  45. 45. Drummond KN, Kramer MS, Suissa S, Lévy-Marchal C, Dell'Aniello S, Sinaiko A, et al. Effects of duration and age at onset of type 1 diabetes on preclinical manifestations of nephropathy. Diabetes. 2003;52(7):1818–24. pmid:12829652
  46. 46. Fong DS, Aiello LP, F L. F, K R.. Diabetic retinopathy. Diabetes Care. 2004;27:2540–53. pmid:15451934
  47. 47. Klein R, Zinman B, Gardiner R, Suissa S, Donnelly SM, Sinaiko AR, et al. The relationship of diabetic retinopathy to preclinical diabetic glomerulopathy lesions in type 1 diabetic patients: the Renin-Angiotensin System Study. Diabetes. 2005;54(2):527–33. pmid:15677511
  48. 48. Penno G, Solini A, Zoppini G, Orsi E, Zerbini G, Trevisan R, et al. Rate and determinants of association between advanced retinopathy and chronic kidney disease in patients with type 2 diabetes: the Renal Insufficiency And Cardiovascular Events (RIACE) Italian multicenter study. Diabetes Care. 2012;35(11):2317–23. pmid:23093684
  49. 49. Stein PK, Domitrovich PP, Hui N, Rautaharju P, Gottdiener J. Sometimes higher heart rate variability is not better heart rate variability: results of graphical and nonlinear analyses. Journal of Cardiovascular Electrophysiology. 2005;16(9):954–9. pmid:16174015
  50. 50. Forbes JM, Cooper ME. Mechanisms of diabetic complications. Physiological Reviews. 2013;93(1):137–88. pmid:23303908
  51. 51. Christensen JH, Toft E, Christensen MS, Schmidt EB. Heart rate variability and plasma lipids in men with and without ischaemic heart disease. Atherosclerosis. 1999;145(1):181–6. pmid:10428309
  52. 52. Jelinek HF, Md Imam H, Al-Aubaidy H, Khandoker AH. Association of cardiovascular risk using non-linear heart rate variability measures with the framingham risk score in a rural population. Front Physiol. 2013;4:186. pmid:23898302
  53. 53. Wang DD, Bakhotmah BA, Hu FB, Alzahrani HA. Prevalence and correlates of diabetic peripheral neuropathy in a Saudi Arabic population: a cross-sectional study. PLoS One. 2014;9(9):e106935. pmid:25184511
  54. 54. Singer W, Spies JM, McArthur J, Low J, Griffin JW, Nickander KK, et al. Prospective evaluation of somatic and autonomic small fibers in selected autonomic neuropathies. Neurology. 2004;62(4):612–8. pmid:14981179
  55. 55. Dabby R, Vaknine H, Gilad R, Djaldetti R, Sadeh M. Evaluation of cutaneous autonomic innervation in idiopathic sensory small-fiber neuropathy. JPNS. 2007;12(2):98–101. pmid:17565534
  56. 56. Tarvainen MP, Laitinen TP, Lipponen JA, Cornforth DJ, Jelinek HF. Cardiac autonomic dysfunction in type 2 diabetes—effect of hyperglycemia and disease duration. Frontiers in endocrinology. 2014;5:130. pmid:25152747
  57. 57. Karmakar CK, Khandoker AH, Jelinek HF, Palaniswami M. Risk stratification of cardiac autonomic neuropathy based on multi-lag Tone-Entropy. Medical & Biological Engineering & Computing. 2013;51(5):537–46.
  58. 58. Lerma C, Infante O, Perez-Grovas H, Jose MV. Poincare plot indexes of heart rate variability capture dynamic adaptations after haemodialysis in chronic renal failure patients. Clinical Physiology and Functional Imaging. 2003;23(2):72–80. pmid:12641600