Metabolism

Shape of the OGTT glucose response curve: relationship with β-cell function and differences by sex, race, and BMI in adults with early type 2 diabetes treated with metformin

Abstract

Introduction The shape of the glucose curve during an oral glucose tolerance test (OGTT) reflects β-cell function in populations without diabetes but has not been as well studied in those with diabetes. A monophasic shape has been associated with higher risk of diabetes, while a biphasic pattern has been associated with lower risk. We sought to determine if phenotypic or metabolic characteristics were associated with glucose response curve shape in adults with type 2 diabetes treated with metformin alone.

Research design and methods This is a cross-sectional analysis of 3108 metformin-treated adults with type 2 diabetes diagnosed <10 years who underwent 2-hour 75 g OGTT at baseline as part of the Glycemia Reduction Approaches in Diabetes: A Comparative Effectiveness Study (GRADE). Insulin sensitivity (homeostasis model of insulin sensitivity, HOMA2-S) and β-cell function (early, late, and total incremental insulin and C peptide responses adjusted for HOMA2-S) were calculated. Glucose curve shape was classified as monophasic, biphasic, or continuous rise.

Results The monophasic profile was the most common (67.8% monophasic, 5.5% biphasic, 26.7% continuous rise). The monophasic subgroup was younger, more likely male and white, and had higher body mass index (BMI), while the continuous rise subgroup was more likely female and African American/black. HOMA2-S and fasting glucose did not differ among the subgroups. The biphasic subgroup had the highest early, late, and total insulin and C peptide responses (all p<0.05 vs monophasic and continuous rise). Compared with the monophasic subgroup, the continuous rise subgroup had similar early insulin (p=0.3) and C peptide (p=0.6) responses but lower late insulin (p<0.001) and total insulin (p=0.008) and C peptide (p<0.001) responses.

Conclusions Based on the large multiethnic GRADE cohort, sex, race, age, and BMI were found to be important determinants of the shape of the glucose response curve. A pattern of a continuously rising glucose at 2 hours reflected reduced β-cell function and may portend increased glycemic failure rates.

Trial registration number NCT01794143.

Significance of this study

What is already known about this subject?

  • The shape of the glucose response curve from an OGTT has been evaluated mostly in persons without diabetes, where it has been shown to be associated with β-cell function and predicts the development of both type 1 and type 2 diabetes.

What are the new findings?

  • The prevalence of a monophasic glucose response curve shape in adults is more common if male, white or obese.

  • The prevalence of a continuously rising glucose response curve at 2 hours is more common if female or African American/black.

  • The continuous rise pattern in adults with type 2 diabetes was associated with lower late insulin and total insulin and C peptide responses to glucose, indicating reduced β-cell function.

How might these results change the focus of research or clinical practice?

  • The shape of the glucose response curve may provide a biomarker to identify and target those patients with type 2 diabetes who are more likely to have either better (biphasic pattern) or poorer (continuous rise pattern) β-cell function for both research and clinical interventions.

Introduction

Type 2 diabetes is a heterogeneous disorder affecting men and women with varying degrees of obesity from all races and ethnicities. Insulin sensitivity and β-cell function are key determinants in the pathophysiology underlying the development of type 2 diabetes and vary from person to person. Glucose tolerance is typically evaluated by assessing the ability to respond to an oral glucose load. In addition to diagnostic glucose cut-points, the shape of the glucose response curve during an oral glucose tolerance test (OGTT) has been noted to fall into patterns that reflect differences in underlying β-cell function and metabolic risk.1–9 A biphasic or multiphasic glucose profile, characterized by a rise, fall, and subsequent rise in glucose after an oral glucose challenge, has been associated with better β-cell function and lower glucose concentrations compared with a monophasic pattern (simple rise and fall) in adults2–5 and youth7–9 with normal and impaired glucose tolerance. In youth with type 2 diabetes, a third pattern of a glucose that has not yet started to decline within the 2 hours following a glucose load was noted, and this pattern, which we have designated as a continuous rise pattern, was associated with poorer β-cell function and greater risk of progression compared with a monophasic or biphasic pattern.10 Only two studies have evaluated glucose response curve shape in populations that included large numbers of adults with type 2 diabetes (a large Chinese population5 and in European women with history of gestational diabetes6), but both studies only identified monophasic or multiphasic patterns.

