|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
CLINICAL STUDY |
Departments of1 Clinical Sciences Malmö, Clinical Obesity2 Community Medicine, Malmö University Hospital MAS, Lund University, Malmö, Sweden S-205023 Skaraborg Institute, Skövde, Sweden S-541304 Lund University Diabetes Centre (LUDC), Malmö, Sweden S-20502
(Correspondence should be addressed to L E Johansson; Email: lovisa.johansson{at}med.lu.se)
| Abstract |
|---|
|
|
|---|
Methods: rs738409 (Met148Ile) and rs2072907 (C to G) were genotyped using TaqMan allelic discrimination assay in a Swedish population-based sample (n=1811). Oral glucose tolerance test (OGTT) with data from three time points (0, 30, and 120 min) were available from individuals under the age of 50 years (n=973).
Results: Both variant alleles were associated with decreased prevalence of obesity (P<0.05); odds ratio 0.75 (0.61–0.93) per carried Ile-allele for rs738409 and 0.80 (0.64–1.00) per carried G-allele for rs2072907. For obesity as a quantitative trait, there was no association in the whole population, but in obese subjects body mass index (BMI; P=0.023) and waist (P=0.0098) were higher in carriers of the Ile-allele. The Ile-carriers also displayed decreased insulin secretion in response to OGTT (30 min insulin; P=0.007, insulinogenic index; P=0.0051) with no significant differences in fasting plasma glucose (P=0.31), β-cell function (disposition index; P=0.17) or homeostasis model of assessment insulin resistance (HOMA-IR; P=0.063). The correlation between BMI and HOMA-IR differed (Met/X versus Ile/Ile, P=0.028), Met-allele carriers were seemingly more insulin resistant at a lower BMI. The rs2072907 variant shows similar results for insulin secretion. The significance of this finding remained after adjusting for age, gender, and level of self-reported leisure-time physical activity.
Conclusion: We confirm the association between PNPLA3 and obesity. In addition, the rs738409 variant was associated with insulin secretion. There seems to be a differential effect of the Ile-allele depending on the degree of obesity, possibly as a consequence of insulin resistance.
| Introduction |
|---|
|
|
|---|
In a recent study by us, we have shown that the wild-type alleles of common variants in the PNPLA3 gene were associated with human obesity in Swedish case–control material (ncase/control=234/234) (6). The risk allele of rs2072907 was also associated with decreased subcutaneous but not visceral adipose tissue PNPLA3 mRNA expression (6). Since adipocytes from carriers of this allele exhibited increased basal lipolysis it was suggested that the main function of adiponutrin may be lipogenic rather than lipolytic, i.e., that decreased expression would result in a net increase in lipolysis due to loss of lipogenesis.
The aim of this study was to replicate previous findings of an association between PNPLA3 variants and obesity and to investigate whether there is an effect of these genetic variants on insulin sensitivity and secretion as forerunners of obesity-associated development of diabetes.
| Subjects and methods |
|---|
|
|
|---|
This investigation is based on data from a population-based study in Vara, a small community in a rural area of southwestern Sweden, as part of the Skaraborg Project (10, 11). A total of 1811 subjects were surveyed (Table 1). Obesity was defined as a body mass index (BMI)
30 kg/m2. A diagnosis of diabetes was confirmed if there was a history of a physician's diagnosis, and in new cases when the fasting plasma glucose value was
7.0 mmol/l (twice) and/or when the 2 h plasma glucose value from an oral glucose tolerance test (OGTT) was
11.1 mmol/l. Self-reported leisure-time physical activity (LTPA) was divided into four different categories, ranging from sedentary to strenuous exercise (11, 12). Per definition low level of LTPA includes subjects with an estimated LTPA of 1 and 2 and high level of LTPA subjects with LTPA of 3 and 4. All subjects underwent an OGTT (75 g glucose) with blood drawn at 0 and 120 min. Blood drawn at 30 min were available for all subjects under the age of 50 years and these data were used for further analysis (n=973, non-diabetic subjects, age 40 (36–45) years, BMI 25.4 (23.1–28.1) kg/m2 and 54% males). Serum insulin was analyzed using ELISA with <0.3 percent cross-reactivity for pro-insulin (DRACO Diagnostics Ltd, Ely, Cambridgeshire, UK) (13). Data are presented as pmol/l but calculations using insulin values are based on µU/ml. Plasma glucose levels were analyzed using standard procedures. Insulin resistance was estimated by homeostasis model of assessment insulin resistance (HOMA-IR), calculated by dividing the product of fasting glucose (mmol/l) and insulin levels (µU/ml) by 22.5 (14). Insulinogenic index ((insulin at 30 min–insulin at baseline)/glucose at 30 min), disposition index (insulinogenic index/HOMA-IR) and area under the curve (AUC) were calculated. The regional ethical review board of Göteborg University approved the study and informed written consent was obtained from the subjects prior to their participation.
|
DNA was extracted from whole blood by using Qiagen MiniPrep at the DNA/RNA Genotyping Lab, SWEGENE Resource Center for Profiling Polygenic Disease, Lund University, Malmö University Hospital, Malmö, Sweden. Two single nucleotide polymorphisms (SNPs) were chosen for genotyping based on a previous study (6). Two SNPs showed strong association with obesity, the rs1010022 and rs2072907. These two SNPs were in high linkage disequilibrium and therefore only the defined tagSNP rs2072907 was chosen for genotyping in this study. SNP rs738409 was chosen because it results in an amino acid switch (Met148Ile) in the important lipolytic and lipogenic patatin domain in the protein and has previously shown marginal association with obesity (6). Genotyping was performed using TaqMan allelic discrimination on the ABI 7900HT from Applied Biosystems (Foster City, CA, USA, Assay by Design). Sequences for rs2072907 forward: CTGGAAGGCAGGTGTAACCA; reverse: CTTTGTGGGCTCCATCTACATATCT; probes: CAGTTGTTATAAACGAACACTA (VIC) and CAGTTGTTATAAACCAACACTA (FAM). Sequences for rs738409 forward: GCTTTCACAGGCCTTGGTATG; reverse: GGAGGGATAAGGCCACTGTAGAA; probes: VIC-TTCCTGCTTCATCCC-MGB and FAM-TTCCTGCTTCATGCC-MGB. Genotyping was performed in 5 µl reaction volumes with 5 ng DNA and Universal TaqMan 2xPCR MasterMix according to the manufacturer's recommendations (Applied Biosystems). Success rate for all genotyping was >99%. For both SNPs, genotyping failed in six subjects. Using the same method, 3–5% of the samples were rerun with a concordance of 100%. There were no deviations from Hardy–Weinberg equilibrium.
Statistical analysis
The
2-test was used to compare the proportions and frequencies. Associations between obesity and categorical variables were analyzed by using binary logistic regression in SPSS v15.0 (Chicago, IL, USA) defining the change in odds ratio (OR) per allele, expressed as OR with 95% confidence interval. Analysis of covariance was used to investigate the effect of genotypes on clinical variables. Age, gender, and LTPA were used as covariates. Regression was analyzed using Pearson's correlation in subgroups based on a dominant model of rs738409 (Met/X versus Ile/Ile) defined by data from Table 3. The slopes of the two regression lines were compared using seemingly unrelated estimation test in STATA. P<0.05 was considered statistically significant. Data were log transformed for normal distribution and presented as median with interquartile range (25th–75th percentile) within brackets or mean±S.E.M. Statistical calculations were performed using Number Cruncher Statistical Systems 2000 software (NCSS, Kaysville, UT, USA) unless stated otherwise.
|
was set at 0.05 and the current sample size of 379 cases and 1432 controls were used the power for SNP rs738409 would be 59 and 90% for rs2072907. | Results |
|---|
|
|
|---|
The study population has been described previously (11). In brief, the distribution of males and females was equal and showed no differences in age (Table 1). The prevalence of diabetes was 6.0%. Males had slightly higher BMI and reported higher LTPA compared with females. There was an increased prevalence of obesity among subjects reporting low LTPA (OR 2.27 (1.71–3.01), P<0.0001). BMI was 26.5 (24.1–30.0) and 25.2 (23.3–27.7) kg/m2 for subjects reporting low and high LTPA respectively (P<0.0001).
The minor allele frequencies in the whole population were 21 and 17% for rs738409 and rs2072907 respectively. There were no differences in age or gender between rs738409 genotype carriers but there was a difference in age concerning rs2072907 (P=0.0067, Table 2).
|
Binary logistic regression suggested a decreased risk of obesity with an OR of 0.75 (0.61–0.93) per Ile-allele of the rs738409 polymorphism when adjusting for age, gender, and LTPA, P=0.009 (adjusting for age and gender only, OR 0.77 (0.63–0.95) per Ile-allele, P=0.016). Similarly, for rs2072907, there was a decreased risk of obesity with an OR of 0.80 (0.64–1.00) per carried G-allele adjusted for age, gender, and LTPA, P=0.049 (adjusting for age and gender only, OR 0.82 (0.66–1.02) per carried G-allele, P=0.069).
There was no significant difference when analyzing obesity as a quantitative trait in the whole population for either polymorphism (Table 2). However, analyzing the traits in obese subjects separately revealed a significant difference. BMI was higher in carriers of the rs738409 Ile-allele after adjustments for age, gender, and LTPA (Met/Met 32.2 (31.1–34.9) n=255, Met/Ile 33.2 (31.5–37.1) n=114 and Ile/Ile 32.2 (31.0–38.2) kg/m2 n=9, P=0.023). The significance of this difference disappeared when adjusting for HOMA-IR (P=0.12). Waist was also associated with the rs738409 genotypes after adjustments for age, gender, and LTPA (Met/Met 105 (98–111) n=254, Met/Ile 107 (100–115) n=114 and Ile/Ile 109 (106–116) cm n=9, P=0.0098) but not when adjusting for HOMA-IR (P=0.050). There were no effects of the rs2072907 polymorphism on quantitative obesity traits in the obese (data not shown).
Genetic variants in the adiponutrin gene are associated with insulin secretion and insulin resistance
Data on insulin levels at 0, 30, and 120 min from an OGTT were available from 973 non-diabetic individuals, all of whom were under the age of 50. With no differences in plasma glucose at these time points, first-phase insulin secretion defined as insulin at 30 min during OGTT was significantly different between Met148Ile genotypes (Table 3). The same relationship was found for AUC between 0 and 30 min (Table 3 and Fig. 1). Insulinogenic index was higher in carriers of the Met-allele but insulin resistance measured as HOMA-IR, and β-cell function measured as disposition index, did not differ significantly (Table 3). Moreover, the slope of the regression line when looking at the correlation between BMI and HOMA-IR was significantly steeper in homozygous carriers of the variant Ile-allele than in carriers of the Met-allele (R=0.67 vs 0.5367, P=0.028, Fig. 2) suggesting that the carriers of the Met-allele are more insulin resistant at a lower BMI. There were no differences between the correlations for disposition index (R=–0.31 vs –0.16, P=0.61, Fig. 2). Results concerning rs2072907 trended towards the same differences (Table 3). Adjusting the analysis regarding rs738409 for SNP rs2072907 did not change the results but performing the analysis vice versa attenuated the findings (data not shown).
|
|
| Discussion |
|---|
|
|
|---|
In a previous study we found that four PNPLA3 SNPs (rs738409, rs2076211, rs2072907, and rs1010022) were associated with obesity, two of them (the tagSNP rs2072907 and rs1010022) remaining significant after correction for multiple testing (6). Contrary to our previous investigation, rs2072907, a non-coding variant located in intron 5, had only minor or no effects on the risk of obesity in the current study. Previous data suggested that the two SNPs were located in two separate haplotype blocks although some linkage was seen between the two (D'=0.87 (0.82–0.92), r2=0.60). Here, the two SNPs show similar linkage disequilibrium (D'=0.90 (0.86–0.92), r2=0.65) (6). We cannot exclude that rs2072907 reflects the effect of another causative variant(s), possibly rs738409, by linkage disequilibrium in the previous study. It is important to bear in mind that there are differences in both design and sample selection between the studies. The subjects in this study were randomly selected from a population while the previous study was based on severely obese patients referred to an obesity clinic at the hospital (BMI 32.5 (31.1–35.7) vs 40.3 (35.5–45.3) kg/m2). Differences in the level of physical activity might also influence the difference between the two studies. We lack these data for the previously investigated patient sample but due to their severe obesity the level of physical activity is expected to be limited.
In this Swedish population-based study we found that the wild-type Met-allele of rs738409 shows association with obesity and measures of increased insulin secretion. The fact that the Met-allele is associated with both increased risk of obesity and increased insulin secretion but not directly with insulin resistance in the whole population may at first seem contradictory. Focusing on the obese subjects separately, the variant Ile-allele rather than the previously defined risk, the Met-allele, was associated with a slightly higher BMI (6). A difference in insulin resistance might have explained these contradictory results but there does not seem to be an immediate effect of rs738409 on HOMA-IR. A possible key to this issue comes from the correlation between BMI and insulin resistance. This analysis revealed that insulin resistance rose more rapidly with increasing BMI in carriers of the Ile/Ile genotype compared with Met-allele carriers. In other words, carriers of the Ile-allele may be more sensitive to developing disease as a consequence of being obese. However, in light of the greater degree of insulin resistance at lower levels of BMI found in carriers of the Met-allele, obesity may have developed as an anabolic consequence of the resulting increase in insulin secretion in these subjects. In fact, a pilot study investigating rs738409 in a small type 2 diabetes case–control material suggests that the Ile/Ile genotype is more frequent in type 2 diabetes cases (data not shown). There are other examples of genetic variants that exhibit these apparently perplexing associations. For example, whereas the Ala-allele of the Pro12Ala polymorphism of the peroxisome proliferator-activated receptor gamma (PPARG) gene is associated with a protection against type 2 diabetes, it is also associated with an increased risk of developing obesity (15, 16). In this case, the common denominator is the increased insulin sensitivity associated with PPARG 12Ala. Both increased insulin secretion and enhanced insulin sensitivity will ultimately lead to the same result, a small but significant weight gain over time. Clearly, the effect of rs738409 on both insulin resistance and type 2 diabetes will require further investigation.
Adipose tissue PNPLA3 mRNA expression is influenced by nutritional status, insulin and insulin sensitivity, degree of obesity, and genetic variants in the PNPLA3 gene (3, 4, 6, 7, 8, 9, 17, 18, 19). The knowledge regarding the role of adiponutrin in the adipocyte or any other cell is still very limited, making it difficult to functionally understand a link between adiponutrin and insulin secretion. Likely, the link resides in insulin resistance rather than a direct interaction between the two in the insulin secreting β-cell. All the same, since fatty acids have been suggested to play a pivotal role both in the production and secretion of insulin (20, 21) and rs738409 is located in the part of the gene coding for the so-called patatin domain, which is responsible for catalyzing the cleavage of fatty acids from triglycerides (22), a direct link cannot be excluded. Hypothetically, the PNPLA3 protein could influence insulin secretion through either its lipolytic or lipogenic properties by altering the load of fatty acids into the circulation and/or in the β-cell itself.
In conclusion, we have presented data confirming the association between a genetic variant in the adiponutrin gene with obesity and further extend the findings by showing an association with insulin secretion in a Swedish population-based sample. Whether this association is a cause or a consequence of the association with obesity remains to be elucidated.
| Declaration of interest |
|---|
|
|
|---|
| Funding |
|---|
| Acknowledgements |
|---|
| References |
|---|
|
|
|---|
Pro12Ala polymorphism is associated with decreased risk of type 2 diabetes. Nature Genetics 2000; 26:76–80.[CrossRef][Web of Science][Medline]16. Masud S & Ye S. The SAS Group. Effect of the peroxisome proliferator activated receptor-gamma gene Pro12Ala variant on body mass index: a meta-analysis. Journal of Medical Genetics 2003; 40:773–780.This article has been cited by other articles:
![]() |
N. A. Noordeen, T. K. Khera, G. Sun, E. R. Longbottom, T. J. Pullen, G. da Silva Xavier, G. A. Rutter, and I. Leclerc Carbohydrate-Responsive Element-Binding Protein (ChREBP) Is a Negative Regulator of ARNT/HIF-1{beta} Gene Expression in Pancreatic Islet {beta}-Cells Diabetes, January 1, 2010; 59(1): 153 - 160. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Kantartzis, A. Peter, F. Machicao, J. Machann, S. Wagner, I. Konigsrainer, A. Konigsrainer, F. Schick, A. Fritsche, H.-U. Haring, et al. Dissociation Between Fatty Liver and Insulin Resistance in Humans Carrying a Variant of the Patatin-Like Phospholipase 3 Gene Diabetes, November 1, 2009; 58(11): 2616 - 2623. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. C. Kienesberger, M. Oberer, A. Lass, and R. Zechner Mammalian patatin domain containing proteins: a family with diverse lipolytic activities involved in multiple biological functions J. Lipid Res., April 1, 2009; 50(Supplement): S63 - S68. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |