Skip to main content
  • Original investigation
  • Open access
  • Published:

Insulin resistance and circadian rhythm of cardiac autonomic modulation



Insulin resistance (IR) has been associated with cardiovascular diseases (CVD). Heart rate variability (HRV), an index of cardiac autonomic modulation (CAM), is also associated with CVD mortality and CVD morbidity. Currently, there are limited data about the impairment of IR on the circadian pattern of CAM. Therefore, we conducted this investigation to exam the association between IR and the circadian oscillations of CAM in a community-dwelling middle-aged sample.


Homeostasis models of IR (HOMA-IR), insulin, and glucose were used to assess IR. CAM was measured by HRV analysis from a 24-hour electrocardiogram. Two stage modeling was used in the analysis. In stage one, for each individual we fit a cosine periodic model based on the 48 segments of HRV data. We obtained three individual-level cosine parameters that quantity the circadian pattern: mean (M), measures the overall average of a HRV index; amplitude (Â), measures the amplitude of the oscillation of a HRV index; and acrophase time (θ), measures the timing of the highest oscillation. At the second stage, we used a random-effects-meta-analysis to summarize the effects of IR variables on the three circadian parameters of HRV indices obtained in stage one of the analysis.


In persons without type diabetes, the multivariate adjusted β (SE) of log HOMA-IR and M variable for HRV were -0.251 (0.093), -0.245 (0.078), -0.19 (0.06), -4.89 (1.76), -3.35 (1.31), and 2.14 (0.995), for log HF, log LF, log VLF, SDNN, RMSSD and HR, respectively (all P < 0.05). None of the IR variables were significantly associated with  or θ of the HRV indices. However, in eight type 2 diabetics, the magnitude of effect due to higher HOMA-IR on M, Â, and θ are much larger.


Elevated IR, among non-diabetics significantly impairs the overall mean levels of CAM. However, the  or θ of CAM were not significantly affected by IR, suggesting that the circadian mechanisms of CAM are not impaired. However, among persons with type 2 diabetes, a group clinically has more severe form of IR, the adverse effects of increased IR on all three HRV circadian parameters are much larger.


Insulin resistance (IR) is a precursor and a characteristic feature of type 2 diabetes [1]. IR is also associated with higher risk of cardiovascular diseases (CVD) [2, 3]. Homeostasis model assessment of IR (HOMA-IR) has been applied to quantify IR in people with or without glucose intolerance, and it has been a reliable tool in the assessment of IR, especially before the clinical diagnosis of type 2 diabetes [4, 5].

Heart rate variability (HRV), an index of cardiac autonomic modulation (CAM) [6] is associated with CVD mortality and CVD morbidity in various populations [713]. Several epidemiological studies have shown that individuals with IR or increased fasting glucose [14] have increased heart rate [1416] and reduced HRV [14, 16, 17]. Galinier et al. demonstrated that patients with hyperinsulinemia or IR had a significant decrease in HRV [18]. In the population-based Atherosclerosis Risk in communities study, a consistent association between IR, metabolic syndrome and impaired CAM has been reported [12, 1923]. However, all of these studies were based on the overall mean levels of HRV.

Several studies have described a circadian pattern of CAM, [2426] which can be quantified with a cosine periodic regression model consisting of three cosine function parameters: mean (M), amplitude (Â), and acrophase (θ). The cosine function parameter M measures the overall average of a HRV index, the  measures the amplitude of the oscillation of a HRV index, and the θ measures the clock time when the highest oscillation (amplitude) is reached. Lack of circadian variation of HRV is associated with increased vulnerability to cardiovascular events [27]. However, no study has examined the quantitative effect of insulin resistance on the three major circadian parameters that quantity the circadian pattern of CAM. Therefore, the objective of this study is to examine the effects of IR on the circadian pattern of CAM in a community-dwelling sample of non-diabetics.


Study population

For this report, we used the data collected for the Air Pollution and Cardiac Risk and its Time Course (APACR) study, which we designed to investigate the mechanisms and the time course of the adverse effects of fine particulate matter (PM2.5) on cardiac electrophysiology, blood coagulation, and systemic inflammation. Recruitment methods for the APACR study have been published elsewhere [28, 29]. All study participants were recruited from communities in central Pennsylvania, primarily from the Harrisburg metropolitan area. The inclusion criteria for the study included nonsmoking adults, ≥ 45 years old, who had not been diagnosed with severe cardiac problems (defined as diagnosed valvular heart disease, congenital heart disease, acute myocardial infarction or stroke within 6 months, or congestive heart failure). Approximately 75% of the individuals who were contacted and who met our inclusion criteria were enrolled in the APACR study. Our targeted sample size was 100 individuals, and we enrolled and examined 102 individuals in the APACR study.

Study participants were examined in the General Clinical Research Center (GCRC) during the morning hours of Day 1, between 8:00 AM and 10:00 AM. After administering the informed consent, the participants completed a health history questionnaire. A trained research nurse measured seated blood pressure three times, height, weight, and drew 50 ml of blood for blood biomarker assays according to the blood sample preparation protocols. A trained investigator connected the PM2.5 and Holter electrocardiography (ECG) recorders. Participants were given an hourly activity log to record special events that occurred during 24 hours, including outdoor activities, exposure to traffic on the street, travel in an automobile, and any physical activities. Before participants were released, they received detailed instructions on how to operate both monitors. The entire session lasted approximately 1 hour. The next morning (Day 2), the participants came back to the GCRC to disconnect the PM2.5 and Holter monitors, return the completed activity log, have another 50 ml of blood drawn, and provide a urine sample. Then, an exercise echocardiography was performed on each participant according to a standardized protocol (Bruce protocol) to measure the participant's ventricular function and structure. The entire Day 2 session lasted approximately 1 hour and 45 minutes. The study protocol was approved by the Penn State University College of Medicine Institutional Review Board (IRB). All participants gave written informed consent prior to their participation in the study. Each participant was reimbursed with $50.00, a breakfast certificate in the hospital cafeteria, and the mileage of the transportation required for participating in the study.