The Glycemia Reduction Approaches in Diabetes: A Comparative Effectiveness Study (GRADE) comprises a large, multiethnic, well-characterized population of mostly middle-aged adults with <10 years duration of diagnosed type 2 diabetes. Its participants are broadly representative of the subset of early type 2 diabetes in the USA, the vast majority of whom are treated with metformin on diagnosis. Participants in GRADE underwent a standard 2-hour OGTT at baseline from which glucose response curves were evaluated. The objectives of the current analyses were threefold: first, determine the prevalence of different glucose response curve shapes in adults with type 2 diabetes treated with metformin alone; second, identify demographic and/or anthropometric factors that associate within glucose response curve subgroups by using classification analysis and multiple linear regression models; and third, determine the association between OGTT glucose response curve shape with estimates of insulin sensitivity and β-cell function. We hypothesized that glucose response curve shape would reflect β-cell function measures in type 2 diabetes and that the biphasic curve shape would be associated with better β-cell function consistent with findings in populations without diabetes.

Subjects, materials, and methods

Study design

We report cross-sectional analyses of baseline data of individuals with type 2 diabetes enrolled in the GRADE. GRADE is a multicenter National Institute of Diabetes and Digestive and Kidney Diseases-funded study conducted in 36 centers across the USA.11 The full GRADE protocol has been published11 and is also available on the GRADE website (https://grade.bsc.gwu.edu/ancillary-study-info). GRADE enrollment occurred between 2013 and 2017. All participants provided written informed consent.

Participants

Participants were eligible if they were ≥30 years of age (≥20 years if Native American) at the time of type 2 diabetes diagnosis, had been diagnosed with type 2 diabetes <10 years before study entry, and were treated with metformin alone. Participants eligible after initial screening underwent a run-in period of 6–14 weeks during which the metformin dose was adjusted, aiming for 2000 mg/day, with a minimum of 1000 mg/day. Participants who were unable to tolerate at least 1000 mg/day were ineligible. Participants were enrolled in the study if their hemoglobin A1c (HbA1c) measured at the end of run-in was between 6.8% and 8.5% (51–69 mmol/mol). A complete list of exclusion criteria is provided separately,11 but included suspected type 1 diabetes and treatment with glucose-lowering medication(s) other than metformin within the past 6 months.

A total of 5047 adults were enrolled in GRADE and completed the baseline visit. We present here a cross-sectional analysis on a subgroup of this cohort with complete baseline data available currently (n=3108). Baseline samples for those randomized to glargine (n=1263), one of the four intervention groups, have not yet been assayed for insulin given the need to validate the assay for measurement of glargine insulin and its metabolites and thus were excluded. Also excluded were 676 participants with incomplete glucose, insulin, or C peptide values.

Study procedures

After at least an 8-hour overnight fast, eligible participants underwent a baseline visit. Information about their medical history and current medication use was obtained along with assessment of body size (waist circumference, weight, and height), blood pressure, demographic characteristics (age, sex, race, and ethnicity), and laboratory measurements. A standardized 75 g OGTT was performed. Metformin was held on the morning of the OGTT. Blood samples were drawn at baseline and then at 30, 60, 90 and 120 min after the start of glucose ingestion. Plasma was separated by centrifugation, aliquoted into cryovials, and frozen prior to shipment on dry ice to the central laboratory (University of Minnesota, Minneapolis, Minnesota) where samples were assayed for glucose, insulin, and C peptide.

Assays

HbA1c was measured in EDTA whole blood on the Tosoh high performance liquid chromatography (HPLC) Glycohemoglobin Analyzer (Tosoh Medics, San Francisco, California). This method is calibrated against the National Glycohemoglobin Standardization Program standards. EDTA plasma samples were used to measure glucose with a hexokinase method on the Roche cobas c501 Chemistry Analyzer (Roche Diagnostics, Indianapolis, Indiana) and insulin and C peptide using a sandwich immunoassay (Roche Diagnostics) on the Roche cobas e601 immunoassay analyzer.

Calculations

Insulin sensitivity

The homeostasis model of insulin sensitivity (HOMA2-S) was calculated using the HOMA2 Calculator V.2.2.3 (Diabetes Trials Unit, University of Oxford, Oxford, UK).12 13

Insulin and C peptide responses to glucose

The insulinogenic index, a measure of the early insulin response, was calculated as the change in insulin divided by the change in glucose from 0 to 30 minutes. A similar C peptide index for the early response was calculated as the change in C peptide divided by the change in glucose from 0 to 30 minutes. For each OGTT, the mean incremental area under the curve (incAUC) for glucose, insulin, and C peptide was calculated as the area under the curve using the trapezoidal rule, divided by 120 minus the fasting value. The ratios of the total incAUC insulin and C peptide responses were divided by the incAUC glucose (incAUCins/incAUCglu and incAUC-CP/incAUCglu) to reflect insulin and C peptide responses to glucose over the OGTT. The late insulin response was calculated as 100×(InsulinAUC60–120/60−I0)/GlucoseAUC60–120/60−G0). An estimate of relative fasting insulin clearance was calculated as fasting C peptide divided by fasting insulin.

