The Peroxisome Poliferator–Activated Receptor-{gamma}2 Pro12Ala Variant

Association With Type 2 Diabetes and Trait Differences

Julie A. Douglas1, Michael R. Erdos2, Richard M. Watanabe3, Andi Braun4, Cristy L. Johnston4, Paul Oeth4, Karen L. Mohlke2, Timo T. Valle5, Christian Ehnholm5, Thomas A. Buchanan6, Richard N. Bergman7, Francis S. Collins2, Michael Boehnke1, and Jaakko Tuomilehto5,8

1 Department of Biostatistics, School of Public Health, University of Michigan, Ann Arbor, Michigan
2 Genetics and Molecular Biology Branch, National Human Genome Research Institute, Bethesda, Maryland
3 Division of Biostatistics, Department of Preventative Medicine, Keck School of Medicine, University of Southern California, Los Angeles
4 Sequenom Inc., San Diego, California
5 Department of Epidemiology and Health Promotion, Diabetes and Genetic Epidemiology Unit, and the Department of Biochemistry, National Public Health Institute, Helsinki, Finland
6 Department of Medicine and the
7 Department of Physiology and Biophysics, Keck School of Medicine, University of Southern California, Los Angeles, California
8 Department of Public Health, University of Helsinki, Helsinki, Finland


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 RESEARCH DESIGN AND METHODS
 REFERENCES
 
Recent studies have identified a common proline-to-alanine substitution (Pro12Ala) in the peroxisome proliferator–activated receptor-{gamma}2 (PPAR-{gamma}2), a nuclear receptor that regulates adipocyte differentiation and possibly insulin sensitivity. The Pro12Ala variant has been associated in some studies with diabetes-related traits and/or protection against type 2 diabetes. We examined this variant in 935 Finnish subjects, including 522 subjects with type 2 diabetes, 193 nondiabetic spouses, and 220 elderly nondiabetic control subjects. The frequency of the Pro12Ala variant was significantly lower in diabetic subjects than in nondiabetic subjects (0.15 vs. 0.21; P = 0.001). We also compared diabetes-related traits between subjects with and without the Pro12Ala variant within subgroups. Among diabetic subjects, the variant was associated with greater weight gain after age 20 years (P = 0.023) and lower triglyceride levels (P = 0.033). Diastolic blood pressure was higher in grossly obese (BMI >40 kg/m2) diabetic subjects with the variant. In nondiabetic spouses, the variant was associated with higher fasting insulin (P = 0.033), systolic blood pressure (P = 0.021), and diastolic blood pressure (P = 0.045). These findings support a role for the PPAR-{gamma}2 Pro12Ala variant in the etiology of type 2 diabetes and the insulin resistance syndrome.

AIRG, acute insulin response to glucose; dBP, diastolic blood pressure; DI, disposition index; FUSION, Finland–United States Investigation of Non–Insulin-Dependent Diabetes Mellitus Genetics; MALDI-TOF, matrix-assisted laser desorption/ionization time-of-flight; OGTT, oral glucose tolerance test; PPAR, peroxisome proliferator–activated receptor; PROBE, primer oligo base extension reaction; sBP, systolic blood pressure; SI, insulin sensitivity


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 RESEARCH DESIGN AND METHODS
 REFERENCES
 
Peroxisome proliferator–activated receptors(PPARs) are members of the nuclear hormone receptor family of transcription factors and are involved in adipocyte differentiation and gene expression. They are also believed to play an important role in type 2 diabetes and diabetes-related traits, including insulin sensitivity and lipid and energy metabolism (1). In fact, studies have shown that ligands for PPAR-{gamma}, including both endogenous ones and those that are synthetic (e.g., thiazolidinedione drugs), stimulate adipogenesis and increase insulin action (2). A common proline-to-alanine substitution at codon 12 (Pro12Ala) of exon B has been inconsistently associated with protection against type 2 diabetes and diabetes-related traits (314). These findings encouraged us to investigate the role of PPAR-{gamma}2 in our sample from the Finnish population. The objective of our study was to examine whether the Pro12Ala variant was associated with type 2 diabetes and to examine the relationship between the Pro12Ala variant and diabetes-related traits among subgroups of diabetic and nondiabetic subjects.