Assessment of insulin resistance: insulin, glucose, and HOMA-IR

During the Day 1 visit, after an overnight fast, blood was drawn from each participant. Fasting glucose and insulin were measured by the Penn State College of Medicine GCRC central laboratory. HOMA-IR was calculated as [fasting insulin (μU/mL) Ø fasting glucose (mmol/L)/22.5], with higher HOMA-IR values indicative of more IR [5].

Continuous Holter ECG recording

Continuous ambulatory ECG

A high-fidelity (sampling frequency 1,000 Hz) 12-lead HScribe Holter System (Mortara Instrument, Inc., Milwaukee, WI) was used to collect the 24-hour Holter beat-to-beat ECG data. The high-fidelity ECG significantly increases the resolution and enhances the accuracy of various wave form measurements. All Holter recordings started between 9:00 AM and 10:00 AM. The Holter ECG data were scanned to a designated computer for offline processing by an experienced investigator using specialized SuperECG software (Mortara Instrument, Inc.). The standardized operation procedures (SOPs) for the APACR study developed by the study investigators were followed rigorously in the data collection, offline processes, HRV analysis and interpretation processes. Briefly, the Holter ECG Data Collection and Analysis Procedures were followed to prepare, hook up, calibrate, and start the Holter digital recorder. After 24 hours of recording, a trained investigator followed the SOP to retrieve and archive the beat-to-beat ECG data for offline processing. The main objective of the offline processing was to verify the Holter-identified ECG waves and to identify and label additional electronic artifacts and arrhythmic beats in the ECG recording. Finally, a single research investigator performed beat-to-beat HRV analysis using the normal beat-to-beat RR interval data.

Cardiac autonomic modulation (CAM) measures

The entire 24-hour normal beat-to-beat RR interval data were divided into 30-minute segments of RR data. Thus, each individual provided 48 segments of 30-minute RR data. The RR data for HRV analysis were processed according to current recommendations [9]. Within each segment, any RR interval <400 ms, >2000 ms, or where the ratio from two adjacent RR intervals was <0.80 or >1.20 were excluded from the HRV analysis. The time- and frequency-domain HRV analysis were performed on the remaining normal RR interval data if the total length of such normal RR intervals was greater than 20 minutes (67% of original data), using the HRV Analysis Software v1.1 [30]. When performing frequency-domain HRV analysis, we used Fast Fourier Transformation (FFT). Briefly, the adjacent RR interval data were interpolated using a piecewise cubic spline approach, with a 2 Hz sampling rate. The FFT was performed on the equidistantly interpolated RR time series. We used a second order polynomial model to remove the slow non-stationary trends of the HRV signal. The following HRV indices were calculated as measures of CAM: standard deviation of all RR intervals (SDNN, ms), square root of the mean of the sum of the squares of differences between adjacent RR intervals (RMSSD, ms), power in the high frequency range (0.15-0.40 Hz, HF), power in the low frequency range (0.04-0.15 Hz, LF), the ratio of LF to HF (LF/HF), and in the very low frequency power (0.00-0.04 Hz, VLF). Following current recommendations [9], we performed logarithmic transformations on HF and LF prior to statistical analysis.


A standardized questionnaire administered on Day 1 to the participants was used to collect the following individual-level information: (1) demographic variables, including age, race, sex, and highest education level; (2) medication uses, including anti-hypertensive medication and glucose lowering medications; and (3) physician diagnosed chronic disease history, including coronary heart disease (e.g., revascularization procedures and myocardial infarction), hypertension, and diabetes. Body weight and standing height were obtained in the morning of Day 1 and were used to calculate body mass index (BMI) as weight (kg)/height (m2). Seated blood pressure was obtained three times after 5 minutes of resting. The averages of the second and third measures of seated systolic and diastolic blood pressures were used to represent a participant's blood pressure levels. Hypertension was defined as systolic blood pressure ≥ 140 mmHg, diastolic blood pressure ≥ 90 mmHg, or physician diagnosed hypertension and currently using anti-hypertensive medication. History of CVD was defined as having a physician diagnosed history of myocardial infarction or coronary arterial revascularization procedures. Type 2 diabetes was defined as fasting glucose > 126 mg/dL, physician diagnosed type 2 diabetes, or currently using glucose lowering medication.

Statistical analysis

From the 102 individuals, we excluded 8 participants with type 2 diabetes. As a result, this report used the data from the remaining 94 individuals. Each individual contributed up to 48 segments of 30-minute RR interval data. After excluding segments where the total length of RR interval data was less than 20 minutes, we analyzed 4404 segments of 30-minute HRV data from 94 individuals.

