Paradoxical Lower Serum Triglyceride Levels and Higher Type 2 Diabetes Mellitus Susceptibility in Obese Individuals with the PNPLA3 148M Variant

Abstract

Background

Obesity is highly associated with elevated serum triglycerides, hepatic steatosis and type 2 diabetes (T2D). The I148M (rs738409) genetic variant of patatin-like phospholipase domain-containing 3 gene (PNPLA3) is known to modulate hepatic triglyceride accumulation, leading to steatosis. No association between PNPLA3 I148M genotype and T2D in Europeans has been reported. Aim of this study is to examine the relationship between PNPLA3 I148M genotypes and serum triglycerides, insulin resistance and T2D susceptibility by testing a gene-environment interaction model with severe obesity.

Methods and Findings

PNPLA3 I148M was genotyped in a large obese cohort, the SOS study (n = 3,473) and in the Go-DARTS (n = 15,448), a T2D case-control study. Metabolic parameters were examined across the PNPLA3 I148M genotypes in participants of the SOS study at baseline and at 2- and 10-year follow up after bariatric surgery or conventional therapy. The associations with metabolic parameters were validated in the Go-DARTS study. Serum triglycerides were found to be lower in the PNPLA3 148M carriers from the SOS study at baseline and from the Go-DARTS T2D cohort. An increased risk for T2D conferred by the 148M allele was found in the SOS study (O.R. 1.09, 95% C.I. 1.01-1.39, P = 0.040) and in severely obese individuals in the Go-DARTS study (O.R. 1.37, 95% C.I. 1.13-1.66, P = 0.001). The 148M allele was no longer associated with insulin resistance or T2D after bariatric surgery in the SOS study and no association with the 148M allele was observed in the less obese (BMI<35) individuals in the Go-DARTS study (P for interaction  = 0.002). This provides evidence for the obesity interaction with I48M allele and T2D risk in a large-scale cross-sectional and a prospective interventional study.

Conclusions

Severely obese individuals carrying the PNPLA3 148M allele have lower serum triglyceride levels, are more insulin resistant and more susceptible to T2D. This study supports the hypothesis that obesity-driven hepatic lipid accumulation may contribute to T2D susceptibility.

About the Authors:

Colin N. A. Palmer

Contributed equally to this work with: Colin N. A. Palmer, Cristina Maglio

Affiliation: Medical Research Institute, Ninewells Hospital and Medical School, University of Dundee, Dundee, Scotland, United Kingdom

Cristina Maglio

Contributed equally to this work with: Colin N. A. Palmer, Cristina Maglio

Affiliation: Department of Molecular and Clinical Medicine and Center for Cardiovascular and Metabolic Research, Sahlgrenska Academy, University of Gothenburg, Gothenburg, Sweden

Carlo Pirazzi

Affiliation: Department of Molecular and Clinical Medicine and Center for Cardiovascular and Metabolic Research, Sahlgrenska Academy, University of Gothenburg, Gothenburg, Sweden

Maria Antonella Burza

Affiliation: Department of Molecular and Clinical Medicine and Center for Cardiovascular and Metabolic Research, Sahlgrenska Academy, University of Gothenburg, Gothenburg, Sweden

Martin Adiels

Affiliation: Department of Molecular and Clinical Medicine and Center for Cardiovascular and Metabolic Research, Sahlgrenska Academy, University of Gothenburg, Gothenburg, Sweden

Lindsay Burch

Affiliation: Medical Research Institute, Ninewells Hospital and Medical School, University of Dundee, Dundee, Scotland, United Kingdom

Louise A. Donnelly

Affiliation: Medical Research Institute, Ninewells Hospital and Medical School, University of Dundee, Dundee, Scotland, United Kingdom

Helen Colhoun

Affiliation: Medical Research Institute, Ninewells Hospital and Medical School, University of Dundee, Dundee, Scotland, United Kingdom

Alexander S. Doney

Affiliation: Medical Research Institute, Ninewells Hospital and Medical School, University of Dundee, Dundee, Scotland, United Kingdom

John F. Dillon

Affiliation: Medical Research Institute, Ninewells Hospital and Medical School, University of Dundee, Dundee, Scotland, United Kingdom

Ewan R. Pearson

Affiliation: Medical Research Institute, Ninewells Hospital and Medical School, University of Dundee, Dundee, Scotland, United Kingdom

Mark McCarthy

Affiliation: Oxford Centre for Diabetes, Endocrinology and Metabolism, University of Oxford, Oxford, United Kingdom

Andrew T. Hattersley

Affiliation: Peninsula NIHR Clinical Research Facility, Peninsula College of Medicine and Dentistry, University of Exeter, Exeter, United Kingdom

Tim Frayling

Affiliation: Peninsula NIHR Clinical Research Facility, Peninsula College of Medicine and Dentistry, University of Exeter, Exeter, United Kingdom

Andrew D. Morris

Affiliation: Medical Research Institute, Ninewells Hospital and Medical School, University of Dundee, Dundee, Scotland, United Kingdom

Markku Peltonen

Affiliation: Chronic Disease Epidemiology and Prevention Unit, Department of Chronic Disease Prevention, National Institute for Health and Welfare, Helsinki, Finland

Per-Arne Svensson

Affiliation: Department of Molecular and Clinical Medicine and Center for Cardiovascular and Metabolic Research, Sahlgrenska Academy, University of Gothenburg, Gothenburg, Sweden

Peter Jacobson

Affiliation: Department of Molecular and Clinical Medicine and Center for Cardiovascular and Metabolic Research, Sahlgrenska Academy, University of Gothenburg, Gothenburg, Sweden

Jan Borén

Affiliation: Department of Molecular and Clinical Medicine and Center for Cardiovascular and Metabolic Research, Sahlgrenska Academy, University of Gothenburg, Gothenburg, Sweden

Lars Sjöström