Glucose response curve classification

Participants were classified into subgroups according to the shape of the glucose response curve over the 120 minutes of the test. For this analysis, we required complete baseline data for glucose, insulin, and C peptide at the 0, 30, 60, 90, and 120 minute time-points (n=3108). A monophasic shape was defined as a simple increase to a peak with a subsequent decrease in glucose. A biphasic shape was defined as a rise in glucose followed by a decrease and a subsequent increase. A continuous rise shape was defined as a profile that increased monotonically with no inflection downward during the 120 min time period. A change in glucose of less than 0.25 mmol/L (4.5 mg/dL) between consecutive time-points was ignored, as done by others.3

Statistical analysis

Participant characteristics were summarized using mean and SD for quantitative variables and counts and column percentages for qualitative variables. Comparisons between pairs of glucose pattern subtypes use the χ² test of independence for qualitative variables and the Student’s t-test with unequal variances for quantitative variables. The overall comparison tests the null hypothesis of equality for all three subtypes and uses the χ² test of independence for qualitative variables and analysis of variance type III F test for quantitative variables. Adjusted p values comparing the proportions of glucose response curve shapes in men and women and between races were obtained from a type III test in a multinomial regression model with race, sex, age, body mass index (BMI), duration of diabetes, estimated glomerular filtration rate (eGFR), and metformin dose as potential determinants.

Extreme outliers were present for the following variables: insulin, HOMA2-S, insulinogenic and C peptide indices, late insulin response, and ratios of incremental insulin and C peptide AUC to incremental glucose AUC. To reduce the influence of outliers on analyses, these variables were winsorized, that is, values above or below specified cut-offs were replaced by the cut-offs. For each variable, the winsorization upper (lower) cut-off was set to the median plus (minus) 8.9 times the distance from the median to the upper (lower) quartiles. For a normally distributed variable, this results in cut-offs 6 SD above and below the mean. The number of winsorized values ranged from 5 (0.2% of values) for insulin at 30 minutes to 26 (0.8%) for the late insulin response.

Differences in metabolic responses between pairs of glucose pattern subtypes were compared using the least square means adjusting for sex, age, race, BMI, diabetes duration, eGFR, and metformin dose. Insulin and C peptide responses were further adjusted for HOMA2-S to assess β-cell function (these are denoted by an asterisk). Adjusted p values were calculated using the Tukey Honest Significant Difference (HSD) approach. The overall p value is based on a type III F test adjusting for the same covariates.

The characteristics of participants with various glucose patterns were summarized by fitting a classification tree model14 with glucose pattern as the response. To examine the effects of both sex and race on the distribution of the glucose patterns, a limited classification tree model was run including only these two covariates (figure 1A). A single classification tree model was also performed including the following covariates: race, sex, BMI as a categorical variable (0–24.9, 25–29.9, 30–34.9, 35–39.9, ≥40), waist to hip ratio, and duration of diabetes (figure 1B). Those with race classified as other, unknown, or Pacific Islander were not included in the classification analysis as the numbers were too small for Pacific Islanders and the groups classified as unknown and other were too varied, including multiracial, that it was hard to draw conclusions (n=2859). An additional four participants were lacking BMI values and thus were excluded from the model that included BMI (n=2855). Classification trees for figure 1 are presented in online supplemental figure 1. In each case, minimal cost-complexity pruning with 20-fold cross-validation was used to obtain the final classification.

Figure 1
Figure 1

Classification analysis of OGTT glucose response curve shapes by sex and race (A) and sex, race, BMI, waist to hip ratio, and duration of diabetes (B). In B, only sex and BMI remained as significant independent classifiers. The percentages within each column are shown on the plot. AA, African American; AI, American Indian; BMI, body mass index; OGTT, oral glucose tolerance test.

All analyses were performed using R V.3.6.015 and p<0.05 was considered statistically significant.

Results

Participant characteristics

Overall, individuals in the cohort were mostly male, middle-aged, obese, and had diabetes for an average of 4 years (table 1). The cohort was racially diverse, with one-third being non-white. Over 90% of participants were taking the maximum dose of metformin (2000 mg/day). Most participants had a monophasic glucose response curve shape (67.8%), while 26.7% demonstrated a continuous rise pattern and only 5.5% a biphasic shape.

Table 1
|
Baseline characteristics by glucose response curve shape

There were several notable differences in baseline characteristics between the three subgroups. Those with a monophasic glucose response curve shape were more likely to be male, younger, and white, have a shorter duration of diabetes, higher BMI and waist circumference, and taking a higher dose of metformin compared with the other two subgroups. Participants with a continuous rise shape were more likely to be female and African American/black. Differences in prevalence by sex remained significant after adjusting for race, age, BMI, duration of diabetes, eGFR, and metformin dose (p<0.001). Differences in prevalence by race also remained significant after adjusting for sex, age, BMI, duration of diabetes, eGFR, and metformin dose (p=0.019). There were no statistically significant differences between the subgroups by duration of diabetes or renal function.

Prevalence of the glucose response curve shape by classification analysis

Classification analysis revealed subsets by sex, race, and BMI (figure 1). When sex and race were combined in the model, six separate subsets emerged, confirming differential distributions by both race and sex (figure 1A). Asian and American Indian men and women had the highest prevalence of the biphasic shape, while African American/black men had the lowest.

To evaluate if obesity changed the classification analysis, BMI (0–24.9, 25–29.9, 30–34.9, 35–39.9, ≥40) was added to a model including race, sex, waist to hip ratio, and duration of diabetes (figure 1B). Sex and BMI ≥30 identified three distinct subsets, while race, waist to hip ratio, and duration of diabetes dropped out. BMI <30 had a greater proportion of biphasic curve shape, while obese men had the largest fraction of the monophasic curve shape.

Metabolic differences by glucose response curve shape

Figure 2 illustrates the mean and 95% CI for glucose, insulin, and C peptide concentrations over time derived from the OGTT for each subgroup. The metabolic data are presented in table 2. The biphasic subgroup had the lowest peak glucose and mean incAUCglu. The 2-hour glucose and HbA1c were highest in the continuous rise subgroup, but peak glucose and mean incAUCglu were the highest in the monophasic subgroup. Fasting glucose did not differ.

Figure 2
Figure 2

OGTT glucose response curve shape for the whole cohort classified as monophasic, biphasic, or continuous rise. Mean glucose (A), insulin (B), and C peptide (C) are depicted at each OGTT time-point separately for participants within each of the subgroups. Pointwise 95% CIs are drawn around the means at each time-point. OGTT, oral glucose tolerance test.

Table 2
|
Comparison of metabolic parameters by glucose response curve shape

Participants with the biphasic response pattern had the highest early, late, and total insulin and C peptide responses to glucose load adjusted for insulin sensitivity. Participants with the monophasic and continuous rise response patterns had similar early insulin and C peptide responses, but those with continuous rise had lower late insulin and total insulin and C peptide responses (table 2 and illustrated in figure 2B,C). Those with the continuous rise pattern also appeared to have a relatively delayed glucose rise that was matched by a concomitant delay in the rise in insulin and C peptide concentrations (figure 2). This delay in glucose rise in the continuous rise subgroup was observed in both men and women (figure 3).

Figure 3
Figure 3

OGTT glucose response curve shapes stratified by sex: monophasic, biphasic, or continuous rise. Mean glucose (A women, D men), insulin (B women, E men), and C peptide (C women, F men) are depicted at each OGTT time-point separately for participants within each of the subgroups. Pointwise 95% CIs are drawn around the means at each time-point. OGTT, oral glucose tolerance test.