A two-stage analysis was performed to assess the relationship between IR measures and the circadian pattern of HRV. At the first stage, for each individual we fit the HRV data based on all available 30-minute segments to a cosine periodic regression model [31] using nonlinear least squares: HRVi(t) = Mi + Ai•cos [2π •(t-θi)/T] + εi, i = 1, ..., 94, in which Mi is the daily average of HRV of the ith subject, Ai is the amplitude of HRV of the ith subject around Mi, t is the time-specific segment order number, T is the total number of 30-minute segments in 24 hours, θi is the acrophase (the lag from the reference time point (9:00 AM) to the time of the zenith of the cosine curve fit to the data of the ith subject), and εi is the error term of the ith subject. One unit of t corresponds to 30 minutes, with 1 indicating 9:00 AM to 9:30 AM, 2 indicating 9:30 AM to 10:00 AM, etc. Thus, from the above described cosine model, we obtained the estimated individual-level cosine periodic regression parameters, namely the M, Â, and θ to quantify the periodicity of the HRV variables. At the second stage, we used random-effects meta-analyses to obtain overall estimates of M, Â and θ, and their 95% confident intervals (CIs) to assess the associations between IR and the three components of the circadian pattern of HRV [32]. Because HOMA-IR and fasting insulin were not normally distributed, we used a logarithmic transformation for these two IR variables before performing the analysis. To visualize the impact of elevated IR on the three circadian parameters, we plotted the population level circadian pattern of HF, and superimposed the circadian pattern predicted HF from one standard deviation increase in HOMA-IR. To confirm the associations between the IR and the means from the random-effects meta-analysis models, we also applied linear mixed-effects models, specifying a first-order autoregressive covariance structure, to assess the associations between HRV variables and HOMA-IR using 24-hour data, analyzing the daytime (9 AM to 9 PM) and nighttime (9 PM to 9 AM next day) data separately. This model will also enable us to examine whether IR affect daytime and nighttime autonomic modulation differently. In this approach we ignored the assumption of the specific cosine form for the data and treated the HRV variables from each 30-minute segment as repeated measures. All analyses were performed using SAS version 9.1 software (SAS Institute Inc., Cary, NC, USA).


Study population characteristics

The main characteristics of the study population are shown in Table 1. From the first 102 APACR study enrollees, we excluded 8 individuals who had a history of type 2 diabetes. The sample size for this analysis was 94 individuals. The mean (SD) age of the entire study population was 56.5 (7.8) years old, with 37% male and 77% white. The prevalence of CVD was 7.4%. In addition, the mean (SD) for BMI (kg/m2), insulin (μU/mL), glucose (mg/dL), and HOMA-IR of the entire cohort were 27.1 (4.9), 6.6 (6.7), 84.4 (11.3), and 1.45 (1.73) respectively. The means (SD) for the HRV indices: log HF, log LF, LF/HF ratio, log VLF, SDNN, RMSSD, and HR were 4.4 (0.9), 5.4 (0.8), 3.8 (2.2), 6.9 (0.6), 61.8 (18.4), 29.6 (14.8), and 76.5 (9.7), respectively.

Table 1 The study population demographic characteristics and summaries of HRV indices.

Circadian pattern of HRV variables

The estimated means of the three cosine periodic regression parameters and their 95% CIs for each of the HRV indices are presented in Table 2. The acrophases of log HF and RMSSD (both are reflective of the vagal modulation) were very similar, 4:00 AM and 4:15 AM, respectively. The acrophases of log LF and SDNN (both are reflective of sympathetic and vagal modulation) were also very similar - the point estimates are around 6:00 AM ± 45 minutes. The acrophase of log VLF was 8:00 AM, LF/HF ratio was 6:00 PM, and that of heart rate was 2:00 PM. The tests of zero amplitude were highly significant (P value < 0.0001) for every HRV variable, suggesting a significant circadian variation.

Table 2 Estimated means of the three cosine periodic regression parameters and their 95% CIs for each HRV index.

Insulin resistance and HRV circadian pattern

The associations between the IR measurements and cosine parameters of the circadian pattern (M, Â, and θ) estimated from the entire sample using random-effects meta-analysis models are presented in Table 3. Two different random-effects meta-analysis regression models were reported: Model 1, unadjusted; Model 2, adjusted for age sex, race, hypertension, and history of CVD. In general, higher levels of IR, especially IR measured by HOMA-IR, was significantly associated with lower M of all HRV indices and higher HR in both unadjusted and multivariable adjusted models, except that LF/HF ratio was not significantly associated with any measure of IR. In this healthy non-diabetes sample, IR was not significantly associated with the  or θ of any HRV indices. Log VLF was also included in Table 3 as a HRV variable that may be reflective of nocturnal sympathetic hyperactivity. The association between IR measures and VLF is very similar to that of HF and LF - elevated insulin resistance is only associated with lower mean levels of VLF, but not with the amplitude or the acrophase.

Table 3 The association between cosine parameters of circadian pattern with insulin resistance measures based on random-effects meta-analysis

The circadian variation of Log HF calculated from the entire population is graphically presented in Figure 1. The predicted circadian pattern of Log HF due to a one SD increase in HOMA-IR is superimposed in Figure 1 to graphically illustrate the impact of elevated IR on the circadian pattern of log HF. It clearly show, similar to the numeric numbers in Table 3, that one SD increase in IR is associated with lower overall mean [β (M)] of Log HF, but not with the other two circadian parameters (Â and θ).

Figure 1
figure 1

Multivariable adjusted (Model 2) cosine periodic model estimates for log HF over the clock time at the population mean and one SD above the mean in log HOMA-IR. Abbreviations: Â, amplitude; HF, high frequency power; HOMA-IR, homeostasis model assessment of insulin resistance; log, logarithm; M, mean; SD, standard deviation; θ, acrophrase time.

We also performed additional analyses examining the associations between HOMA-IR and HRV cosine parameters of circadian pattern in 8 individuals with physician-diagnosed type 2 diabetes. These analyses were preformed in an unadjusted model because small sample size. The results are presented in Table 4. Compared to that regression coefficients in Table 3 Model 1, the magnitude of effect due to higher Log HOMA-IR on the mean (M), amplitude (Â), and acrophase time (θ), are much larger among persons with type 2 diabetes, a group that clinically have more severe form of insulin resistance. Because of the small sample size, most of the statistical tests of the regression coefficients in Table 4 did not reach traditional significant level (P value < 0.05).

