Electrocardiographic changes with the onset of diabetes and the impact of aerobic exercise training in the Zucker Diabetic Fatty (ZDF) rat

Background Early markers of diabetic autonomic neuropathy (DAN) in an electrocardiogram (ECG) include elevated R wave amplitudes, widening of QTc intervals and decreased heart rate variability (HRV). The severity of DAN has a direct relationship with mortality risk. Aerobic exercise training is a common recommendation for the delay and possible reversal of cardiac dysfunction. Limited research exists on ECG measures for the evaluation of aerobic exercise training in Zucker Diabetic Fatty (ZDF) rat, a model of type 2 diabetes. The objective of this study was to assess whether aerobic exercise training may attenuate diabetes induced ECG changes. Methods Male ZDF (obese fa/fa) and control Zucker (lean fa/+) rats were assigned to 4 groups: sedentary control (SC), sedentary diabetic (SD), exercised control (EC) and exercised diabetic (ED). The exercised groups began 7 weeks of treadmill training after the development of diabetes in the ED group. Baseline (prior to the training) and termination measurements included body weight, heart weight, blood glucose and glycated hemoglobin levels and ECG parameters. One way repeated measures ANOVA (group) analyzed within and between subject differences and interactions. Pearson coefficients and descriptive statistics described variable relationships and animal characteristics. Results Diabetes caused crucial changes in R wave amplitudes (p < 0.001), heart rate variability (p < 0.01), QT intervals (p < 0.001) and QTc intervals (p < 0.001). R wave amplitude augmentation in SD rats from baseline to termination was ameliorated by exercise, resulting in R wave amplitude changes in ED animals similar to control rats. Aerobic exercise training neither attenuated QT or QTc interval prolongation nor restored decreases in HRV in diabetic rats. Conclusion This study revealed alterations in R wave amplitudes, HRV, QT and QTc intervals in ZDF rats. Of these changes, aerobic exercise training was able to correct R wave amplitude changes. In addition, exercise has beneficial effect in this diabetic rat model in regards to ECG correlates of left ventricular mass.


Background
Cardiovascular disease (CVD) accounts for the majority of deaths for people with type 2 diabetes mellitus. CVD is a broad term which includes any condition causing pathological changes in blood vessels, cardiac muscle or valves and cardiac rhythm. The electrocardiogram (ECG) offers a quick, non-invasive clinical and research screen for the early detection of CVD.
Electrocardiographic changes in raw and corrected QT intervals and R wave amplitudes are early indicators of evolving CVD and increased cardiovascular risk. Prolonged QT and QT c intervals are considered reliable predictors of heart disease and fatal ventricular arrhythmias [1][2][3]. A positive linear relationship exists between QT c interval prolongation and diabetic cardiac autonomic neuropathy (DAN) severity in diabetic population [4]. Heart rate variability (HRV), one indicator of DAN, decreases with diabetes which indicates increased mortality risk [5]. QT and QT c interval abnormalities reflect changes in cardiac architecture. A positive correlation between QT or QT c interval prolongation and left ventricular (LV) mass has been reported [6][7][8]. LV hypertrophy presents as exaggerated R wave amplitudes on ECG recordings. Elevated R wave amplitudes are an independent risk factor for cardiovascular events [9]. LV hypertrophy and QT interval alterations coupled with decreased cardiac function are commonly observed with diabetes related CVD [10].
Non pharmacological interventions for CVD focus primarily on lifestyle changes with physical activity as the primary focus and a risk reduction strategy. Physical activity reduces QT c interval prolongation and cardiac dysfunction in healthy subjects [11,12]. Exercise lowers heart rate and increases HRV in healthy and diseased populations [13,14]. Physical activity can serve as potent prescription in the delay and attenuation of the CVD complications for persons with type 2 diabetes but additional comparative studies are needed regarding the cardiac response to exercise under diabetic conditions at various time points of the disease.
The Zucker Diabetic Fatty (ZDF) rat is a model of type 2 diabetes. The ZDF rat develops hyperglycemia and hyperlipidemia by week 8 and overt diabetes by week 12. The progression mimics the obesity-related insulin resistance and inflammation seen in humans [15,16]. The ZDF rat is commonly used to investigate prevention of diabetes; however, research related to the diabetic heart disease including ECG studies with this animal model is limited. We hypothesized that ECG changes occur in ZDF rats early in the disease process and aerobic exercise training is able to alleviate the changes. We detected changes in ECG parameters that were partially corrected by exercise training. Our findings add to the characterization of the ZDF model for studying type 2 diabetes effects on the heart and explore the benefits of an early exercise intervention in the presentation and progression of diabetes related CVD.

