© 2003 by the American Diabetes Association, Inc.
Analysis of the Relationship Between the Pro12Ala Variant in the PPAR-
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
(PPAR-
) gene variants on the response rate to therapy with the thiazolidinedione (TZD) pioglitazone, because in vitro studies have suggested that genetic variants of the PPAR-
gene may influence the drug efficacy of TZD.
RESEARCH DESIGN AND METHODSA total of 131 patients were treated in an open-label, randomized, multicenter study with pioglitazone (45 mg o.d.) during a course of
26 weeks. Response to the pioglitazone therapy was defined by either a >20% decrease in fasting plasma glucose or a >15% decrease in HbA1c values after 26 weeks of pioglitazone treatment. We evaluated the association between the PPAR-
genotype and the response rate to pioglitazone treatment.
RESULTSThe Pro12Ala and the Pro12Pro variants in the PPAR-
gene are not associated with the response rate to pioglitazone treatment in patients with type 2 diabetes. However, we identified initial fasting plasma glucose level >11.0 mmol/l, HbA1c value >9.0%, BMI >32 kg/m2, and fasting C-peptide concentrations at baseline >2.5 pmol/l as predominant confounding factors for the responder frequency to pioglitazone treatment.
CONCLUSIONSThe Pro12Ala variant in the PPAR-
gene does not affect the therapy efficacy of pioglitazone, suggesting that the drug-treatment response is independent from pharmacogenetic effects between PPAR-
and its ligand pioglitazone. Whether the Ala12Ala genotype plays a role in the response rate to TZD therapy remains to be determined.
Abbreviations: DGGE, denaturing gradient gel electrophoresis PPAR-
, peroxisome proliferator-activated receptor-
TZD, thiazolidinedione
| INTRODUCTION |
|---|
|
|
|---|
Therefore, heterogeneity at the molecular level provides an important explanation for differential individual drug response and might also explain the nonresponder rate or an unexpectedly high response to drug treatment.
The interaction of thiazolidinedione (TZD) with the peroxisome proliferator-activated receptor-
(PPAR-
) represents another model of pharmacogenomics. TZD has been shown to decrease plasma glucose concentrations in patients with type 2 diabetes (911). Clinical studies showed that
1025% of the patients treated with TZD did not achieve a 15% reduction in fasting plasma glucose or did not convert from impaired glucose tolerance to normal glucose tolerance (12). However, the molecular reasons for the different responses to TZD therapy have not been determined, and differences in the PPAR-
genotype may modify the response to TZD treatment. Therefore, we investigated whether polymorphisms in the PPAR-
gene are a cause for different responses or therapy resistance in 131 patients with type 2 diabetes treated for at least 26 weeks with the TZD pioglitazone. We compared the frequencies of the previously described Pro115Gln (13), Pro12Ala (14), Val290Met, Pro467Leu (15), and one silent polymorphism (16) between a group of patients responsive to the therapy with pioglitazone, defined by either a >20% decrease in fasting plasma glucose concentrations or a >15% decrease in HbA1c values after 26 weeks of TZD treatment, and a group of patients resistant to the therapy. Moreover, with the aim of defining the best responders to TZD therapy on the molecular level, we evaluated the association between the PPAR-
genotype and a high response to pioglitazone treatment in a subgroup of individuals who achieved a >40% decrease in fasting plasma glucose or HbA1c values with pioglitazone treatment.
| RESEARCH DESIGN AND METHODS |
|---|
|
|
|---|
7.5 and
11.5% and fasting serum glucose
7.8 and
14.0 mmol/l were chosen according to the following inclusion criteria: 1) men or women aged
35 or
75 years; 2) BMI
25 or
35 kg/m2; 3) written informed consent given; and 4) women who were postmenopausal, hysterectomized, surgically sterilized, or using appropriate contraceptive methods. In addition, patients with type 1 diabetes, current exposure to antidiabetic medication or necessary concomitant therapy with other antidiabetic drugs, any history of ketoacidosis, any disease with malabsorption, acute or chronic pancreatitis, chronic gastrointestinal diseases with significant digestion and resorption disturbances, conditions that may be aggravated by increased gas in the intestine, any history of unstable angina pectoris or any other unstable heart disease, heart failure of New York Heart Association (NYHA) stage II-IV, history of any malignant disease within the last 10 years, any infection with HIV, significant anemia of any etiology or any other clinically relevant hematological disease, any serious preexisting impairment of liver, kidney, or bone marrow function, >3 years of phenacetin use, donation of blood during the study, any history of alcohol or drug abuse, pregnancy/lactation, as well as use of medications with a laxative or constipating effect or affecting gastrointestinal motility and/or absorption were excluded from the study. From the initially enrolled 532 patients who initially met the inclusion criteria, 264 patients were withdrawn from the study because at least one of the exclusion/withdrawal criteria was fulfilled at some point during the study. A total of 76 patients were withdrawn before entering the treatment arms after the run-in period at the baseline visit. From the remaining 188 subjects, 97 patients were withdrawn from the pioglitazone study arm and 91 subjects were withdrawn from the acarbose arm.
Responders and nonresponders to the therapy with TZD were defined as follows: responders, >15% decrease in HbA1c levels and/or >20% decrease in fasting blood glucose levels after 12 or 26 weeks of pioglitazone treatment as compared with baseline; and nonresponders, HbA1c >11.5% or fasting blood glucose >14.0 mmol/l during the study course and <15% decrease in HbA1c levels and/or <20% decrease in fasting blood glucose levels after 12 or 26 weeks of pioglitazone treatment as compared with baseline. In addition, the occurrence of clinical symptoms of diabetes after beginning pioglitazone therapy and unexplained excretion of ketones in the urine was considered nonresponse.
Clinical laboratory
Blood samples were collected at all visits after an overnight fast to determine fasting blood glucose, HbA1c, fasting C-peptide, fasting insulin, total cholesterol, HDL cholesterol, LDL cholesterol, VLDL cholesterol, triglycerides, and routine chemistry and hematology standard laboratory parameters. Fasting blood glucose, total cholesterol, and triglyceride concentrations were measured with enzymatic assays. C-peptide and insulin plasma levels were measured using immunoassays, the HbA1c values were determined with a cation exchange/high-performance liquid chromatography method (reference values 3.46.0%), and the HDL, LDL, and VLDL cholesterol concentrations were measured using lipoprotein electrophoresis. Albuminuria was evaluated at baseline and 26 weeks after initiation of treatment. Analysis of all blood samples was performed using standard methods by a central laboratory (Dr. W. März, University of Freiburg, Germany).
Preparation of genomic DNA and screening for the PPAR-
mutations
Genomic DNA was isolated from human leukocyte nuclei isolated from whole blood by proteinase K digestion using QIAamp DNA Blood Mini Kit (Qiagen, Valencia, CA). PCR amplification of the segments with the PPAR-
mutations was performed in a total volume of 50 µl containing PCR buffer (10x Optiperform Buffer III: 500 mmol/l potassium hydroxide, pH 9.2 at 25°C; 160 mmol/l (NH4)2SO4; 0.1% Tween 20; and 2 mmol/l MgCl2), 5 µl cDNA with 1 unit of CombiPool DNA Polymerase Mix in vitro, 0.5 mmol/l dNTPs, and 1 µmol/l sense and antisense primers. The PCR conditions for the Pro12Ala mutation, the Pro467Leu mutation, and the silent polymorphism CAC478CAT were as follows: denaturation at 94°C for 3 min, followed by 34 cycles of denaturation at 94°C for 30 s, annealing at 56°C for 30 s, extension at 72°C, and a final extension at 72°C for 7 min. For the Pro115Gln and the Val290Met mutation, PCR was performed as described above except for the number of cycles (36 cycles) and the annealing temperature (54°C). The selection of all primers was based on the human genomic PPAR-
2 sequence. All primer pairs span introns. The following primer pairs were selected: Pro12Ala: forward primer: 5'-GGC TCC TAA TAG GAC AGT GCC-3', reverse primer: 5'-TAC CCT TAC ATA AAT GCC CCC-3'; Pro115Gln: forward primer: 5'-TTG CCC TGT TGC CTT TTT AG-3', reverse primer: 5-CAC CTA AAA AAG GGG TTC TGC-3'; Val290Met: forward primer: 5'-AGT CAT CCA CGT TTT CCC TG -3', reverse primer: 5'-TCC AAA ATT CTT TTG GCC AC -3'; and Pro467Leu and CAC478CAT: forward primer: 5-TTT TTT GAC TGA ACC CCC TG-3', reverse primer: 5'-TGG AAG AAG GGA AAT GTT GG-3. As a screening method for the detection of the PPAR-
2 gene mutations, denaturing gradient gel electrophoresis (DGGE) was used. For DGGE, one GC-clamped primer of each described primer pair was used. For mutation detection by DGGE, heteroduplices between wild-type and mutant alleles were generated and the optimal DGGE conditions were established for each DNA. The determined optimal DGGE conditions were as follows: Pro12Ala, Pro115Gln, Pro467Leu, 6% polyacrylamide gel containing a linear gradient of 2040%, 3-h running time; and Val290Met, 6% polyacrylamide gel containing a linear gradient of 2060%, 6 h running time. Confirmation of the detected Pro12Ala and Ala12Ala variants and the detection of the silent polymorphism CAC478CAT was performed using direct sequencing (ABI Prism Sequencer System).
Statistical analysis
Data are shown as the mean ± SD. All calculations and statistics were performed using SPSS for Windows (release 7.5; SPSS, Chicago, IL). The differences between the groups were tested using one-way ANOVA. If the P value was <0.05, the groups were compared using the appropriate test (Students t test for unpaired samples or
2 test). A P value <0.05 was considered significant. Correlations between variables were tested using Spearmans correlation test. A correlation coefficient (r) of P <0.05 was considered significant. To obtain the odds ratios for a possible specific influence of the PPAR-
genotype on BMI, C-peptide, and plasma lipid levels, multivariate logistic regression analysis was performed.
| RESULTS |
|---|
|
|
|---|
genotype. The clinical and biochemical characteristics of the 131 patients who fulfilled the complete course of pioglitazone treatment are summarized in Table 1. The frequencies of the Pro12Ala and Ala12Ala variant of the PPAR-
gene were not significantly different between the acarbose- and pioglitazone-treated groups. Furthermore, the PPAR-
allele and genotype distribution was not different between the acarbose- and pioglitazone-treated early-withdrawn patients and was not different from the distribution in the 268 patients who completed the study. From the remaining 188 subjects, 97 patients were withdrawn from the pioglitazone study arm and 91 subjects were withdrawn from the acarbose arm. Of the 97 withdrawn patients from the pioglitazone study, 88 patients were withdrawn due to HbA1c values >11.5% or fasting plasma glucose >14.0 mmol/l. The PPAR-
genotype distribution for these patients was as follows: Pro12Pro 84.9% (75 patients), Pro12Ala 11.7% (10 patients), and Ala12Ala 3.4% (3 patients). The distribution of the different PPAR-
genotypes was not different from the PPAR-
genotype distribution in the 131 patients who completed the study. In the 131 patients who completed the study, the allelic frequency of the Pro12Ala polymorphism was 12.2% (10 men, 6 women); the frequency of the Ala12Ala variant was 3.8% (four men, one woman) in the pioglitazone treatment group, and in the acarbose-treated group, the frequency was 12.8% for the Pro12Ala variant and 2.7% for the Ala12Ala variant. Therefore, we excluded a PPAR-
genotype selection bias in the pioglitazone treatment group. Only two patients (1.5%) with the Pro12Pro genotype were detected with an additional silent polymorphism C478T in exon 6. The previously described Pro115Gln (13), Pro467Leu, and Val290Met mutations (15) were not detected in any patient.
|
gene and the response to pioglitazone treatment. There were no significant differences between the responder frequency to pioglitazone treatment between the Pro12Pro and the Pro12Ala, variant. However, because only five carriers of the Ala12Ala variant were detected, the statistical power in this study is not sufficient to detect significant differences in the response rate as compared with the other PPAR-
genotypes in position 12.
|
variants in position 12 at baseline (Table 2). However, after 6 months of pioglitazone treatment, plasma triglyceride concentrations were significantly decreased in patients with the Pro12Pro and Pro12Ala genotype (P < 0.05) and HDL cholesterol concentrations were significantly increased in the four men who were Ala12Ala carriers (P < 0.05) (Table 2). At baseline, plasma C-peptide concentrations were not significantly different between the different PPAR-
genotypes (Table 1). After 6 months of pioglitazone treatment, C-peptide levels decreased in 81.5% of the patients (P < 0.05). However, a different response in function of the PPAR-
genotype was not detectable using significantly decreased C-peptide concentrations as definition of the responders.
|
genotype in position 12 and the response rate to pioglitazone treatment using the different criteria, BMI, C-peptide, total, HDL, LDL, and VLDL cholesterol levels, and serum triglyceride concentration. However, only the Ala12Ala variant was significantly related to BMI (Table 3). There was no evidence that the PPAR-
genotype in position 12 is associated with the response rate to pioglitazone treatment or plasma lipid parameters (Table 3).
|
genotype in position 12 has no significant effect on the responder rate to pioglitazone treatment (Table 4).
|
| CONCLUSIONS |
|---|
|
|
|---|
was identified as the target for TZD (18). Moreover, different genetic variants of the PPAR-
gene were shown to cause a different drug efficacy in vitro (19). Therefore, variants in the PPAR-
gene could cause differences in the efficiency of TZD therapy in the clinical application and might provide a molecular definition of either very good responders or nonresponders to TZD treatment. Therefore, we investigated whether known variants in the PPAR-
gene are associated with either resistance or high response to the therapy with the TZD pioglitazone.
The most common variant in the PPAR-
gene, the Pro12Ala variant, occurs with a frequency of
1215% (14,16,2023). Other mutations are very rare (13,15,24) and therefore unlikely to be a major cause for therapy modification by TZD. Because of its prevalence, the Pro12Ala variant is the first candidate for the possible relationship between PPAR-
genotype and response to TZD treatment. In our study of patients with type 2 diabetes, the frequency of the Pro12Ala and Ala12Ala polymorphism in the PPAR-
2 gene was similar to that reported in previous studies investigating obese patients with type 2 diabetes (14,16,2022,25). The genotype frequencies did not deviate from those expected on the basis of the Hardy-Weinberg equilibrium, although subjects with the Pro12Ala and Ala12Ala genotype were somewhat more frequent than expected based on the allele frequency (P = 0.08). A prediction of the expected phenotype of the Pro12Ala carriers is difficult, because the Pro12Ala polymorphism has been variably associated with either increased (21,26), decreased (16), or unchanged (20,22,27) BMI and improved insulin sensitivity (16,23,25). Moreover, it was recently suggested that the Pro12Ala variant may influence the susceptibility for obesity (27) and that this variant might be a genetic marker indicating the risk of obesity in children persisting into adolescence (28). In contrast to Japanese men (29) but in accordance with results from previous studies among Caucasians (21,30), the homozygous Ala12 variant was associated with increased BMI. In accordance with the normal weight or slightly overweight subjects (16,25), the overweight Ala12Ala carriers in our study had C-peptide concentrations that were indistinguishable from the other genotypes in position 12. Because only five carriers of the Ala12Ala were detected in our study, the statistical power to identify significant differences between the investigated parameters was far below 80%. Therefore, we were not able to establish a statistically significant association of this genotype with the response rate to TZD therapy, although this variant may be important in the response rate to pioglitazone therapy. More than 1,000 patients would be needed for detection of a possible association between the Ala12Ala variant and the response rate to pioglitazone treatment according to calculation based on allele frequency, effect size, and population variance, as suggested in a recent meta-analyis of the Pro12Ala variant in the PPAR-
2 gene (31).
We did not find any subject with the recently reported Pro115Gln (13), Pro467Leu, or Val290Met mutations (15). Therefore, no conclusion regarding a potential effect of these mutations on the response rate to TZD therapy can be drawn from our data. A more frequent silent polymorphism CAC478CAT in exon 6 of the PPAR-
2 gene (16) was also not associated with the therapeutic effect of pioglitazone in our study group.
To define the responders to pioglitazone therapy, we used a decrease of >20% in fasting blood glucose concentration as compared with baseline or a >15% decrease in HbA1c values after 3 or 6 months as compared with the baseline of pioglitazone therapy. The definition of response rate (>15% decrease in HbA1c values and/or >20% decrease in fasting blood glucose) was chosen according to calculations that a significant decrease in fasting blood glucose and HbA1c must be at least >2 SD of the variance of the method and at least >2 SD of the day-to-day differences in glucose and HbA1c measurements. Using these criteria, the nonresponder rate in our study population was much higher (48.1%) as compared with the reported
2030% nonresponder frequency in previous studies in patients with type 2 diabetes or impaired glucose tolerance (11,12,17,32,33). In these studies, the criteria for the nonresponse rate was defined as decreasing HbA1c by 0.51% compared with the initial glycemic status and/or decreasing fasting plasma glucose by
10%. However, using these criteria would result in an overestimation of the response rate in our study group, because a subgroup of patients with not statistically significant changes in higher initial HbA1c values and fasting plasma glucose concentrations after pioglitazone treatment would be considered responders. This aspect is supported by our finding that the nonresponder rate is significantly decreased in patients with higher initial fasting blood glucose, higher initial C-peptide levels, and insufficient glycemic control (HbA1c >9.0% at baseline). However, the therapeutic target, fasting plasma glucose concentration <6.0 mmol/l and HbA1c value <6.5% was only achieved by <30% of the treated patients, because this threshold is also dependent on other confounding factors.
To test the hypothesis that a different definition of the response rate might reveal a potential relationship between the Pro12Ala and Pro12Pro variant and the response rate to pioglitazone therapy, we analyzed our data using the following definitions of response to pioglitazone therapy: decreasing fasting plasma glucose and/or HbA1c value by >10, 20, or 30% as well as decreasing HbA1c value to <7% or decreasing fasting plasma glucose to <6.0 mmol/l to meet clinically relevant target values. Independent of the definition of response rate to pioglitazone treatment, there was no association between the Pro12Ala and the Pro12Pro variant in the PPAR-
gene and the response rate to pioglitazone therapy.
In our population, we show that the major confounding factors for TZD treatment efficacy are high fasting plasma glucose (>11.0 mmol/l) at baseline, HbA1c value >9.0%, BMI >33 kg/m2, and hyperinsulinemia (C-peptide concentrations at baseline >2.5 pmol/l) rather than the PPAR-
genotype. These findings are in accordance with the previous observation that the glucose-lowering effects of TZD are dependent on the degree and duration of the decompensation of glucose metabolism (9,11,33,34). Therefore, comparisons of response rates to TZD should be based on the same initial glycemic status. Because there is evidence for a PPAR-
gene-nutrient interaction (35), we cannot exclude an effect of individual differences in the dietary regimen on our study results. In our study population, we excluded an effect of several influence factors, such as differences in the initial glycemic control (fasting plasma glucose, HbA1c), age, BMI, and duration of diabetes on the impact of a potential association between PPAR-
genotype and the response rate to pioglitazone treatment by introducing matched subgroups for these factors. However, we could not find evidence of an association between the pioglitazone therapy response and the PPAR-
genotype in any of these matched groups. Furthermore, the multivariate regression analysis could not identify the Pro12Ala or the Ala12Ala genotype as independent modulating factors of the pioglitazone therapy effect. Because C-peptide concentrations were significantly lower in 81.5% of the patients after pioglitazone treatment, a significant decrease in fasting hyperinsulinemia, as determined by fasting C-peptide levels, might be a better parameter for testing the therapy efficacy of TZD in patients with type 2 diabetes. Moreover, it was recently shown that pioglitazone treatment has beneficial effects on serum lipid concentrations in patients with type 2 diabetes (36,37). In accordance with these studies, we found decreased serum triglyceride concentrations and increased HDL cholesterol concentrations, independent of the genotype of the PPAR-
gene in position 12, suggesting that these variants do not modify the beneficial effects of pioglitazone treatment on the lipid profile in our patients.
In conclusion, the Pro12Ala and the Pro12Pro variants in the PPAR-
gene are not associated with the response rate to pioglitazone treatment in patients with type 2 diabetes. Further studies are necessary to investigate the potential importance of the Ala12Ala variant in the response to thiazolidinedione treatment. We identified initial fasting plasma glucose >11.0 mmol/l, HbA1c value >9.0%, BMI >33 kg/m2, and fasting C-peptide concentrations at baseline >2.5 pmol/l, but not the PPAR-
genotype in position 12, as predominant confounding factors for the responder frequency to pioglitazone treatment.
| Acknowledgments |
|---|
| Footnotes |
|---|
Received for publication 22 May 2002 and accepted in revised form 15 November 2002.
A table elsewhere in this issue shows conventional and Système International (SI) units and conversion factors for many substances.
| References |
|---|
|
|
|---|
2, a regulator of adipocyte differentiation. N Engl J Med 339: 953959, 1998
(PPAR-
) gene in diabetic Caucasians: identification of a Pro12Ala PPAR
2 missense mutation. Biochem Biophys Res Commun 240: 270276, 1997
associated with severe insulin resistance, diabetes mellitus and hypertension. Nature 402: 880883, 1999[Medline]
2 associated with decreased receptor activity, lower body mass index and improved insulin sensitivity. Nat Genet 20: 284287, 1998[Medline]
2 in 3T3L1 cells. Biochem Biophys Res Commun 264: 9399, 1999[Medline]
2 gene in lean and obese subjects. Eur J Endocrinol 141: 9092, 1999[Abstract]
2 gene with obesity in two Caucasian populations. Diabetes 47: 18061808, 1998[Medline]
2 (PPAR-
2): divergent modulating effects on body mass index in obese and lean Caucasian men. Diabetologia 42: 892895, 1999[Medline]
2 amino acid polymorphism Pro12Ala is prevalent in offspring of type II diabetic patients and is associated to increased insulin sensitivity in a subgroup of obese subjects. Diabetologia 42: 758762, 1999[Medline]
and obesity. Diabetes Care 23: 126127, 2000
2 gene (Pro12Ala) affects metabolic responses to weight loss and subsequent weight gain. Diabetes 50: 21722176, 2001
2 gene on adiposity, fat distribution, and insulin sensitivity in Japanese men. Biochem Biophys Res Commun 251: 195198, 1998[Medline]
gene are associated with severe overweight among obese women. J Clin Endocrinol Metab 84: 37083712, 1999
locus. Diabetes 50: 686689, 2001This article has been cited by other articles:
![]() |
J. C. Florez, K. A. Jablonski, M. W. Sun, N. Bayley, S. E. Kahn, H. Shamoon, R. F. Hamman, W. C. Knowler, D. M. Nathan, D. Altshuler, et al. Effects of the Type 2 Diabetes-Associated PPARG P12A Polymorphism on Progression to Diabetes and Response to Troglitazone J. Clin. Endocrinol. Metab., April 1, 2007; 92(4): 1502 - 1509. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. K. Wolford, K. A. Yeatts, S. K. Dhanjal, M. H. Black, A. H. Xiang, T. A. Buchanan, and R. M. Watanabe Sequence Variation in PPARG May Underlie Differential Response to Troglitazone Diabetes, November 1, 2005; 54(11): 3319 - 3325. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Snitker, R. M. Watanabe, I. Ani, A. H. Xiang, A. Marroquin, C. Ochoa, J. Goico, A. R. Shuldiner, and T. A. Buchanan Changes in Insulin Sensitivity in Response to Troglitazone Do Not Differ Between Subjects With and Without the Common, Functional Pro12Ala Peroxisome Proliferator-Activated Receptor-{gamma}2 Gene Variant: Results from the Troglitazone in Prevention of Diabetes (TRIPOD) study Diabetes Care, June 1, 2004; 27(6): 1365 - 1368. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Sesti, M. A. Marini, M. Cardellini, A. Sciacqua, S. Frontoni, F. Andreozzi, C. Irace, D. Lauro, A. Gnasso, M. Federici, et al. The Arg972 Variant in Insulin Receptor Substrate-1 Is Associated With an Increased Risk of Secondary Failure to Sulfonylurea in Patients With Type 2 Diabetes Diabetes Care, June 1, 2004; 27(6): 1394 - 1398. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Diabetes | Diabetes Care | Clinical Diabetes | Diabetes Spectrum |