Table 4 The association between cosine parameters of circadian pattern with Log HOMA-IR based on random-effects meta-analysis of 8 physician diagnosed type 2 diabetes individuals.

Ignoring the specific cosine form for the data and treating HRV variables from each 30-minute segment as repeated measures, the regression coefficients, SE, and P value relating HOMA-IR and HRV indices according to daytime (9 AM to 9 PM) and nighttime (9 PM to 9 AM next day) are presented in Table 5. In general, these results confirmed that higher IR is associated with lower HRV indices, and the patterns of association are similar daytime and night time.

Table 5 The regression coefficients (β) estimates and standard error (SE) of Log-HOMA-IR on HRV indices using 24-hour data, daytime (9 AM to 9 PM) and nighttime (9 PM to 9 AM next day) data separately from linear mixed-effects models analysis.


The estimated means and their 95% CIs of the three cosine periodic regression parameters presented in Table 2 clearly show a distinctive 24-hour based circadian rhythm of HRV indices. For instance, using Log HF as an example, a recognized HRV marker of parasympathetic modulation, had an overall mean of 4.36 ms2 and reached the acrophase (highest amplitude of  = 0.60 ms2) around 4:00 AM. Furthermore, the tests of zero amplitude were highly significant (P value < 0.0001) for every HRV variable, suggesting a significant circadian variation. The mean (M), amplitude (Â), and acrophase time (θ) together define the circadian pattern of HRV, which is reflective of the balance of two opposing autonomic modulation branches, namely parasympathetic and sympathetic modulations.

Past epidemiological studies performed in different populations investigating the relationship between IR and the overall mean levels of HRV indices have shown that elevated IR has a significant adverse effect on the overall mean levels of CAM [17, 33]. Pikkujämsä et al. demonstrated that subjects with IR syndrome have reduced HRV indices [17]. In addition, Perciaccante et al. confirmed that impaired autonomic activity was present in IR patients [33]. On the other hand, other studies have shown that individuals without type 2 diabetes or offspring of type 2 diabetic persons, even when their glucose was in the normal range, had lower HRV indices [14, 16, 34]. Stein et al. found that in non-diabetic individuals with normal fasting glucose levels, higher glucose levels [14] and higher HOMA-IR [16] were associated with faster heart rate and lower HRV. In addition, Fiorentini et al. showed that higher HOMA-IR was associated with autonomic impairment among the offspring of individuals with type 2 diabetes, even when they had no history of type 2 diabetes, normal glucose tolerance test, and no history of hypertension [34]. Similar to these previous findings, in this community-dwelling sample of healthy non-smokers without type 2 diabetes, we found that a higher level of IR is significantly associated with lower mean levels of HRV. The pattern of these associations remained significant after adjustment for age, sex, race, hypertension, and history of CVD. This conclusion is also supported in the analysis treating HRV variables from each 30-minute segment as repeated measures, and stratified according to daytime (9 AM to 9 PM) and nighttime (9 PM to 9 AM next day). Moreover, the daytime and night time stratified model suggest a similar daytime and nighttime association between IR and CAM.

To our knowledge, none of the previously published studies investigated the effects of IR on the other two parameters of HRV circadian pattern, the amplitude and the acrophase. In this study, we used a two-stage analytic approach to first obtain individual-level cosine function parameters (M, Â, and θ) from 24 hours of HRV data, and then used a random-effects meta-analysis regression model to investigate the impact of IR on each of the circadian cosine parameters. This approach allows us to estimate the impact of IR, not only on overall mean levels of HRV variables as in previously published studies, but it also enabled us to examine the impact of IR on the amplitude and acrophase time of HRV. In our data obtained from obviously healthy non-diabetics and non-smokers, IR was not associated with impaired circadian variation of HRV as suggested by lack of association between amplitude and acrophase time of HRV index and IR.

In humans, circadian rhythms in blood pressure [35, 36], and in heart rate have been intensively studied, mainly due to the increased cardiovascular death reported during the morning hours [37, 38]. There are a large number of different biological circadian rhythms, including numerous aspects of cardiovascular functions [39]. These are generated by an endogenous oscillator, composed of a central clock which resides in the hypothalamic suprachiasmic nucleus (SCN), and peripheral clock in peripheral tissues [40]. These two have an important role in the timing and organization/coordination of sleep with other physiological rhythms, such as the balance of autonomic modulation [41]. The central clock is composed of self-sustaining single cell circadian oscillators which produce circadian signals [40, 42]. These oscillators are not precisely 24 hours, so they adjust to the external light/dark signal from the retinohypothalamic tract [40, 42]. These signals give information about the 24-hour day pattern, the relative proportion of light/dark, and the amplitude of light intensity that tracks the seasons [43]. Most of the peripheral tissue cells have similar oscillations like the SCN [44].

Various animal studies have examined the impact of disruption clock genes [4548]. Turek et al. showed that homozygous Clock mutant mice had a loss of sense mutation, lost circadian rhythmicity, became obese, and developed hyperglycemia among other health problems [32]. Rudic et al. showed that mutations in Bmal and Clock modified circadian variation in glucose and triglycerides and influenced the development of IR [46]. Marcheva et al. found that phase of oscillation of islet genes involved in growth, glucose metabolism and insulin signaling is delayed in circadian mutant mice, and both Clock and Bmal1 mutants show reduce insulin secretion, and impaired glucose tolerance [47].