Animals
Male Zucker Diabetic Fatty (fa/fa) rats of 11 weeks were utilized for the study. Male lean, age-matched Zucker (fa/+) rats (both from Charles River Laboratory, Saint Louis, MO) served as non-diabetic controls. The animals were allowed food and water ad libitum and were placed on a 12:12 light-dark cycle. As per vendor's recommendations, the animals were fed Purina 5008 diet during the entire study for the development of a disease process resembling type 2 diabetes and its complications. All animal procedures were performed according to the University of Kansas Medical Center Institutional Animal Care and Use Committee guidelines and an approved Animal Care and Use Protocol.

Measurements
Body weights and blood glucose levels were measured weekly on all animals. Blood glucose levels were measured from rat tail using Accu-Check Active meter (Roche Diagnostics, Indianapolis, IN). Glycated hemoglobin (HbA1c) levels were measured at the end of the experiment using antibody-based A1cNow meter (Metrika, Sunnyvale, CA). When rats had blood glucose or HbA1c levels higher than detectable by the method used, we used the highest detectable value (600 mmol/L or 13%, respectively) for statistical purposes. Animals were killed within 36 hours of the last exercise training episode.

Aerobic exercise training
The rats started a treadmill exercise program at 12 weeks of age, immediately after the onset of diabetes, and continued exercising for 7 weeks. Four groups of rats were used: sedentary control (SC, n = 12), sedentary diabetic (SD, n = 10), exercised control (EC, n = 10) and exercised diabetic (ED, n = 12). This training protocol has been published by our group previously for a rat model of type 1 diabetes [17] and was adapted for the obese diabetic rats as they were unable to perform at that intensity level. During the first week of training the animals ran at 10 m/min with time increased from 10 min per day to 40 min per day at the end of the week. The progression allowed the rats to acclimate to the treadmill. Starting at week two and until the completion of the training session the rats ran at 15 m/min, for 40 minutes, 5 days per week. In order to accommodate for the disease progression in the diabetic rats, any animals showing signs of fatigue were allowed breaks of a few minutes until they were able to continue, for a total run time of 40 min per day. All rats assigned to the exercise groups completed the exercise training protocol.

Electrocardiogram (ECG) assessment
Animals received ketamine (60 mg/kg) and xylazine (7 mg/kg) prior to the resting ECG recording. ECG leads I, II, III, aVR, aVL, aVF were recorded with surface electrodes (ADInstruments, Colorado Springs, CO). Measurements were collected at baseline, prior to training and after 7 weeks of exercise training. The mean value for each rat was obtained from four values consisting of four consecutive cardiac cycles using LabChart software (ADInstruments, Colorado Springs, CO). Corrected QT (QT C ) was calculated with mean values and the Bazett's Formula (QT c = QT Interval/√ (RR interval) [18]. The heart rate (bpm) for each animal was calculated by dividing 60 by the mean RR interval. Heart rate variability (HRV) was measured as the standard deviation of the RR intervals.

Statistical analysis
Descriptive statistics was performed on animals' means for each group. One way repeated measures ANOVA (group) analyzed within and between subject differences and interactions. Single time point measurements or change scores were completed with one way ANOVA (group) with Least Significant Difference (LSD) post-hoc analysis. Pearson correlations were utilized for relationship values. Partial eta squared values are reported for the proportion of total variability attributed to a factor. Statistics were conducted with PASW Version 17 software (SPSS Inc, Chicago, IL, USA). Significance was measured at p < 0.05. Results are presented as means ± standard errors. The effect size of baseline body weights was large with Cohen's d = 2.8 and power was greater than 90% with sample sizes of 10-12 per group.

