Abstract

Aim. The present study was designed to evaluate the heart rate variability (HRV) in nonalcoholic fatty liver disease (NAFLD) and to assess the effect of grade of NAFLD and diabetic status on HRV. Methods. This cross-sectional study included 75 subjects (25 NAFLD without diabetes, 25 NAFLD with diabetes, and 25 controls). Measurements included anthropometry, body composition analysis, estimation of plasma glucose, serum lipids, hsCRP, and serum insulin. HRV analysis was performed in both time and frequency domains. Results. The time and frequency domain indices of overall variability (SDNN, total power) were significantly lower in NAFLD with diabetes as compared to the controls. However, the LF : HF ratio did not differ among the three groups. The variables related to obesity, lipid profile, and glucose metabolism were also higher in NAFLD with diabetes and those with Grade II NAFLD without diabetes, as compared to controls. Multivariate stepwise regression analysis showed a negative correlation between HRV and total cholesterol and fat percentage. Conclusion. The grade of NAFLD as well as diabetic status contributes to the decrease in the cardiovascular autonomic function, with diabetic status rather than grade of NAFLD playing a critical role. Serum lipids and adiposity may also contribute to cardiac autonomic dysfunction.

1. Introduction

Nonalcoholic fatty liver disease (NAFLD) is a clinicopathological condition characterised by lipid deposition in the liver and is a common cause of liver dysfunction. The prevalence of NAFLD in Asian population ranges from 9% to 45% [1] and is estimated to be about 30% in western population [2]. Apart from an increased risk of liver-related morbidity and mortality, patients of NAFLD also have higher cardiovascular risk [3, 4] especially when present along with type 2 diabetes (T2DM) [4, 5].

A balanced autonomic output to liver is important for maintenance of circadian rhythms of liver metabolic enzymes and glucose level [6]. Imbalance in autonomic function has been proposed as a component in pathogenesis of NAFLD [7]. The autonomic dysfunction has been shown to be higher in the NAFLD. The autonomic symptom burden assessed by orthostatic grading scale was higher in nondiabetic NAFLD [8] and the sudomotor dysfunction was also higher in the NAFLD after accounting for all confounders [9].

A recent systematic review of 27 studies showed an association of NAFLD with subclinical atherosclerosis independent of traditional risk factors and metabolic syndrome [10]. Similar independent association of NAFLD with subclinical atherosclerosis has been observed in Asian Indians also [11]. In other clinical settings, it has been shown that arterial stiffness is inversely related to heart rate variability [12, 13]. However, there are only a few studies in which heart rate variability has been evaluated in NAFLD [9, 14, 15]. Decrease in indices of HRV has been reported by Liu et al. with decrease in parasympathetic component after adjustment of covariates [14]. A higher LF : HF ratio (low frequency : high frequency) of HRV was reported in patients () with NAFLD along with a lower baroreflex sensitivity [15].

Cardiovascular autonomic neuropathy is a common complication of diabetes [16]. NAFLD commonly coexists with diabetes. Both may act synergistically leading to varied clinical outcomes [17]. The independent contribution of diabetes and grade of NAFLD in development of autonomic dysfunction is not known. The present study was designed to evaluate cardiovascular autonomic status in the patients of NAFLD with or without diabetes as compared to control population of individuals without either diabetes or NAFLD. In the present study we evaluated the association of indices of HRV with anthropometric variables, lipid profile, and diabetic status in patients with and without NAFLD.

2. Methods

This was a cross-sectional study approved by the Institute Ethics Committee of the All India Institute of Medical Sciences (AIIMS), New Delhi. The subjects were recruited from the medicine outpatient department. Informed written consent was obtained after explanation of the purpose, type, and duration of the study.

2.1. Study Subjects

A total of 75 subjects were recruited for the study. The subjects were grouped into NAFLD patients without diabetes (, 13 females and 12 males), NAFLD patients with diabetes of <5 years of duration (, 17 females and 8 males), and healthy controls with similar age distribution (, 9 females and 16 males). Subjects with alcohol consumption >140 gm/week, any secondary cause of fatty liver, hereditary disorders, any severe acute or chronic illness, seropositivity for hepatitis B, hepatitis C, and HIV, established coronary heart disease, pregnancy, and lactation were excluded. Diagnosis of NAFLD was made using ultrasonography performed by an experienced sonologist who was blinded to the clinical data of the patients, using a 3.5 MHz convex transducer by subcostal and intercostal approach (Volusion, GE).

2.2. Measurements
2.2.1. Anthropometry and Body Composition