HOMA2-S did not differ by glucose curve shape (table 2). Note that all insulin and C peptide responses to glucose were adjusted for insulin sensitivity (HOMA2-S) given the well-known inverse relationship between responses and sensitivity, thereby providing an estimate of β-cell function (table 2, adjustment for HOMA2-S denoted by asterisk).

Conclusions

This study presents novel data on several aspects related to the shape of the glucose response curve during an OGTT. First, it is the first large study of its kind in people with type 2 diabetes, all on metformin only, that includes both men and women over a broad range of ethnicities and BMI. Second, the analysis highlights novel associations in the prevalence of glucose response curve patterns by sex, race, age, and BMI. Third, it investigates a new continuous rise glucose pattern, which has only previously been evaluated in youth in the Treatment Options for Type 2 Diabetes in Adolescents and Youth (TODAY) study,10 and more recently in youth and adults with impaired glucose tolerance or newly diagnosed type 2 diabetes in the Restoring Insulin Secretion (RISE) study.16 Fourth, the study highlights the differences in β-cell function linked with each glucose response curve pattern, with a biphasic shape, although uncommon, demonstrating better β-cell function and the continuous rise subgroup showing reduced measures of late β-cell function relative to those with a monophasic shape.

Differences by sex in the prevalence of the monophasic and continuous rise subgroups were present that were independent of race and BMI using both a least squares and classification approach. Previous studies found that the biphasic or multiphasic pattern was more common in women (3-hour OGTT analysis in a large Chinese adult population with and without diabetes5 and a European population without diabetes3) or did not differ by sex.2 4 However, none of these studies included an analysis of a continuously rising glucose pattern. In GRADE, we found that more women exhibited a continuous rise shape, which has lower β-cell function, findings that are opposite to that in populations without diabetes, where the biphasic or multiphasic shape, associated with better β-cell function, was more common in women. The reason for a sex difference in the prevalence of the monophasic and continuous rise glucose curve shapes among GRADE participants is not clear but warrants further investigation.

Race also impacted the prevalence of the glucose curve subgroups but was in part due to differences in BMI. When BMI category was added to the classification analysis, race dropped out, but race remained independently associated with glucose response curve shape when adjusting for BMI in the least squares analysis. The most significant finding was a higher prevalence of the continuous rise shape in African American/black. Racial differences in curve shape distribution have not been reported for adults with type 2 diabetes and there were no differences in glucose response curve shape by race in youth with type 2 diabetes in the TODAY study.10 If the continuous rise pattern in adults follows the more aggressive course observed in youth in the TODAY study,10 the greater prevalence of this pattern in African American/black could portend a higher risk of progression to glycemic failure. African American/black individuals are disproportionally affected by type 2 diabetes, with a prevalence in the USA of 16.4% compared with 11.9% in non-Hispanic white.17 They also are at higher risk of developing type 2 diabetes, with the Coronary Artery Risk Development in Young Adults study finding a more than twofold higher risk of developing type 2 diabetes, although the disparity was no longer significant after adjusting for known modifiable risk factors, with biological factors having the greatest impact.18

Differences in both age and diabetes duration were also noted in the prevalence of the monophasic and continuous rise subgroups, with the monophasic subgroup being older and having a longer duration of diabetes. This is perhaps not surprising. A systematic review of the literature indicates that β-cell function declines with aging, consistent with the increase in incidence of type 2 diabetes in the aging population.19

The continuous rise curve shape has not yet been described in a large cohort of adults with type 2 diabetes. The two large studies that evaluated glucose response curve shape in adults with type 2 diabetes5 6 based classification on a 3-hour OGTT rather than the standard 2-hour test. Delays or deficiencies in insulin secretion and response that result in a continuous rise shape on a 2-hour test need a longer sampling time frame to capture the eventual fall in glucose. In essence, the classification of monophasic versus continuous rise glucose curve shapes reflects differences in timing. Thus, classification by curve shape alone during a 3-hour test is likely to miss important differences in timing between monophasic and continuous rise subgroups identified on a 2-hour test.

In GRADE, the early insulin and C peptide responses to glucose did not differ between the continuous rise and monophasic subgroups, but the continuous rise subgroup was characterized by lower late insulin and total insulin and C peptide responses, suggesting greater loss of insulin secretory ability. Peak glucose also differed between these two subgroups, being higher in the monophasic subgroup. Glucose acts not only to stimulate insulin secretion directly but can prime the β-cell to enhance subsequent insulin secretion.20 21 The higher responses in the monophasic subgroup may reflect an effect of higher peak glucose concentrations to potentiate late-phase insulin secretion or β-cells that are more responsive to the effects of potentiation. Consistent with our results, the TODAY study in youth with type 2 diabetes10 and the RISE study comparing youth and adults with impaired glucose tolerance (n=249 adults) or newly diagnosed type 2 diabetes (n=104 adults)16 both found that the continuous rise pattern was associated with the lowest β-cell function.

The continuous rise subgroup appears to have unique characteristics. The profiles for insulin and C peptide depicted in figure 1 suggest a ‘sluggish’ glucose response in those with a continuously rising glucose response curve. The rise in the glucose concentration, which acts as the stimulus for insulin secretion, also appears to be somewhat delayed, and the insulin and C peptide curves mimic the glucose response curve. While not tested in GRADE, one possible explanation for this delay could be that gastric emptying or glucose absorption is delayed. This hypothesis would need to be tested with studies designed to directly measure rates of glucose absorption.

The continuous rise glucose curve shape may identify a subgroup of individuals at higher risk of progression or who may respond differentially to specific treatment approaches. In youth with type 2 diabetes in the TODAY study, a pattern of a continuously rising glucose was associated with a higher HbA1c and lower late insulin and total insulin and C peptide responses, and presaged increased glycemic failure rates and accelerated deterioration of β-cell function.10 In RISE, progression from impaired glucose tolerance to diabetes or change in glycemic outcomes at 12 months did not differ by baseline glucose response curve categorization but may have been limited by the short follow-up time.16 Interestingly, progression from a biphasic to a monophasic and then ultimately to a continuous rise pattern along with corresponding decreasing β-cell function was also recently reported to occur during progression to type 1 diabetes in the TrialNet Pathway to Prevention study.22 Once GRADE is completed, we will be able to test a similar hypothesis in adults to determine if the continuous rise pattern identifies those at highest risk of progression and whether the biphasic pattern identifies those with more stable disease.

The biphasic glucose response pattern was uncommon, being present in only 5.5% of participants. These data are consistent with findings from a large population study in China in which only 3% of those with a diagnosis of type 2 diabetes had a multiphasic pattern based on a 3-hour OGTT.5 In populations without diabetes, biphasic or multiphasic curve shapes correlate with lower glucose values, higher insulin sensitivity and β-cell function,2–4 7 9 and lower risk of progression to diabetes23 and are present in a higher percentage of the population (20%–30% in adults). In GRADE, HbA1c and fasting glucose values were similar with the monophasic and continuous rise subgroups, but estimates of β-cell function were higher. The presence of the biphasic pattern within a small proportion of the diabetes population suggests that this may identify a subgroup that may differ in terms of pathophysiology and/or indicate a lower risk of progression.

The strengths of this study include the large multiethnic cohort of men and women with relatively early type 2 diabetes (<10 years) and use of OGTT to assess β-cell function. As all participants were treated with metformin alone for their diabetes, there was no confounding by use of other glucose-lowering medications. The study is generalizable to a large portion of the US population that has type 2 diabetes treated with metformin alone. Limitations include the cross-sectional analysis and the lack of more precise measures of insulin sensitivity and β-cell function.

In conclusion, this is the first large multiethnic cohort of adults with early type 2 diabetes to identify differences in OGTT glucose curve shape distribution by sex, race, and BMI and metabolic differences between the two most common curve shapes, monophasic and continuous rise. The continuous rise curve shape exhibited more advanced β-cell dysfunction and higher HbA1c compared with the monophasic subgroup and may serve as a biomarker for more advanced disease. Further mechanistic studies are needed to evaluate the metabolic processes underlying these differences. In the future, once GRADE is completed, analysis of GRADE outcome data will allow us to determine if these patterns predict differential progression and/or response to the four glucose-lowering therapies being studied.