Animal characteristics
A summary of animal characteristics are reported in Table 1 and weekly body weights are plotted in Figure 1. A significant difference in body weight between diabetic and control animals was observed at baseline and termination time points, F(3,40) = 19.37; p < 0.001. The mean difference between baseline and termination body weights indicated that control animals gained approximately 14% of their baseline body weights compared to a 5% gain for diabetic animals, p < 0.001. Diabetic animals outweighed their control counterparts by 9-12% at termination even with the discrepancy in weight gain. Although ED rats weighed less than SD animals at the termination, the ED rats were significantly heavier than SC and EC animals, p = 0.016 and p < 0.001 respectively. However, aging accounted for 64% of the variance in body weights as evident by average weight gains of 16-57 grams in all animal groups, F(1,39) = 71.29; p < 0.001. A moderate relationship between body weight and blood glucose levels existed at the baseline r(44) = 0.56; p < 0.001 and termination r(44) = 0.41; p < 0.01. The main effect of the group on blood glucose levels was modified by aging, F(1,39) = 3.02; p < 0.04. Although all animals had an increase in blood glucose levels, ED animals had only a 10% change compared to the 21% change of SD rats. HbA1C levels were elevated in diabetic groups compared to control groups, F(3,40) = 99.27; p < 0.001. Exercise did not attenuate HbA1C levels in ED animals, p = 0.92. Hyperglycemia is a factor in left ventricular hypertrophy development and progression but a correlation between A1C levels and heart weight/body weight ratios was not observed, r(44) = -0.17 p = 0.27. Heart weight/ body weight ratios were similar for the animal groups, F(3,40) = 1.39; p = 0.26.
ECG wave amplitudes R wave amplitude R wave amplitudes were similar for groups at baseline except for the SD groups, F(3,40) = 15.16, p < 0.001. Therefore, statistical analyses used the change value for R wave amplitude to account for the difference at baseline. Gains in R wave amplitudes from baseline to termination of the experiment were only observed in the SD animals as reported in Figure 2. The SD rats had a 17% increase in R wave amplitude, suggesting left ventricular hypertrophy. A reduction in R wave amplitude was found in ED animals at termination. ED rats had change values similar to SC and EC animals, F(3,40) = 4.13, p = 0.84 and p = 0.87.

T wave amplitude
Due to the significant differences between diabetic groups at baseline, we analyzed the change of T wave amplitudes from baseline to termination ( Figure 3). T wave amplitudes changes were similar in the four groups, F(3,40) = 1.81, p = 0.16.

ECG intervals RR intervals
A significant shortening of RR intervals and increased heart rates was observed in SD group at baseline when compared to the SC group, F(1,38) = 8.83, p < 0.01 (Table 2). At termination, only a trend toward tachycardia was observed SD animals compared to SC animals, F(3,36) = 2.79, p = 0.06. Exercise did not reverse the heart rate pattern for diabetic animals, p = 0.92. Table 2 shows that all animal groups had various levels of RR interval shortening at termination, but statistically an aging effect was not observed, F(1,38) = 0.34, p = 0.56. Decreased HRV was observed in all diabetic animals, sedentary and exercised compared to control animals F (1,38) = 1662.23, p < 0.01. HRV was not improved with exercise, p = 0.91 as calculated with the termination RR intervals of the ED animals.

QRS intervals
Widening of QRS intervals, a sign of abnormal intraventricular conduction, was found at baseline in diabetic animals, F(3,40) = 3.72, p < 0.01 (Table 2). An 18% difference existed in the duration of QRS intervals between the SC and SD animals. At termination, QRS intervals increased 6-9% for SD animals in comparison to 27-29% increases in SC animals. These changes implied that aging interacted with the main effect of the group, F(1,41) = 16.78, p < 0.001 and accounted for 29% of the variability of QRS intervals at termination. Therefore, the difference between SC and SD animals at baseline was lost at termination, F(3,40) = 0.57, p = 0.64. An impact of exercise on QRS intervals was not observed in the EC and ED animals.

QT and QT c intervals
SD animals presented with QT intervals (a measure of ventricular repolarization) 31% wider than SC animals at baseline, F(3,40) = 9.37, p < 0.001 (Table 2). After 7 weeks of exercise, QT intervals in ED animals remained widened, F(3,40) = 14.13, p = 0.85. At termination, the difference between SC and SD animals decreased to 21%. Aging accounted for 20% of the variability in termination QT intervals, but the group effect of diabetes accounted for 51% of the variability. An interaction between aging and group factors was not significant, and QT intervals were not affected by exercise in control or diabetic animals, p = 0.91 and p = 0.61. QT c intervals (QT intervals corrected with RR intervals) were analyzed for an evaluation independent of heart rate. SD animals displayed significant widening of QT c intervals compared to SC animals, F(1,42) = 10.58, p < 0.001 (Table 2). Exercise did not attenuate the widening of QT c intervals in ED animals, p = 0.66. Compared to QT intervals, aging accounted for 85% of the variability in QT c intervals and 43% was due to group effects at termination. The interaction between time and group was significant, F(3,40) = 4.54, p < 0.01.

