Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Evaluation of Fasting State-/Oral Glucose Tolerance Test-Derived Measures of Insulin Release for the Detection of Genetically Impaired β-Cell Function

  • Silke A. Herzberg-Schäfer ,

    Contributed equally to this work with: Silke A. Herzberg-Schäfer, Harald Staiger

    Affiliation Division of Endocrinology, Diabetology, Angiology, Nephrology, and Clinical Chemistry, Department of Internal Medicine, University Hospital Tübingen and Eberhard Karls University Tübingen, Member of the German Center for Diabetes Research (DZD), Tübingen, Germany

  • Harald Staiger ,

    Contributed equally to this work with: Silke A. Herzberg-Schäfer, Harald Staiger

    Affiliation Division of Endocrinology, Diabetology, Angiology, Nephrology, and Clinical Chemistry, Department of Internal Medicine, University Hospital Tübingen and Eberhard Karls University Tübingen, Member of the German Center for Diabetes Research (DZD), Tübingen, Germany

  • Martin Heni,

    Affiliation Division of Endocrinology, Diabetology, Angiology, Nephrology, and Clinical Chemistry, Department of Internal Medicine, University Hospital Tübingen and Eberhard Karls University Tübingen, Member of the German Center for Diabetes Research (DZD), Tübingen, Germany

  • Caroline Ketterer,

    Affiliation Division of Endocrinology, Diabetology, Angiology, Nephrology, and Clinical Chemistry, Department of Internal Medicine, University Hospital Tübingen and Eberhard Karls University Tübingen, Member of the German Center for Diabetes Research (DZD), Tübingen, Germany

  • Martina Guthoff,

    Affiliation Division of Endocrinology, Diabetology, Angiology, Nephrology, and Clinical Chemistry, Department of Internal Medicine, University Hospital Tübingen and Eberhard Karls University Tübingen, Member of the German Center for Diabetes Research (DZD), Tübingen, Germany

  • Konstantinos Kantartzis,

    Affiliation Division of Endocrinology, Diabetology, Angiology, Nephrology, and Clinical Chemistry, Department of Internal Medicine, University Hospital Tübingen and Eberhard Karls University Tübingen, Member of the German Center for Diabetes Research (DZD), Tübingen, Germany

  • Fausto Machicao,

    Affiliation Division of Endocrinology, Diabetology, Angiology, Nephrology, and Clinical Chemistry, Department of Internal Medicine, University Hospital Tübingen and Eberhard Karls University Tübingen, Member of the German Center for Diabetes Research (DZD), Tübingen, Germany

  • Norbert Stefan,

    Affiliation Division of Endocrinology, Diabetology, Angiology, Nephrology, and Clinical Chemistry, Department of Internal Medicine, University Hospital Tübingen and Eberhard Karls University Tübingen, Member of the German Center for Diabetes Research (DZD), Tübingen, Germany

  • Hans-Ulrich Häring ,

    hans-ulrich.haering@med.uni-tuebingen.de

    Affiliation Division of Endocrinology, Diabetology, Angiology, Nephrology, and Clinical Chemistry, Department of Internal Medicine, University Hospital Tübingen and Eberhard Karls University Tübingen, Member of the German Center for Diabetes Research (DZD), Tübingen, Germany

  • Andreas Fritsche

    Affiliations Division of Endocrinology, Diabetology, Angiology, Nephrology, and Clinical Chemistry, Department of Internal Medicine, University Hospital Tübingen and Eberhard Karls University Tübingen, Member of the German Center for Diabetes Research (DZD), Tübingen, Germany, Division of Nutritional and Preventive Medicine, Department of Internal Medicine, University Hospital Tübingen, Tübingen, Germany

Abstract

Background

To date, fasting state- and different oral glucose tolerance test (OGTT)-derived measures are used to estimate insulin release with reasonable effort in large human cohorts required, e.g., for genetic studies. Here, we evaluated twelve common (or recently introduced) fasting state-/OGTT-derived indices for their suitability to detect genetically determined β-cell dysfunction.

Methodology/Principal Findings

A cohort of 1364 White European individuals at increased risk for type 2 diabetes was characterized by OGTT with glucose, insulin, and C-peptide measurements and genotyped for single nucleotide polymorphisms (SNPs) known to affect glucose- and incretin-stimulated insulin secretion. One fasting state- and eleven OGTT-derived indices were calculated and statistically evaluated. After adjustment for confounding variables, all tested SNPs were significantly associated with at least two insulin secretion measures (p≤0.05). The indices were ranked according to their associations' statistical power, and the ranks an index obtained for its associations with all the tested SNPs (or a subset) were summed up resulting in a final ranking. This approach revealed area under the curve (AUC)Insulin(0-30)/AUCGlucose(0-30) as the best-ranked index to detect SNP-dependent differences in insulin release. Moreover, AUCInsulin(0-30)/AUCGlucose(0-30), corrected insulin response (CIR), AUCC-Peptide(0-30)/AUCGlucose(0-30), AUCC-Peptide(0-120)/AUCGlucose(0-120), two different formulas for the incremental insulin response from 0–30 min, i.e., the insulinogenic indices (IGI)2 and IGI1, and insulin 30 min were significantly higher-ranked than homeostasis model assessment of β-cell function (HOMA-B; p<0.05). AUCC-Peptide(0-120)/AUCGlucose(0-120) was best-ranked for the detection of SNPs involved in incretin-stimulated insulin secretion. In all analyses, HOMA-β displayed the highest rank sums and, thus, scored last.

Conclusions/Significance

With AUCInsulin(0-30)/AUCGlucose(0-30), CIR, AUCC-Peptide(0-30)/AUCGlucose(0-30), AUCC-Peptide(0-120)/AUCGlucose(0-120), IGI2, IGI1, and insulin 30 min, dynamic measures of insulin secretion based on early insulin and C-peptide responses to oral glucose represent measures which are more appropriate to assess genetically determined β-cell dysfunction than fasting measures, i.e., HOMA-B. Genes predominantly influencing the incretin axis may possibly be best detected by AUCC-Peptide(0-120)/AUCGlucose(0-120).

Introduction

Recently, genome-wide association (GWA) scans in tens of thousands of human cases and controls using high-density single nucleotide polymorphism (SNP) arrays and subsequent meta-analyses of these data provided important new insights into the genetic architecture of complex diseases [1]. In the course of these studies, a series of nearly 20 novel type 2 diabetes risk loci were identified. In smaller but extensively and thoroughly phenotyped cohorts, many of the diabetogenic alleles were shown to affect β-cell function [2]. Despite this recent scientific progress, a shortcoming of the genetic findings up to now is that the sum of all reported common GWA-derived risk alleles only marginally improves the prediction of future type 2 diabetes, when combined with established clinical parameters, and only explains about 6% of the heritability of the disease [3]. Thus, it is anticipated that many more loci remain to be discovered that act in an additive or even synergistic manner to increase the type 2 diabetes risk. Amongst others, the following strategies are currently discussed to find them: (i) use of SNP arrays of higher density, (ii) assessment of rare variants, and (iii) realization of GWA analyses using quantitative traits known to be crucially involved in the pathogenesis of type 2 diabetes, such as insulin secretion and insulin sensitivity [2]; [3].

One possibility to identify more loci affecting β-cell function is to determine insulin release in cohorts large enough to allow reliable genetic analyses. To estimate insulin release in such cohorts with reasonable effort, i.e., at low expenses in time and costs, fasting state- and several different oral glucose tolerance test (OGTT)-derived indices calculated from plasma insulin, C-peptide, and glucose concentrations are available [4][10]. However, which of these indices are best-suited to detect genetically determined β-cell dysfunction is unknown. Therefore, we evaluated, in this study, fasting state- (homeostasis model assessment of β-cell function [HOMA-B]) and OGTT-derived indices (insulin and C-peptide concentrations at 30 min of OGTT, insulinogenic indices [IGIs], area under the curve [AUC]Insulin(0-30)/AUCGlucose(0-30), AUCC-Peptide(0-30)/AUCGlucose(0-30), AUCInsulin(0-120)/AUCGlucose(0-120), AUCC-Peptide(0-120)/AUCGlucose(0-120), oral disposition index [DI oral], corrected insulin response [CIR], and first-phase insulin secretion) for their suitability to detect altered insulin release due to confirmed type 2 diabetes risk SNPs convincingly described to affect specific aspects of β-cell function, such as glucose-stimulated insulin secretion (GSIS), incretin-stimulated insulin secretion (ISIS), or incretin release. For this investigation, we included the type 2 diabetes risk loci/SNPs MTNR1B rs10830963, HHEX rs7923837, CDKAL1 rs7754840, TCF7L2 rs7903146, WFS1 rs10010131, and KCNQ1 rs151290.

Materials and Methods

Ethics statement

From all participants, informed written consent to the study was obtained, and the Ethics Committee of the Medical Faculty of the University of Tübingen approved the study protocol.

Subjects

A cohort of 1364 White individuals was recruited from the ongoing Tübingen family study for type 2 diabetes (TÜF) that currently encompasses ∼2000 participants at increased risk for type 2 diabetes (non-diabetic individuals from Southern Germany with family history of type 2 diabetes or diagnosis of impaired fasting glycaemia) [11]. More than 99.5% of the TÜF participants are of European ancestry. All participants underwent the standard procedures of the study protocol including medical history and physical examination, assessment of smoking status and alcohol consumption habits, routine blood tests, and an OGTT. Selection of the present study cohort was based on the absence of newly diagnosed diabetes and the availability of complete sets of clinical and genotype data. Moreover, the participants were not taking any medication known to affect glucose tolerance or insulin secretion. The subject characteristics are given in Table 1. From this cohort, a subset of 274 individuals additionally underwent a frequently sampled intravenous glucose tolerance test (IVGTT).

thumbnail
Table 1. Subject characteristics including the fasting state- and OGTT-derived indices of insulin release (N = 1364).

https://doi.org/10.1371/journal.pone.0014194.t001

OGTT and IVGTT

A standard 75-g OGTT was performed after a ten-hour overnight fast, and venous blood samples were drawn at time-points zero, 30, 60, 90, and 120 min for the determination of plasma glucose, insulin, and C-peptide concentrations. In those individuals who agreed to undergo the IVGTT, baseline samples (-10, -5, and 0 min) were collected before a glucose dose of 0.3 g/kg body weight was given. Blood samples for the measurement of plasma glucose and insulin were obtained at 2, 4, 6, 8, 10, 20, 30, 40, 50, and 60 min.

Laboratory measurements

Plasma glucose was determined using a bedside glucose analyzer (glucose oxidase method, Yellow Springs Instruments, Yellow Springs, OH, USA). Plasma insulin and C-peptide concentrations were measured by commercial chemiluminescence assays for ADVIA Centaur (Siemens Medical Solutions, Fernwald, Germany) according to the manufacturer's instructions.

Selection of loci/SNPs

From each confirmed type 2 diabetes risk locus previously reported to affect specific aspects of β-cell function, we selected one representative SNP based on the availability of genotype data and on the robustness of the SNP's β-cell effect in our cohort. As loci/SNPs associated with GSIS, we selected MTNR1B rs10830963 [12]; [13], HHEX rs7923837 [14]; [15], and CDKAL1 rs7754840 [16]; [17]. As loci/SNPs predominantly associated with ISIS/incretin release, we selected TCF7L2 rs7903146 [18]; [19], WFS1 rs10010131 [20], and KCNQ1 rs151290 [21]. All SNPs were genotyped in the whole cohort in the course of earlier studies [12]; [15]; [17]; [18]; [20]; [21] using TaqMan assays (Applied Biosystems, Foster City, CA, USA) and passed the quality controls. Details on this as well as on minor allele frequencies, genotyping success rates, and Hardy-Weinberg equilibrium are reported in the aforementioned references.

Calculations

Insulin secretion derived from the fasting state was calculated as HOMA-B: 20·I0/(G0–3.5) with I0 = fasting insulin in µU/ml and G0 = fasting glucose in mmol/l [5]. All other insulin secretion indices were derived from the OGTT with insulin and C-peptide concentrations given in pmol/l, and glucose concentration given in mmol/l. AUCs of insulin, C-peptide, and glucose concentrations during the entire 120 min of the OGTT were calculated according to the trapezoid method as: 0.5·(0.5·c0+c30+c60+c90+0.5·c120) with c = concentration. AUCInsulin(0-30)/AUCGlucose(0-30) was calculated as: (I0+I30)/(G0+G30) [9]. AUCC-Peptide(0-30)/AUCGlucose(0-30) was calculated analogously. IGI1 was calculated as: (I30–I0)/(G30–G0) [10]. IGI2 was calculated as: (I30–I0)/G30 [6]. DI oral was calculated as: IGI1/I0 [8]. CIR was calculated as: 100·I30/[G30·(G30–3.89)] [4]. First-phase insulin secretion was calculated as: 1283+1.829·I30–138.7·G30+3.772·I0 [7]. Insulin sensitivity derived from the OGTT was estimated as proposed by Matsuda and DeFronzo [22]: 10000/(G0·I0·Gmean·Imean)½. Fasting insulin clearance was calculated as CP0/I0 with CP0 = fasting C-peptide, insulin clearance during the OGTT was calculated as AUCC-Peptide(0-120)/AUCInsulin(0-120). Acute insulin response (AIR) derived from the IVGTT was used as gold standard for the assessment of insulin secretion and calculated as: 0.5·(0.5·I0+I2+I4+I6+I8+0.5·I10).

Statistical analyses

Prior to analysis, all continuous data were loge-transformed in order to approximate normal distribution. Multiple linear regression analysis was performed using the least-squares method. In the regression models, the insulin secretion parameter was chosen as dependent variable, the SNP genotype (additive inheritance model) as independent variable, and gender, age, BMI, and OGTT-derived insulin sensitivity as confounding variables. In addition, the SNP genotype (additive inheritance model) was tested as dependent variable and the insulin secretion parameter as independent variable with inclusion of the aforementioned confounders in the models. Since the critical confounding variables age, BMI, and OGTT-derived insulin sensitivity did not achieve normal distribution even after applying the ladder of powers (probably due to the inclusion/exclusion criteria of our study), we additionally performed linear regression models including these parameters as nominal variables after stratification into quartiles. A p-value ≤0.05 was considered statistically significant. Multiple linear regression analyses, post hoc power calculations [statistical power (1-β) and least significant number (lsn; i.e., the sample size expected to be needed to achieve statistical significance) to detect the effect size given by the default settings (square root of the sum of squares for the hypothesis divided by N)], and Wilcoxon rank sum tests were performed using the statistical software package JMP 4.0 (SAS Institute, Cary, NC, USA).

Results

Seventeen subjects with calculated negative values in single insulin secretion measures were excluded from all analyses resulting in a final cohort of 1347 individuals. Since the phenotype (insulin release) is determined by the genotype, we started our analyses using the insulin secretion index as dependent variable and the SNP genotype (additive inheritance model) as independent variable. As expected, all tested loci/SNPs were significantly associated with at least two of the indices after adjustment for the confounding variables gender, age, BMI, and OGTT-derived insulin sensitivity (p≤0.05, Table 2; additional statistical data given in Table S1). Most secretion indices identified three, four, or five of the six tested loci/SNPs to be significantly associated with insulin release, whereas HOMA-B detected MTNR1B rs10830963 only. Inclusion of the confounding parameters age, BMI, and OGTT-derived insulin sensitivity as nominal variables (after stratification into quartiles) in the linear regression models resulted in very similar statistical data (Table S2). After adjustment for gender, age, and BMI, none of the SNPs showed significant association with insulin clearance either in the fasting state (p>0.1) or during the OGTT (p≥0.06).

thumbnail
Table 2. Statistical data of the SNPs' associations with indices of insulin release and ranking of the indices according to their lsn.

https://doi.org/10.1371/journal.pone.0014194.t002

To evaluate which indices are most appropriate to detect genetically determined differences in insulin release, we first calculated the post hoc least significant numbers for all associations and converted them into ranks with indices that displayed the lowest least significant number being the best-ranked (Table 2). Then, we summed up the ranks of each insulin secretion index obtained for all the tested SNPs and ranked the indices according to their rank sums (Table 3). Using this approach, AUCInsulin(0-30)/AUCGlucose(0-30) was identified as the best-ranked index (Table 3). Moreover, AUCInsulin(0-30)/AUCGlucose(0-30), CIR, AUCC-Peptide(0-30)/AUCGlucose(0-30), AUCC-Peptide(0-120)/AUCGlucose(0-120), IGI2, IGI1, and insulin 30 min, but not C-peptide 30 min, first-phase insulin secretion, AUCInsulin(0-120)/AUCGlucose(0-120), or DI oral, were significantly higher-ranked than HOMA-B (p<0.05; Table 3). To avoid over-adjustment of DI oral, a secretion parameter already normalised for a rough estimate of insulin sensitivity (i.e., fasting insulin), this parameter was also tested in the absence of additional adjustment for OGTT-derived insulin sensitivity. This analysis resulted in a somewhat higher rank sum (57) that, however, had no impact on this index' overall rank (rank 11). When summing up the ranks of the indices obtained for the three loci/SNPs affecting GSIS, i.e., MTNR1B rs10830963, HHEX rs7923837, and CDKAL1 rs7754840, AUCInsulin(0-30)/AUCGlucose(0-30) again turned out to be the highest-ranked index (Table 3). Notably, when summing up the ranks obtained for the three loci/SNPs predominantly affecting ISIS, i.e., TCF7L2 rs7903146, WFS1 rs10010131, and KCNQ1 rs151290, AUCC-Peptide(0-120)/AUCGlucose(0-120) was the best-ranked index (Table 3). In the GSIS and ISIS subgroups, statistical analysis of the rankings was inappropriate due to the small sample sizes. In all rankings, HOMA-β displayed the highest rank sums and, thus, represented the lowest-ranked index. Assessing the SNP genotype (additive inheritance model) as dependent variable and the insulin secretion parameter as independent variable with inclusion of the aforementioned confounders in the multiple regression models yielded very similar rankings (Tables S3 and S4).

thumbnail
Table 3. Ranking of the indices of insulin release according to their rank sums.

https://doi.org/10.1371/journal.pone.0014194.t003

Interestingly, the indices that performed best in all these analyses, i.e., AUCInsulin(0-30)/AUCGlucose(0-30) and CIR, also revealed the best correlations with IVGTT-derived AIR (both r = 0.76), and HOMA-B, the lowest ranked index, showed the weakest correlation with AIR (r = 0.64, N = 274; Table 4).

thumbnail
Table 4. Association of the fasting- and OGTT-derived indices of insulin release with IVGTT-derived AIR (N = 274).

https://doi.org/10.1371/journal.pone.0014194.t004

Discussion

In this study, we intended to identify, among twelve fasting state- and common (or recently introduced) OGTT-derived measures feasible for genetic studies in large cohorts, the indices best-suited to detect genetically determined alterations of insulin release. Since the suitability of the indices for detection of altered β-cell function may depend on the SNPs' pathomechanisms, we additionally analysed the SNPs affecting GSIS separately from those affecting the incretin axis (ISIS or incretin release). It was not the primary aim of this study to evaluate the performance of the fasting state- and OGTT-derived estimates of insulin secretion by comparing them with gold standard measures derived from laborious and expensive methods, such as IVGTT or hyperglycemic clamp.

Using summation of the ranks derived from post hoc least significant numbers, we show here that AUCInsulin(0-30)/AUCGlucose(0-30), a recently proposed index validated against first-phase insulin release in a frequently sampled IVGTT [9], represents the best-ranked index for the detection of SNP effects on overall insulin release as well as on GSIS. By contrast, loci/SNPs affecting the incretin axis may be better captured by AUCC-Peptide(0-120)/AUCGlucose(0-120). One explanation for this divergent result is that plasma concentrations of incretins, as compared to plasma glucose, do not rapidly decline after having reached their maximum during the first 60 min of the OGTT, but remain elevated until the end of the protocol [18]; [21]. Thus, OGTT-induced levels of incretins, compared to glucose, may exert more prolonged effects on the β-cell which are best assessed using indices covering the entire OGTT period. The observation that AUCInsulin(0-120)/AUCGlucose(0-120) is not among the best-suited indices to detect alterations of the incretin axis may be due to the shorter circulating half-life of insulin as compared to C-peptide [23].

AUCInsulin(0-30)/AUCGlucose(0-30), CIR, AUCC-Peptide(0-30)/AUCGlucose(0-30), AUCC-Peptide(0-120)/AUCGlucose(0-120), IGI2, IGI1, and insulin 30 min significantly outperformed HOMA-B in the detection of genetically determined differences in overall insulin release, and thus are clearly superior to HOMA-B in this regard. Since the ranks of C-peptide 30 min, first-phase insulin secretion, AUCInsulin(0-120)/AUCGlucose(0-120), and DI oral were statistically indistinguishable from that of HOMA-B – and HOMA-B displayed the lowest ranks in all analyses –, these indices may not be recommended for genetic studies aimed at the identification of novel loci/SNPs affecting β-cell function. The validity of the OGTT-derived index AUCInsulin(0-30)/AUCGlucose(0-30) as a preferable proxy for the assessment of β-cell function in large genetic studies is underscored by its strong correlation with IVGTT-derived AIR. In these latter analyses with IVGTT-derived AIR as gold standard, HOMA-B again revealed the weakest correlation and, thus, was confirmed to be less useful for the detection of impaired β-cell function.

A limitation of our study is that the results were generated in a single study population and, thus, clearly need replication in other comparably genotyped and phenotyped cohorts of similar or larger sample size. Furthermore, the ranking of the insulin secretion indices may also depend on the ethnicity. Since our study cohort was nearly exclusively comprised of White European subjects, similar analyses in other ethnicities would be interesting.

Finally, we conclude that, according to our data, AUCInsulin(0-30)/AUCGlucose(0-30), along with CIR, AUCC-Peptide(0-30)/AUCGlucose(0-30), AUCC-Peptide(0-120)/AUCGlucose(0-120), IGI2, IGI1, and insulin 30 min, represents an appropriate surrogate parameter to assess genetically determined β-cell dysfunction. HOMA-B, DI oral, AUCInsulin(0-120)/AUCGlucose(0-120), first-phase insulin secretion, and C-peptide 30 min, however, are of limited informative value for genetic studies on β-cell function in humans. The influence of genes on ISIS or incretin release may possibly be better detected by calculating AUCC-Peptide/AUCGlucose. These findings, if replicated in comparably sized and phenotyped cohorts, should facilitate the identification of novel loci/SNPs affecting insulin release in large cohorts metabolically characterized by OGTT with glucose, insulin, and C-peptide measurements.

Supporting Information

Table S1.

Additional statistical data of the SNPs' associations with indices of insulin release. Given are the estimate and the standard deviation of the minor allele's effect. Seventeen subjects with calculated negative values were excluded (N = 1347). Prior to multiple linear regression analysis, all continuous variables were loge-transformed to approximate normal distribution. In the multiple linear regression models, the insulin secretion parameter was chosen as dependent variable, the SNP genotype (additive inheritance model) as independent variable and gender, age, BMI, and OGTT-derived insulin sensitivity as confounding variables. AUC - area under the curve; BMI - body mass index; CIR - cleared insulin response; DI - disposition index; HOMA-B - homeostasis model assessment of beta-cell function; IGI - insulinogenic index; SD - standard deviation; SNP - single nucleotide polymorphism.

https://doi.org/10.1371/journal.pone.0014194.s001

(0.07 MB DOC)

Table S2.

Statistical data of the SNPs' associations with indices of insulin release using the covariates age, BMI, and OGTT-derived insulin sensitivity as nominal variables (stratified in quartiles). Given are the p-value, estimate and the standard deviation of the minor allele's effect. Seventeen subjects with calculated negative values were excluded (N = 1347). In the multiple linear regression models, the insulin secretion parameter was chosen as dependent variable, the SNP genotype (additive inheritance model) as independent variable and gender, age, BMI, and OGTT-derived insulin sensitivity as confounding variables. AUC - area under the curve; BMI - body mass index; CIR - cleared insulin response; DI - disposition index; HOMA-B - homeostasis model assessment of beta-cell function; IGI - insulinogenic index; SD - standard deviation; SNP - single nucleotide polymorphism.

https://doi.org/10.1371/journal.pone.0014194.s002

(0.07 MB DOC)

Table S3.

Statistical data of the SNPs' associations with indices of insulin release using the genotype as dependent variable. Seventeen subjects with calculated negative values were excluded (N = 1347). Prior to multiple linear regression analysis, all continuous variables were loge-transformed to approximate normal distribution. In the multiple linear regression models, the SNP genotype (additive inheritance model) was chosen as dependent variable, the insulin secretion parameter as independent variable and gender, age, BMI, and OGTT-derived insulin sensitivity as confounding variables. AUC - area under the curve; BMI - body mass index; CIR - cleared insulin response; DI - disposition index; HOMA-B - homeostasis model assessment of beta-cell function; IGI - insulinogenic index; lsn - least significant number (sample size expected to be needed to achieve statistical significance); SNP - single nucleotide polymorphism.

https://doi.org/10.1371/journal.pone.0014194.s003

(0.08 MB DOC)

Table S4.

Ranking of the indices of insulin release according to their rank sums derived from the statistics presented in Supplemental Table S3 (genotype as dependent variable). * Significantly different from HOMA-B (p<0.05; Wilcoxon rank sum test). AUC - area under the curve; CIR - cleared insulin response; DI - disposition index; GSIS - glucose-stimulated insulin secretion; HOMA-B - homeostasis model assessment of beta-cell function; IGI - insulinogenic index; ISIS - incretin-stimulated insulin secretion/incretin production; lsn - least significant number; SNP - single nucleotide polymorphism.

https://doi.org/10.1371/journal.pone.0014194.s004

(0.06 MB DOC)

Acknowledgments

We thank all study participants for their cooperation. We gratefully acknowledge the excellent technical assistance of Anna Bury, Heike Luz, Alke Guirguis, Melanie Weisser, and Roman-Georg Werner.

Author Contributions

Conceived and designed the experiments: NS HUH AF. Performed the experiments: SAHS MH CK MG KK FM. Analyzed the data: SAHS HS. Contributed reagents/materials/analysis tools: HUH. Wrote the paper: SAHS HS. Contributed to discussion: SAHS HS FM NS HUH AF. Reviewed/edited manuscript: NS HUH AF.

References

  1. 1. McCarthy MI, Abecasis GR, Cardon LR, Goldstein DB, Little J, et al. (2008) Genome-wide association studies for complex traits: consensus, uncertainty and challenges. Nat Rev Genet 9: 356–369.
  2. 2. Staiger H, Machicao F, Fritsche A, Haring HU (2009) Pathomechanisms of type 2 diabetes genes. Endocr Rev 30: 557–585.
  3. 3. Manolio TA, Collins FS, Cox NJ, Goldstein DB, Hindorff LA, et al. (2009) Finding the missing heritability of complex diseases. Nature 461: 747–753.
  4. 4. Sluiter WJ, Erkelens DW, Reitsma WD, Doorenbos H (1976) Glucose tolerance and insulin release, a mathematical approach I. Assay of the beta-cell response after oral glucose loading. Diabetes 25: 241–244.
  5. 5. Matthews DR, Hosker JP, Rudenski AS, Naylor BA, Treacher DF, et al. (1985) Homeostasis model assessment: insulin resistance and beta-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia 28: 412–419.
  6. 6. Wareham NJ, Phillips DI, Byrne CD, Hales CN (1995) The 30 minute insulin incremental response in an oral glucose tolerance test as a measure of insulin secretion. Diabet Med 12: 931.
  7. 7. Stumvoll M, Mitrakou A, Pimenta W, Jenssen T, Yki-Jarvinen H, et al. (2000) Use of the oral glucose tolerance test to assess insulin release and insulin sensitivity. Diabetes Care 23: 295–301.
  8. 8. Utzschneider KM, Prigeon RL, Faulenbach MV, Tong J, Carr DB, et al. (2009) Oral disposition index predicts the development of future diabetes above and beyond fasting and 2-h glucose levels. Diabetes Care 32: 335–341.
  9. 9. Stancakova A, Javorsky M, Kuulasmaa T, Haffner SM, Kuusisto J, et al. (2009) Changes in insulin sensitivity and insulin release in relation to glycemia and glucose tolerance in 6,414 Finnish men. Diabetes 58: 1212–1221.
  10. 10. Kadowaki T, Miyake Y, Hagura R, Akanuma Y, Kajinuma H, et al. (1984) Risk factors for worsening to diabetes in subjects with impaired glucose tolerance. Diabetologia 26: 44–49.
  11. 11. Stefan N, Machicao F, Staiger H, Machann J, Schick F, et al. (2005) Polymorphisms in the gene encoding adiponectin receptor 1 are associated with insulin resistance and high liver fat. Diabetologia 48: 2282–2291.
  12. 12. Staiger H, Machicao F, Schafer SA, Kirchhoff K, Kantartzis K, et al. (2008) Polymorphisms within the novel type 2 diabetes risk locus MTNR1B determine beta-cell function. PLoS ONE 3: e3962.
  13. 13. Sparso T, Bonnefond A, Andersson E, Bouatia-Naji N, Holmkvist J, et al. (2009) G-allele of intronic rs10830963 in MTNR1B confers increased risk of impaired fasting glycemia and type 2 diabetes through an impaired glucose-stimulated insulin release: studies involving 19,605 Europeans. Diabetes 58: 1450–1456.
  14. 14. Grarup N, Rose CS, Andersson EA, Andersen G, Nielsen AL, et al. (2007) Studies of association of variants near the HHEX, CDKN2A/B, and IGF2BP2 genes with type 2 diabetes and impaired insulin release in 10,705 Danish subjects: validation and extension of genome-wide association studies. Diabetes 56: 3105–3111.
  15. 15. Staiger H, Machicao F, Stefan N, Tschritter O, Thamer C, et al. (2007) Polymorphisms within novel risk loci for type 2 diabetes determine beta-cell function. PLoS ONE 2: e832.
  16. 16. Stancakova A, Pihlajamaki J, Kuusisto J, Stefan N, Fritsche A, et al. (2008) Single-nucleotide polymorphism rs7754840 of CDKAL1 is associated with impaired insulin secretion in nondiabetic offspring of type 2 diabetic subjects and in a large sample of men with normal glucose tolerance. J Clin Endocrinol Metab 93: 1924–1930.
  17. 17. Groenewoud MJ, Dekker JM, Fritsche A, Reiling E, Nijpels G, et al. (2008) Variants of CDKAL1 and IGF2BP2 affect first-phase insulin secretion during hyperglycaemic clamps. Diabetologia 51: 1659–1663.
  18. 18. Schafer SA, Tschritter O, Machicao F, Thamer C, Stefan N, et al. (2007) Impaired glucagon-like peptide-1-induced insulin secretion in carriers of transcription factor 7-like 2 (TCF7L2) gene polymorphisms. Diabetologia 50: 2443–2450.
  19. 19. Lyssenko V, Lupi R, Marchetti P, Del Guerra S, Orho-Melander M, et al. (2007) Mechanisms by which common variants in the TCF7L2 gene increase risk of type 2 diabetes. J Clin Invest 117: 2155–2163.
  20. 20. Schafer SA, Mussig K, Staiger H, Machicao F, Stefan N, et al. (2009) A common genetic variant in WFS1 determines impaired glucagon-like peptide-1-induced insulin secretion. Diabetologia 52: 1075–1082.
  21. 21. Mussig K, Staiger H, Machicao F, Kirchhoff K, Guthoff M, et al. (2009) Association of type 2 diabetes candidate polymorphisms in KCNQ1 with incretin and insulin secretion. Diabetes 58: 1715–1720.
  22. 22. Matsuda M, DeFronzo RA (1999) Insulin sensitivity indices obtained from oral glucose tolerance testing: comparison with the euglycemic insulin clamp. Diabetes Care 22: 1462–1470.
  23. 23. Polonsky KS, Rubenstein AH (1984) C-peptide as a measure of the secretion and hepatic extraction of insulin. Pitfalls and limitations. Diabetes 33: 486–494.