Abstract
Objective:
To report the impact of continuous glucose monitoring (CGM) on glycemic variability (GV) indices, factors predictive of change, and to correlate variability with conventional markers of glycemia.
Methods:
Data from the JDRF study of CGM in participants with type 1 diabetes were used. Participants were randomized to CGM or self-monitored blood glucose (SMBG). GV indices at baseline, at 26 weeks in both groups, and at 52 weeks in the control group were analyzed. The associations of demographic and clinical factors with change in GV indices from baseline to 26 weeks were evaluated.
Results:
Baseline data were available for 448 subjects. GV indices were all outside normative ranges (P < 0.001). Intercorrelation between GV indices was common and, apart from coefficient of variation (CV), low blood glucose index (LBGI), and percentage of glycemic risk assessment diabetes equation score attributable to hypoglycemia (%GRADEhypoglycemia), all indices correlate positively with HbA1c. There was strong correlation between time spent in hypoglycemia, and CV, LBGI, and %GRADEhypoglycemia, but not with HbA1c. A significant reduction in all GV indices, except lability index and mean absolute glucose change per unit time (MAG), was demonstrated in the intervention group at 26 weeks compared with the control group. Baseline factors predicting a change in GV with CGM include baseline HbA1c, baseline GV, frequency of daily SMBG, and insulin pump use.
Conclusions:
CGM reduces most GV indices compared with SMBG in people with type 1 diabetes. The strong correlation between time spent in hypoglycemia and CV, LBGI, and %GRADEhypoglycemia highlights the value of these metrics in assessing hypoglycemia as an adjunct to HbA1c in the overall assessment of glycemia.
Introduction
C
Despite evidence supporting a role for glycemic variability (GV), 6 –9 its contribution to diabetes-related complications remains controversial and its clinical value as an additional measure to HbA1c remains unclear. 10,11 A recent systematic review has suggested that GV is associated with an increased risk of microvascular complications in type 2 diabetes only and that the relationship between GV and vascular complications in type 1 diabetes is less clear. 12
This analysis evaluates several indices of GV in a large cohort with type 1 diabetes with reference to ranges previously defined in people without diabetes and explores the primary hypothesis that the use of real-time continuous glucose monitoring (CGM) reduces GV. We have also explored the extent to which different GV measures positively correlate with HbA1c and time spent in hypoglycemia, and whether baseline factors predict GV. Our aim is to identify GV measures that could be a useful adjunct to HbA1c for a better assessment of the overall glycemic status of people with type 1 diabetes.
Methods
Data source
Data from the JDRF CGM study were used.
13
–15
The publicly accessible data were obtained from the Jaeb Center for Health Research (
The JDRF study protocol and clinical characteristics of enrolled subjects have been previously described in detail. 13 –15 In summary, the study is a 6-month randomized, parallel group, efficacy, and safety study designed to evaluate the impact of CGM on glycemic control in children and adults with type 1 diabetes. Participants underwent blinded CGM for 1 week before randomization to either standard self-monitored blood glucose (SMBG) (control) or use of unblinded CGM as a supplement to SMBG (intervention). All trial participants were provided with written instructions on how to use the data provided by CGM and capillary blood glucose meters to make real-time adjustments to insulin, and on the use of computer software to retrospectively review the glucose data to alter future insulin doses. Patients using CGM received additional instructions for modifying insulin doses and treatment of hypoglycemia on the basis of the glucose trend. Changes were made to diabetes management, as needed, for all participants during scheduled contacts. The randomized trial was followed by a 6-month extension study in which unblinded CGM was continued in the CGM group, but with unblinding to CGM data in the control group. CGM profiles were obtained at 26 weeks postrandomization from both control and intervention groups.
Measures of GV
We computed 13 measures of GV using EasyGV (v8.8.2.R2) software. Evaluated GV measures included standard deviation (SD), coefficient of variation (CV), mean amplitude of glycemic excursions (MAGE), continuous overall net glycemic action (CONGA), mean of daily differences (MODD), lability index (LI), and mean absolute glucose change per unit time (MAG), glycemic risk assessment diabetes equation (GRADE), M-value, average daily risk range (ADRR), J-Index, low blood glucose index (LBGI), and high blood glucose index (HBGI). GRADE score is also reported as %GRADEhypoglycemia, %GRADEeuglycemia, and %GRADEhyperglycemia representing percentages of GRADE scores attributable to glucose values <3.9 mmol/L, and between 3.9–7.8 mmol/L and >7.8 mmol/L, respectively. EasyGV is a Microsoft Excel workbook that has a number of options to define the sampling interval, CONGA length, LI interval, reference value of M-value, and whether SMBG or CGM is used for MAGE calculation. 16 A description of the GV measures, formulae used for their calculations, and a critical review of their limitations is described elsewhere. 6,17
Analysis design
GV indices were assessed in comparison to previously published GV reference ranges in people without diabetes. 6 To evaluate the effect of CGM on GV, between-group differences in GV indices were evaluated at 26 weeks and glucose profile data collected at 52 weeks from subjects in the control group who crossed over to CGM (unblinded control) were compared with data collected at 26 weeks in the same group (blinded control). The relationship between baseline GV indices and HbA1c was explored. Finally, baseline predictors of change in GV indices from baseline to 26 weeks in the CGM group were evaluated.
Statistical analysis
Data were examined for normality; non-normally distributed variables were logarithmically transformed with the use of geometric mean, SD, and 95% confidence interval (CI) for descriptive statistics of baseline data (n = 448).
An unpaired t-test was used to compare the effect of CGM on the change of GV indices from baseline to 26 weeks in the control and intervention arms, while a paired t-test was used for analysis of within-group differences in each of those arms. To study the effect on GV of unmasking CGM in the control group, a paired t-test was used to compare GV indices in the 207 participants in the control group with unblinded CGM at 52 weeks with their blinded GV indices at 26 weeks. Spearman correlation was used to examine the relationship between HbA1c and GV indices at baseline.
The associations of baseline demographic and clinical factors with change in GV indices from baseline to 26 weeks were evaluated in the CGM group (n = 231) using regression analysis. The analysis was constructed using the following predictor variables: age, gender, race, education level of caregiver, insulin modality (pump or multiple daily injections), frequency of daily self-reported blood glucose monitoring, occurrence of one or more episodes of severe hypoglycemia in last 6 months, diabetes duration, baseline HbA1c, and baseline GV. The analysis was performed with change in GV measures and continuous variables expressed as z scores. Categorical variables were included as dummy variables. Baseline factors with P ≤ 0.2 in the univariate analysis were carried forward to multivariable analysis.
Data are presented as means (SD), unless otherwise stated. Statistical tests were two tailed and for descriptive and exploratory analyses, a significance level of P < 0.05 was adopted. Where significance levels were consistently <0.001, t values were reported to aid evaluation of relative differences in magnitude and scatter of the differences observed. For exploring our primary hypothesis that the use of real-time CGM reduces GV, 13 different measures of glycemic variability were tested. Therefore, to test our primary hypothesis, we have adopted a significance level of P < 0.004 (i.e., 0.05/13). Statistical tests were performed using SPSS 21.0 for Mac (SPSS Inc., Chicago, IL).
Results
Participants
CGM profile data were available for 448 subjects at baseline (54.9% women, 94.4% white race), following exclusion of three subjects due to missing HbA1c data at baseline. Mean age was 25.1 years (SD 15.8) and mean diabetes duration 13.6 years (SD 11.7). Two hundred thirty-one subjects were randomized to the CGM group and 217 subjects to the control group. At 26 weeks, CGM profile data were available for all the 231 subjects in the CGM group and 214 subjects in the control group (blinded CGM). At 52 weeks, data were available for 207 subjects in the control group (unblinded CGM). Baseline characteristics were similar between the two groups (Table 1). There was no statistically significant difference in HbA1c and the evaluated GV measures between the two groups at baseline.
Data are presented as number (%), mean (SD), or median (IQR).
CGM, continuous glucose monitoring; IQR, interquartile range; MDI, multiple daily injections; SD, standard deviation; SMBG, self-monitored blood glucose.
Measures of GV and glycemic control at baseline
Reference means and 95% CIs for several measures of GV have been previously described by analyzing CGM profiles of 70 subjects without diabetes. 6 As shown in Table 2, the mean values of GV indices at baseline in the 448 subjects with type 1 diabetes were all appreciably higher than the reference means, with complete nonoverlap of 95% CIs throughout (P < 0.001).
Difference between measures of GV and glycemic control in the two groups is statistically significant (P < 0.001). Measures are expressed in mmol/L.
CI, confidence interval; GV, glycemic variability.
Intercorrelation between measures of GV and glycemic control at baseline and correlation with HbA1c
Statistically significant intercorrelation between measures of GV and HbA1c was common. The GV index, ADRR, showed the strongest intercorrelation with other measures with r > 0.7 in 9 out of 13 intercorrelations. By contrast, CV, LBGI, and %GRADEhypoglycemia showed the weakest intercorrelation. Similarly, all these measures, apart from CV, LBGI, and %GRADEhypoglycemia, correlated significantly, but moderately (r 0.35–0.66), with HbA1c as shown in Table 3. The relationship between percent of time spent in hypoglycemia (glucose level <3.9, <3.3, and <2.8 mmol/L), CV, LBGI, %GRADEhypoglycemia, and HbA1c was analyzed using Spearman correlation. CV, LBGI, and %GRADEhypoglycemia correlated strongly with time spent in hypoglycemia (<3.9, 3.3, and 2.8 mmol/L) (r > 0.88, P < 0.001). No relationship between HbA1c and time spent in hypoglycemia was seen (Table 4).
Correlation is significant at the 0.05 level (two tailed).
Correlation is significant at the 0.01 level (two tailed).
Correlation is significant at the 0.01 level (two tailed).
Correlation is significant at the 0.05 level (two tailed).
Effect of CGM on GV measures in the intervention group
Analysis of between-group differences demonstrated a significant difference between the intervention and control group in all measures of GV, with the exception of CV, LI, and MAG (Table 5). At 26 weeks, there was a significant reduction in all measures of GV, with the exception of LI and MAG, from baseline in the intervention group. There was a significant reduction in HbA1c in the intervention group from 58 mmol/mol (7.4%) to 55 mmol/mol (7.2%) (P < 0.001). Correspondingly, the achieved relative reductions in M-value, LBGI, and GRADE were 25.7%, 24.9%, and 16.5%, respectively (P < 0.001) (Table 5). In contrast, there was no statistically significant reduction in any of the measures of GV and glycemic control in the control group at 26 weeks compared to baseline.
The table also shows the absolute and relative change in the mean values in the intervention group. Measures are expressed in mmol/L.
Between-group difference for change in measures of GV and quality of glycemic control at end of 26 weeks from baseline.
Effect of unmasking of CGM in the control group
The effect of unmasking CGM in the control group at 26 weeks was evaluated by comparing GV measures at 52 weeks to those at 26 weeks (immediately before unmasking CGM). Despite the nonsignificant change in HbA1c from 26 to 52 weeks, there was a significant reduction in all measures of GV, with the exception of LI and MAG, from baseline (Supplementary Table S1; Supplementary Data are available at
Factors predictive of response in the CGM group
Univariate analysis showed that baseline GV was a significant predictor of change in all measures of GV at 26 weeks in the intervention group, with higher baseline GV associated with greater reduction in GV at 26 weeks. Similarly, baseline HbA1c was also a significant predictor of change in GV in the majority of the evaluated GV measures (all except LBGI, CV, LI, and MAG), with higher baseline HbA1c associated with greater reduction in GV measures at 26 weeks. However, multivariable analysis showed that higher HbA1c at baseline was associated with less of a reduction of GV measures at 26 weeks (Supplementary Table S2). In the multivariable analysis, treatment with insulin pump predicted a reduction in LBGI and M-value, while frequent use of SMBG predicted a reduction in LI. Other variables, including education level of caregiver, did not predict change in any of the evaluated GV measures.
Discussion
The increased availability of CGM provides a wealth of data and enables assessment of GV. Several measures of GV have been described. 18 –23 These can be broadly subdivided into measures based on glucose distribution (e.g., SD, CV, MAGE, CONGA, MODD, LI, and MAG), and measures based on risk and quality of glycemic control that are also sensitive to GV (e.g., GRADE, M-value, ADRR, J-Index, LBGI, and HBGI). 24 There are several challenges related to GV measurements and interpretation, 24 including correlation with mean glycemia, making it difficult to evaluate if the change in GV measures following an intervention is related to change in mean glycemia, GV, or both. These challenges are further complicated by the plethora of GV measures proposed in the literature with the lack of a “gold standard” measure. 18 –23
The data reported here are from the largest available CGM data set, with evaluation of measures of GV, the effect of real-time CGM on these measures, and predictors of changes to variability. The analysis demonstrates the magnitude of variability in people with type 1 diabetes and the impact that real-time CGM has on GV measures. Measures of dispersion (SD and CV), glucose risk indices (LGBI, HGBI, ADRR, GRADE), measures of unknown significance (CONGA, MAGE, J-index, M-value), and day to day change (MODD) all fell significantly in the intervention group over 26 weeks, and (with the exception of CV) fell significantly more than in the control group. Moreover, over a further 26 weeks, in the control group, unblinding resulted in significant falls in these measures of GV. However, two measures of variability over time (LI and MAG) were entirely unchanged. This may suggest that the use of real-time glucose trends and alarms addresses glucose dispersion and risk while the underlying temporal glucose variability remains unchanged. It could also reflect the impact of interventions based on CGM use such as more frequent insulin correction boluses or carbohydrate ingestion, resulting in more high-frequency (rapid) glucose fluctuations that LI and MAG are more sensitive to.
The impact of CGM on GV measures is consistent with previous studies. 25 –27 However, these studies have been for a shorter duration 25,26 or in limited numbers of participants. 27 Similarly, the correlation between GV measures reported in this analysis confirms previous findings. A correlation analysis of GV measures in 48 subjects with type 1 diabetes showed similar relationships to our study between MAGE and other measures of GV, with a significant correlation of 0.74 between MAGE and SD and a correlation of 0.58 between SD and CONGA1. 28 Rodbard also reported a strong correlation between SD and MAGE (r = 0.89), SD and CONGA1 (r = 0.71), SD and MODD (r = 0.81), and between MAGE and MODD (r = 0.74). 29 Compared to our analysis, other studies also showed similar relationships between HbA1c and HBGI (r = 0.63) 30 and between HbA1c and MAGE (r = 0.49). 31 In another retrospective analysis of 72-h CGM data in 815 outpatients (48 with type 1 diabetes and 767 with type 2 diabetes), correlation between MAG and SD was stronger (r = 0.88 compared to 0.66 in our analysis), as was the relationship between MAG and other GV measures. 32
Common intercorrelation between evaluated GV measures is demonstrated alongside a moderate correlation between most of the evaluated GV measures and HbA1c (Table 3). This suggests that these measures convey similar information, reflecting mean glycemia, as well as information on glucose variability.
LBGI and %GRADEhypoglycemia, which are sensitive to hypoglycemia alone, correlated poorly with other GV measures and not at all with HbA1c, supporting the hypothesis that HbA1c has a limited role in reflecting or predicting the risk of hypoglycemia 20 and suggesting that these two metrics may offer additional information when assessing CGM data. The relationship between time spent in hypoglycemia, LBGI, and %GRADEhypoglycemia may suggest that direct assessment of time spent in hypoglycemia, routinely reported from CGM data, is sufficient. However, these two metrics provide a continuous scale for hypoglycemia with a different risk score assigned to glucose values rather than the categorical information provided by time spent below a glucose threshold, which may not consider severity, with all values assigned an equal weight. 33 Similarly, CV correlated with time spent in hypoglycemia (Table 4) and not with HbA1c or mean glucose (Table 3). CV has been previously proposed as the best parameter to characterize GV since it is corrected for the mean and avoids the dependency of SD and other measures of GV on mean glucose or HbA1c. 34 The weaker impact of CGM on CV demonstrated in the intervention group supports the role of CV in providing distinct information from the other GV measures, and from HbA1c.
There was a complex of associations between HbA1c and GV at baseline and subsequent change in GV in the CGM group. In univariate analyses, higher baseline HbA1c was associated with greater GV (except for CV and LGBI), and higher baseline HbA1c and GV were associated with greater reductions in GV over 26 weeks. However, multivariable analysis showed that the association between higher baseline HbA1c and greater reduction in GV was secondary to the association between baseline GV and reduction in GV. In fact, when variation in baseline GV was taken into account, the association between baseline HbA1c and subsequent reduction in GV changed direction: higher baseline HbA1c was associated with an increase in the majority of GV measures over 26 weeks, indicating that among those with high baseline HbA1c, the subsequent reduction in GV was not as great as would have been expected from their baseline GV.
Limitations of this analysis include the reliability of some of the evaluated GV measures. ADRR requires data collected typically over 1 month with a minimum of 14 days and with typical frequency of 3–5 glucose measurements per day, whereas CGM data used to compute ADRR in control group at 26 weeks in this analysis were collected for only 1 week.
Another limitation is the possible different outcomes when calculating GV measures using different automated GV calculators. Several methods and automated calculators have been described for calculating MAGE. 35 –38 EasyGV uses a modified method (MAGE-CGM) for MAGE calculation. The MAGE-CGM formula selects a peak or trough based on direction of change (rising or falling) of the preceding and succeeding data points. It also contains a 15-min lag window for the direction of change based on the lag between interstitial fluid glucose measurement and plasma glucose concentrations. It also contains an algorithm that eliminates short-term fluctuations related to sensor inaccuracies. 6 The correlation between MAGE calculated using EasyGV and MAGE calculated based on the methods of Fritzsche or Baghurst was reported to be of 0.87. 39 Comparison between EasyGV (v8.8.2.R2) and a validated automated GV calculator showed a correlation of 0.76 in calculating MODD but poor correlation in calculating CONGA1. 40 EasyGV has now been updated to address an issue with the MODD calculation but the reason for discrepancies in CONGA1 between calculators remains unclear and is the subject of ongoing investigation.
It is also important to note that the normative range of GV measures referred to in our article was based on the analysis of 55,000 data points from only 72 h of CGM obtained from 70 subjects with fasting plasma glucose of <120 mg/dL. Although this excludes subjects with diabetes, it is possible that some subjects had impaired glucose tolerance. This could explain why SD in this cohort, reported to be 1.5, is larger compared to SD reported in other studies in subjects with normal glucose tolerance (ranged 0.5–0.9). However, other GV measures were similar to those reported from analysis of other data sets. 28,31,41 –43 The accuracy of the CGM systems used in the study, particularly in the hypoglycemic range, is a further limitation and may affect GV measurements. Since the CGM study was conducted, the accuracy of available systems has improved significantly and the large volume of data analyzed here mitigates much of this problem.
Conclusion
In summary, we report the largest study to date of the long-term impact of CGM on GV in people with type 1 diabetes. There was a significant reduction in the evaluated measures of GV with CGM use, with the exception of LI and MAG, suggesting that CGM reduces mean glycemia, glucose dispersion, and risk but not glucose fluctuation. The strong correlation between time spent in hypoglycemia, CV, LBGI, and %GRADEhypoglycemia, but not with HbA1c, highlights the value of these metrics in assessing hypoglycemia and as a useful adjunct to HbA1c in overall assessment of glycemic status in people with type 1 diabetes. A large-scale longitudinal intervention study is required to evaluate the HbA1c-independent role of GV in the development of diabetes-related vascular complications and this will additionally require the definition of an agreed “gold-standard” GV measure.
Footnotes
Authors' Contributions
A.H.L. analyzed the data and wrote the manuscript; I.F.G. supervised data analysis and reviewed/edited the manuscript; D.G.J. contributed to the discussion and reviewed/edited the manuscript; N.S.O. processed original glucose profile data and reviewed/edited the manuscript.
Author Disclosure Statement
N.S.O. has received honoraria for consultancy from Abbott Diabetes Care and Roche.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