PR intervals
Baseline measurements revealed significant prolongation of PR intervals in SD animals compared to SC animals, F(1,39) = 5.40, p < 0.01 (Table 2). At termination, the difference between groups was not observed, p = 0.34. Exercise did not impact PR intervals in EC and ED animals.

Discussion
CVD risk is increased up to four-fold in people with diabetes compared with their nondiabetic counterparts [19]. Researchers are aggressively trying to identify early detection methods and explore the factors contributing to diabetes related heart disease. The ZDF rat model is routinely used to investigate physiological and molecular hypotheses regarding diabetes and its related complications. Reports are available indicating cardiac dysfunction in the ZDF rat established primarily with hemodynamic or echocardiographic measurements or through experiments on the isolated heart [20][21][22][23]. However, limited information exists about the ECG changes that occur in the ZDF rat and its response to physical activity, specifically to aerobic exercise that is commonly recommended for those at risk or diagnosed with diabetes related heart disease [24,25]. Our project aimed to address this gap in the current literature. Although we did not measure autonomic function directly, ECG indicators of autonomic dysfunction were observed in the ZDF rats at baseline with HRV alterations, tachycardia and QT interval prolongation. Autonomic innervations control HRV and cardiac function through a delicate balance of sympathetic and parasympathetic responses. Diabetes stimulates the sympathetic nervous system initially but prolonged exposure to hyperglycemia and elevated catecholamine levels cause a decrease in adrenergic receptors [26]. As diabetes creates a sympathetic predominance, it also produces a corresponding parasympathetic imbalance by denervation of the vagus nerve. Subclinical symptoms of DAN, primarily decreased HRV, are believed to appear in humans within one year of type 2 diabetes diagnosis and clinical presentations may not emerge until years into the diagnosis [27]. Bergstrom et al identified symptoms of DAN in type 1 diabetic patients with durations of diabetes as low as two months [28]. In our animal study, diminished HRV was noted in ECG recordings as early as within one week of diabetes onset for the ZDF.
Clinical DAN symptoms include resting tachycardia, exercise intolerance and orthostatic hypotension and heart rate syndromes [29]. Tachycardia was observed in diabetic animals at baseline, within a week of diabetes  onset. Tachycardia is also a common finding in humans with uncontrolled diabetes [30]. Similarly, our animals were not treated for their hyperglycemia as evident by blood glucose and HbA1C levels. The role of tachycardia is controversial because research indicates that it may be a diabetic complication or a causative factor of diabetes. Nagaya et al argue that elevated resting heart rates and systolic blood pressure increase the risk for type 2 diabetes [31]. Another work by the same group indicated that prolonged QT c intervals were also an independent risk factor for the development of diabetes [32]. Our analysis indicates that widening of QT c intervals was present with the onset of diabetes at twelve weeks of age in ZDF rats. Thus abnormalities of ventricular repolarization are present at an early stage of diabetes in this model. An earlier time study would be useful to verify the presence of tachycardia and QT c interval prolongation during a pre-diabetic state. Of interest, several ECG parameters in our study showed an impact of aging, in agreement with findings by others [33] indicating cardiac and renal changes in the ZDF rat model with aging, up to 36 months. P wave amplitudes were unaffected with the onset of diabetes. Future studies need to investigate the relationship between the dispersion of the P wave and the onset of diabetes. Obesity is commonly linked to diabetes and is reported to increase the dispersion of the P wave [34]. PR intervals were widened in diabetic animals at baseline, but normalized at termination. PR intervals are commonly associated with atrial fibrillation [35,36]; however our ECG assessments of diabetic animals did not indicate atrial dysfunction. Diabetes and exercise appeared to have no effect on T wave amplitudes. T wave amplitudes were skewed at baseline and therefore change scores were used for analysis. However, the morphology of the T wave has shown to have prognostic value for CVD, not amplitude changes [37,38].
After seven weeks of aerobic exercise, ED animals showed R wave amplitudes comparable to control animals. Only SD animals had an increase in R wave amplitude. Hyperglycemia is associated with left ventricular hypertrophy in type 2 diabetes [39]. The increase in R wave amplitude may indicate a change in left ventricular mass. With diabetes, pathological hypertrophy results from myocardial damage and fibrosis [40][41][42]. Fibrosis affects the filling and contractility of the ventricles. Subsequently, cardiac dysfunction presents as decreased activity tolerance, ejection fraction, cardiac output and heart failure [43]. However, other factors can alter R wave amplitudes, including electrical axis deviations, altered electrodes position and differences in chest wall thickness. Ideally, the ECG should be correlated with an echocardiography. Unfortunately, we were unable to use echocardiography approach, thus we recognize this as a limitation of this study. In analyzing other indications of heart hypertrophy, we found no difference in the heart weight/body weight ratio in our study. Alternatively, the LV weight/tibial length ratio has been validated as an index of cardiac hypertrophy in mature rats [44]. Due to the fact that the onset of diabetes in the ZDF rat occurs early in their age, we were restricted to using relatively young animals that were in the latter stages of their growth phase. Thus, normalizing heart size to the changing tibial bone length during growth may not have provided an accurate estimation in our study. Darmellah et al reported that normalization of the heart weight per body weight or tibial length resulted in similar measurements of cardiac hypertrophy in Goto-Kakizaki animals, another rat model of type 2 diabetes [45].
Aerobic exercise did not impact the hyperglycemia. Exercise is postulated to improve glucose uptake and decrease lipid accumulation in persons with controlled diabetes but the protective mechanism of exercise is lost if hyperglycemia persist [46]. A comparative study of fenobirate and metformin validated the role of lipid oxidation in the development and progression of diabetes related heart disease, with fenobirate decreasing triglycerides content and fibrosis in diabetic myocardium [47]. The switch in myocardial substrate from glucose to fatty acids has been shown to result in systolic and diastolic dysfunction in the ZDF model [48]. Exercise training has also been suggested to improve microcirculation through enhanced endothelial function [49] through normalizing glycemic levels. However, the severity of diabetes will determine if the body can adapt to the demands of exercise or regional flow has already been compromised beyond recovery [50]. Microcirculatory disturbances or small vessel disease may lead to declines in myocardial blood flow which could influence ECG parameters [51].
A 31% difference in QT intervals was reported between SC and SD animals at baseline. We hypothesize that cardiac remodeling was already in the process in the SD rats when we were taking the baseline measures, prior to any training, due to animals already displaying hyperglycemia and obesity associated with the genotype, in contrast to the SC rat. QTc intervals did not respond to exercise, but the chronic tachycardia shortened the intervals as a compensatory effect. Commonly, exercise may cause a decrease in heart rate and increased ventricular relaxation which presents as longer QT intervals. A study investigating the effects of a seven month endurance training program in dogs revealed an increase in QT intervals [52]. An acute resistance exercise resulted in a similar effect on QT c intervals [53]. Since our obese, diabetic animals were showing early signs of autonomic disturbance, exercise tolerance was lowered and animals required frequent rest breaks during our training program. The mode and duration of exercise might not have been sufficient for QT c interval adaptation. In a future study we will investigate longer durations of exercise training to see whether it may return diabetic animals to normal heart rates or restore HRV. Pagkalos et al reported improvements in cardiac autonomic function with six months of aerobic exercise training [54]. Another alternative is the evaluation of HRV during post-exercise recovery. Training may not affect resting HRV, but the benefits may be evident with the post-exercise recovery. This conclusion is supported by a study investigating cardiac autonomic function in women with and without diabetes [55].

Conclusions
In summary, our investigation proved that ECG alterations do occur with diabetes in the ZDF rat. The alterations include prolongation of the QT c interval and tachycardia which constitute important electrophysiologic alterations in this animal model of diabetes. These modifications coupled with high R wave amplitude illustrate the early cardiac anatomic and electrophysiologic alterations in this diabetic model. After seven weeks of exercise training, R wave amplitude changes in diabetes were similar to control animals from baseline to termination. However, aging may have an impact on several ECG parameters and the ZDF model showed changes in atrial and ventricular conduction possibly due to an interaction of aging and group effects. Future studies are needed to investigate ECG changes in the ZDF model before the onset of diabetes which will provide additional information about the use of QT c intervals and HRV in the early detection of DAN.