Affiliation: Department of Molecular and Clinical Medicine and Center for Cardiovascular and Metabolic Research, Sahlgrenska Academy, University of Gothenburg, Gothenburg, Sweden

Lena M. S. Carlsson

Affiliation: Department of Molecular and Clinical Medicine and Center for Cardiovascular and Metabolic Research, Sahlgrenska Academy, University of Gothenburg, Gothenburg, Sweden

Stefano Romeo

* E-mail: [email protected]

Affiliations Department of Molecular and Clinical Medicine and Center for Cardiovascular and Metabolic Research, Sahlgrenska Academy, University of Gothenburg, Gothenburg, Sweden, Department of Clinical and Experimental Medicine, University of Catanzaro, Catanzaro, Italy

Introduction

Hepatic steatosis, serum triglycerides, insulin resistance and type 2 diabetes are tightly associated [1]–[5]. The most widely replicated genetic variant associated with hepatic steatosis is an isoleucine to methionine substitution at position 148 (I148M, rs738409) in the patatin-like phospholipase domain-containing 3 gene (PNPLA3) [6]–[11]. To date, no association between PNPLA3 I148M genotype and type 2 diabetes in Europeans has been reported [9], [12]–[14]. However, a recent study showed that hepatic Pnpla3 overexpression leads to changes in the glucose and triglycerides homeostasis in mouse models of obesity and diabetes [15]. Moreover, genetic reports in Europeans indicate that obesity modifies signals arising from the genome, showing that the PNPLA3 I148M genotype association with increased alanine transferase (ALT) is greatly exacerbated in obesity [12], [16], [17].

To test a gene-environment interaction model with severe obesity in determining metabolic traits, we examined the relationship between PNPLA3 I148M genotypes, serum triglycerides, insulin resistance and type 2 diabetes susceptibility in a large obese cohort, the Swedish Obese Subjects (SOS) study, before and after bariatric surgery. We also validated our findings in the Genetics of Diabetes Audit and Research Tayside Scotland (Go-DARTS) study, a type 2 diabetes case-control study.

Methods

Study design

This genetic study involved two independent European cohorts: the SOS and the Go-DARTS studies (Figure 1). The SOS is a prospective controlled intervention trial that compares the long term effects of bariatric surgery and conventional care in obese subjects (Figure 1A). The Go-DARTS is a cross sectional population case control study for type 2 diabetes comprising normal weight to morbidly obese individuals (Figure 1B).

[Figure omitted. See PDF.]

Figure 1. Study design.

PNPLA3 I148M was genotyped in the Swedish Obese Subjects (SOS) study (A) and Genetics of Diabetes Audit and Research Tayside Scotland (Go-DARTS) study (B) participants (n = 3,473 and 15,448 respectively).

https://doi.org/10.1371/journal.pone.0039362.g001

The hypothesis of an association between PNPLA3 I148M genotype and serum triglycerides, insulin resistance (assessed by homeostasis model assessment for insulin resistance, HOMA-IR) and type 2 diabetes was tested in obese participants of the SOS study before any treatment (baseline, Figure 1A). Subsequently, the associations found at baseline were investigated at 2- and 10-year follow up by analyzing separately the control (individuals who did not have substantial changes in body weight), and the surgery group (individuals who underwent bariatric surgery) to test if the associations were preserved at different time points and if they were affected by severe weight loss.

Next, the Go-DARTS study (Figure 1B) participants were examined to confirm the SOS study findings in an independent study cohort. Specifically, the PNPLA3 I148M association with serum triglycerides and type 2 diabetes risk in severely obese individuals (Body-mass index ≥35) was tested.

Ethic statement