The mean and standard deviation of selected trait values are given in Table 1 by subgroup. The entire sample of 935 subjects consisted of 636 Pro/Pro subjects, 271 Pro/Ala subjects, and 28 Ala/Ala subjects. The allele frequency of the Pro12Ala variant in the PPAR-{gamma}2 gene was 0.15 among diabetic subjects, 0.19 among spousal control subjects, and 0.22 among elderly control subjects (Table 2), which mirrors the continuum of diabetes susceptibility across these subgroups. The frequency of the variant was significantly lower in diabetic subjects than in elderly control subjects ({chi}2 = 11.72, df = 1, P < 0.0007) and marginally lower than in spousal control subjects (P = 0.069). Comparison with combined spousal and elderly control subjects gave a significant association result ({chi}2 = 10.60, df = 1, P = 0.001). A second independent sample of 263 Finnish diabetic subjects in our study subsequently confirmed the variant frequency of 0.15 in the original 522 diabetic subjects (data not shown). The observed genotype data were consistent with Hardy-Weinberg equilibrium.


View this table:
[in this window]
[in a new window]
 
TABLE 1 Characteristics of the subjects by clinical subgroup

 

View this table:
[in this window]
[in a new window]
 
TABLE 2 Frequency of the PPAR-{gamma}2 Pro12Ala variant by clinical subgroup

 
Results for the quantitative traits were less compelling. Genotype-specific means for all traits for diabetic subjects, elderly control subjects, and spousal control subjects, respectively, are available in an online appendix (ATables 13) at www.diabetes.org/diabetes/appendix.asp. Table 3 shows the significant trait differences between subjects with and without the Pro12Ala variant by subgroup. In diabetic subjects, the presence of the variant was associated with greater weight change after 20 years of age (22.2 ± 14.0 vs. 19.5 ± 13.0 kg) and lower serum triglyceride levels (2.29 ± 1.65 vs. 2.68 ± 2.21 mmol/l). Both results were significant after adjustment for sex, age, and (for triglyceride levels) BMI (P = 0.023 and 0.033, respectively). There was a significant interaction (P = 0.038) between the variant and BMI for diastolic blood pressure (dBP); the variant was associated with higher dBP only among grossly obese diabetic subjects (Table 4). A similar trend was also observed for systolic blood pressure (sBP), although the interaction was not statistically significant (P = 0.299).


View this table:
[in this window]
[in a new window]
 
TABLe 3 Significant results by clinical subgroup and presence/absence of the PPAR-{gamma}2 Pro12Ala variant

 

View this table:
[in this window]
[in a new window]
 
TABLE 4 dBP and sBP in diabetic subjects by presence/absence of the PPAR-{gamma}2 Pro12Ala variant and BMI

 
Among elderly control subjects, only maximum lifetime weight was significantly associated with the Pro12Ala variant (P = 0.049), and it was no longer significant after adjustment for sex (P = 0.191) (Table 3). Among nondiabetic spousal control subjects, the variant was significantly associated with higher sBP (151.4 ± 25.5 vs. 142.5 ± 18.9 mm Hg) (P = 0.014). When both sBP and dBP and fasting serum insulin were adjusted for sex, age, and BMI, differences between spousal control subjects with and without the Pro12Ala variant remained and/or became significant (P = 0.021, 0.045, and 0.033, respectively) (Table 3). All three traits were significantly higher among subjects with the variant.

In our analysis, we found a significantly lower frequency of the Pro12Ala variant of the PPAR-{gamma}2 gene in diabetic subjects than in nondiabetic subjects. The directionality of these highly significant (P = 0.001) findings is consistent with results from the studies of Deeb et al. (3), Mancini et al. (4), and Altshuler et al. (14), although the difference in allele frequencies in the second study failed to reach statistical significance. Coupled with these studies and the biological importance of PPAR{gamma}, these findings suggest a link between the Pro12Ala variant of the PPAR-{gamma}2 gene and the pathogenesis of type 2 diabetes. The increased frequency of the variant in nondiabetic subjects would seem to suggest that the Pro12Ala variant confers some protective effect against diabetes.

Despite the increased frequency of the Pro12Ala variant among elderly control subjects, we failed to find any significant trait associations within this subgroup. Instead, we observed weak but significant associations between the Pro12Ala variant and traits characteristic of the insulin resistance syndrome in both diabetic and nondiabetic subjects. For example, greater weight gain was associated with the Pro12Ala variant in diabetic subjects, whereas higher fasting insulin, sBP, and dBP were associated with the variant in nondiabetic spouses. It should be emphasized that, in contrast to the spousal control subjects, the elderly control subjects represent a quite distinct subgroup of nondiabetic subjects who are unlikely to ever develop type 2 diabetes. As such, they are unlikely to carry the cluster of susceptibility genes that may interact with variants in PPAR-{gamma}2 to result in the insulin resistance syndrome phenotype. The spousal control subjects are somewhat younger and remain at risk for developing type 2 diabetes during their lifetime.

Alterations in functional characteristics of the PPAR-{gamma}2 gene induced by the Ala isoform may be partly responsible for the manifestation of some characteristics of the insulin resistance syndrome. Deeb et al. (3) identified lowered transactivation capacity and reduced stimulation of PPAR-{gamma} target genes as a potential molecular mechanism underlying the association of the Pro12Ala variant with lower BMI and increased insulin sensitivity, a hypothesis consistent with their observations in Finnish subjects. Although this hypothesis may appear to be at odds with (or at least not supported by) our trait findings, several points should be clarified. First, the middle-aged subjects in the study by Deeb et al. (3) were much younger and leaner than our nondiabetic spouses and elderly control subjects. Second, although the elderly subjects from both studies were better matched, we could not parallel their genotype-based analysis because of insufficient numbers of Pro12Ala homozygotes. If the Pro/Ala and Ala/Ala subjects from their study had been pooled, it is unlikely that they would have observed significant trait differences because trends within their elderly subjects were inconsistent (e.g., fasting insulin was highest for Pro/Ala heterozygotes).

Consistent with at least one report of a differential effect of the PPAR-{gamma}2 Pro12Ala variant in the lean and obese states (12), we also found an interaction between BMI and the variant for dBP in the diabetic subjects. Among severely obese subjects, those with the Pro12Ala variant had substantially higher blood pressure. Higher values for sBP and dBP were also associated with the variant among nondiabetic spousal control subjects, though there was no evidence for an interaction between BMI and the Pro12Ala variant. These associations are of interest, given the recent report by Barroso et al. (13) of three type 2 diabetic subjects with early-onset hypertension and polymorphisms in the PPAR-{gamma}2 gene, suggesting that this receptor is important in both blood pressure and glucose homeostasis.

In summary, we found that the Pro12Ala variant of the PPAR-{gamma}2 gene was associated with protection against type 2 diabetes in Finnish subjects, a finding consistent with several reports in the literature (3,4,14). Because we only screened for this particular variant, we cannot exclude the role of other PPAR-{gamma}2 variants or variants in nearby genes, possibly in linkage disequilibrium with the Pro12Ala variant. Further studies, including functional analyses, will be required to fully understand the role of this gene in type 2 diabetes. Our data suggest that the PPAR-{gamma}2 Pro12Ala variant has variable effects among subgroups of individuals with different levels of diabetes risk.


    RESEARCH DESIGN AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 RESEARCH DESIGN AND METHODS
 REFERENCES
 
The Finland–United States Investigation of Non–Insulin-Dependent Diabetes Mellitus Genetics (FUSION) Study is an international collaborative effort to map and clone genes predisposing to type 2 diabetes and related traits in Finnish subjects. The FUSION study design and family material have been described previously (15). For the present investigation, our sample included 522 unrelated subjects with type 2 diabetes, 193 nondiabetic spouses of a diabetic subject or his/her affected sibling, and 220 unrelated elderly nondiabetic control subjects. Diabetes was diagnosed by World Health Organization (16) criteria. Spouses had a single normal oral glucose tolerance test (OGTT). Elderly control subjects had normal glucose tolerance at ages 65 and 70 years.

A total of 14 traits were analyzed on all subjects: BMI, waist circumference, waist-to-hip ratio, current weight, maximum lifetime weight, fasting plasma glucose, fasting serum insulin, total cholesterol, HDL cholesterol, HDL ratio (HDL cholesterol/total cholesterol), LDL cholesterol, triglycerides, sBP, and dBP. Values for sBP and dBP were each determined as the mean of two measurements. Seven additional traits were ascertained on diabetic subjects: weight at 20 years of age, change in weight after 20 years of age, maximum lifetime weight change after 20 years of age, age at diagnosis of diabetes, diabetes duration, age at which insulin treatment started (if applicable), and fasting plasma C-peptide concentrations. In addition, glucose and insulin concentrations 2 h after OGTT were analyzed in nondiabetic subjects, whereas the insulin sensitivity index (SI), the glucose effectiveness index, the acute insulin response to glucose (AIRG), and the disposition index (DI) were analyzed (DI = SI x AIRG) only in the nondiabetic spouses; the latter analyses used tolbutamide-modified frequently sampled intravenous glucose tolerance tests and minimal model analysis (17). Glucose, insulin, C-peptide, and lipid concentrations were assayed using standard methods (15).

Genotyping by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry.
The PPAR-{gamma}2 Pro12Ala variant was analyzed by matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) mass spectrometry. A 69-bp fragment containing the Pro12Ala variant site was amplified by polymerase chain reaction from 20 ng genomic DNA using 25 pmol forward primer 5'-GCTGTTATGGGTGAAACTCTG, 2 pmol of a universal sequence-tailed reverse primer 5'-AGCGGATAACAATTTCACACAGGCAGTGTATCAGTGAAGGAATCG, and 10 pmol of a biotinylated universal primer 5'-biotin-AGCGGATAACAATTTCACACAGG under standard reaction conditions (Fig. 1A). After 15 min of denaturation at 95°C, 55 cycles (5 s at 95°C, 20 s at 53°C, and 30 s at 72°C) were performed. To recover the single-stranded DNA template, the product was immobilized on streptavidin-coated magnetic beads (Dynal, Great Neck, NY), washed with 10 mmol/l Tris-HCl at pH 8.0, denatured in 50 µl 0.1 mol/l NaOH, and washed again with 10 mmol/l Tris-HCl.



View larger version (18K):
[in this window]
[in a new window]
 
FIG. 1. Genotype analysis by MALDI-TOF spectrometry. A: PROBE reaction. The region containing the PPAR{gamma}2 Pro12Ala (CCA->GCA) variant is amplified with a biotinylated primer to enable purification of the single-stranded template. Next, the PROBE primer anneals to the template and is extended. When the single nucleotide polymorphism (SNP) is C, the probe is extended by one nucleotide, dideoxy-CTP. When the SNP is G, the probe is extended by two nucleotides, deoxy-GTP, and dideoxy-CTP. B: Mass spectrometry profiles of primer extension products. Peaks at 6,989.6 and 7,318.8 Da correspond to the mass of the probe primer extended by one or two nucleotides, respectively. Genotypes of the spectra are 1) CC, 2) CG, and 3) GG. The mass of the unextended PROBE primer is indicated at 6,716.4 Da, but in these examples, none is detected.

 
The primer oligo base extension reaction (PROBE) was performed by the addition of 20 pmol extension primer 5'-TCTGGGAGATTCTCCTATTGAC under conditions similar to those previously described (18). The extension reaction products were applied to a SpectroChip (Sequenom, San Diego, CA) prespotted with a matrix of 3-hydroxypicolinic acid using a Spectrojet piezoelectric nanoliter dispensing system (19). A modified Bruker Biflex III MALDI-TOF mass spectrometer (DNA MassArray; Sequenom) was used to determine genotypes by the appearance of peaks corresponding to the expected extension product masses (Fig. 1B).

Statistical analyses.
Associations of the Pro12Ala variant of the PPAR-{gamma}2 gene between diabetic subjects and both nondiabetic spouses and elderly control subjects were examined by {chi}2 tests of independence. Trait differences within diabetic, elderly control, or spousal control subgroups were examined by analysis of variance. Initially, we tested whether trait means differed significantly among subjects with the Pro/Pro, Pro/Ala, and Ala/Ala genotypes. Due to the small number of individuals with the Ala/Ala genotype, we subsequently tested whether the trait means differed between subjects with and without the Pro12Ala variant (Pro/Ala and Ala/Ala versus Pro/Pro). All analyses were performed with and without adjustment for covariates, including sex, age, and BMI. Preselected interactions between the variant and sex or BMI were also tested. Standard regression diagnostics were computed to examine the adequacy of model assumptions, and traits were transformed to approximate normality when necessary. P values <0.05 were considered statistically significant. No adjustments for multiple comparisons were made. We excluded from the analyses any subject who, on the day of their examinations, took medications that could influence the trait of interest. We also excluded subjects whose diabetic status was uncertain and those with a first-degree relative with type 1 diabetes.


    ACKNOWLEDGMENTS
 
The FUSION study is made possible by intramural funds from the National Human Genome Research Institute (Project number OH95-C-N030), by grants from the Finnish Academy (38387 and 46558), and by National Institiutes of Health grants HG00040 (J.A.D.), HG00376 (M.B.), DK09525 (R.M.W.), DK27619, and DK29867 (R.N.B.). Currently, J.A.D. is supported by a University of Michigan Rackham Predoctoral Fellowship, and R.M.W. is supported by a Career Development Award from the American Diabetes Association.

We wish to thank all of the subjects for their invaluable contribution to the FUSION study. We also gratefully acknowledge Peter Chines for his exceptional work in pre-paring the data. Family studies were approved by institutional review boards at the National Institutes of Health (assurance number SPA S-5737-05) and at the National Public Health Institute in Helsinki, Finland.


    FOOTNOTES
 
Address correspondence and reprint requests to Michael Boehnke, University of Michigan, Department of Biostatistics, 1420 Washington Heights, Ann Arbor, Michigan 48109-2029. E-mail: boehnke{at}umich.edu.

Received for publication 12 June 2000 and accepted in revised form 29 December 2000.

J.A.D. and M.R.E. contributed equally to this work.

Additional information can be found in an online appendix at www.diabetes.org/diabetes/appendix.asp.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 RESEARCH DESIGN AND METHODS
 REFERENCES
 

  1. Auwerx J: PPAR{gamma}, the ultimate thrifty gene. Diabetologia 42:1033–1049, 1999[Medline]
  2. Spiegelman BM: PPAR-{gamma}: adipogenic regulator and thiazolidinedione receptor (Review). Diabetes 47:507–514, 1998[Abstract]
  3. Deeb SS, Fajas L, Nemoto M, Pihlajamaki J, Mykkanen L, Kuusisto J, Laakso M, Fujimoto W, Auwerx J: A Pro12Ala substitution in PPAR{gamma}2 associated with decreased receptor activity, lower body mass index and improved insulin sensitivity. Nat Genet 20:284–287, 1998[Medline]
  4. Mancini FP, Vaccaro O, Sabatino L, Tufano A, Rivellese AA, Riccardi G, Colantuoni V: Pro12Ala substitution in the peroxisome proliferator-activated receptor-{gamma}2 is not associated with type 2 diabetes. Diabetes 48:1466–1468, 1999[Abstract]
  5. Ringel J, Engeli S, Distler A, Sharma AM: Pro12Ala missense mutation of the peroxisome proliferator activated receptor {gamma} and diabetes mellitus. Biochem Biophys Res Commun 254:450–453, 1999[Medline]
  6. Clement K, Hercberg S, Passinge B, Galan P, Varroud-Vial M, Shuldiner AR, Beamer BA, Charpentier G, Guy-Grand B, Froguel P, Vaisse C: The Pro115Gln and Pro12Ala PPAR {gamma} gene mutations in obesity and type 2 diabetes. Int J Obes 24:391–393, 2000
  7. Meirhaeghe A, Fajas L, Helbecque N, Cottel D, Auwerx J, Deeb SS, Amouyel P: Impact of the peroxisome proliferator activated receptor {gamma}2 Pro12Ala polymorphism on adiposity, lipids, and non-insulin-dependent diabetes mellitus. Int J Obes 24:195–199, 2000
  8. Ek J, Urhammer SA, Sorensen TI, Andersen T, Auwerx J, Pedersen O: Homozygosity of the Pro12Ala variant of the peroxisome proliferation-activated receptor-{gamma}2 (PPAR-{gamma}2): divergent modulating effects on body mass index in obese and lean Caucasian men. Diabetologia 42:892–895, 1999[Medline]
  9. Beamer BA, Yen CJ, Andersen RE, Muller D, Elahi D, Cheskin LJ, Andres R, Roth J, Shuldiner AR: Association of the Pro12Ala variant in the peroxisome proliferator–activated receptor-{gamma}2 gene with obesity in two Caucasian populations. Diabetes 47:1806–1808, 1998[Free Full Text]
  10. Koch M, Rett K, Maerker E, Volk A, Haist K, Deninger M, Renn W, Haring HU: The PPAR{gamma}2 amino acid polymorphism Pro 12 Ala is prevalent in offspring of type II diabetic patients and is associated to increased insulin sensitivity in a subgroup of obese subjects. Diabetologia 42:758–762, 1999[Medline]
  11. Cole SA, Mitchell BD, Hseuh W, Pineda P, Beamer BA, Shuldiner AR, Comuzzie AG, Blangero J, Hixson JE: The Pro12Ala variant in peroxisome proliferator-activated receptor-{gamma}2 (PPAR-{gamma}2) is associated with measures of obesity in Mexican Americans. Int J Obes 24:522–524, 2000
  12. Ristow M, Muller-Wieland D, Pfeiffer A, Krone W, Kahn CR: Obesity associated with a mutation in a genetic regulator of adipocyte differentiation. N Engl J Med 339:953–959, 1998[Abstract/Free Full Text]
  13. Barroso I, Gurnell M, Crowley VE, Agostini M, Schwabe JW, Soos MA, Maslen GL, Williams TD, Lewis H, Schafer AJ, Chatterjee VK, O’Rahilly S: Dominant negative mutations in human PPAR{gamma} associated with severe insulin resistance, diabetes mellitus and hypertension. Nature 402:880–883, 1999[Medline]
  14. Altshuler D, Hirschhorn JN, Klannemark M, Lindgren CM, Vohl M, Nemesh J, Lane CR, Schaffner SF, Bolk S, Brewer C, Tuomi T, Gaudet D, Hudson TJ, Daly M, Groop L, Lander ES: The common PPAR{gamma} Pro12Ala polymorphism is associated with decreased risk of type 2 diabetes. Nat Genet 26:76–80, 2000[Medline]
  15. Valle T, Tuomilehto J, Bergman RN, Ghosh S, Hauser ER, Eriksson J, Nylund SJ, Kohtamaki K, Toivanen L, Vidgren G, Tuomilehto-Wolf E, Ehnholm C, Blaschak J, Langefeld CD, Watanabe RM, Magnuson V, Ally DS, Hagopian WA, Ross E, Buchanan TA, Collins F, Boehnke M: Mapping genes for NIDDM: design of the Finland-United States Investigation of NIDDM (FUSION) Genetics Study. Diabetes Care 21:949–958, 1998[Abstract]
  16. World Health Organization: Diabetes Mellitus: Report of a WHO Study Group . Geneva, World Health Org., 1985 (Tech. Rep. Ser. no. 727)
  17. Bergman RN: Lilly Lecture 1989: Toward physiological understanding of glucose tolerance: minimal-model approach. Diabetes 38:1512–1527, 1989[Abstract]
  18. Braun A, Little DP, Köster H: Detecting CFTR gene mutations by using primer oligo base extension and mass spectrometry. Clin Chem 43:1151–1158, 1997[Abstract/Free Full Text]
  19. Little DP, Cornish TJ, O’Donnell MJ, Braun A, Cotter RJ, Köster H: MALDI on a chip: analysis of arrays of low-femtomole to subfemtomole quantities of synthetic oligonucleotides and DNA diagnostic products dispensed by a piezoelectric pipette. Anal Chem 69:4540–4546, 1997