© 2001 by the American Diabetes Association, Inc.
Interactions Between Variants in the ß3-Adrenergic Receptor and Peroxisome ProliferatorActivated Receptor-
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
2. We hypothesized that these single gene variants may mark mutations that act through convergent pathways to produce synergistic effects on obesity. RESEARCH DESIGN AND METHODSThe sample included 453 subjects from 10 large Mexican-American families participating in the population-based San Antonio Family Heart Study. The effects of each gene variant singly and jointly were estimated as fixed effects using the measured genotype approach framework. Analyses were conditioned on the pedigree structures to account for the correlations among family members. Statistical significance was evaluated by the likelihood ratio test with adjustment for age, sex, and diabetes status.
RESULTSThe allele frequencies for the ADRß3 Trp64Arg and PPAR
2 Pro12Ala variants were 18 and 12%, respectively. The ADRß3 variant was not significantly associated with any of the obesity-related traits, but subjects with the PPAR
2 variant (n = 98) had significantly higher levels of fasting insulin (P = 0.03), leptin (P = 0.009), and waist circumference (P = 0.03) than those without. Subjects with both gene variants (n = 32) had significantly higher BMI, insulin, and leptin levels than those with only the PPAR
2 variant (n = 66) (P for interaction: 0.04, 0.02, and 0.01 for BMI, fasting insulin, and leptin, respectively).
CONCLUSIONSOur results suggest that epistatic models with genes that have modest individual effects may be useful in understanding the genetic underpinnings of typical obesity in humans.
Abbreviations: ADRß3, ß3-adrenergic receptor BAT, brown adipose tissue FFA, free fatty acid PPAR, peroxisome proliferatoractivated receptor TZD, thiazolidinedione
| INTRODUCTION |
|---|
|
|
|---|
In our previous studies of Mexican-Americans, we have implicated a role for variants in two obesity candidate genes, the ß3-adrenergic receptor (ADRß3) and the peroxisome proliferatoractivated receptor-
2 (PPAR
2), in the etiology of obesity (36). Both variants are associated with amino acid substitutions that may result in altered function of their respective gene products (79). ADRß3 is expressed in adipocytes and mediates the rate of lipolysis in response to catecholamines, whereas PPAR
2 is a nuclear receptor that acts as a transcription factor to regulate adipocyte differentiation and insulin sensitivity. Among our study participants, we observed an association between the Trp64Arg variant of ADRß3 and obesity (as have many others [10,11]), although the association in our population could be detected only by matching siblings for a previously detected obesity quantitative trait locus on chromosome 2 (4). We have subsequently detected an association between the Pro12Ala variant of PPAR
2 and obesity (6), which has also been reported by others (1214).
The ADRß3 Arg and PPAR
2 Ala variants are reasonably common in diverse populations, and thus, the frequency of individuals who harbor both variants is appreciable. Because both of these gene products influence fat accumulation, we evaluated evidence for interaction between these two variants to determine whether subjects having variants in both ADRß3 and PPAR
2 are at increased risk for obesity beyond the expected risk based on the additive effects associated with each individual variant. In addition, we have estimated the combined contribution of these two variants (both independent and interaction effects) to the total variation of obesity-related traits.
| RESEARCH DESIGN AND METHODS |
|---|
|
|
|---|
16 years of participating probands who were ascertained without regard to the probands medical condition. This report is based on 453 subjects comprising the initial set of 10 families who were genotyped for both the ADRß3 Trp64Arg and the PPAR
2 Pro12Ala variants. Informed consent was obtained from all participants, and all protocols were approved by the institutional review board. The examination procedure included measurements of blood pressure, anthropometry, and an oral glucose tolerance test. Serum glucose concentrations were measured using an Abbott V/P Analyzer (Abbott Laboratories, Abbott Park IL), and serum concentrations of insulin and leptin were measured using commercial radioimmunoassay kits (insulin, Diagnostic Products, Los Angeles, CA; leptin, Linco, St. Louis, MO), as previously described (15,16). The coefficient of variation for fasting insulin and leptin were 6.5 and 7.5%, respectively. Diabetes was diagnosed according to World Health Organization criteria (17). BMI was calculated as weight (kilograms) divided by height (meters squared). Related obesity measures included waist circumference, measured to the nearest centimeter, and total body fat (or fat mass, in kilograms), assessed by bioimpedance. In addition, we accessed smoking habits, physical activity by 7-day recall, and daily dietary caloric intake by food frequency questionnaire (15).
Genotyping
DNA was isolated from lymphocytes using standard methods. The Trp64Arg variant of the ADRß3 and the Pro12Ala variant of PPAR
2 were inferred from genotypes generated using polymerase chain reactionbased assays as previously described (4,12).
Statistical analysis
Mean levels of obesity traits were estimated according to ADRß3 and PPAR
2 genotypes. To account for the relatedness among family members, we used the measured genotype approach (18), in which we estimated the likelihood of specific genetic models given the pedigree structure. For example, we compared the likelihood of a full model, which allowed for genotypic-specific means, to that of a nested model in which genotypic means were restricted to be equal to each other. Parameter estimates were obtained by maximum likelihood methods. The significance of the association was tested by likelihood ratio tests, which compare the difference in the likelihoods between the full and the nested models. Two times the difference between the logarithm of the likelihoods of the two models is distributed asymptotically as a
2 statistic with degrees of freedom equal to the difference in the numbers of parameters in the models being compared. Within each model, we simultaneously estimated the effects of age, sex, and diabetes status (yes/no). In addition, we allowed for a residual heritability of the trait. Fasting insulin and leptin values were transformed by their natural logarithms (ln) to reduce skewness. Analyses of insulin levels were performed using data from nondiabetic subjects only. All analyses were conducted using the SOLAR program (19).
We tested for interaction between the ADRß3 Arg and the PPAR
2 Ala variants by comparing the likelihood of a full model in which we parameterized the genetic effect into three components: an independent effect of having the ADRß3 variant, an independent effect of having the PPAR
2 variant, and a joint effect of having both variants (ADRß3 x PPAR
2) with the likelihood of a nested model in which we parameterized the genetic effect into two independent componentsan effect of having the ADRß3 variant and an effect of having the PPAR
2 variant. We simultaneously estimated the effects of age, sex, and diabetes status in each model and allowed for a residual heritable component of the trait.
We estimated the proportion of the total trait variance accounted for by these variants (both their main effects and interaction effects) by calculating the reduction in the total trait variance (after adjusting for age, sex, and diabetes status only) that occurs after also estimating the effects of these two variants and their interaction and then expressing this reduction as a proportion of the total trait variance. In other words, we estimated (
2tot
2adj)/
2tot, where
2tot equals the total trait variance (after adjusting for age, sex, and diabetes status only) and
2adj equals the trait variance after further adjustment for genetic effects.
| RESULTS |
|---|
|
|
|---|
2 Ala variants were 18 and 12%, respectively. The genotype frequencies did not deviate from those expected on the basis of the Hardy-Weinberg equilibrium. Mean ages of men and women were 39.3 and 38.3 years, respectively, whereas mean BMI was 28.7 kg/m2 for men and 30.4 kg/m2 for women (P < 0.05). Women had significantly higher mean fat mass and leptin levels than men. There were no significant sex differences in waist circumference or insulin levels.
|
2 revealed that subjects with the Ala variant were significantly more obese (i.e., had larger waist circumference [P = 0.03] and higher leptin [P = 0.009] and insulin levels [P = 0.03]) than subjects homozygous for the more common Pro variant. The association of the PPAR
2 Ala variant with obesity has been previously reported (6).
We next evaluated the effect of both variants on the variation of obesity-related traits simultaneously. Table 2 shows the adjusted mean levels of the obesity-related traits according to subjects combined ARDß3 and PPAR
2 genotypes. Only eight combinations are shown because none of the subjects were homozygous for both variants. There were 234 subjects who did not have either gene variant, 121 subjects who had the ADRß3 Arg variant only, 66 subjects who had the PPAR
2 Ala variant only, and 32 subjects who had both variants. In general, mean levels of most traits (especially BMI, insulin, and leptin) tended to be higher among those with the PPAR
2 Ala variant than in those without and higher still in those having both variants.
|
2 only, or in both genes. Figure 1 shows the effect of each variant combination compared with the baseline group in which subjects did not have a variant in either gene. Mean BMI and mean concentrations of insulin and leptin were all higher in subjects with variants in both genes than in subjects having variants in only one of the two genes. The formal test for interaction revealed that the model in which we allowed for an interaction effect between these two variants (i.e., the model with three genetic effect parameters: ADRß3 Arg variant + PPAR
2 Ala variant + ADRß3 Arg variant x PPAR
2 Ala variant) provided a significantly better fit to the data (i.e., higher likelihood) than the model in which we did not allow for an interaction (i.e., the model with only two independent genetic effect parameters: ADRß3 Arg variant + PPAR
2 Ala variant). The age-, sex-, and diabetes-adjusted significance levels for the interaction effects were 0.04 for BMI, 0.02 for insulin, 0.01 for leptin, 0.26 for waist circumference, and >0.90 for fat mass. Further adjustment for the effects of smoking, level of physical activity, and total daily caloric intake altered these results very little (P = 0.08 for BMI, P = 0.007 for insulin, and P = 0.02 for leptin), and the results were not altered significantly when the 69 subjects with diabetes were removed from the analysis.
|
| CONCLUSIONS |
|---|
|
|
|---|
2 may not be uniform for all carriers of the Ala allele but may be restricted primarily to the group of subjects who also carry the Arg variant of ADRß3. In fact, when we eliminated the 32 subjects carrying variants in both genes from the analysis, the association between the Ala variant and obesity was no longer statistically significant (adjusted P values: 0.28 for BMI, 0.22 for insulin, and 0.13 for leptin; data not shown). These results raise the intriguing speculation that the effect of the PPAR
2 Ala variant on obesity might be expressed only in the presence of some other factors (such as the ADRß3 variant). Conversely, this speculation also suggests a possible explanation for the failure of many association studies to detect an effect of the ADRß3 Trp64Arg variant on obesitynamely that these studies may have contained an insufficient number of subjects with both factors (i.e., variants in both genes).
Emerging data suggest that PPAR
2 is regulated by a complex set of factors, and possibly, the obese state itself may modify the expression of this gene (20). Expression of PPAR
2 mRNA in adipose tissue is increased in obese humans but is downregulated while following a low-calorie diet (21). A potential role for the Ala variant in PPAR
2 expression has been postulated (20), and epidemiological data provide some support. For example, Ek et al. (13) observed that in obese subjects, individuals homozygous for the Ala variant had significantly higher BMIs than other subjects, whereas in lean subjects, those homozygous for the Ala variant had significantly lower mean BMI. Several additional studies have also reported an association between the Ala variant and higher BMI in obese individuals (12,14), and an association is also observed in our study of a relatively obese population (6). On the other hand, Deeb et al. (7) studied a relatively leaner Finnish population and found that subjects with the PPAR
2 Ala variant were less obese and had greater insulin sensitivity than subjects lacking the variant.
We did not observe an association between type 2 diabetes and the PPAR
2 Ala variant in this Mexican-American population (6), although the power to detect a small variant effect in our sample may be limited. Somewhat paradoxically, however, several studies (7,22,23), including a recent meta-analysis (23), have reported the PPAR
2 Ala variant to be associated with a decreased risk of type 2 diabetes. It is intriguing to speculate how the same mutation could be associated with both obesity and reduced risk of diabetes. Although epidemiological studies show a strong association between obesity and type 2 diabetes, it is important to realize that only 30% of obese subjects are diabetic, suggesting that obesity is not sufficient to cause diabetes. Indeed, some investigators (24) have even argued that not all types of obesity increase the risk of type 2 diabetes. Consistent with this view is the observation that enhanced insulin sensitivity, which is associated with a decreased risk of diabetes, is a predictor of weight gain (25). An increased number of adipocytes may not lead to increased risk of diabetes, and in fact, in mice with lipoatrophic diabetes, hyperglycemia was reversed after a surgical transplantation of adipose tissue (26). Adipocyte volume, on the other hand, is a predictor of type 2 diabetes, independent of body fat and body fat distribution, as shown in Pima Indians (27). These observations suggest that adipocyte tissue may promote antidiabetic effects via metabolic or endocrine mechanisms, such as increasing uptake of glucose, free fatty acid (FFA), or hormonal secretion to increase insulin sensitivity. However, if the adipocytes ability to differentiate is compromised, its antidiabetic effect may be lost once it has reached its maximum capacity to store energy (e.g., glucose and FFA).
The mechanism(s) by which mutations in both PPAR
2 and ADRß3 might interact to enhance obesity and insulin resistance is not yet proven. However, because both genes play important regulatory roles in adipocytes, there are several plausible mechanisms. PPAR
is a nuclear hormone receptor that stimulates transcription of multiple genes necessary for adipogenesis and insulin signaling, whereas the ADRß3 is a membrane receptor that increases intracellular cAMP concentrations and stimulates lipolysis. Singly, functional studies indicate that the Trp64Arg substitution in ADRß3 and the Pro12Ala substitution in PPAR
2 result in decreased activity of their respective proteins (79,28,29). The most direct piece of evidence supporting a direct interaction between ADRß3 and PPAR
2 expression comes from the recent work of Bakopanos and Silva (30), who observed that thiazolidinediones (TZDs), insulin sensitizers that exert their effects through the specific binding and activation of PPAR
s (31), reduce ADRß3 expression at the transcriptional level. Furthermore, this effect was evident only on adipose cell lines containing PPAR
2. In a second experiment, the inhibition effect by TZDs was also observed in a cell line that lacked PPAR
but which was transfected by PPAR
2, and the inhibition increased in a dose-dependent manner with the amount of PPAR
2 cDNA introduced in these cells.
Another study has indicated that adrenergic stimulation of rat brown adipose tissue (BAT) through cold exposure results in a marked increase in the expression of PPAR
(32). Because ADRß3 is the major adrenergic receptor subtype in rat BAT, these findings suggest a direct interaction between ADRß3 activation and PPAR
expression. Because the Trp64Arg substitution in ADRß3 leads to a decrease in intracellular signal transduction (28), it is possible that subjects with the Trp64Arg variant may have lower levels of expression of PPAR
, which may compound defects in PPAR
action caused by the Pro12Ala PPAR
2 substitution. Indeed, Gros et al. (33) demonstrated that forced expression of human ADRß3 in CHO/K1 fibroblasts resulted in a marked increase in PPAR
expression. The wild-type form (Trp64 ADRß3) increased PPAR
expression, whereas the mutant form (Arg64 ADRß3) did not. Together, these findings strongly suggest that the presence of both variants could result in significantly lesser PPAR
activity. The simplest model of PPAR
action (consistent with in vitro studies using PPAR
agonists) would suggest that lower PPAR
activity in vivo should associate with lesser obesity and lesser insulin sensitivity. However, it was shown recently that mice heterozygous for a knockout of PPAR
have surprisingly greater insulin sensitivity than the wild-type mice (34). This paradox makes it clear that the mechanisms whereby PPAR
influences obesity and insulin sensitivity in vivo are poorly understood and likely to be complex. Additional studies in humans will be required to clarify these and perhaps other mechanisms whereby these two genes may interact.
In summary, we observed that the PPAR
2 Ala variant, but not the ADRß3 Arg variant, was significantly associated with obesity. However, subjects with variants in both genes were more obese than subjects with the PPAR
2 Ala variant only, and, in fact, it appears that only those subjects carrying both the PPAR
2 Ala and ADRß3 Arg variants have increased obesity. In terms of population significance, the effects of these two variants on obesity in this population are small. Together, they accounted for 23% of the total variation in BMI, insulin, and leptin, with the interaction accounting for about half of their total effects. Of course, we evaluated the effects of only two specific gene variants in this study. If there are other functional variants in these genes or in their regulatory regions, then the effects of these genes on obesity will be larger. More importantly, as there is an interaction between these two gene variants, it is likely that there are interactions among variants of other genes. Although these analyses are to some extent exploratory, they support the importance of these variants in obesity and reinforce the idea that gene effects will be underestimated if only their individual effects are considered. Incorporating interactions, whether gene by gene or gene by environment, may enhance our ability to detect genetic effects and may contribute to our understanding of the etiology of obesity.
| ACKNOWLEDGMENTS |
|---|
We are deeply grateful for the cooperation of the families participating in the San Antonio Family Heart Study.
A summary of this study was previously published as an abstract in Am J Human Genet 65 (Suppl. 4):A14, 1999.
| FOOTNOTES |
|---|
Received for publication 30 October 2000 and accepted in revised form 2 January 2001.
A table elsewhere in this issue shows conventional and Système International (SI) units and conversion factors for many substances.
| References |
|---|
|
|
|---|
2 (PPAR-
2) is associated with measures of obesity in Mexican Americans. Int J Obes 24:522524, 2000
2 associated with decreased receptor activity, lower BMI and improved insulin sensitivity. Nat Genet 20:284287, 1998[Medline]
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]
gene are associated with severe overweight among obese women. J Clin Endocrinol Metab 84:37083712, 1999
, the ultimate thrifty gene. Diabetologia 42:10331049, 1999[Medline]
2 on thiazolidinedione-induced adipogenesis. Biochem Biophys Res Commun 268:178182, 2000[Medline]
deficiency. J Clin Invest 105:287292, 2000[Medline]This article has been cited by other articles:
![]() |
M. H. Black, T. E. Fingerlin, H. Allayee, W. Zhang, A. H. Xiang, E. Trigo, J. Hartiala, A. B. Lehtinen, S. M. Haffner, R. N. Bergman, et al. Evidence of Interaction Between PPARG2 and HNF4A Contributing to Variation in Insulin Sensitivity in Mexican Americans Diabetes, April 1, 2008; 57(4): 1048 - 1056. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. Yang, T. Kelly, and J. He Genetic Epidemiology of Obesity Epidemiol. Rev., June 12, 2007; (2007) mxm004v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Soriguer, S. Morcillo, F. Cardona, G. Rojo-Martinez, M. de la Cruz Almaraz, M. de la Soledad Ruiz de Adana, G. Olveira, F. Tinahones, and I. Esteva Pro12Ala Polymorphism of the PPARG2 Gene Is Associated with Type 2 Diabetes Mellitus and Peripheral Insulin Sensitivity in a Population with a High Intake of Oleic Acid J. Nutr., September 1, 2006; 136(9): 2325 - 2330. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Vanttinen, P. Nuutila, J. Pihlajamaki, K. Hallsten, K. A. Virtanen, R. Lautamaki, P. Peltoniemi, J. Kemppainen, T. Takala, A. P. M. Viljanen, et al. The Effect of the Ala12 Allele of the Peroxisome Proliferator-Activated Receptor-{gamma}2 Gene on Skeletal Muscle Glucose Uptake Depends on Obesity: A Positron Emission Tomography Study J. Clin. Endocrinol. Metab., July 1, 2005; 90(7): 4249 - 4254. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. K. Hansen, E.-M. D. Nielsen, J. Ek, G. Andersen, C. Glumer, B. Carstensen, P. Mouritzen, T. Drivsholm, K. Borch-Johnsen, T. Jorgensen, et al. Analysis of Separate and Combined Effects of Common Variation in KCNJ11 and PPARG on Risk of Type 2 Diabetes J. Clin. Endocrinol. Metab., June 1, 2005; 90(6): 3629 - 3637. [Abstract] [Full Text] [PDF] |
||||
![]() |
S Masud and S Ye Effect of the peroxisome proliferator activated receptor-{gamma} gene Pro12Ala variant on body mass index: a meta-analysis J. Med. Genet., October 1, 2003; 40(10): 773 - 780. [Full Text] [PDF] |
||||
![]() |
M. Stumvoll and H. Haring The Peroxisome Proliferator-Activated Receptor-{gamma}2 Pro12Ala Polymorphism Diabetes, August 1, 2002; 51(8): 2341 - 2347. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Diabetes | Diabetes Care | Clinical Diabetes | Diabetes Spectrum |