Written Informed consent has been obtained by all study participants. All clinical investigations have been conducted according to the principles expressed in the Declaration of Helsinki. The SOS study protocol was approved the following Swedish ethics committees: Regional Institutional Review Board of Gothenburg University, Regional Institutional Review Board of Linköping University, Regional Institutional Review Board of Lund University, Regional Institutional Review Board of Karolinska Intitute, in Stockholm, Regional Institutional Review Board of Umeå University, Regional Institutional Review Board of Örebro University and Regional Institutional Review Board of Uppsala University. SOS trial has been registered in the ClinicalTrials.gov registry (NCT01479452, http://clinicaltrials.gov/ct2/show/NCT01479452?term  =  NCT01479452&rank  = 1). The Go-DARTS study has been approved by the Tayside Medical ethics Committee at the Ninewells Hospital and Medical School, Dundee, UK.

The SOS study

The SOS study has been previously described [18], [19]. Briefly, the study enrolled 4,047 obese individuals in Sweden between September 1987 and January 2001. A total of 2,010 individuals constituted the bariatric surgery arm and a matched control group of 2,037 individuals was enrolled based on 18 matching variables. Inclusion and exclusion criteria have been presented previously [18]. A total of 3,585 DNA samples from SOS study participants (SOS Version 1.0) was available for analyses. PNPLA3 was successfully genotyped in 3,473 individuals: 1,719 from the control group and 1,754 from the surgery group (success rate: 97%).

The surgery patients and the conventionally treated controls both started the study 4 weeks before the date of bariatric surgery (termed baseline). Biochemical parameters were measured at baseline and after 2 and 10 years in both the control and the surgery group. All blood samples were obtained in fasting conditions [18], [19]. Type 2 diabetes (n = 521 at baseline) was classified as fasting blood glucose ≥6.1 mmol/L (corresponding to fasting plasma glucose ≥7.0 mmol/L or 126 mg/dL) and/or therapy with glucose-lowering medications. At each examination, weight, height and blood pressure were measured.

The Go-DARTS study

The Go-DARTS study is a Wellcome Trust-sponsored ongoing cohort study examining the genetic factors that contribute to type 2 diabetes and related conditions [20]–[22]. Recruitment was performed in the Tayside region of Scotland. All participants are of European ancestry. Between 1997 and 2009, Go-DARTS recruited 9,000 individuals with type 2 diabetes and 8,187 individuals without type 2 diabetes. PNPLA3 I148M genotyping was available for a total of 7,691 and 7,757 individuals in the case and control group, respectively.

Retrospective and follow up data are ascertained from the Tayside region electronic medical record. Biochemical parameters were measured at the study visit. ALT levels were available from routine clinical database.

Genotyping of the PNPLA3 I148M genetic variant (rs738409)

Fluorogenic 5′-nucleotidase (Taqman) allelic discrimination assay for the PNPLA3 rs738409 sequence variant was used to genotype all participants, as previously described [12], [16], with a success rate of 97%. To ensure for the quality and consistency of genotyping, three positive controls (one for each PNPLA3 rs738409 I148M genotype) which were previously genotyped by direct Sanger sequencing were included in each plate of the allelic discrimination assay. Furthermore, each DNA sample was run in triplicate and genotype calls with a 100% consistency were included in the analyses.

Statistical analyses

Variables are shown as means ± standard deviations. HOMA-IR was calculated only in non-diabetic participants as: [fasting insulin mIU/L x fasting glucose mmol/L]/22.5 [23]. Genotype and allele as well as categorical variable distributions across the genotype classes were compared using χ2 test. Linear regression analysis was used to assess the effect of the three different PNPLA3 genotypes on continuous dependent variables and P values were adjusted for age, gender and body-mass index (BMI). P values were similar after adjustment of triglycerides for HOMA-IR and ALT or after adjustment of HOMA-IR for ALT and triglycerides. Changes at 2- and 10- year follow up have been calculated as follows: [(follow up value – baseline value)/ baseline value] ×100. Binary logistic regression analysis was used to calculate allelic odds ratio and 95% confidence interval for insulin resistance, type 2 diabetes and combined risk; age, gender, and BMI were included in the model as covariates.

Non-normally distributed variables were log transformed before entering the statistical models. Statistical analyses were carried out using the Statistical Package for Social Sciences (SPSS, version 18.0.0, Inc. Chicago, IL, USA). P values <0.05 were considered significant.

Results

PNPLA3 148M allele is associated with metabolic traits in obese subjects in the SOS study before bariatric surgery

To assess the effect of the PNPLA3 I148M genotype on serum lipids and insulin resistance in obese subjects, DNA from SOS participants was genotyped (Figure 1A). Clinical characteristics of participants at baseline (n = 3,473; combined control and surgery groups) are shown in Table S1. Genotype and allele frequencies in these subjects were in Hardy-Weinberg equilibrium (Table S2) and they were not different from those previously observed in Europeans [24]. The clinical characteristics at baseline stratified for PNPLA3 I148M genotypes are presented in Table 1. All participants were severely obese at baseline (mean BMI 41±5, Table S1). As a positive control for the genetic analyses ALT were examined. Increased serum levels of ALT were found in carriers of the 148M allele at baseline (P<0.001; Table 1). Triglycerides were found to be lower (P<0.001; Table 1) and HOMA-IR (P = 0.004; Table 1) to be higher in the PNPLA3 148M carriers at baseline.

[Figure omitted. See PDF.]

Table 1. Baseline Clinical Characteristics of SOS Study Participants Stratified by PNPLA3 I148M Genotype.

https://doi.org/10.1371/journal.pone.0039362.t001

Next the risk of being insulin resistant was examined in non diabetic individuals at baseline. The individuals were divided into insulin resistant and sensitive according to the median of the HOMA-IR (in men 3.9 and in women 2.9 U). The PNPLA3 148M variant was associated with a higher risk of being insulin resistant (Odds ratio, O.R.: 1.10, 95% Confidence interval, C.I. 1.02-1.35, P = 0.038; Panel A, Table S3). Next, the risk of type 2 diabetes was examined and the PNPLA3 148M allele was also associated with a significantly higher risk of type 2 diabetes at baseline after adjusting for age, gender, and BMI (O.R. 1.09, 95% C.I. 1.01-1.39, P = 0.040, Panel B, Table S3). The combined risk of developing insulin resistance or diabetes was consistently increased (O.R. 1.11, 95% C.I. 1.03-1.37, P = 0.021, Panel C, Table S3).

Associations between PNPLA3 148M allele and metabolic traits are abolished by surgically induced weight loss

The clinical characteristics at 2- and 10- year follow up in the control and the surgery group separately, stratified for PNPLA3 I148M genotypes are shown in table S4 and S5, respectively. Weight changes in the surgery and the control group from baseline to 10-year follow up have been previously described [19]. Briefly, there was almost no weight loss in the control group at the 2- and 10-year follow up, whereas the mean BMI in the surgery group was reduced by 24% at 2- year and 17% at 10- year follow up. Serum triglycerides were lower in PNPLA3 148M carriers in the control at 2- and 10-year follow up (P<0.001 in both, Figure 2A) but not in the surgery group (Figure 2B). HOMA-IR was greater in carriers of the 148M allele in the control group at the 2-year (P = 0.001) and 10-year (P = 0.045) follow up (Figure 2C), but no differences were found across PNPLA3 genotypes in the surgery group (Figure 2D). In the 2- and 10- year follow up there was no difference in the PNPLA3 I148M genotype distribution in participants with and without type 2 diabetes in both the control and surgery group (Table S4 and S5).

[Figure omitted. See PDF.]

Figure 2. Serum triglycerides and homeostasis model assessment for insulin resistance (HOMA-IR) at 2- and 10-year follow-up.

Serum triglycerides (panel A and B) levels and HOMA-IR (panel C and D) at 2- and 10- year follow up (FU) in the control (A, C) and the surgery (B, D) group across the PNPLA3 genotypes. Values are means and standard deviations. * HOMA-IR values are shown only in non-diabetic individuals.

https://doi.org/10.1371/journal.pone.0039362.g002

As a positive control for the genetic analyses, ALT levels were examined. Increased serum levels of ALT were found in carriers of the 148M allele in the control group at the 2-year (P<0.001) and 10-year (P = 0.004) follow up (Table S4). In the surgery group no differences in ALT levels were found across PNPLA3 genotypes (Table S5) at 2- year follow up whereas at 10-year follow up ALT increased across the PNPLA3 genotypes (P = 0.002; table S5).

Moreover, 2- and 10- year changes in metabolic parameters were tested across the genotype classes. In the surgery group, the PNPLA3 148M allele was associated with a higher reduction in HOMA-IR and ALT levels and with a lower reduction in triglyceride levels at 2- and 10-year follow up (Table S6). No associations were observed in the individuals without surgery.

Association of PNPLA3 148M Allele with Metabolic Traits in the Go-DARTS Study

The effect of the PNPLA3 148M allele on metabolic traits was also examined in the Go-DARTS study, a large case (n = 7,691) control (n = 7,757) study of type 2 diabetes based in the Tayside region of Scotland (for clinical characteristics Table S7). Genotype and allele frequencies were in Hardy-Weinberg Equilibrium (Table S8). Lower levels of fasting serum triglycerides coupled with an increase in ALT levels (both P<0.001) were observed in 148M diabetic carriers but not in the control group (Table 2). Subsequently, the Go-DARTS diabetic participants with BMI ≥35 (n = 1,842) and <35 (n = 5,849) were analyzed separately. The association with ALT levels was present in both groups (data not shown), but the association with serum triglycerides was only found in the severely obese (BMI≥35) group (II: 225±131, IM: 208±132, MM: 191±131 mg/dL, P = 0.003) and not in those with BMI<35 (II: 190±122, IM: 184±117, MM: 184±120 mg/dL, P = 0.219), although the interaction term did not reach statistical significance (P = 0.18). After stratifying the Go-DARTS control group by BMI, no associations between serum triglyceride levels and PNPLA3 were found in individuals either with BMI≥35 or <35. However, it should be pointed out that the severely obese group without type 2 diabetes was small (n = 456).

[Figure omitted. See PDF.]

Table 2. Clinical Characteristics of Go-DARTS Study Participants Stratified by PNPLA3 I148M Genotype.

https://doi.org/10.1371/journal.pone.0039362.t002

Moreover, the PNPLA3 148M allele was associated with a marginally increased risk of developing type 2 diabetes in the overall study (O.R. 1.04, 95% C.I. 0.98-1.11, P = 0.104, table 3A). After stratifying the cohort according to BMI, the 148M allele was found to confer a significantly increased risk in those with BMI ≥35 (n = 1,842 in the type 2 diabetes group and n = 456 in the control group; O.R. 1.37, 95% C.I. 1.13-1.66, P = 0.001, table 3B) but not in those with BMI <35 (n = 5,849 in the type 2 diabetes group and n = 7,301 in the control group; O.R. 1.01, 95% C.I. 0.94-1.07, P = 0.855, table 3B). A highly significant statistical interaction term between PNPLA3 genotype and severe obesity status was observed in determining risk of type 2 diabetes (P = 0.002).

[Figure omitted. See PDF.]

Table 3. Multivariate regression analysis of Risk of Type 2 Diabetes Mellitus in the Go-DARTS study.

https://doi.org/10.1371/journal.pone.0039362.t003

Discussion

This is the first study showing that the PNPLA3 148M allele is associated with insulin resistance and increased type 2 diabetes risk specifically in severely obese individuals despite relatively lower serum triglycerides. The results of this study are supported by extensive prior research that has established the robust gene/environment interaction between the PNPLA3 148M allele and obesity in determining individual susceptibility to hepatocyte lipid accumulation and damage [10], [12], [16], [24]. This interaction is clearly evident in our study, as the obese PNPLA3 148M carriers in the SOS study had higher ALT levels before, but not after, weight loss caused by bariatric surgery. We have now shown that the association of the PNPLA3 148M variant with other phenotypes is also dependent on increased body mass per se or on factors related to it.

Lower serum triglycerides and higher insulin resistance in 148M carriers were found before but not after surgery. This difference in insulin resistance was seen in the individuals without diabetes but was also reflected in an observed 15% increased risk for type 2 diabetes in the pre-surgery SOS study participants for each PNPLA3 148M allele. The absence of the association with type 2 diabetes in the SOS control group at 2- and 10- year follow up is likely due to a lack of statistical power given the reduced sample size.

Furthermore, we examined if PNPLA3 I148M variant affects changes in insulin resistance, triglyceride and transaminase levels. In the surgery group, the PNPLA3 148M allele was associated with a higher reduction in HOMA-IR and ALT levels and with a lower reduction in triglyceride levels at 2- and 10-year follow up. No associations were present in the control group. This result reflects the association of PNPLA3 with HOMA-IR, triglyceride and ALT levels before but not after weight loss.

Similar observations were made in the Go-DARTS study, where the association with triglycerides was dependant on obesity, with the lean group showing no association. The association with diabetes was also consistent between the two studies. While the overall increased risk of the PNPLA3 148M allele for type 2 diabetes in the Go-DARTS study was only 4%, in the severely obese a 37% increased risk per 148M allele for type 2 diabetes was found. Again, no association with type 2 diabetes was seen in those with a lower BMI.

The mean BMI of the Go-DARTS diabetic cohort is similar to that of the SOS surgery group in the follow up. However, no effect of the PNPLA3 genotype on serum parameters was observed in the SOS surgery group. Surgery may lead to a great improvement of the metabolic profile due to profound reductions in calorific intake and thus blunt the effect of the PNPLA3 genotype. Indeed, despite the similar BMI, the SOS surgery group shows an overall healthier metabolic profile compared to the Go-DARTS diabetic cohort.

The gene/obesity interaction, where the association is only seen in a fairly extreme subgroup, may explain why the PNPLA3 148M allele has not been widely associated with serum triglyceride levels, insulin resistance and type 2 diabetes in other studies [6], [24], [25]. In many population-based or healthy volunteer studies, the prevalence of severe obesity would be rather limited. Moreover, a low statistical power may explain the lack of significant associations in previous studies performed on obese subjects [10], [12], [16]. To date, two well-powered genetic studies have shown an association of the PNPLA3 148M allele with lower fasting serum triglycerides; this association was specifically present in overweight/obese individuals and no association was found in the other non-obese cohorts [24], [26]. Furthermore, a recent genetic report showed an association of the 148M allele with insulin resistance in individuals of Asian descent with NAFLD where the authors raised the question regarding the possibility of an association of the 148M allele and diabetes mellitus [27], however this study was confounded by the fact that the cases were ascertained by hepatic steatosis status, therefore the causal relationships between the genotype and insulin resistance could not be established.

The cut off point of a BMI of 35 used in the analysis of the Go-DARTS cohort was selected to match the severe obesity in the SOS study, but further analysis showed that this was an optimal split point with the associations with serum triglycerides and type 2 diabetes being limited to individuals with severe obesity. In fact, none of these associations were apparent in the overweight or moderately obese groups (25.1–34.9, data not shown), although the association between 148M and ALT levels was only absent in individuals with a BMI <25 in the Go-DARTS cohort, presumably reflecting the more pronounced gene environment interaction of this variant and hepatocyte damage. This notion is supported by the fact that only the association with ALT levels was apparent upon the 10 year weight rebound in the SOS surgery group.

The novelty of our analysis resides not only in the PNPLA3 148M allele being associated with type 2 diabetes but also simultaneously with lower serum triglycerides and increased insulin resistance. Furthermore, it supports a recent study in which serum triglyceride raising alleles were not found associated with increased risk for type 2 diabetes [28]. In line with this, the GCKR gene sequence variant (rs780094) has been found consistently associated with increased serum triglycerides, decreased insulin resistance and protection from type 2 diabetes [13], [29], [30]. Moreover, results of this study also suggest the possibility of finding new genetic loci associated with type 2 diabetes, insulin resistance or fasting triglycerides levels uncovered by obesity which would serve as a stressor acting on the genetic background. In vivo genetically modified murine models often reveal glucose or lipid metabolism changes after obesity is induced by genetic manipulation or high fat diet [31].

Although no effect on lipid or glucose metabolism was found in the Pnpla3 knock out mouse model [32], [33], a recent study of mice overexpressing hepatic Pnpla3 demonstrated marked changes in serum triglycerides and glucose [15]. Fasting serum triglycerides are a reflection of the basal hepatic VLDL efflux as opposed to postprandial serum triglycerides which consist mainly of chylomicrons [34]. PNPLA3 is a membrane bound protein highly expressed in hepatocytes [35]. It is tightly associated with the endoplasmic reticulum and lipid droplets [36] where nascent VLDL are enriched in triglycerides. This protein shows an in vitro hydrolytic lipase activity [36] on triglycerides which may be responsible for the nascent VLDL lipidation. The PNPLA3 148M mutant protein fails to hydrolyze triglycerides in vitro [36]. The loss of function mutation could account for a reduction of the VLDL lipidation and a subsequent impaired hepatic triglycerides efflux resulting in lower serum triglycerides and liver fat retention. Under these circumstances liver fat accumulation may lead to insulin resistance and possibly to type 2 diabetes in obese individuals. Further physiological studies are warranted to understand the exact role of PNPLA3 in mediating the observed phenotypes.

A strength of this study is the evaluation of two large, well-powered and phenotyped European study cohorts and the use of a surgical intervention to prospectively modify the gene/obesity interaction within individuals. A limitation of this study is the lack of imaging data regarding the accumulation of lipids in the livers of the obese subjects; however the role of PNPLA3 148M in modulating hepatic triglyceride accumulation and the corresponding relationships with altered ALT values is well established by many previous studies. Another limitation of this study is the use of HOMA-IR to assess insulin resistance; studies on selected severe obese individuals using euglycaemic clamp to measure insulin resistance are needed to confirm the association of the 148M allele and insulin resistance.

In conclusion, this study demonstrates a replicated gene by obesity interaction in determining type 2 diabetes risk, and provides insight into the causal role of obesity-driven hepatic lipid accumulation and type 2 diabetes.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

Clinical Characteristics of SOS Study Participants at Baseline.

https://doi.org/10.1371/journal.pone.0039362.s001

(DOC)

Table S2.

Genotype and Allele Frequencies of the PNPLA3 I148M Sequence Variant (rs738409) in the SOS Study Participants.

https://doi.org/10.1371/journal.pone.0039362.s002

(DOC)

Table S3.

Multivariate regression analysis of insulin resistance, type 2 diabetes and combined risk in the SOS study.

https://doi.org/10.1371/journal.pone.0039362.s003

(DOC)

Table S4.

Clinical Characteristics of SOS Study Control Group Stratified by PNPLA3 I148M Genotype at 2- and 10-Year Follow Up.

https://doi.org/10.1371/journal.pone.0039362.s004

(DOC)

Table S5.

Clinical Characteristics of SOS Study Surgery Group Stratified by PNPLA3 I148M Genotype at 2- and 10-Year Follow Up.

https://doi.org/10.1371/journal.pone.0039362.s005

(DOC)

Table S6.

Two- and ten-year changes in HOMA-IR, triglyceride and ALT values in the surgery and the control group from the SOS study.

https://doi.org/10.1371/journal.pone.0039362.s006

(DOC)

Table S7.

Clinical Characteristics of Go-DARTS Study Participants.

https://doi.org/10.1371/journal.pone.0039362.s007

(DOC)

Table S8.

Genotype and Allele Frequencies of the PNPLA3 I148M Sequence Variant in the Type 2 Diabetes and Control Groups from the Go-DARTS Study.

https://doi.org/10.1371/journal.pone.0039362.s008

(DOC)

Acknowledgments

We acknowledge the support of the Health Informatics Centre, University of Dundee for managing and supplying the anonymised data and NHS Tayside, the original data owner. We are grateful to all the participants who took part in the Go-DARTS study, to the general practitioners, to the Scottish School of Primary Care for their help in recruiting the participants, and to the whole team, which includes interviewers, computer and laboratory technicians, clerical workers, research scientists, volunteers, managers, receptionists, and nurses.

Author Contributions

Conceived and designed the experiments: CNAP SR CM. Performed the experiments: CM CP MAB. Analyzed the data: CNAP SR CM CP MAB. Contributed reagents/materials/analysis tools: CNAP LMSC JB SR. Wrote the paper: CNAP SR CM CP MAB MA LB LAD HC ASD JFD ERP MM ATH TF ADM MP PAS PJ JB LS LMSC.

Citation: Palmer CNA, Maglio C, Pirazzi C, Burza MA, Adiels M, Burch L, et al. (2012) Paradoxical Lower Serum Triglyceride Levels and Higher Type 2 Diabetes Mellitus Susceptibility in Obese Individuals with the PNPLA3 148M Variant. PLoS ONE7(6): e39362. https://doi.org/10.1371/journal.pone.0039362

References

1. Kotronen A, Seppala-Lindroos A, Bergholm R, Yki-Jarvinen H (2008) Tissue specificity of insulin resistance in humans: fat in the liver rather than muscle is associated with features of the metabolic syndrome. Diabetologia 51: 130–138.A. KotronenA. Seppala-LindroosR. BergholmH. Yki-Jarvinen2008Tissue specificity of insulin resistance in humans: fat in the liver rather than muscle is associated with features of the metabolic syndrome.Diabetologia51130138

2. Kotronen A, Juurinen L, Tiikkainen M, Vehkavaara S, Yki-Jarvinen H (2008) Increased liver fat, impaired insulin clearance, and hepatic and adipose tissue insulin resistance in type 2 diabetes. Gastroenterology 135: 122–130.A. KotronenL. JuurinenM. TiikkainenS. VehkavaaraH. Yki-Jarvinen2008Increased liver fat, impaired insulin clearance, and hepatic and adipose tissue insulin resistance in type 2 diabetes.Gastroenterology135122130

3. Lewis GF, Carpentier A, Adeli K, Giacca A (2002) Disordered fat storage and mobilization in the pathogenesis of insulin resistance and type 2 diabetes. Endocr Rev 23: 201–229.GF LewisA. CarpentierK. AdeliA. Giacca2002Disordered fat storage and mobilization in the pathogenesis of insulin resistance and type 2 diabetes.Endocr Rev23201229

4. Adiels M, Westerbacka J, Soro-Paavonen A, Hakkinen AM, Vehkavaara S, et al. (2007) Acute suppression of VLDL1 secretion rate by insulin is associated with hepatic fat content and insulin resistance. Diabetologia 50: 2356–2365.M. AdielsJ. WesterbackaA. Soro-PaavonenAM HakkinenS. Vehkavaara2007Acute suppression of VLDL1 secretion rate by insulin is associated with hepatic fat content and insulin resistance.Diabetologia5023562365

5. Adiels M, Olofsson SO, Taskinen MR, Boren J (2008) Overproduction of very low-density lipoproteins is the hallmark of the dyslipidemia in the metabolic syndrome. Arterioscler Thromb Vasc Biol 28: 1225–1236.M. AdielsSO OlofssonMR TaskinenJ. Boren2008Overproduction of very low-density lipoproteins is the hallmark of the dyslipidemia in the metabolic syndrome.Arterioscler Thromb Vasc Biol2812251236

6. Romeo S, Kozlitina J, Xing C, Pertsemlidis A, Cox D, et al. (2008) Genetic variation in PNPLA3 confers susceptibility to nonalcoholic fatty liver disease. Nat Genet 40: 1461–1465.S. RomeoJ. KozlitinaC. XingA. PertsemlidisD. Cox2008Genetic variation in PNPLA3 confers susceptibility to nonalcoholic fatty liver disease.Nat Genet4014611465

7. Browning JD, Cohen JC, Hobbs HH (2010) Patatin-like phospholipase domain-containing 3 and the pathogenesis and progression of pediatric nonalcoholic fatty liver disease. Hepatology 52: 1189–1192.JD BrowningJC CohenHH Hobbs2010Patatin-like phospholipase domain-containing 3 and the pathogenesis and progression of pediatric nonalcoholic fatty liver disease.Hepatology5211891192

8. Kantartzis K, Peter A, Machicao F, Machann J, Wagner S, et al. (2009) Dissociation between fatty liver and insulin resistance in humans carrying a variant of the patatin-like phospholipase 3 gene. Diabetes 58: 2616–2623.K. KantartzisA. PeterF. MachicaoJ. MachannS. Wagner2009Dissociation between fatty liver and insulin resistance in humans carrying a variant of the patatin-like phospholipase 3 gene.Diabetes5826162623

9. Kotronen A, Johansson LE, Johansson LM, Roos C, Westerbacka J, et al. (2009) A common variant in PNPLA3, which encodes adiponutrin, is associated with liver fat content in humans. Diabetologia 52: 1056–1060.A. KotronenLE JohanssonLM JohanssonC. RoosJ. Westerbacka2009A common variant in PNPLA3, which encodes adiponutrin, is associated with liver fat content in humans.Diabetologia5210561060

10. Valenti L, Alisi A, Galmozzi E, Bartuli A, Del Menico B, et al. (2010) I148M patatin-like phospholipase domain-containing 3 gene variant and severity of pediatric nonalcoholic fatty liver disease. Hepatology 52: 1274–1280.L. ValentiA. AlisiE. GalmozziA. BartuliB. Del Menico2010I148M patatin-like phospholipase domain-containing 3 gene variant and severity of pediatric nonalcoholic fatty liver disease.Hepatology5212741280

11. Goran MI, Walker R, Le KA, Mahurkar S, Vikman S, et al. (2010) Effects of PNPLA3 on liver fat and metabolic profile in Hispanic children and adolescents. Diabetes 59: 3127–3130.MI GoranR. WalkerKA LeS. MahurkarS. Vikman2010Effects of PNPLA3 on liver fat and metabolic profile in Hispanic children and adolescents.Diabetes5931273130

12. Romeo S, Sentinelli F, Dash S, Yeo GS, Savage DB, et al. (2010) Morbid obesity exposes the association between PNPLA3 I148M (rs738409) and indices of hepatic injury in individuals of European descent. Int J Obes (Lond) 34: 190–194.S. RomeoF. SentinelliS. DashGS YeoDB Savage2010Morbid obesity exposes the association between PNPLA3 I148M (rs738409) and indices of hepatic injury in individuals of European descent.Int J Obes (Lond)34190194

13. Dupuis J, Langenberg C, Prokopenko I, Saxena R, Soranzo N, et al. (2010) New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk. Nat Genet 42: 105–116.J. DupuisC. LangenbergI. ProkopenkoR. SaxenaN. Soranzo2010New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk.Nat Genet42105116

14. Voight BF, Scott LJ, Steinthorsdottir V, Morris AP, Dina C, et al. (2010) Twelve type 2 diabetes susceptibility loci identified through large-scale association analysis. Nat Genet 42: 579–589.BF VoightLJ ScottV. SteinthorsdottirAP MorrisC. Dina2010Twelve type 2 diabetes susceptibility loci identified through large-scale association analysis.Nat Genet42579589

15. Qiao A, Liang J, Ke Y, Li C, Cui Y, et al. (2011) Mouse PNPLA3 influences systemic lipid and glucose homeostasis. Hepatology 54: 509–521.A. QiaoJ. LiangY. KeC. LiY. Cui2011Mouse PNPLA3 influences systemic lipid and glucose homeostasis.Hepatology54509521

16. Romeo S, Sentinelli F, Cambuli VM, Incani M, Congiu T, et al. (2010) The 148M allele of the PNPLA3 gene is associated with indices of liver damage early in life. J Hepatol 53: 335–338.S. RomeoF. SentinelliVM CambuliM. IncaniT. Congiu2010The 148M allele of the PNPLA3 gene is associated with indices of liver damage early in life.J Hepatol53335338

17. Miraglia Del Giudice E, Grandone A, Cirillo G, Santoro N, Amato A, et al. (2011) The Association of PNPLA3 Variants with Liver Enzymes in Childhood Obesity Is Driven by the Interaction with Abdominal Fat. PLoS One 6: e27933.E. Miraglia Del GiudiceA. GrandoneG. CirilloN. SantoroA. Amato2011The Association of PNPLA3 Variants with Liver Enzymes in Childhood Obesity Is Driven by the Interaction with Abdominal Fat.PLoS One6e27933

18. Sjostrom L, Larsson B, Backman L, Bengtsson C, Bouchard C, et al. (1992) Swedish obese subjects (SOS). Recruitment for an intervention study and a selected description of the obese state. Int J Obes Relat Metab Disord 16: 465–479.L. SjostromB. LarssonL. BackmanC. BengtssonC. Bouchard1992Swedish obese subjects (SOS). Recruitment for an intervention study and a selected description of the obese state.Int J Obes Relat Metab Disord16465479

19. Sjostrom L, Narbro K, Sjostrom CD, Karason K, Larsson B, et al. (2007) Effects of bariatric surgery on mortality in Swedish obese subjects. N Engl J Med 357: 741–752.L. SjostromK. NarbroCD SjostromK. KarasonB. Larsson2007Effects of bariatric surgery on mortality in Swedish obese subjects.N Engl J Med357741752

20. Doney AS, Lee S, Leese GP, Morris AD, Palmer CN (2005) Increased cardiovascular morbidity and mortality in type 2 diabetes is associated with the glutathione S transferase theta-null genotype: a Go-DARTS study. Circulation 111: 2927–2934.AS DoneyS. LeeGP LeeseAD MorrisCN Palmer2005Increased cardiovascular morbidity and mortality in type 2 diabetes is associated with the glutathione S transferase theta-null genotype: a Go-DARTS study.Circulation11129272934

21. Zhou K, Bellenguez C, Spencer CC, Bennett AJ, Coleman RL, et al. (2011) Common variants near ATM are associated with glycemic response to metformin in type 2 diabetes. Nat Genet 43: 117–120.K. ZhouC. BellenguezCC SpencerAJ BennettRL Coleman2011Common variants near ATM are associated with glycemic response to metformin in type 2 diabetes.Nat Genet43117120

22. Palmer CN, Kimber CH, Doney AS, Proia AS, Morris AD, et al. (2010) Combined effect of inflammatory gene polymorphisms and the risk of ischemic stroke in a prospective cohort of subjects with type 2 diabetes: a Go-DARTS study. Diabetes 59: 2945–2948.CN PalmerCH KimberAS DoneyAS ProiaAD Morris2010Combined effect of inflammatory gene polymorphisms and the risk of ischemic stroke in a prospective cohort of subjects with type 2 diabetes: a Go-DARTS study.Diabetes5929452948

23. Matthews DR, Hosker JP, Rudenski AS, Naylor BA, Treacher DF, et al. (1985) Homeostasis model assessment: insulin resistance and beta-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia 28: 412–419.DR MatthewsJP HoskerAS RudenskiBA NaylorDF Treacher1985Homeostasis model assessment: insulin resistance and beta-cell function from fasting plasma glucose and insulin concentrations in man.Diabetologia28412419

24. Kollerits B, Coassin S, Beckmann ND, Teumer A, Kiechl S, et al. (2009) Genetic evidence for a role of adiponutrin in the metabolism of apolipoprotein B-containing lipoproteins. Hum Mol Genet 18: 4669–4676.B. KolleritsS. CoassinND BeckmannA. TeumerS. Kiechl2009Genetic evidence for a role of adiponutrin in the metabolism of apolipoprotein B-containing lipoproteins.Hum Mol Genet1846694676

25. Speliotes EK, Butler JL, Palmer CD, Voight BF, Hirschhorn JN, et al. (2010) PNPLA3 variants specifically confer increased risk for histologic nonalcoholic fatty liver disease but not metabolic disease. Hepatology 52: 904–912.EK SpeliotesJL ButlerCD PalmerBF VoightJN Hirschhorn2010PNPLA3 variants specifically confer increased risk for histologic nonalcoholic fatty liver disease but not metabolic disease.Hepatology52904912

26. Krawczyk M, Gruenhage F, Mahler M, Tirziu S, Acalovschi M, et al. (2011) The common adiponutrin variant I148M does not confer gallstone risk but affects fasting glucose and triglyceride levels. J Physiol Pharmacol 62: 369–375.M. KrawczykF. GruenhageM. MahlerS. TirziuM. Acalovschi2011The common adiponutrin variant I148M does not confer gallstone risk but affects fasting glucose and triglyceride levels.J Physiol Pharmacol62369375

27. Wang CW, Lin HY, Shin SJ, Yu ML, Lin ZY, et al. (2011) The PNPLA3 I148M polymorphism is associated with insulin resistance and nonalcoholic fatty liver disease in a normoglycaemic population. Liver Int 31: 1326–1331.CW WangHY LinSJ ShinML YuZY Lin2011The PNPLA3 I148M polymorphism is associated with insulin resistance and nonalcoholic fatty liver disease in a normoglycaemic population.Liver Int3113261331

28. De Silva NM, Freathy RM, Palmer TM, Donnelly LA, Luan J, et al. (2011) Mendelian randomization studies do not support a role for raised circulating triglyceride levels influencing type 2 diabetes, glucose levels, or insulin resistance. Diabetes 60: 1008–1018.NM De SilvaRM FreathyTM PalmerLA DonnellyJ. Luan2011Mendelian randomization studies do not support a role for raised circulating triglyceride levels influencing type 2 diabetes, glucose levels, or insulin resistance.Diabetes6010081018

29. Saxena R, Voight BF, Lyssenko V, Burtt NP, de Bakker PI, et al. (2007) Genome-wide association analysis identifies loci for type 2 diabetes and triglyceride levels. Science 316: 1331–1336.R. SaxenaBF VoightV. LyssenkoNP BurttPI de Bakker2007Genome-wide association analysis identifies loci for type 2 diabetes and triglyceride levels.Science31613311336

30. Sparsø T, Andersen G, Nielsen T, Burgdorf KS, Gjesing AP, et al. (2008) The GCKR rs780094 polymorphism is associated with elevated fasting serum triacylglycerol, reduced fasting and OGTT-related insulinaemia, and reduced risk of type 2 diabetes. Diabetologia 51: 70–75.T. SparsøG. AndersenT. NielsenKS BurgdorfAP Gjesing2008The GCKR rs780094 polymorphism is associated with elevated fasting serum triacylglycerol, reduced fasting and OGTT-related insulinaemia, and reduced risk of type 2 diabetes.Diabetologia517075

31. Postic C, Girard J (2008) Contribution of de novo fatty acid synthesis to hepatic steatosis and insulin resistance: lessons from genetically engineered mice. J Clin Invest 118: 829–838.C. PosticJ. Girard2008Contribution of de novo fatty acid synthesis to hepatic steatosis and insulin resistance: lessons from genetically engineered mice.J Clin Invest118829838

32. Chen W, Chang B, Li L, Chan L (2010) Patatin-like phospholipase domain-containing 3/adiponutrin deficiency in mice is not associated with fatty liver disease. Hepatology 52: 1134–1142.W. ChenB. ChangL. LiL. Chan2010Patatin-like phospholipase domain-containing 3/adiponutrin deficiency in mice is not associated with fatty liver disease.Hepatology5211341142

33. Basantani MK, Sitnick MT, Cai L, Brenner DS, Gardner NP, et al. (2011) Pnpla3/Adiponutrin deficiency in mice does not contribute to fatty liver disease or metabolic syndrome. J Lipid Res 52: 318–329.MK BasantaniMT SitnickL. CaiDS BrennerNP Gardner2011Pnpla3/Adiponutrin deficiency in mice does not contribute to fatty liver disease or metabolic syndrome.J Lipid Res52318329

34. Gotto AM Jr, Pownall HJ, Havel RJ (1986) Introduction to the plasma lipoproteins. Methods Enzymol 128: 3–41.AM Gotto JrHJ PownallRJ Havel1986Introduction to the plasma lipoproteins.Methods Enzymol128341

35. Huang Y, He S, Li JZ, Seo YK, Osborne TF, et al. (2010) A feed-forward loop amplifies nutritional regulation of PNPLA3. Proc Natl Acad Sci U S A 107: 7892–7897.Y. HuangS. HeJZ LiYK SeoTF Osborne2010A feed-forward loop amplifies nutritional regulation of PNPLA3.Proc Natl Acad Sci U S A10778927897

36. He S, McPhaul C, Li JZ, Garuti R, Kinch L, et al. (2010) A sequence variation (I148M) in PNPLA3 associated with nonalcoholic fatty liver disease disrupts triglyceride hydrolysis. J Biol Chem 285: 6706–6715.S. HeC. McPhaulJZ LiR. GarutiL. Kinch2010A sequence variation (I148M) in PNPLA3 associated with nonalcoholic fatty liver disease disrupts triglyceride hydrolysis.J Biol Chem28567066715