A detailed clinical history and clinical examination was done. All the subjects underwent anthropometric measurements and body composition analysis by bioelectrical impedance method as described earlier [11]. Height and weight were recorded to the nearest 0.1 cm and nearest 0.1 kg with stadiometer and electronic scale. Body Mass Index (BMI) was calculated by weight (kg)/height (m2). Waist circumference (WC) was measured midway between the iliac crest and the lower most margin of the ribs. Hip circumference (HC) was measured at the maximum circumference of buttocks. Body fat and percentage body fat were estimated by foot-to-foot bioelectrical impedance technique (TANITA, Japan).

2.2.2. Biochemical Measurements

Venous blood samples were obtained after an overnight fast of at least 10 hours. Biochemical measurements included liver function tests, fasting blood glucose (FBG) and postprandial blood glucose (PPBG) levels, and fasting lipid profile using standard methods as described earlier [11]. Serum levels of hsCRP were measured using a commercially available reagent kit based on the principle of solid phase enzyme-linked immune sorbent assay (Biocheck, Inc., Foster City, USA). Fasting insulin was measured using a commercially available reagent kit based on the principle of electrochemiluminescence (Liaison® Insulin, DiaSorin Inc., USA). The inter- and intra-assay coefficient of variation was <5%. The HOMA-IR was calculated by using the following formula: HOMA-IR = fasting insulin (µU/mL) fasting blood glucose (mmol/L)/22.5.

2.2.3. Assessment of Heart Rate Variability

The cardiovascular autonomic status was estimated in the Autonomic Function Testing Laboratory (AFT Lab), Department of Physiology, AIIMS, New Delhi. The subject was requested to lie down on the couch where electrodes for lead II ECG acquisition by Labchart Pro 7® (AD Instruments, Australia) were attached. The recording was performed in a controlled ambient temperature of 22 to 25°C. After an initial period of rest of 5 minutes, a 5-minute lead II ECG (sampling rate 1 KHz) was recorded for later offline analysis. The subject was instructed to avoid movement during data acquisition to prevent artefacts in the recording.

Offline analysis of ECG was done using Labchart Pro 7 (AD Instruments, Australia) and HRV was computed using Hemolab software (version 8.5, Harald Stauss Scientific, Iowa). For time domain analysis the beat-to-beat heart rate series was computed from ECG using Labchart Pro and was then exported as a .txt file. The file format was converted to .asc as required by the Hemolab software for batch processing. The time domain indices of HRV computed were SDNN (standard deviation of all NN intervals), SDSD (standard deviation of differences between adjacent NN intervals), and pNN50% (NN50 count divided by the total number of all NN intervals). For frequency domain analysis RR intervals were computed from ECG using Labchart Pro and were then exported as a.txt file. The file format was converted to .asc and then to .tsa for batch processing. This was followed by spline interpolation at sampling frequency of 5 Hz followed by the power spectral density (PSD) calculation using Fast Fourier Transformation (FFT). The measurement of low frequency (LF: 0.04–0.15 Hz) and high frequency (HF: 0.15–0.40 Hz) power components was made in absolute values of power (ms2). The powers were further normalized to account for changes in the total power of the HRV.

2.3. Statistical Analysis

The data was analysed using SPSS software (Version 22, IBM). The normality of distribution was estimated with Shapiro-Wilk test. Quantitative data is expressed as mean ± SD for normally distributed data, and median with interquartile range (IQR) is used to express skewed data. For parametric data, one-way ANOVA followed by Tukey’s test was used and for nonparametric data Kruskal Wallis H test followed by Dunn’s comparison was done. was considered to be statistically significant. Two-way ANOVA was employed to evaluate the relative contribution of the diabetic status and grade of NAFLD. Correlation analysis was done with Spearman’s test followed by stepwise multivariate analysis.

3. Results

The anthropometric measures, biochemical investigation, and lead II ECG recordings were measured in all the subjects. The data of 7 subjects in the NAFLD without diabetes group was excluded from HRV due to presence of artefacts that prevented computation.

3.1. Anthropometric Features and Biochemical Profile

Table 1 shows the general, anthropometric, and biochemical parameters of the subjects. Patients with NAFLD and diabetes were older (42.9 ± 7.6 yrs) than patients with NAFLD without diabetes (41.8 ± 7.3 yrs) and controls (37.9 ± 8.4 yrs). The BMI, fat mass, and fat percentage were higher in NAFLD with diabetes group as well as NAFLD without diabetes group as compared to the controls but there was no difference between NAFLD with and without diabetes. Value of WHR was higher only in NAFLD without diabetes as compared to the controls. The fat percentage was the only measure that was higher in the NAFLD with diabetes as compared to the NAFLD without diabetes.