Researchers have been concerned over the disruption of these circadian rhythms in the human clock, which regulate the metabolism in time with day length [49, 50]. For example, glucose and insulin have been shown to exhibit diurnal variation in humans [50]. The normal physiological variation of these variables between day and night could be lost (with an increase over the 24 hours) contributing to the development of CVD [50] or IR.

Earlier studies have investigated the circadian rhythm of HRV in humans [15, 51, 52]. Huikuri et al. observed that the circadian rhythm of HRV (measured from a 24-hour ECG) in healthy subjects had a maximum occurring early in the morning before waking up, reflecting higher parasympathetic tone, but decreased abruptly during the hours after waking [51]. Malpa et al. found that HRV variation over 24 hours in alcoholics and persons with diabetes was significantly reduced compared to normal controls [25]. Huikuri et al. demonstrated that there was a circadian pattern of HRV in both survivors of cardiac arrest and controls, with acrophase occurring in the morning hours before arousal, followed by a decrease in HRV after arousal; however there was low vagal modulation in survivors of cardiac arrest in the morning [52].

To our knowledge, the effects of IR on the circadian rhythm (amplitude and acrophase time) of HRV have not been reported in the literature. The aim of this study was to investigate the association between IR and the circadian pattern of CAM assessed by three circadian cosine parameters of HRV. In our data, higher levels of IR, as assessed by HOMA-IR, fasting insulin, or fasting glucose, were significantly associated with lower mean levels of HRV. However, IR in this obviously healthy non-diabetic and nonsmoking population sample was not associated with impaired circadian variation of HRV indices, i.e., lower amplitude and acrophase time of any HRV index. This lack of significant effects on the circadian pattern of HRV variation, indicates that the burden of IR in our healthy study population is not large enough to exert its effects yet. We speculate that with prolonged burden of elevated subclinical IR, possibly long before full manifestation of clinical diabetes, all three circadian parameters of CAM may be adversely affected. This is supported by the analysis of the associations between Log HOMA-IR and HRV cosine parameters of circadian pattern in 8 individuals with physician-diagnosed type 2 diabetes (Table 4). Specifically, the adverse effects of increased IR on all three HRV circadian parameters, namely the mean (M), amplitude (Â), and acrophase time (θ), are much larger among persons with type 2 diabetes, a group that clinically have more severe form of IR.

Future longitudinal studies are needed to investigate this hypothesis. From a practical stand point, the results of lack of effects on the circadian pattern of HRV also suggest that early interventions to enhance insulin sensitivity, such as physical activity, may be helpful to preserve the circadian pattern of CAM.

A few limitations of this study should be recognized. First, our cross-sectional data limited our ability to establish a temporal relationship. Thus, our results support the need for a longitudinal study to investigate the complex relationships between IR and CAM assessed by circadian variation of HRV indices. Second, we only collected data from nonsmokers who did not have physician diagnosed heart disease in the 6 months prior to participating in our study. Therefore, our findings cannot be generalized to smokers or individuals with recent acute cardiac events. However, there are several strengths worth mentioning. First, the characteristics of this population were very similar to those in the general population. We analyzed several variables that are characteristics of IR, including insulin, glucose, and HOMA-IR. These variables, especially HOMA-IR, are considered reliable quantitative measures of IR before clinical diabetes [4, 5]. Second, distinct from previous publications, which primarily analyzed short time HRV variables, [17, 53] we utilized cosine periodic functions to fit individual-level models to the 24 hours of HRV data to estimate the individual-level mean, as well as amplitude and crescent time, of each HRV variable. This enabled us to examine not only the impact of various independent variables, such as HOMA-IR, on the mean levels of HRV variables, but also allowed us to examine the impact of such variables on the other circadian cosine parameters of the HRV indices.


