Highly informative yet simple protocols to assess insulin secretion and action would considerably enhance the quality of epidemiological and large-scale clinical trials. In an effort to develop such protocols, a 5-h, 11-sample oral glucose tolerance test (OGTT) was performed in 100 individuals and a 7-h, 21-sample meal in another 100. Plasma glucose, insulin, and C-peptide concentrations were measured. We show that virtually the same minimal model assessment of β-cell responsivity (dynamic [Φd] and static [Φs]), insulin sensitivity (Si), and disposition index (DI) can be obtained with a reduced seven-sample 2-h protocol: Φd, reduced versus full: 871.50 vs. 873.32, r = 0.98 in OGTT and 494.88 vs. 477.99 10−9, r = 0.91 in meal; Φs: 42.36 vs. 44.35, r = 0.88 in OGTT and 35.31 vs. 35.37 10−9 min−1, r = 0.90 in meal; Si: 24.33 vs. 22.77 10−5 dl · kg−1 · min−1 per pmol/l, r = 0.89 in OGTT and 19.03 vs. 19.77 10−5 dl · kg−1 · min−1 per pmol/l, r = 0.85 in meal; and DI: 1,282.26 vs. 1,273.23, r = 0.84 in OGTT and 726.92 vs. 776.97 10−14 dl · kg−1 · min−2 per pmol/l, r = 0.84 in meal. This reduced protocol will facilitate the study of insulin secretion and action under physiological conditions in nondiabetic humans.
The oral meal or oral glucose tolerance test (OGTT) glucose and C-peptide minimal models can simultaneously measure insulin action, β-cell function, and the rate of meal glucose appearance (1–5). These models have been subjected to various validation strategies, including multiple tracer protocols (6) and euglycemic and hyperglycemic clamp studies (7, 8) and have been used in pathophysiological studies (9). All of the above studies have relied on a 300-min OGTT and 420-min meal protocols in which plasma glucose, insulin, and C-peptide concentrations have been sampled 11 and 21 times, respectively (Fig. 1A). Since this large number of samples and prolonged period of study is not feasible for large population-based studies, we sought to determine whether a seven-sample protocol performed over 120 min could provide estimates of insulin secretion and action equivalent to those obtained with the longer and more complex studies. To do this, we analyzed data from 200 studies consisting of 100 5-h, 11-sample OGTTs and 100 7-h, 21-sample meal studies performed in subjects with various degrees of glucose tolerance. We report that a seven-sample protocol performed over 120 min following ingestion of either a mixed meal or 75 g glucose enables accurate assessment of both insulin secretion and action in nondiabetic humans.
RESEARCH DESIGN AND METHODS
One hundred nondiabetic subjects underwent an OGTT performed at the Division of Endocrinology, Metabolism, and Lipid Research, Washington University School of Medicine, St. Louis, Missouri, and 100 nondiabetic subjects underwent a mixed meal test performed at the Division of Endocrinology, Diabetes, Metabolism, and Nutrition, Department of Internal Medicine, Mayo Clinic and Foundation, Rochester, Minnesota. All subjects provided informed consent.
The OGTT consisted of oral administration of 75 g glucose at time 0; as detailed in (2), blood samples were collected at time 0, 10, 20, 30, 60, 90, 120, 150, 180, 240, and 300 min, and plasma glucose, insulin, and C-peptide concentrations were measured. Subject characteristics are reported in Table 1 (right column).
The mixed meal (10 kcal/kg, 45% carbohydrate, 15% protein, and 40% fat) contained 1 ± 0.02 g/kg glucose. As detailed in (9), plasma samples were collected at 0, 5, 10, 15, 20, 30, 40, 50, 60, 75, 90, 120, 150, 180, 210, 240, 260, 280, 300, 360, and 420 min, and plasma glucose, insulin, and C-peptide concentrations were measured. Subject characteristics are reported in Table 1 (left column).
Oral glucose minimal model.
Insulin sensitivity (Si; dl · kg−1 · min−1 per pmol/l) was estimated from plasma glucose and insulin concentrations measured during the test by using the oral glucose minimal model (4, 6). The model is shown in Fig. 2A. Si measures the overall effect of insulin to stimulate glucose disposal and inhibit glucose production. The model also reconstructs the rate of appearance (Ra; mg · kg−1 · min−1) in plasma of ingested glucose.
Oral C-peptide minimal model.
β-Cell responsivity indexes were estimated from plasma glucose and C-peptide concentrations measured during the test by using the oral C-peptide minimal model (2, 3, 5). The model is shown in Fig. 2B. Insulin secretion is made up of two components. The dynamic component is likely to represent secretion of promptly releasable insulin and is proportional to the rate of increase of glucose concentration through a parameter, Φd (10−9), that defines the dynamic responsivity index. The static component derives from provision of new insulin to the releasable pool and is characterized by a static index, Φs (10−9 min−1), and by a delay time constant (T; min). The meaning of Φs and T can be made clear with reference to a response to an above-basal step increase of glucose: provision tends with time constant T toward a steady state that is linearly related to the glucose step size through parameter Φs.
Φd and Φs can also be expressed in relation to insulin sensitivity through the dynamic (DId) and static (DIs) disposition indexes, i.e., DId = Φd × Si and DIs = Φs × Si. Since from Φd and Φs one can also calculate a single overall β-cell responsivity, Φ (10−9 min−1), a single overall DI can also be derived as DI = Φ × Si. The model also reconstructs insulin secretion (SR; pmol/min) and its dynamic (SRd) and static (SRs) components.
Model identification.
The oral glucose and C-peptide minimal models of Fig. 2 were numerically identified in both full and reduced protocol OGTT and meal studies, as detailed in (2, 4, 6). Two comments are in order. The first concerns the fraction of ingested glucose that is absorbed (area under the curve [AUC] of Ra from time 0 to 7 h divided by the dose). This value is fixed a priori to 0.90 in the full protocol identification of the glucose minimal model (6). In the reduced protocol identification, it is fixed at the same value by extrapolating Ra from 2 to 7 h with a decay constant of 60 and 125 min (obtained from full protocol studies) in OGTT and meal studies, respectively. The second concerns the glucose threshold above which secretion occurs in the C-peptide minimal model. This was fixed to the basal glucose value instead of estimated, since the 2-h data does not have enough information to allow a reliable estimation of the threshold. This approach introduces no appreciable bias, since the estimated threshold is within a few percent of the basal glucose values (2, 5).
Statistical analysis.
Data are presented as means ± SE. Two-sample comparisons were done by Wilcoxon signed-rank test (significance level set to 5%). Pearson’s correlation was used to evaluate univariate correlation. Bland-Altman plot was used to represent the agreement of the methods.
RESULTS
Plasma concentrations.
Average plasma glucose, insulin, and C-peptide concentration in OGTT and meal studies in the basal state are reported in Table 1 together with their SE and range of variability, while their average profiles during OGTT and meal are shown in Fig. 3 (left and right panel, respectively) together with their range of variability (gray area). Range of variability in overall OGTT and meal studies was at basal 4.31–6.94 mmol/l, 9.00–222.60 pmol/l, and 110–2,158 pmol/l and at peak 6.77–16.26 mmol/l, 132–2,028 pmol/l, and 1,130–6,710 pmol/l for glucose, insulin, and C-peptide concentrations, respectively.
Glucose minimal model: Ra and Si.
The model fit in OGTT and meal studies is shown in Fig. 4A and B. Model predictions during reduced protocol are almost superimposable on those of the full protocol in OGTT, while in meal one can note small differences that are essentially due to the higher number of samples present in the first 2 h in meal studies.
The average Ra of ingested glucose during OGTT and meal is shown in Fig. 5A and B. Ra of the reduced protocol well describes the full protocol Ra up to 120 min in both oral tests. Si estimates are shown in Fig. 6 (first panel, left and right, respectively) for the two tests. The reduced protocol Si values are not significantly different and well correlated to those obtained with the full protocol: (reduced versus full) 24.33 ± 1.92 vs. 22.77 ± 1.45 10−5 dl · kg−1 · min−1 per pmol/l, r = 0.89 in OGTT and 19.03 ± 1.38 vs. 19.77 ± 1.20 10−5 dl · kg−1 · min−1 per pmol/l, r = 0.85 in meal. For both tests, the regression line is not different from the identity line (both slopes and zero intercept). The good agreement between the two tests is also evident from Bland-Altman plots, which show that the differences between full and reduced protocol estimates are not related to the size of the measurement (Fig. 6, right upper panel). Precision of Si estimate (expressed as coefficient of variation) was decreased in the reduced protocol with respect to full protocol (18 and 15% vs. 6 and 5% in OGTT and meal, respectively).
C-peptide minimal model: insulin secretion and β-cell responsivity.
The model fit in OGTT and meal studies is shown in Fig. 4C and D. The average insulin secretion (SR) during OGTT and meal is shown in Fig. 5C and D, where its dynamic (SRd) and static (SRs) components of SR are also shown. The reduced and full protocol reconstructed SR, SRd, and SRs are virtually superimposable during both OGTT and meal. Dynamic (Φd) and static (Φs) β-cell responsivity indexes are shown in Fig. 6 for the two tests (second and third panel). The reduced protocol values of Φd and Φs are not significantly different and well correlated to those obtained with the full protocol: Φd is (reduced versus full) 871.50 ± 45.80 vs. 873.32 ± 47.34 10−9, r = 0.98 and 494.88 ± 25.16 vs. 477.99 ± 24.69 10−9, r = 0.91 and Φs is 42.36 ± 1.57 vs. 44.35 ± 1.87 10−9 min−1, r = 0.88 and 35.31 ± 1.13 vs. 35.37 ± 1.12 10−9 min−1, r = 0.90 in OGTT and meal, respectively. The delay constant T is also not significantly different and correlated in the reduced and full protocol: 10.59 ± 0.46 vs. 9.98 ± 0.67 min, r = 0.60 and 12.46 ± 0.77 vs. 13.71 ± 0.83 min, r = 0.54 in OGTT and meal, respectively. Finally, the single overall β-cell responsivity index Φ is 54.57 ± 2.03 vs. 56.43 ± 2.53 10−9 min−1, r = 0.85 and 39.74 ± 1.21 vs. 40.02 ± 1.27 10−9 min−1, r = 0.87 in OGTT and meal, respectively (Fig. 6, lower panel). For both tests, the regression line is not different from the identity line (both slopes and zero intercept). Bland-Altman plots show that the differences between full and reduced protocol estimates are not related to the size of the measurement (Fig. 6, on the right of each panel). Precision of Φd, Φs, Φ, and T estimates (expressed as coefficient of variation) in the reduced protocol was similar to that of the full protocol in both OGTT (21 vs. 24%, 11 vs. 15%, 7 vs. 15%, and 43 vs. 40%, respectively) and meal (30 vs. 28%, 10 vs. 9%, 8 vs. 9%, and 41 vs. 31%, respectively).
DI.
The reduced and full protocol DIs, DId, DIs, and DI, are shown in Fig. 7. They are not significantly different and well correlated. DId is (reduced versus full) 18,477.58 ± 1,432.84 vs. 17,845.13 ± 1,214.89 10−14 dl · kg−1 · min−1 per pmol/l, r = 0.85 and 9,363.59 ± 846.94 vs. 9,427.57 ± 757.15 10−14 dl · kg−1 · min−1 per pmol/l, r = 0.89; DIs is 980.12 ± 84.4 vs. 1,010.00 ± 85.06 10−14 dl · kg−1 · min−2 per pmol/l, r = 0.91 and 636.17 ± 47.40 vs. 674.26 ± 41.90 10−14 dl · kg−1 · min−2 per pmol/l, r = 0.84; and DI is 1,282.26 ± 110.38 vs. 1,273.23 ± 105.49 10−14 dl · kg−1 · min−2 per pmol/l, r = 0.84 and 726.92 ± 55.48 vs. 776.97 ± 51.27 10−14 dl · kg−1 · min−2 per pmol/l, r = 0.84 in OGTT and meal, respectively. For both tests, regression line is not different from the identity line (both slopes and zero intercept). Bland-Altman plots show that the differences between full and reduced protocol estimates are not related to the size of the measurement (Fig. 7, on the right of each panel).
DISCUSSION
There is an increasing demand of highly informative yet relatively simple protocols to assess insulin action and β-cell function in epidemiological and large-scale clinical trials. We have focused on the oral glucose route of delivery, i.e., OGTT or meal in nondiabetic individuals with a large spectrum of glucose tolerance (Table 1, Fig. 3). We have shown that indexes of insulin action and β-cell function can be reliably estimated using the glucose and C-peptide oral minimal models by measuring plasma glucose, insulin, and C-peptide concentrations at 0, 10, 20, 30, 60, 90, and 120 min after glucose or meal ingestion and interpreting these measurements.
The oral route of glucose delivery is clearly more physiological than intravenous glucose injection or continuous infusion of insulin during an hyperinsulinemic clamp. However, measuring insulin action following ingestion of glucose or a mixed meal is more difficult than after intravenous glucose injection. This is because the systemic rate of appearance of exogenous glucose following intravenous glucose injection equals the administered dose, whereas it must be estimated with a model of the glucose system following glucose ingestion. In an effort to avoid this limitation, a new oral glucose minimal model that enables measurement of insulin sensitivity has been developed and validated against multitracer (6) and euglycemic-hyperinsulinemic clamp (7) protocols. In addition, use of a C-peptide minimal model (2, 5) that has been validated both against hyperglycemic clamp (8) and intravenous glucose tolerance test (9) protocols allows concurrent measurement of insulin secretion. This enables determination of whether β-cell function is appropriate for the prevailing level of insulin action (10, 11).
The database used in the present studies consisted of 100 OGTT and 100 meal tolerance tests performed in individuals who had a wide range of glucose tolerance. All 200 individuals underwent what we refer to as a full oral protocol, i.e., a 5-h OGTT with 11 samples or a 7-h meal with 21 samples (2, 9). Plasma measurements of glucose, insulin, and C-peptide interpreted with the full glucose and C-peptide oral minimal models provided reference values to which the shortened protocol were compared. As shown in Figs. 5–7, the full and reduced protocol exogenous glucose appearance rates and indexes of insulin secretion and action were highly correlated during both the OGTT and meal.
The proposed 2-h duration of the reduced protocol is the same as that of a standard OGTT. It also retains the typical 0- and 120-min samples with the addition of samples at 10, 20, 30, 60, and 90 min. Of note, the samples at 10 and 20 min are required for an accurate and precise estimation of the dynamic β-cell responsivity, Φd, which by definition relies on the data where glucose is increasing. When these sample times are included, the ability of the reduced protocol to reconstruct the secretion indexes obtained during the full protocol is remarkable. As is evident in Fig. 6, Φ, Φd, Φs, and T, calculated with the reduced protocol, are virtually identical to those calculated with the full protocol. Of interest, the correlation of T was somewhat lower, likely due to lower precision of its estimation during the reduced protocol. Since insulin sesitivity is also highly correlated, the reduced and full protocols yielded virtually identical estimates of the DIs DI, DId, and DIs, thereby enabling assessment as to whether insulin secretion was appropriate for the degree of insulin resistance
We believe the ability to obtain a virtually identical insulin secretion and action portrait by using only the first 2-h portion of a 5-h OGTT or 7-h meal is due to the predictive power inherent in a model of system. In other words, the two oral minimal models, with modest additional knowledge, i.e., the extrapolation of the rate of appearance of ingested glucose, Ra, beyond the 2 h for the glucose model, only need the first 2 h of information to provide an accurate picture in a nondiabetic population. To better appreciate the predictive power of the oral minimal model method, it is of interest to contrast it with an AUC approach. This is possible because, e.g., dynamic (Φd) and static (Φs) β-cell responsivity indexes are also amenable to an AUC interpretation. It is easy to show (2, 5) that Φd, besides being a parameter of the C-peptide oral minimal model, also represents the AUC of SRd per unit increase of glucose concentration; similarly, Φs is the AUC of above-basal SRs per AUC of above-basal glucose concentration. We are thus in the position to compare the reduced with the full protocol AUC values of Φd and Φs calculated directly from the data. Of interest, when calculated in this manner, AUC Φs values of the reduced protocol are in this case significantly higher that those of the full protocol: 50.19 ± 1.97 vs. 44.34 ± 1.87 (P < 0.0001) and 39.36 ± 1.29 vs. 35.37 ± 1.12 min−1 (P < 0.0001) in OGTT and meal, respectively. Also, correlation deteriorated (r = 0.70) in both OGTT and meal with respect to C-peptide model Φs.
In conclusion, we have shown that seven samples obtained during the first 2 h after ingestion of either 75 g glucose or a mixed meal enables accurate assessment of both insulin secretion and action, thereby facilitating the conduct of large-scale epidemiologic studies and clinical trials. Whether a similar approach can be used in individuals with impaired insulin secretion (e.g., diabetes) awaits further study.
Article Information
This study was supported by the National Insitutes for Health Grants EB-01975 and DK29953-23.