The ESR was higher in NAFLD with diabetes while hsCRP was higher in NAFLD without diabetes as compared to controls. All the indices of liver function (bilirubin, SGOT, SGPT, and ALP) were not significantly different among the three groups. As expected, the values of FBG, PPBG, and HOMA-IR were significantly higher in NAFLD with diabetes group as compared to the other two groups. Serum levels of TC and LDL were higher in NAFLD with diabetes group as compared to the control group.

3.2. Heart Rate Variability (HRV)

The time and frequency domain indices of HRV are shown in Table 2. The overall variability (time domain: SDNN and frequency domain: total power) was lower in NAFLD with diabetes group as compared to the control group. Even though median values of the overall variability of HRV were lower in NAFLD group without diabetes as compared to controls, the difference was not statistically significant. Similarly, the overall variability of NAFLD with diabetes was not statistically different from NAFLD without diabetes even though the median values were lower. The low frequency and high frequency component of HRV and the LF : HF ratio were similar in all three groups. In time domain analysis, the indices of parasympathetic component of HRV (SDSD and pNN50) were lower in NAFLD with diabetes group as compared to controls as well as NAFLD without diabetes.

The proportional distribution of Grade I and Grade II was not different in the NAFLD groups with or without diabetes. The subjects were regrouped depending upon the grades of NAFLD irrespective of their diabetic status and statistical analysis was done (Table 3). The distribution of diabetics in NAFLD Grade I and Grade II was not statistically different. The indices of overall variability (SDNN and total power) and parasympathetic component (SDSD, pNN50) were lower in Grade II as compared to controls. Even though the median values of the same parameters were lower in Grade I when compared with controls, it was not statistically significant.

In order to delineate the relative contribution of the diabetic status and grade of NAFLD, two-way ANOVA was done on all the NAFLD subjects. The presence of diabetes was found to be significantly associated with decrease in SDSD and total power (Table 4).

3.3. Association between Indices of HRV and Anthropometric and Biochemical Parameters

The result of Spearman’s correlation performed between the selected indices of HRV (SDNN, SDSD, and total power) and anthropometric and biochemical parameters is shown in Table 5. No significant association was found between indices of HRV, parameters of inflammation (ESR, hsCRP), liver function (serum bilirubin, SGOT, SGPT, and ALP), or insulin resistance (HOMA-IR, fasting insulin). The rest of the parameters were negatively correlated with the indices of HRV.

A further multivariate stepwise regression analysis was performed between indices of HRV as dependent variable and selected independent variables (with significant association in univariate analysis). Table 6 shows the two models. Of all the parameters that were negatively associated with HRV, only total cholesterol (TC) and fat percentage significantly associated with SDNN and total power after multivariate analysis.

4. Discussion

The study was done to assess the heart rate variability in patients of NAFLD with or without diabetes. Since diabetes is known to independently cause autonomic dysfunction, the patients of nonalcoholic fatty liver disease were grouped into NAFLD without diabetes and NAFLD with diabetes.

4.1. Anthropometric Features and Biochemical Profile

As expected, the variables related to obesity, lipid profile, and glucose metabolism were higher in NAFLD with diabetes and those with Grade II NAFLD as compared to controls [1820]. A higher ESR in NAFLD with diabetes and higher hsCRP in NAFLD without diabetes indicate ongoing inflammatory process. Similar high hsCRP was reported by Nigam et al. [21] and Ajmal et al. [22], even after adjusting for covariates.

4.2. Heart Rate Variability (HRV)

The time domain and frequency domain indices of overall variability of HRV were lower in the NAFLD with diabetes group as compared to the control. Even though the median values of these parameters were lower in NAFLD without diabetes, the difference was not significant. The indices of parasympathetic component of HRV (SDNN, pNN50) were lower in NAFLD with diabetes as compared to controls as well as NAFLD without diabetes. When the analysis was done on the basis of grades of NAFLD irrespective of the diabetic status, the indices of overall HRV were found to be lower in Grade II NAFLD as compared to controls. Though the median values of the theses indices were lower in Grade I NAFLD they were not statistically significant. This indicates an independent contribution of diabetic status as well as grade of NAFLD in development of autonomic dysfunction. Within the patients of NAFLD, the diabetic status had significant main effect in decreasing HRV rather than the grade.

Similar decrease in indices of overall HRV (Ln SDNN) and parasympathetic components (Ln RMSSD) was noted by Liu et al. [23]. In the present study, no change in normalized low frequency, high frequency, or their ratio was found. Liu et al., however, have reported lower low frequency and high frequency components of HRV [14, 23]. This difference is likely due to reporting absolute values rather than normalizing the values to total power of the HRV as was done in the present study. The data of the present study suggests that decrease in total power of HRV results from a decrease in both sympathetic and parasympathetic components. In the presence of diabetes, the decrease in the parasympathetic component is more [24]. Jakovljevic et al. had reported higher values of LF : HF ratio though it had not been compared with any control group. In the present study LF : HF ratio was similar in NAFLD with or without diabetes as compared to controls. This difference could be due to inclusion of diabetics with less 5 years of history as well as lower BMI of the study subjects in the present study as compared to the report by Jakovljevic et al. [15].

In univariate analysis, the decrease in indices of HRV was significantly and negatively correlated with variables of lipid profile, obesity, and diabetic status. Interestingly, the decrease in HRV was not associated either with HOMA-IR or with fasting insulin but with fasting blood and postprandial glucose level. These observations are in tune with the known role blood glucose level and glycation play in the pathogenesis of diabetic autonomic neuropathy [25]. However, in multivariate stepwise regression analysis, indices of HRV were negatively associated with total cholesterol and fat percentage only. The importance of total serum cholesterol in the development of diabetic autonomic neuropathy has also been reported earlier [26]. These observations are similar to a report by Pimenta et al. where the autonomic control in NAFLD (measured by heart rate recovery after maximum graded exercise test) correlated with body composition and body fat [27]. In a recent review, importance of hepatic accumulation of fat in a setting of obesity for the development of hepatic/peripheral insulin resistance has been proposed [20]. NAFLD and diabetes then independently lead to the development of cardiac autonomic dysfunction.

The observation of importance of diabetic status in decrease in overall HRV is in line with observations of increase prevalence of cardiovascular disease in NAFLD with diabetes [2830]. A recent systematic review of 27 studies showed an association of NAFLD with subclinical atherosclerosis independent of traditional risk factors and metabolic syndrome [10]. In other clinical settings, it has been shown that arterial stiffness is inversely related to heart rate variability [12, 13]. It is probable that increase in stiffness of arteries due to subclinical atherosclerosis leads to decrease in the transducer function of the baroreceptors. A lower baroreflex sensitivity in NAFLD has been reported in a single study [15]. A decrease in baroreflex function will decrease the heart rate variability of both sympathetic and parasympathetic components.

5. Conclusion

The result of the present study indicates that the grade of NAFLD as well as diabetic status contributes to the decrease in the cardiovascular autonomic function with a decrease in overall variability but an unchanged sympathovagal balance. It also shows that once NAFLD is developed a further decrease is more likely due to diabetes rather than further increase in the grade of NAFLD and the important role of dyslipidemia and obesity in the cardiac autonomic dysfunction represented by heart rate variability.

Abbreviations

ANOVA:Analysis of variance
ALP:Alkaline Phosphatase
ALT:Alanine Transaminase
AST:Aspartate Transaminase
BMI:Body Mass Index
ESR:Erythrocyte Sedimentation Rate
FBG:Fasting blood glucose
FFT:Fast Fourier Transformation
HC:Hip circumference
HDL:High Density Lipoprotein
HF:High frequency (0.15–0.40 Hz)
HOMA-IR:Homeostatic Model Assessment for Insulin Resistance
HRVHeart rate variability
hsCRP:High-Sensitivity C Reactive Protein
IQR:Interquartile range.
Reactive Protein
IU:International Unit
LDL:Low Density Lipoprotein
LF:Low frequency (0.04–0.15 Hz)
NAFLD:Nonalcoholic fatty liver disease
nu:Normalized unit
pNN50%:NN50 count divided by the total number of all NN intervals
PPBG:Postprandial blood glucose
PSD:Power spectral density
SDNN:Standard deviation of all NN intervals
SDSD:Standard deviation of differences between adjacent NN intervals
SGOT:Serum Glutamic Oxaloacetic Transaminase
SGPT:Serum Glutamate Pyruvate Transaminase
T2DM:Type 2 diabetes mellitus
TC:Total cholesterol
TG:Triglycerides
VLDL:Very Low Density Lipoprotein
WC:Waist circumference
WHR:Waist Hip Ratio.

Additional Points

Strengths and Limitations. The present study is first study on effect of grade and diabetic status of NAFLD on heart rate variability. The exclusion of subjects with duration longer than 5 years of diabetes provides a suitable group to elucidate the contribution of diabetic status. A limitation of the study is that it is cross-sectional in design with number not sufficient to have subgroup analysis based on Grade I with or without diabetes and Grade II with or without diabetes. A prospective design with long term follow-up can further establish the observation of the present study.

Competing Interests

Authors declare no conflict of interests.

Acknowledgments

The authors are thankful to Dr. Surya Prakash for his assistance in performing anthropometric and biochemical measurements. The cooperation of subjects who agreed to participate in this study is highly appreciated.