In this cross-sectional study of community dwelling non-smokers, after excluding the few individuals with a history of diabetes, subclinical IR only exhibited significant adverse effects on the mean levels of HRV. The subclinical IR in this sample did not show an adverse impact on the amplitude and the acrophase time of the circadian cosine parameters. It might be anticipated that with fully manifested diabetes, the consequence of prolonged IR, the entire circadian variation of CAM will be adversely affected. From a practical stand point, these data suggest the potential for early intervention to improve insulin sensitivity and to preserve the circadian pattern of CAM, which is protective from cardiovascular adverse events.


  1. Lindmark S, Lonn L, Wiklund U, Tufvesson M, Olsson T, Eriksson JW: Dysregulation of the autonomic nervous system can be a link between visceral adiposity and insulin resistance. Obes Res. 2005, 13: 717-728. 10.1038/oby.2005.81.

    Article  CAS  PubMed  Google Scholar 

  2. Pyorala K, Savolainen E, Kaukola S, Haapakoski J: Plasma insulin as coronary heart disease risk factor: relationship to other risk factors and predictive value during 9 1/2-year follow-up of the Helsinki Policemen Study population. Acta Med Scand Suppl. 1985, 701: 38-52.

    CAS  PubMed  Google Scholar 

  3. Casassus P, Fontbonne A, Thibult N, Ducimetiere P, Richard JL, Claude JR, Warnet JM, Rosselin G, Eschwege E: Upper-body fat distribution: a hyperinsulinemia-independent predictor of coronary heart disease mortality. The Paris Prospective Study. Arterioscler Thromb. 1992, 12: 1387-1392.

    Article  CAS  PubMed  Google Scholar 

  4. Bonora E, Targher G, Alberiche M, Bonadonna RC, Saggiani F, Zenere MB, Monauni T, Muggeo M: Homeostasis model assessment closely mirrors the glucose clamp technique in the assessment of insulin sensitivity: studies in subjects with various degrees of glucose tolerance and insulin sensitivity. Diabetes Care. 2000, 23: 57-63. 10.2337/diacare.23.1.57.

    Article  CAS  PubMed  Google Scholar 

  5. Bonora E, Formentini G, Calcaterra F, Lombardi S, Marini F, Zenari L, Saggiani F, Poli M, Perbellini S, Raffaelli A, Cacciatori V, Santi L, Targher G, Bonadonna R, Muggeo M: HOMA-estimated insulin resistance is an independent predictor of cardiovascular disease in type 2 diabetic subjects: prospective data from the Verona Diabetes Complications Study. Diabetes Care. 2002, 25: 1135-1141. 10.2337/diacare.25.7.1135.

    Article  PubMed  Google Scholar 

  6. Kleiger RE, Stein PK, Bigger JT: Heart rate variability: measurement and clinical utility. Ann Noninvasive Electrocardiol. 2005, 10: 88-101. 10.1111/j.1542-474X.2005.10101.x.

    Article  PubMed  Google Scholar 

  7. Kleiger RE, Miller JP, Bigger JT, Moss AJ: Decreased heart rate variability and its association with increased mortality after acute myocardial infarction. Am J Cardiol. 1987, 59: 256-262. 10.1016/0002-9149(87)90795-8.

    Article  CAS  PubMed  Google Scholar 

  8. Zuanetti G, Neilson JM, Latini R, Santoro E, Maggioni AP, Ewing DJ: Prognostic significance of heart rate variability in post-myocardial infarction patients in the fibrinolytic era. The GISSI-2 results. Gruppo Italiano per lo Studio della Sopravvivenza nell' Infarto Miocardico. Circulation. 1996, 94: 432-436.

    Article  CAS  PubMed  Google Scholar 

  9. Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology: Heart rate variability: standards of measurement, physiological interpretation and clinical use. Circulation. 1996, 93: 1043-1065.

    Article  Google Scholar 

  10. Tsuji H, Larson MG, Venditti FJ, Manders ES, Evans JC, Feldman CL, Levy D: Impact of reduced heart rate variability on risk for cardiac events. The Framingham Heart Study. Circulation. 1996, 94 (11): 2850-2855.

    Article  CAS  PubMed  Google Scholar 

  11. Stein PK, Barzilay JI, Chaves PH, Domitrovich PP, Gottdiener JS: Heart rate variability and its changes over 5 years in older adults. Age Ageing. 2009, 38 (2): 212-218. 10.1093/ageing/afn292.

    Article  PubMed Central  PubMed  Google Scholar 

  12. Liao D, Barnes RW, Chambless LE, Simpson RJ, Sorlie P, Heiss G: Population based study of heart rate variability and prevalent myocardial infarction -- The ARIC Study. J Electrocardiol. 1996, 29: 189-198. 10.1016/S0022-0736(96)80082-3.

    Article  CAS  PubMed  Google Scholar 

  13. Dekker JM, Crow RS, Folsom AR, Hannan PJ, Liao D, Swenne CA, Schouten EG: Low heart rate variability in a two minute rhythm strip predicts risk of coronary heart disease and mortality from several causes -- The ARIC Study. Circulation. 2000, 102: 1239-1244.

    Article  CAS  PubMed  Google Scholar 

  14. Stein PK, Barzilay JI, Domitrovich PP, Chaves PM, Gottdiener JS, Heckbert SR, Kronmal RA: The relationship of heart rate and heart rate variability to non-diabetic fasting glucose levels and the metabolic syndrome: the Cardiovascular Health Study. Diabet Med. 2007, 24: 855-863. 10.1111/j.1464-5491.2007.02163.x.

    Article  CAS  PubMed  Google Scholar 

  15. Troisi RJ, Weiss ST, Parker DR, Sparrow D, Young JB, Landsberg L: Relation of obesity and diet to sympathetic nervous system activity. Hypertension. 1991, 17: 669-77.

    Article  CAS  PubMed  Google Scholar 

  16. Stein PK, Barzilay JI, Chaves PH, Traber J, Domitrovich PP, Heckbert SR, Gottdiener JS: Higher levels of inflammation factors and greater insulin resistance are independently associated with higher heart rate and lower heart rate variability in normoglycemic older individuals: the Cardiovascular Health Study. J Am Geriatr Soc. 2008, 56: 315-321. 10.1111/j.1532-5415.2007.01564.x.

    Article  PubMed  Google Scholar 

  17. Pikkujamsa SM, Huikuri HV, Airaksinen KE, Rantala AO, Kauma H, Lilja M, Savolainen MJ, Kesaniemi YA: Heart rate variability and baroreflex sensitivity in hypertensive subjects with and without metabolic features of insulin resistance syndrome. Am J Hypertens. 1998, 11: 523-531. 10.1016/S0895-7061(98)00035-1.

    Article  CAS  PubMed  Google Scholar 

  18. Galinier M, Fourcade J, Ley N, Boveda S, Solera S, Solera ML, Massabuau P, Elhabaj S, Fauvel JM, Valdiguie P, Bounhoure JP: Hyperinsulinism, heart rate variability and circadian variation of arterial pressure in obese hypertensive patients. Arch Mal Coeur Vaiss. 1999, 92: 1105-1109.

    CAS  PubMed  Google Scholar 

  19. Liao D, Cai J, Brancati F, Crow R, Barnes RW, Tyroler HA, Heiss G: Association of vagal tone with serum insulin, glucose and diabetes mellitus -- The ARIC Study. Diabetes Res Clin Pract. 1995, 30: 211-221. 10.1016/0168-8227(95)01190-0.

    Article  CAS  PubMed  Google Scholar 

  20. Liao D, Cai J, Rosamond W, Barnes R, Hutchinson R, Whitsel E, Rautaharju P, Heiss G: Cardiac autonomic function and incident CHD: a population based case-cohort study -- The ARIC Study. Am J Epidemiol. 1997, 145: 696-706.

    Article  CAS  PubMed  Google Scholar 

  21. Liao D, Sloan RP, Cascio WE, Folsom AR, Liese AD, Evans GW, Cai J, Sharrett AR: The multiple metabolic syndrome is associated with lower heart rate variability -- The ARIC Study. Diabetes Care. 1998, 21: 2116-2122. 10.2337/diacare.21.12.2116.

    Article  CAS  PubMed  Google Scholar 

  22. Liao D, Carnethon MR, Evans GW, Cascio WE, Heiss G: Lower Heart Rate Variability Is Associated with the Development of Coronary Heart Disease in Individuals with Diabetes - The ARIC Study. Diabetes. 2002, 51: 3524-3531. 10.2337/diabetes.51.12.3524.

    Article  CAS  PubMed  Google Scholar 

  23. Carnethon MR, Golden SH, Folsom AR, Haskell W, Liao D: Prospective investigation of autonomic nervous system function and the development of type 2 diabetes: the Atherosclerosis Risk in Communities study, 1987-1998. Circulation. 2003, 107: 2190-5. 10.1161/01.CIR.0000066324.74807.95.

    Article  PubMed  Google Scholar 

  24. Furlan R, Guzzetti S, Crivellaro W, Dassi S, Tinelli M, Baselli G, Cerutti S, Lombardi F, Pagani M, Malliani A: Continuous 24-hour assessment of the neural regulation of systemic arterial pressure and RR variabilities in ambulant subjects. Circulation. 1990, 81: 537-547.

    Article  CAS  PubMed  Google Scholar 

  25. Malpas SC, Purdie GL: Circadian variation of heart rate variability. Cardiovasc Res. 1990, 24: 210-213. 10.1093/cvr/24.3.210.

    Article  CAS  PubMed  Google Scholar 

  26. Lombardi F, Sandrone G, Mortara A, La Rovere MT, Colombo E, Guzzetti S, Malliani A: Circadian variation of spectral indices of heart rate variability after myocardial infarction. Am Heart J. 1992, 123: 1521-1529. 10.1016/0002-8703(92)90804-5.

    Article  CAS  PubMed  Google Scholar 

  27. Kardelen F, Akcurin G, Ertug H, Akcurin S, Bircan I: Heart rate variability and circadian variations in type 1 diabetes mellitus. Pediatr Diabetes. 2006, 7: 45-50. 10.1111/j.1399-543X.2006.00141.x.

    Article  PubMed  Google Scholar 

  28. Liao D, Shaffer ML, Rodriguez-Colon S, He F, Li X, Wolbrette DL, Yanosky J, Cascio WE: Acute adverse effects of fine particulate air pollution on ventricular repolarization. Environ Health Perspect. 2010, 118: 1010-1015. 10.1289/ehp.0901648.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  29. He F, Shaffer M, Li X, Rodriguez-Colon S, Wolbrette D, Williams R, Cascio W, Liao D: Individual-level PM(2.5) exposure and the time course of impaired heart rate variability: The APACR study. J Expo Sci Environ Epidemiol. []

  30. Niskanen JP, Tarvainen MP, Ranta-Aho PO, Karjalainen PA: Software for advanced HRV analysis. Comput Methods Programs Biomed. 2004, 76 (1): 73-81. 10.1016/j.cmpb.2004.03.004.

    Article  PubMed  Google Scholar 

  31. D'Negri CE, Marelich L, Vigo D, Acunzo RS, Girotti LA, Cardinali DP, Siri LN: Circadian periodicity of heart rate variability in hospitalized angor patients. Clin Auton Res. 2005, 15: 223-32.

    Article  PubMed  Google Scholar 

  32. DerSimonian R, Laird N: Meta-analysis in clinical trials. Control Clin Trials. 1986, 7: 177-88. 10.1016/0197-2456(86)90046-2.

    Article  CAS  PubMed  Google Scholar 

  33. Perciaccante A, Fiorentini A, Paris A, Serra P, Tubani L: Circadian rhythm of the autonomic nervous system in insulin resistant subjects with normoglycemia, impaired fasting glycemia, impaired glucose tolerance, type 2 diabetes mellitus. BMC Cardiovasc Disord. 2006, 6: 19-10.1186/1471-2261-6-19.

    Article  PubMed Central  PubMed  Google Scholar 

  34. Fiorentini A, Perciaccante A, Paris A, Serra P, Tubani L: Circadian rhythm of autonomic activity in non diabetic offsprings of type 2 diabetic patients. Cardiovasc Diabetol. 2005, 4: 15-10.1186/1475-2840-4-15.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Felício JS, de Souza ACCB, Kohlmann N, Kohlmann O, Ribeiro AB, Zanella MT: Nocturnal blood pressure falls as predictor of diabetic nephropathy in hypertensive patients with type 2 diabetes. Cardiovascular Diabetol. 2010, 9: 36.

    Article  Google Scholar 

  36. Gupta AK, Cornelissen G, Greenway FL, Dhoopati V, Halberg F, Johnson WD: Abnormalities in circadian blood pressure variability and endothelial function: pragmatic markers for adverse cardiometabolic profiles in asymptomatic obese adults. Cardiovascular Diabetol. 2010, 9: 58-10.1186/1475-2840-9-58.

    Article  Google Scholar 

  37. Cygankiewicz I, Krzysztof Wranicz J, Bolinska H, Zaslonka J, Zareba W: Circadian changes in heart rate turbulence parameters. J Electrocardiol. 2004, 37: 297-303. 10.1016/j.jelectrocard.2004.07.008.

    Article  PubMed  Google Scholar 

  38. Tofler GH, Gebara OC, Mittleman MA, Taylor P, Siegel W, Venditti FJ, Rasmussen CA, Muller JE: Morning peak in ventricular tachyarrhythmias detected by time of implantable cardioverter/defibrillator therapy. The CPI Investigators. Circulation. 1995, 92: 1203-1208.

    Article  CAS  PubMed  Google Scholar 

  39. Prasai MJ, George JT, Scott EM: Molecular clocks, type 2 diabetes and cardiovascular disease. Diab Vasc Dis Res. 2008, 5: 89-95. 10.3132/dvdr.2008.015.

    Article  PubMed  Google Scholar 

  40. Dunlap JC: Molecular bases for circadian clocks. Cell. 1999, 96: 271-290. 10.1016/S0092-8674(00)80566-8.

    Article  CAS  PubMed  Google Scholar 

  41. Glotzbach SF, Edgar DM, Boeddiker M, Ariagno RL: Biological rhythmicity in normal infants during the first 3 months of life. Pediatrics. 1994, 94: 482-488.

    CAS  PubMed  Google Scholar 

  42. Dunlap JC: Physiology. Running a clock requires quality time together. Science. 2006, 311 (5758): 184-186. 10.1126/science.1122839.

    Article  CAS  PubMed  Google Scholar 

  43. Lincoln GA, Andersson H, Loudon A: Clock genes in calendar cells as the basis of annual timekeeping in mammals--a unifying hypothesis. J Endocrinol. 2003, 179: 1-13. 10.1677/joe.0.1790001.

    Article  CAS  PubMed  Google Scholar 

  44. Ptitsyn AA, Zvonic S, Conrad SA, Scott LK, Mynatt RL, Gimble JM: Circadian clocks are resounding in peripheral tissues. PLoS Comput Biol. 2006, 2 (3): e16-10.1371/journal.pcbi.0020016.

    Article  PubMed Central  PubMed  Google Scholar 

  45. Turek FW, Joshu C, Kohsaka A, Lin E, Ivanova G, McDearmon E, Laposky A, Losee-Olson S, Easton A, Jensen DR, Eckel RH, Takahashi JS, Bass J: Obesity and metabolic syndrome in circadian Clock mutant mice. Science. 2005, 308: 1043-1045. 10.1126/science.1108750.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  46. Rudic RD, McNamara P, Curtis A, Boston RC, Panda S, Hogenesch JB, Fitzgerald GA: BMAL1 and CLOCK, two essential components of the circadian clock, are involved in glucose homeostasis. PLoS Biol. 2004, 2: e377-10.1371/journal.pbio.0020377.

    Article  PubMed Central  PubMed  Google Scholar 

  47. Marcheva B, Ramsey KM, Buhr ED, Kobayashi Y, Su H, Ko CH, Ivanova G, Omura C, Mo S, Vitaterna MH, Lopez JP, Philipson LH, Bradfield CA, Crosby SD, JeBailey L, Wang X, Takahashi JS, Bass J: Disruption of the clock components CLOCK and BMAL1 leads to hypoinsulinaemia and diabetes. Nature. 2010, 466: 627-631. 10.1038/nature09253.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  48. Duez H, Staels B: Nuclear receptors linking circadian rhythms and cardiometabolic control. Arterioscler Thromb Vasc Biol. 2010, 30: 1529-34. 10.1161/ATVBAHA.110.209098.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. Staels B: When the Clock stops ticking, metabolic syndrome explodes. Nat Med. 2006, 12: 54-5. 10.1038/nm0106-54.

    Article  CAS  PubMed  Google Scholar 

  50. Walters J, Skene D, Hampton SM, Ferns GA: Biological rhythms, endothelial health and cardiovascular disease. Med Sci Monit. 2003, 9 (1): RA1-8.

    PubMed  Google Scholar 

  51. Huikuri HV, Kessler KM, Terracall E, Castellanos A, Linnaluoto MK, Myerburg RJ: Reproducibility and circadian rhythm of heart rate variability in healthy subjects. Am J Cardiol. 1990, 65: 391-393. 10.1016/0002-9149(90)90308-N.

    Article  CAS  PubMed  Google Scholar 

  52. Huikuri HV, Linnaluoto MK, Seppanen T, Airaksinen KE, Kessler KM, Takkunen JT, Myerburg RJ: Circadian rhythm of heart rate variability in survivors of cardiac arrest. Am J Cardiol. 1992, 70: 610-5. 10.1016/0002-9149(92)90200-I.

    Article  CAS  PubMed  Google Scholar 

  53. Flanagan DE, Vaile JC, Petley GW, Moore VM, Godsland IF, Cockington RA, Robinson JS, Phillips DI: The autonomic control of heart rate and insulin resistance in young adults. J Clin Endocrinol Metab. 1999, 84: 1263-1267. 10.1210/jc.84.4.1263.

    CAS  PubMed  Google Scholar 

Download references


The Air Pollution and Cardiac Risk and its Time Course (APACR) study was funded by a grant from the National Institute of Environmental Health Sciences (1 R01 ES014010), National Institutes of Health. The authors wish to thank the participants of the APACR study.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Duanping Liao.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

All authors listed on the manuscript participated in the design and coordination of the study and made substantial contribution to the intellectual content of the project to be included as authors. They also read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Rodríguez-Colón, S.M., Li, X., Shaffer, M.L. et al. Insulin resistance and circadian rhythm of cardiac autonomic modulation. Cardiovasc Diabetol 9, 85 (2010).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: