Impact of Two Common Polymorphisms in the PPARγ Gene on Glucose Tolerance and Plasma Insulin Profiles in Monozygotic and Dizygotic Twins
Thrifty Genotype, Thrifty Phenotype, or Both?
Thrifty Genotype, Thrifty Phenotype, or Both?
The Pro12Ala polymorphism in the PPARγ2 gene has been associated with reduced risk of type 2 diabetes and insulin resistance. Recently, an association between dizygotic twinning and PPARγ gene polymorphisms has been proposed. We investigated the phenotypic appearance of the two polymorphisms (Pro12Ala and exon 6 C→T) in PPARγ among elderly twins (207 monozygotic [MZ] and 342 dizygotic [DZ]) and evaluated whether they could explain previously reported differences in plasma glucose and insulin profiles among MZ and DZ twins. We demonstrated a significant impact of the Pro12Ala polymorphism on glucose tolerance, diabetic status, homeostasis model assessment for insulin resistance, and plasma insulin profiles in twins. No impact of the silent exon 6 polymorphism on glucose homeostasis or plasma insulin profiles was found. Independent of the polymorphisms, we observed a significant impact of zygosity status per se on the plasma insulin profile after oral glucose ingestion, with the MZ twins being more hyperinsulinemic, indicating insulin resistance, than the DZ twins. Nonsignificantly higher glucose concentrations were observed among MZ compared with DZ twins. We demonstrated an association between the Ala allele and reduced risk of diabetes and insulin resistance in twins. However, the differences in metabolic profiles among MZ and DZ twins were not explained by differences in frequencies of the genetic variants and may be due to intrauterine environmental factors operating in twins independent of genotype. Accordingly, our study simultaneously supports a role for both the intrauterine environment (thrifty phenotype) and for genetics (thrifty genotype) in the etiology of insulin resistance and perhaps glucose intolerance in twins.
The association between birth anthropometry and glucose intolerance and/or type 2 diabetes has consistently been demonstrated in several epidemiological studies among singletons (1). However, it has been proposed that the association may be due to a genotype leading to both low birth weight and subsequent development of diabetes in adult life (2). Among a homogenous population-based twin sample, we have recently reported elevated levels of plasma glucose and insulin after oral glucose ingestion among monozygotic (MZ) compared with dizygotic (DZ) twins (3), supporting a role of the intrauterine environment on glucose homeostasis given the more adverse intrauterine conditions characterizing MZ pregnancies (4). Although other studies have not replicated these findings in younger twin populations (5,6), we have recently demonstrated a lower in vivo insulin action, measured by the euglycemic hyperinsulinemic clamp technique, among elderly MZ compared with DZ twins (7).
A recent study reported a higher frequency of the silent C→T substitution in exon 6 of PPARγ among DZ relative to MZ twins of different ethnic origins (8). The authors concluded that this silent mutation may be linked to DZ twinning and involved in intrauterine survival of DZ twins. Although the frequency of the Pro12Ala polymorphism was similar among MZ and DZ twins, some evidence of transmission disequilibrium was reported for the combination of the two different polymorphisms from the parents of the DZ twins. Another recent study found no linkage between PPARγ and dizygotic twinning. The study, however, did not report allele or genotype frequencies among MZ and DZ twins (9).
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. The common Pro12Ala polymorphism of PPARγ2, initially identified by Yen et al. (10), has inconsistently been associated with protection against type 2 diabetes and associated diabetic traits, including insulin resistance (11–18). One study has reported an association between the C→T substitution in exon 6 and a reduced risk of cardiovascular disease (19). No studies have directly addressed the phenotypic appearance according to plasma glucose and insulin profiles of the silent C→T substitution in exon 6.
Given the important role of PPARγ in glucose homeostasis, it is tempting to speculate that differences in PPARγ polymorphisms, and not the intrauterine environment, may explain our previous findings of increased plasma glucose and insulin profiles after oral glucose ingestion in MZ compared with DZ twins. Thus, in the present study we report frequencies of the Pro12Ala and C→T polymorphism among a homogeneous population-based sample of Danish Caucasian MZ and DZ twins. Furthermore, we report the phenotypic appearance according to plasma glucose and insulin profiles of the polymorphisms, and in particular, we have determined whether any differences in genotype frequencies can explain our findings of different metabolic profiles among MZ and DZ twins (3).
The twin sample was population-based and ascertained through the Danish twin register as described previously (20). A total of 303 (125 MZ and 178 DZ) twin pairs participated in the clinical examination, which involved a standard oral glucose tolerance test (OGTT) and anthropometric measurements. Zygosity classification was established using the similarity method (1). Genotyping of the twins revealed a misclassification of three pairs (genotypic discordant MZ pairs), who were excluded from further analysis because of uncertain zygosity classification. DNA was not available from 46 twins, leaving a total of 554 twins in the present study; 210 MZ (90 pairs and 30 single twins) and 344 DZ (168 pairs; 8 single twins).
The study was approved by the regional ethical committees and was conducted according to the principles of the Helsinki Declaration.
RESEARCH DESIGN AND METHODS
The clinical examination, including an OGTT and anthropometric measures, was performed as previously described (1,20). Plasma glucose and insulin concentrations were measured and analyzed as described (1,20). Incremental glucose and insulin areas under the curves (AUCs) were calculated using the trapezoidal method. Measures of insulin resistance were calculated from fasting plasma glucose and insulin concentrations using homestasis model assessment (HOMA) (21). Type 2 diabetes and impaired glucose tolerance were defined according to the 1998 World Health Organization criteria.
The genetic analyses were performed on genomic DNA isolated from human leukocytes. The Pro12Ala variant was genotyped as previously described (22). The exon 6 C→T silent histidine variant was genotyped with primers 5′-TGT GAA GCC CAT TGA AGA CA-3′ and 5′-GAG CGG GTG AAG ACT CAT GT-3′ (derived from EMBL no. AB005526) and carried out in a volume of 25 μl containing 100 ng genomic DNA, 1× PCR buffer, 0.2 μmol/l of each primer, 0.2 mmol/l dNTP, 1.4 units AmpliTaq Gold polymerase (Applied Biosystems), and 3 mmol/l MgCl2. The cycling program was a denaturation step at 95°C for 15 min followed by 40 cycles of 94°C for 30 s, 53°C for 30 s, and 72°C for 60 s, with a final elongation step at 72°C for 9 min using a GeneAmp 9600 thermal cycler (Perkin Elmer). The PCR product was digested overnight at 37°C with two units of NlaIII (New England Biolabs), and fragments were separated on a 4% agarose gel. The concordance rates between the PCR–restriction fragment–length polymorphism assay used for genotyping and sequencing or repeat genotyping was 100%.
Fischer’s exact test was applied to compare allele frequencies and genotype distributions in subgroups of twins. Only one random twin from each twin pair was included in the analysis. The test for Hardy-Weinberg equilibrium was performed as described previously (23).
Since MZ twins share their entire genome and DZ twins share half of their genome, one may question the supposition that the twins are independent observations. Accordingly, we have performed an ANOVA for each phenotypic parameter in which we adjusted for the inter–twin-pair relationship. The full ANOVA model includes a random effect term for twin-pair membership and fixed effect terms for PPARγ and zygosity and their interaction, i.e., different effects of PPARγ in MZ and DZ twins. The interaction effect was not significant, and the tests of zygosity and genotype were made in the model without interaction terms. In addition, it was investigated by paired comparisons whether the Pro12Ala and the Ala12Ala genotypes could be pooled. The comparisons showed no significant differences, and therefore these two genotypes were combined. The statistical analyses were performed with PROC MIXED of the SAS/STAT system (Version 8.2, SAS Institute ). Initial descriptive data analysis showed that all variables except waist-to-hip ratio exhibited markedly right-skewed distribution. Therefore, a transformation was considered necessary for application of the mixed ANOVA, and a logarithmic transformation was found to give a much better approximation to a symmetric Gaussian distribution. The means are reported as the antilog of the means of the log-transformed values, and SE was calculated using the standard error of the log-transformed mean. Comparisons of metabolic and anthropometric phenotype data within the genotypically discordant DZ twin pairs were performed using the parametric paired t test for normally distributed data and the Wilcoxon’s signed-rank test for skewed data. A P value ≤0.05 was considered significant. Due to the fact that a meta-analysis of all published papers strongly supports that the Pro12Ala polymorphism (22) represents a significant type 2 diabetes and insulin resistance gene, we have found it reasonable not to correct for multiple comparisons when analyzing the effect of this polymorphism on glucose homeostasis.
Genotype and allele frequencies according to zygosity.
The distribution of both polymorphisms was in Hardy-Weinberg equilibrium among both MZ and DZ twins. The allele and genotype frequencies were similar in both MZ and DZ twins for each polymorphism (Table 1). However, the genotype distribution of the exon 6 polymorphism was nearly significantly different between MZ and DZ twins (P = 0.06) (Table 1). The degree of allele sharing among the DZ twins (phenotype, DZ twin) was 0.91 (SD 0.20) compared with an expected value of 0.90 (P = 0.94).
Genotype and allele frequencies according to glucose tolerance.
There was a near significant difference in the genotype distribution according to glucose tolerance status, with a higher proportion of the Pro/Pro genotype among twins with type 2 diabetes and impaired glucose tolerance (IGT) compared with the group of normal glucose tolerant twins (P = 0.07). Twins with type 2 diabetes and the combined group of twins with type 2 diabetes and IGT had a significantly lower frequency of the minor Ala allele compared with twins with normal glucose tolerance (NGT). No differences in allele frequency of the minor T- allele of the exon 6 polymorphism according to glucose tolerance were observed (Table 2).
Phenotypic appearances of the two variants in the PPARγ gene among MZ and DZ twins.
Both MZ and DZ carriers of the Pro and Ala allele had similar anthropometry. The Ala allele was associated with a significantly lower glucose concentration after 120 min and AUC for glucose compared with the Pro allele. Furthermore, the fasting and 120-min plasma insulin concentrations together with the AUC for insulin were significantly lower among carriers of the Ala allele compared with the Pro allele carriers. Finally, the carriers of the Ala allele were significantly less insulin resistent (HOMA index) compared with the carriers of the Pro allele (Table 3). The exon 6 polymorphism was not associated with measures of anthropometry, plasma glucose, or insulin concentrations (P values between 0.2 and 0.9) among MZ and DZ, respectively (data not shown).
Metabolic profiles according to zygosity.
Independent of genotype, the MZ twins had significantly higher 30-min plasma insulin concentrations and AUC for insulin compared with DZ twins (Table 3). Furthermore, there was a tendency toward higher 30-min plasma glucose concentrations (P = 0.09) and AUC for glucose (P = 0.10) among MZ twins compared with DZ twins.
Twins discordant for the variant genotypes.
We identified 27 dizygotic twin pairs discordant for the Pro12Ala polymorphism (i.e., Pro/Pro versus Pro/Ala). Twins homozygous for the Pro allele had significantly higher 120-min plasma insulin concentrations (P = 0.007) and near significantly higher 120-min plasma glucose concentrations (P = 0.06) and AUC for plasma insulin (P = 0.08) (Table 4). Among the 27 dizygotic twin pairs discordant for the exon 6 C→T polymorphism (i.e., C/C versus C/T), no differences in phenotypes were observed (data not shown).
Several recent studies reported an association between the Ala allele of the codon 12 polymorphism and reduced risk of diabetes (12–15,22), decreased insulin resistance (13,16–18,22), and low (16) or high BMI (24). However, not all studies have been able to replicate these associations (25–31). Our finding of a significantly higher frequency of the Ala allele among twins with NGT provides further evidence in favor of an association between the Ala allele and a reduced risk of type 2 diabetes and/or IGT. This is illustrated by the lower plasma glucose and insulin concentrations during OGTT together with lower HOMA insulin resistance among carriers of the Ala allele compared with the noncarriers. Since measures of obesity were similar in carriers of the different genotypes, obesity does not appear to be a confounding factor in the observed metabolic differences between the different genotypes. Finally, the lower 120-min plasma glucose and insulin concentrations and AUC for insulin in DZ carriers of the Ala allele, compared with their discordant cotwin carrying the Pro allele, further supports this finding.
In the present study, we found no association between the exon 6 C→T polymorphism and variables related to type 2 diabetes or insulin resistance.
The allele frequencies of the exon 6 T-allele (MZ 12.6% and DZ 9.1%) among the present twin population are not in line with the findings of Busjahn et al. (8) (MZ 9% and DZ 19%). Moreover, the degree of allele sharing among the DZ twins did not differ from the expected value for allele sharing. Therefore, our data does not provide any major support for the idea that the silent C→T substitution in exon 6 of the PPARγ gene may be linked to DZ twinning. The present study population, including 550 subjects, is a homogenous population-based sample of Danish Caucasian twins, while the population in the study by Busjahn et al. (8) included 606 subjects from three different European countries. Among these subjects 100 DZ twins were from Finland, a population that exhibits a higher allele frequency of the Pro12Ala polymorphism than other European populations (15). If the exon 6 variant cosegregates with the Pro12Ala variant, the discrepancies in these results could be explained by differences in ethnicity. However, the associations in the study by Busjahn et al. (8) were in fact even more pronounced in the non-Finish populations, which indicates that other ethnic differences or unknown factors may explain the different results.
We observed similar allele frequencies of both PPARγ polymorphisms among MZ and DZ twins. These findings alone indicate that the metabolic differences among MZ and DZ twins cannot be explained by differences in the PPARγ genotypes. When comparing anthropometry and measures of OGTT plasma glucose and insulin concentrations between MZ and DZ twins, MZ twins had significantly higher 30-min and AUCs for plasma insulin compared with DZ twins, independent of genotype. The higher concentrations of 30-min and AUCs for plasma glucose during the OGTT among MZ compared with DZ twins did not reach the level of statistical significance. This can be explained by the smaller number of twins included in the present study and/or the fact that we performed corrections for intrapair relationships. Although the possibility remains that other unknown genotypes may play a role, we still believe that the differences in metabolic profiles between MZ and DZ twins might be explained by differences in their respective intrauterine environment, in accordance with the thrifty phenotype hypothesis (1).
In conclusion, the differences in plasma insulin profiles during OGTT among MZ and DZ twins, the latter being less insulin resistant, cannot be explained by differences in frequencies of the PPARγ variants. The differences may therefore, as previously reported, be due to intrauterine environmental factors independent of genotype. However, with the present study we provide additional evidence in favor of an association between the Pro12Ala variant and reduced risk of diabetes and reduced insulin resistance. Therefore, our study simultaneously provides evidence for a role of both the intrauterine environment and the thrifty phenotype hypothesis and for genetics and the thrifty genotype hypothesis in the etiology of insulin resistance and perhaps glucose intolerance in twins.
This study was supported by grants from the NOVO Foundation, Clinical Research Institute, Odense University, and The Danish Diabetes Association.
Address correspondence and reprint requests to Dr. Pernille Poulsen, Department of Endocrinology, Odense University Hospital, Sdr. Boulevard, 5000 Odense, Denmark. E-mail:; or Dr. Allan Vaag, Steno Diabetes Center, Niels Steensensvej, 2820 Gentofte, Denmark. E-mail: .
Received for publication 28 August 2001 and accepted in revised form 7 October 2002.
P.P. and G.A. contributed equally to this work.
AUC, area under the curve; DZ, dizygotic; HOMA, homeostasis model assessment; IGT, impaired glucose tolerance; MZ, monozygotic; NGT, normal glucose tolerance; OGTT, oral glucose tolerance test; PPAR, peroxisome proliferator–activated receptor.