- Research Paper
- Open Access
The PNPLA3 Ile148Met interacts with overweight and dietary intakes on fasting triglyceride levels
Genes & Nutritionvolume 9, Article number: 388 (2014)
The Erratum to this article has been published in Genes & Nutrition 2016 11:11
The Ile148Met (rs738409, G-allele) in the patatin-like phospholipase domain-containing protein 3 gene (PNPLA3) associates with liver fat content and may lead to loss-of-function (hydrolysis) or gain-of-function (CoA-dependent lysophosphatidic acid acyltransferase) defects. PNPLA3 is up-regulated by dietary carbohydrates, and interactions between rs738409 and carbohydrates, and sugar and ω6:ω3-polyunsaturated fatty acid (PUFA) ratio on hepatic fat accumulation have been reported. We examined interaction between rs738409 and overweight, and between rs738409 and dietary intakes (carbohydrates, sucrose and ω6:ω3-PUFA ratio), on fasting triglyceride levels. From the Malmo Diet and Cancer Study-Cardiovascular Cohort, 4,827 individuals without diabetes aged 58 ± 6 years, 2,346 with BMI ≤ 25 kg/m2 and 2,478 with BMI > 25 kg/m2, were included in cross-sectional analyses. Dietary data were collected by a modified diet history method. Overweight modified the association between rs738409 and fasting triglyceride levels (P interaction = 0.003). G-allele associated with lower triglycerides only among overweight individuals (P = 0.01). Nominally, significant interaction on triglyceride levels was observed between rs738409 and sucrose among normal-weight individuals (P interaction = 0.03). G-allele associated with lower triglycerides among overweight individuals in the lowest tertiles of carbohydrate and ω6:ω3-PUFA ratio (P = 0.04 and P = 0.001) and with higher triglycerides among normal-weight individuals in the highest tertile of sucrose (P = 0.001). We conclude that overweight and dietary sucrose may modify the association between rs738409 and fasting triglyceride levels.
A common missense variation Ile148Met (rs738409, G-allele), in the patatin-like phospholipase domain-containing protein 3 gene (PNPLA3), also named adiponutrin, was identified associated with liver fat content and non-alcoholic fatty liver disease (NAFLD) in a genome-wide association study (GWAS) in 2008 (Romeo et al. 2008). This association between the 148M (G-allele) and fatty liver and histological disease severity has since then been widely replicated in both adults and children (Kotronen et al. 2009; Kantartzis et al. 2009; Goran et al. 2010; Hyysalo et al. 2011; Santoro et al. 2010; Valenti et al. 2010; Del Giudice et al. 2011; Romeo et al. 2010a, b). Two recent studies observed that the PNPLA3 rs738409 GG-genotype associated with lower serum triglyceride and cholesterol levels among obese (Palmer et al. 2012) or glucose-intolerant individuals (Krarup et al. 2012) and was suggested to be a result of an obesity-driven reduced triglyceride hydrolysis and intrahepatic fat accumulation (Pirazzi et al. 2012). Indeed, PNPLA3 has been shown to have a hydrolytic activity against the three major glycerolipids (Huang et al. 2011; Pignitore et al. 2013), and the PNPLA3 148M was also reported to be a loss-of-function (hydrolysis) variant with a markedly decreased V max for glycerolipids (Huang et al. 2011). Loss of the hydrolysis function of PNPLA3 by the 148M variant has been indicated in several studies: analysis of VLDL kinetics after a bolus infusion of stable isotopes in overweight/obese men revealed that the relative secretion of VLDL1 particles was lower, the rate of hydrolysis of [3H]—triglycerides during lipid depletion was decreased, and the cellular accumulation of labeled triglycerides was enhanced in the presence of excess free fatty acids (FFAs) by the 148M variant (Pirazzi et al. 2012; He et al. 2010; Perttila et al. 2012).
In addition to lipase activity, both human and mice PNPLA3 were recently demonstrated to exhibit an acyl-CoA-dependent lysophosphatidic acid acyltransferase (LPAAT) activity, which promotes cellular lipid synthesis by converting LPA to phosphatidic acid (Kumari et al. 2012). High sucrose diet up-regulates liver PNPLA3 of wild-type mice, which concomitantly exhibited increased LPAAT activity. In line with this, human PNPLA3-148M had increased LPAAT activity, and gain of this function was proposed to provide a plausible biochemical mechanism for increased hepatic fat accumulation (Kumari et al. 2012).
PNPLA3 is highly responsive to changes in energy balance. It is down-regulated by fasting and dramatically up-regulated by re-feeding (Baulande et al. 2001; Liu et al. 2004; Johansson et al. 2006; Hoekstra et al. 2010). Carbohydrate (CHO) feeding up-regulates PNPLA3 in murine adipose tissue, and liver and glucose stimulate PNPLA3 expression in human hepatocytes via carbohydrate response element-binding protein (ChREBP) (Perttila et al. 2012; Dubuquoy et al. 2011).
Dietary CHOs and especially sucrose intake has been proposed as a key player in liver fat deposition due to the high lipogenic potential of fructose, which stimulates hepatic de novo lipogenesis (Polson and Thompson 2003; Le et al. 2009; Maersk et al. 2012; Abdemalek et al. 2010; Sevastinova et al. 2012). In a study of 154 overweight Hispanic children, the hepatic fat content correlated positively with reported CHO and sugar intakes among children homozygous for the risk allele (GG-genotype), but not among C-allele carriers (Davis et al. 2010).
In addition to CHOs, PNPLA3 is up-regulated by certain fatty acids that inhibit its degradation (Huang et al. 2010). Moreover, dietary fat and in particular diets rich in ω-6 and low in ω-3 polyunsaturated fatty acids (PUFAs) have been suggested to contribute to metabolic defects in NAFLD (Cortez-Pinto et al. 2006). In line with this, a study in 127 children and adolescents observed a positive correlation between the ratio of ω-6 and ω-3 PUFA intakes and liver fat content among GG-genotype carriers (Santoro et al. 2012). Finally, the presence of the 148M allele in children treated with docosahexaenoic acid (22:6 ω-3) worsened the beneficial reduction induced by this fatty acid on hepatic steatosis (Nobili et al. 2013), and in obese adolescents, the oxidized fatty acids derived from linoleic acid associated with markers of liver injury only in PNPLA3 risk-genotype carriers (Santoro et al. 2014).
In this study, we examine interaction between PNPLA3 rs738409 and overweight on fasting blood triglyceride and serum alanine aminotransferase (ALT) levels and challenge the question if rs738409 interacts with dietary intakes of CHO, sucrose or ω-6:ω-3 PUFA ratio, on fasting triglyceride levels in adult Swedes participating the Malmö Diet and Cancer Study-Cardiovascular Cohort (MDCS-CC).
Subjects and methods
Study population and data collection
The Malmö Diet and Cancer Study (MDCS) is a Swedish population-based prospective cohort. During 1991–1996, all women born between 1923–1950 and men born between 1923–1945, living in the city of Malmö, were invited to participate (Berglund et al. 1993; Manjer et al. 2001). Limited Swedish language skills and mental incapacity were the only exclusion criteria. The participants filled out the questionnaires covering socioeconomic, lifestyle and dietary factors, recorded meals, and underwent a diet history interview. A total of 28,098 (40 % of the eligible persons) completed all baseline examinations (Manjer et al. 2001). From this cohort, 6,103 individuals with a baseline examination between 1991 and 1994 were randomly selected to participate in a cardiovascular cohort (MDC-CC), of whom 5,543 underwent blood sampling under standardized fasting conditions.
After an overnight fast, blood samples were drawn for the determination of serum lipids, serum insulin and whole blood glucose. Blood glucose was determined by a routine hexokinase method. Triglycerides and total cholesterol were determined on a DAX 48 automatic analyzer with use of reagents and calibrators from the supplier of the instrument (Bayer AB, Goteborg, Sweden). High-density lipoprotein (HDL) cholesterol was determined by the same procedure as used for total cholesterol but after precipitation of low-density lipoprotein (LDL) and very low-density lipoprotein (VLDL) with dextran–sulfate. LDL cholesterol was calculated from the values for triglycerides, total cholesterol and HDL cholesterol according to the Friedewald formula: LDL = total cholesterol – HDL − (triglycerides/2.2). Homeostasis model assessment for insulin resistance (HOMA-IR) was used as a measure of insulin resistance and was calculated with the following formula: insulin × blood glucose/22.5 (Nilsson et al. 2007).
From MDC-CC population (5,543), we excluded 213 individuals with prevalent type 2 diabetes, based on self-reported diabetes diagnosis or self-reported anti-diabetic regimen. After exclusion, we were left with 5,330 individuals, of whom we had information on diet, PNPLA3 rs738409 genotype, BMI and fasting triglycerides for 4,827 individuals.
Of all MDCS participants, a subgroup of 7,198 individuals participated in another study, the Malmö Preventive Project (MPP) (Berglund et al. 1996), in which the ALT measurements were performed using standard clinical assay. MPP is a prospective cohort of southern Sweden that was set up during 1974–1992 within the Department of Medicine at Malmö University Hospital, Sweden. ALT levels were used to analyze association with rs738409, but not for interaction with diet as the ALT measurements were from an earlier time as compared to the MDCS diet collection.
The ethical committee at Lund University approved MDCS and MPP, and all the participants have given their written informed consent.
Dietary intakes were collected at baseline using a combination of a 7-day menu book for recording of meals that vary from day to day (lunch and dinner meals), cold beverages and nutrient supplements, and a 168-item diet questionnaire (covering food regularly consumed during the past years) for assessment of consumption frequencies and portion sizes of regularly eaten foods that were not covered by the menu book. Finally, a 45-min interview completed the dietary assessment. The mean daily intake of foods was calculated based on frequency and portion sizes estimates from the questionnaire and menu book. The food intake was converted to energy and nutrient intakes using the MDC nutrient database where the majority of the nutrient information comes from PC-KOST2-93 from the National Food Administration in Uppsala, Sweden. The MDC method is described in detail elsewhere (Callmer et al. 1993; Wirfalt et al. 2001).
The routines for coding dietary data were slightly altered in September 1994 to shorten the interview time. This change did not reveal any major influence on the ranking of individuals (Callmer et al. 1993). The relative validity of the MDC method was evaluated in the Malmo Food study 1984–1985 (Wirfalt et al. 2001). The Pearson correlation coefficients between the reference method and the MDC method, adjusted for total energy, were in men and women 0.66/0.70 for CHO, 0.60/0.74 for the variable sugar, 0.23/0.68 for linoleic acid (ω-6), 0.55/0.44 for arachidonic acid (ω-6), 0.22/0.58 for α-linolenic acid (ω-3), 0.24/0.38 for eicosapentaenoic acid (ω-3), 0.37/0.40 for docosapentaenoic acid (ω-3), and 0.20/0.27 for docosahexaenoic acid (ω-3) (Riboli et al. 1997).
We used the following variables in this study: total energy (kcal) (including energy from fat, CHOs, protein, alcohol and fiber), CHO (E%), sucrose (E%), fat (E%), SFA (E%), MUFA (E%), PUFA (E%), P:S ratio, ω-3 PUFA, ω-6 PUFA, ω-6:ω-3 PUFA ratio, and fiber (g/1,000 kcal). The nutrient densities were calculated by dividing total nutrient intakes by non-alcohol energy intake. Tertiles were used as exposure categories of macronutrients. Individuals with potentially inaccurate reports of energy (N = 932) were identified as having a ratio of energy intake to the basal metabolic rate outside the 95 % CI limits of the physical activity level calculated for each individual as total energy expenditure. This procedure is described in detail elsewhere (Mattisson et al. 2005).
Other variables used as potential confounders
Information on age was obtained from the personal identification number. Body mass index (BMI; kg/m2) was calculated from direct measurement of weight and height. Leisure time physical activity was assessed by asking the participants to estimate the number of minutes per week they spent on 17 different activities. The duration was multiplied by an intensity factor to create a physical activity score that was divided into tertiles and categorized as low, medium and high. Participants were classified as current smokers, ex-smokers or never smokers. Alcohol intake was classified into four categories based on grams of alcohol consumed per day: zero-consumers, low (<15 g/day in women or <20 g/day in men), medium (15–30 g/day in women or 20–40 g/day in men) and high (>30 g/day in women or >40 g/day in men). The education variables were created by classifying participants according to their highest educational level (≤8y, 9–10y, 11–13y, and university degree). Season was defined as season of diet data collection: winter (December–February), spring (March–May), summer (June–August) and fall (September–November). Diet assessment method version was defined as data collection before or after the change of coding routines in September 1994. Missing values for the variables were treated as separate categories.
PNPLA3 rs738409 was genotyped using Taqman PCR method (Applied Biosystems, Foster City, CA USA) according to manufacturer’s instructions. ABI Prism Sequence Detection Systems ABI 7900HT (Applied Systems) was used for post-PCR allelic discrimination by measuring allele-specific fluorescence. Genotypes were successfully determined for 96.5 % of the individuals, and the minor allele frequency was 21 %. The genotype distribution did not deviate from Hardy–Weinberg equilibrium (P = 0.26), and concordance rate of repeated genotyping of all individuals was 99.9 %.
IBM SPSS Statistics version 20.0 was used for all statistical analyses. Statistical significance was set at P < 0.05. We examined baseline characteristics across the PNPLA3 genotypes for continuous variables using the general linear model (GLM) adjusted for age and sex when applicable. All continuous variables except age were Ln transformed to achieve normal distribution when testing for trend across genotypes. Associations between PNPLA3 genotype and serum triglyceride levels in tertiles of dietary intakes were evaluated using an additive model (genotypes coded as 0, 1 or 2 indicating the number of risk G-alleles) using GLM, adjusting for age and sex in a basic model and then with additional adjustment for BMI. Interactions between PNPLA3 genotype and tertiles of dietary intake on serum triglyceride levels were analyzed by introducing a multiplicative factor of genotype and diet tertiles (treated as continuous variables) in addition to genotype and diet to the models. The multivariate analyses were adjusted for age (continuous) and total energy intake (continuous) and the following categorical variables: sex, diet assessment method version, season, leisure time physical activity, smoking, alcohol intake and education. These covariates were identified from the literature or indicated potential confounding in the MDC cohort, due to their associations with dietary intakes. The multivariate analysis was also performed with additional adjustments for BMI (continuous) and finally also for fiber intakes. We used 25 kg/m2 as a limit for overweight following the statements by WHO. As earlier papers reported association between PNPLA3 I148M and lower triglyceride levels in obese or insulin resistant individuals, i.e., 30 kg/m2 in Palmer et al. (2012) and a mean BMI of 28.3 kg/m2 in Krarup et al. (2012), we performed additional interaction analyses between PNPLA3 genotype and severe overweight (BMI > 27.5 kg/m2) and obesity (BMI > 30 kg/m2) on triglycerides. Interaction analyses with dietary intakes were performed separately in normal-weight and overweight individuals defining normal weight as BMI ≤ 25 kg/m2 and overweight as BMI > 25 kg/m2).
Interaction between PNPLA3 rs738409 and obesity status on fasting triglyceride and ALT levels
The clinical characteristics and dietary intakes by the PNPLA3 rs738409 genotypes and by overweight status (normal weight defined as BMI ≤ 25 kg/m2 and overweight as BMI > 25 kg/m2) are presented in (Table 1). Among all individuals without stratification by overweight status, the rs738409 G-allele associated significantly only with higher ALT levels (P = 0.005). In contrast to results in normal-weight individuals, where none of the clinical variables associated significantly with rs738409, we observed associations between the G-allele and lower triglyceride-, VLDL-cholesterol- and LDL-cholesterol levels (P = 0.01, P = 0.002 and P = 0.04), and with higher ALT levels (P = 0.001) in overweight individuals. Overweight (25 kg/m2 ≥ BMI > 25 kg/m2) modified significantly the associations between rs738409 and fasting triglyceride and ALT levels (P interaction = 0.003 and 0.03, respectively). The interactions were similar with severe overweight (BMI > 27.5 kg/m2) and with obesity (BMI > 30 kg/m2) (P interaction = 0.004 and 0.028, for fasting triglyceride levels, and P interaction = 0.002 and 0.004 for ALT levels, respectively).
Total energy intake and intake levels of CHO, sucrose or fiber did not significantly differ between the different rs738409 genotype carriers. The intake levels of ω-3 PUFAs, ω-6 PUFAs and the ratio of them (ω-6:ω-3 PUFA) did not differ significantly among the different genotype carriers in the whole study group nor among overweight individuals. However, in the normal-weight individuals (BMI ≤ 25 kg/m2), G-allele associated with somewhat higher intake levels of fiber and ω-6:ω-3 PUFA ratio (Table 1).
Interaction analyses between PNPLA3 rs738409 and dietary intake levels of CHO, sucrose and ω-6:ω-3 PUFA on fasting triglyceride levels
The rs738409 G-allele associated with lower triglyceride levels in the lowest intake tertile of CHO (P = 0.02) in the whole study group, and this association was not seen among normal-weight individuals, but was strong among overweight individuals (Table 2). Among the overweight individuals, the association between the G-allele and lower triglycerides was only significant in the lowest intake tertiles of CHO and ω-6:ω-3 PUFA ratio (P = 0.04, P = 0.001), and similar tendency was observed in the lowest tertile of sucrose (P = 0.06) (Fig. 1). Among normal-weight individuals, the G-allele associated with higher triglycerides among those in the highest sucrose intake tertile (P = 0.001), whereas no associations were observed at the lower intakes. Statistical tests for interaction indicated a tendency for interaction between PNPLA3 rs738409 and CHO intake (P = 0.07) and sucrose intake (P = 0.07) on triglyceride levels in the whole study population. The interaction between rs738409 and sucrose intake on triglyceride levels was nominally significant among normal-weight individuals (P = 0.03), while among overweight individuals, a tendency for interaction was observed with ω-6:ω-3 PUFA intake ratio (P = 0.08). In sensitivity analyses, after excluding 932 (19.3 %) individuals with potential inaccurate energy reporting, the interactions between PNPLA3 rs738409 and CHO and sucrose intakes on triglyceride levels were somewhat strengthened in all individuals (P = 0.01 and P = 0.06), while the interaction with CHO intake was somewhat strengthened and with sucrose intake somewhat weakened among normal-weight individuals (P = 0.05 for both). The tendency for interaction with ω-6:ω-3 PUFA intake ratio in overweight individuals was weakened (P = 0.20), and all other observations remained virtually unchanged (see Supplementary table 2).
As we observed a tendency for interaction with the ω-6:ω-3 PUFA intake ratio among overweight individuals and as different dietary fats show a strong correlation with each other, we performed secondary interaction analyses with other fat intake variables (intakes of total fat, SFA, MUFA, PUFA, P:S ratio, ω-3 PUFA and ω-6 PUFA). However, we did not observe any significant interactions (see Supplementary table 3). Adjustments for BMI or for fiber intake did not markedly change the results. Additionally, our observations were similar in women and men (data not shown). The Pearson correlation coefficients between energy-adjusted dietary intakes are shown in Supplementary table 4.
We investigated the interaction between the PNPLA3 rs738409 (Ile148Met) and overweight and dietary intakes of CHO, sucrose and ω-6:ω-3 PUFA ratio on fasting triglyceride levels. Our results indicate that the association between rs738409 and fasting triglycerides, as well as the directionality of this association, is dependent on overweight and may be modified by dietary intakes of CHO, sucrose and ω-6:ω-3 PUFA ratio. Importantly, our results suggest, using fasting triglyceride levels as a marker for liver fat accumulation, that the increased hepatic fat accumulation by 148Met may be a result of both a loss of function of the hydrolysis activity and/or a gain of function of the lipid synthase (LPAAT) activity. Moreover, although we lack the data on liver fat content, our results may suggest that overweight and dietary intakes may play a role on which of the defects of PNPLA3 148Met-mutated protein influences the phenotype.
The proportion of total variation in liver fat content attributed to rs738409 genotypes is among the strongest ever reported for common variants that associate with multifactorial phenotypes. Polymorphisms in three other genes (GCKR, NCAN and LYPLAL1) have been associated with liver fat accumulation and NAFLD (Speliotes et al. 2011), and the NAFLD risk variants of GCKR and NCAN associate additionally with fasting blood triglyceride levels (Speliotes et al. 2011). However, the NAFLD risk allele of GCKR associates with elevated, while that of NCAN with decreased triglyceride levels (Speliotes et al. 2011). The notion that a variant increasing the risk for hepatic fat accumulation can associate with higher or lower circulating triglyceride levels is also demonstrated by our results, suggesting that PNPLA3 148Met is capable of contributing to both lower or higher circulating triglyceride levels.
Although the PNPLA3 knock-out mice do not have any phenotype (Basantani et al. 2011; Chen et al. 2010), administration of the human mutant PNPLA3 148Met to the livers of wild-type mice increased their susceptibility of fatty liver (Huang et al. 2011). In the light of these findings and our results, it can be speculated that an increased LPAAT activity of the PNPLA3 148M produces a phenotype. In line with this, over-expression of hepatic PNPLA3 of mice markedly increased circulating triglyceride levels (Qiao et al. 2010). More recently, studies of transgenic mice over-expressing human PNPLA3 148Met provided evidence for both increased triglycerides synthesis and decreased hydrolysis of hepatic triglycerides in these animals (Li et al. 2012). Further, mice over-expressing the PNPLA3 148Met developed liver steatosis on high sucrose diet, which suggests that the PNPLA3 148Met may be involved in the hepatic de novo lipogenesis (Li et al. 2012).
Dietary sucrose up-regulates PNPLA3 in both transcriptional and posttranslational level while the PNPLA3 is virtually absent in fasted animals (Li et al. 2012; Huang et al. 2010). A high ω-6:ω-3 PUFA intake ratio has been shown to stabilize the protein, further contributing to high levels of PNPLA3 protein (Huang et al. 2010). As also high FFA levels stabilize PNPLA3, one may speculate that obese individuals, in whom FFAs in circulation are generally elevated, may have higher levels of liver PNPLA3. Without data on liver fat levels in our study, it is not easy to interpret the results of our study. We can only speculate that the reduced hydrolysis function may be reflected among overweight carriers of 148Met as lower circulating triglyceride levels, in particular, when dietary intakes of CHOs, sucrose or ω-6:ω-3 PUFA ratio were low. In such environment, the level of PNPLA3 protein can be expected to be very low, and thus, the loss of the lipase activity could be the dominating functional defect, as the importance of a partial loss of function can be expected to increase and that of any gain-of-function defect to decrease by decreasing protein levels and vice versa. In line with this, we further speculate that the reverse enzymatic turning into an increased LPAAT activity by 148Met may be reflected as higher triglyceride levels among normal-weight individuals when the level of PNPLA3 protein is expected to be high, i.e., when dietary intakes of sucrose and ω-6:ω-3 PUFA ratio were high. Two earlier studies have reported interaction between dietary intake of sugar or CHOs (Davis et al. 2010) or ω-6:ω-3 PUFA ratio and liver fat content, both in obese children or adolescents. In the first study of 154 overweight Hispanic children, a positive correlation was observed between hepatic fat content and reported CHO and sugar intake among children homozygous for 148Met (GG-genotype), but not among the other genotype carriers. Similarly, the second study of 127 children and adolescents of mainly Caucasian origin observed that the ω-6:ω-3 PUFA intake ratio correlated positively with liver fat content only among GG-genotype carriers (Santoro et al. 2012). The main finding of these studies was that the positive correlation between high intake of sugar or high ω-6:ω-3 PUFA ratio and liver fat content was restricted to individuals homozygous for 148Met (GG-genotype). Further evidence for that high ω-6: ω-3 PUFA ratio may be particularly important for liver injury among the PNPLA3 G-allele carriers was recently provided by a study of 80 obese adolescents in which oxidized fatty acids derived from linoleic acid associated with markers of liver injury only in risk-genotype carriers (Santoro et al. 2014).
Strengths of our study include the large sample size and high-quality dietary data, based on a 7-day diet diary and large questionnaire, with information on many potential confounding factors and the ability to identify inaccurate reporters of energy intake. Still, the main limitation of our study is the lack of data on hepatic fat content. Another limitation is the lack of data on fatty acid concentrations, like those of linoleic- and alpha linolenic acids, which would have been very informative for the interpretation of the data. Further, we only had ALT levels analyzed in connection to another study almost two decades before MDCS and diet data collection, and therefore, interaction analyses between diet and ALT would not have been appropriate. On the other hand, association analyses between serum ALT levels and PNPLA3 genotypes according to overweight status strengthened our observation of that the association between 148Met and lower triglyceride levels, putatively reflecting the hydrolysis defect, was restricted to overweight individuals and aided at the interpretation of the results of our study without information on liver fat levels. Although nutrient habits may differ between populations, our findings from this population-based study may be applicable to other westernized countries, because the nutrient intake levels could by no means be considered as extreme and are comparable with those from Swedish national surveys. Finally, we need to emphasize that we did not correct our analyses for multiple comparisons, and therefore, the nominally significant or borderline significant interactions with the dietary intakes in this study must be interpreted with caution and need to be replicated in other large studies with dietary data of good quality.
In conclusion, the findings of our epidemiological study indicate that overweight modifies the association between PNPLA3 Ile148Met on circulating triglyceride levels and suggest that dietary intakes of carbohydrates and in particular sucrose, and possibly ω-6:ω-3 PUFA ratio, may additionally modify this association. Our results are in line with the hypothesis that overweight and high sucrose intake facilitate liver fat accumulation in 148Met carriers and may suggest the defective hydrolysis function of PNPLA3-148Met to be more important in the overweight state, while the defective increase in LPAAT activity may be a more important defect when dietary intake of sucrose is high. It is thus important to challenge the clinical question of if the excess liver fat accumulation in 148Met carriers can be eliminated by weight control and reduction in dietary sucrose intake by intervention studies.
Patatin-like phospholipase domain-containing protein 3
Genome-wide association study
Malmö Diet and Cancer Study
Non-alcoholic fatty liver disease
Polyunsaturated fatty acid
Lysophosphatidic acid acyltransferase
Abdemalek FM, Suzuki A, Guy C, Unalp AA, Colvin R, Johnson JR, Diehl MA (2010) Increased fructose consumption is associated with fibrosis severity in patients with NAFLD. Hepatology 51(6):1961–1971
Basantani MK, Sitnick MT, Cai L, Brenner DS, Gardner NP, Li JZ, Schoiswohl G, Yang K, Kumari M, Gross RW, Zechner R, Kershaw EE (2011) PNPLA3/Adiponutrin deficiency in mice does not contribute to fatty liver disease or metabolic syndrome. J Lipid Res 52:318–329
Baulande S, Lasnier F, Lucas M, Pairault J (2001) Adiponutrin, a transmembrane protein corresponding to a novel dietary and obesity linked mRNA specifically expressed in the adipose lineage. J Biol Chem 276:33336–33344
Berglund G, Elmstahl S, Janzon L, Larsson A (1993) The Malmo Diet and Cancer study. Design and feasibility. J Intern Med 233:45–51
Berglund G, Eriksson KF, Israelsson B, Kjellstrom T, Lingarde F, Mattiasson I, Nilsson JA, Stavenow L (1996) Cardiovascular risk groups and mortality in an urban swedish male population: the Malmo Preventive Project. J Intern Med 239(6):489–497
Callmer E, Riboli E, Saracci R, Akesson B, Lindgarde F (1993) Dietary assessment methods evaluated in the Malmo food study. J Intern Med 233(1):53–57
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
Cortez-Pinto H, Jesus L, Barros H, Lopes C, Moura MC, Camilo ME (2006) How different is the dietary pattern in non-alcoholic steatohepatitis patients? Clin Nutr 25:816–823
Davis JN, Le KA, Walker RW, Vikman S, Metz DS, Weigensberg JM, Allayee H, Goran MI (2010) Increased hepatic fat in overweight Hispanic youth influenced by interaction between genetic variation in PNPLA3 and high dietary carbohydrates and sugar consumption. Am J Clin Nutr 92:1522–1527
Del Giudice ME, Grandone A, Cirillo G, Santoro N, Amato A, Brienza C, Savarese P, Marzuillo P, Perrone L (2011) The association of PNPLA3 variants with liver enzymes in childhood obesity is driven by the interaction with abdominal fat. PLoS ONE 6(11):e27933
Dubuquoy C, Robichon C, Lasnier F, Langlois C, Dugail I, Fabienne F, Girard J, Burnol FA, Postic C, Moldes M (2011) Distinct regulation of adiponutrin/PNPLA3 gene expression by the transcription factor ChREBP and SREBP1c in mouse and in human hepatocytes. J Hepatol 55:145–153
Goran MI, Walker R, Le KA, Mahurkar S, Vikman S, Davis JN, Metz-Spruijt D, Weigensberg MJ, Allayee H (2010) Effects of PNPLA3 on liver fat and metabolic profile in Hispanic children and adolescents. Diabetes 59:3127–3130
He S, McPaul C, Li JZ, Garuti R, Kinch L, Grishin NV, Cohen JC, Hobbs HH (2010) A sequence variation (I148M) in PNPLA3 associated with nonalcoholic fatty liver disease disrupts triglyceride hydrolysis. JBC 285(9):6706–6715
Hoekstra M, Li Z, Kruijt JK, Van Eck M, Van Berkel TJC, Kuiper J (2010) The expression level of non-alcoholic fatty liver disease-related gene PNPLA3 in hepatocytes is highly influenced by hepatic lipid status. J Hepatol 52:244–251
Huang Y, He S, Li JZ, Seo KY, Osborne FT, Cohen JC, Hobbs HH (2010) A feed forward loop amplifies nutritional regulation of PNPLA3. Proc Natl Acad Sci USA 107(17):7892–7897
Huang Y, Cohen JC, Hobbs HH (2011) Expression and characterization of a PNPLA3 protein isoform (I148M) associated with Nonalcoholic fatty liver disease. JBC 286(43):37085–37093
Hyysalo J, Stojkovic I, Kotronen A, Hakkarainen A, Sevastinova K, Makkonen J, Lundbom N, Rissanen A, Krauss MR, Melander O, Orho-Melander M, Yki-Järvinen H (2011) Genetic variation in PNPLA3 but not APOC3 influences liver fat in NAFLD. J Gastroenterol Hepatol 27:951–956
Johansson LE, Hoffstedt J, Parikh H, Carlsson E, Wabitsch M, Bondenson AG, Hedenbro J, Tornqvist H, Groop L, Ridderstråle M (2006) Variation in the adiponutrin gene influences its expression and associates with obesity. Diabetes 55:826–833
Kantartzis K, Peter A, Machicao F, Machann J, Wagner S, Konigsrainer A, Schick F, Fritsche A, Haring HU, Stefan N (2009) Dissociation between fatty liver and insulin resistance in humans carrying a variant of the patatin-like phospholipase 3 gene. Diabetes 58:2616–2623
Kotronen A, Johansson LE, Johansson LM, Roos C, Westerbacka J, Hamsten A, Bergholm R, Arkkila P, Arola J, Kiviluoto T, Fisher RM, Ehrenborg E, Orho-Melander M, Ridderstråle M, Groop L, Yki-Järvinen H (2009) A common variant in PNPLA3, which encodes adiponutrin, is associated with liver fat content in humans. Diabetologia 52:1056–1060
Krarup TN, Grarup N, Banasik K, Friedrichsen M, Farch K, Sandholt HC, Jorgensen T, Poulsen P, Witte RD, Vaag A, Sørensen T, Pedersen O, Hansen T (2012) The PNPLA3 rs738409 G-allele associates with reduced fasting serum triglyceride and serum cholesterol in Danes with impaired glucose regulation. PLoS ONE 7(7):e40376
Kumari M, Schoiswohl G, Chitraju C, Paar M, Cornaciu I, Rangrez AY, Wongsiriroj N, Nagy HM, Ivanova PT, Scott SA, Knittelfelder O, Rechberger GN, Gruenberger RB, Eder S, Brown HA, Haemmerle G, Oberer M, Lass A, Kershaw EE, Zimmermann R, Zechner R (2012) Adiponutrin functions as a nutritionally regulated lysophosphatidic acid acyltransferase. Cell Metab 15:691–702
Le KA, Ith M, Kreis R, Faeh D, Bortolotti M, Tran C, Boesch C, Tappy L (2009) Fructose overconsumption causes dyslipidemia and ectopic lipid deposition in healthy subjects with and without family history of type 2 diabetes. Am J Clin Nutr 89:1760–1765
Li JZ, Huang Y, Karaman R, Ivanova PT, Brown HA, Roddy T, Castro-Perez J, Cohen JC, Hobbs HH (2012) Chronic overexpression of PNPLA3I148M in mouse liver causes hepatic steatosis. J Clin Invest 122(11):4130–4144
Liu YM, Moldes M, Bastard JP, Bruckert E, Viguerie N, Hainque B, Basdevant A, Langin D, Pairault J, Clement K (2004) Adiponutrin: a new gene regulated by energy balance in human adipose tissue. J Clin Endocrinol Metab 89(6):2684–2689
Maersk M, Belza A, Jorgensen HS, Ringgaard S, Chabanova E, Thomsen H, Pedersen SB, Astrup A, Richelsen B (2012) Sucrose-sweetened beverages increase fat storage in the liver, muscle and visceral fat depot: a 6-mo randomized intervention study. Am J Clin Nutr 95:283–289
Manjer J, Carlsson S, Elmstahl S, Gullberg B, Janzon L, Lindström M, Mattisson I, Berglund G (2001) The Malmö Diet and Cancer study: representativity, cancer incidence and mortality in participants and non-participants. Eur J Cancer Prev 10(6):489–499
Mattisson I, Wirfalt E, Aronsson CA, Wallstrom P, Berglund G (2005) Misreporting of energy: prevalence, characteristics of misreporters and influence on observed risk estimates in the Malmo Diet and Cancer cohort. Br J Nutr 94(5):832–842
Nilsson PM, Engstrom G, Hedblad B (2007) The metabolic syndrome and incidence of cardiovascular disease in non-diabetic subjects-a population-based study comparing three different definitions. Diabet Med 24(5):464–472
Nobili V, Bedogni G, Donati B, Alisi A, Valenti L (2013) The I148M variant of PNPLA3 reduces the response to docosahexaenoic acid in children with non-alcoholic fatty liver disease. J Med Food 00:1–4
Palmer CNA, Maglio C, Pirazzi C, Burza AM, Adiels M, Burch L, Donnelly AL, Colhoun H, Doney SA, Dillon FJ, Pearson ER, McCarthy M, Hattersley AT, Frayling T, Morris AD, Peltonen M, Svensson PA, Jacobson P, Boren J, Sjöström L, Carlsson LMS, Romeo S (2012) Paradoxical lower serum triglyceride levels and higher type 2 diabetes mellitus susceptibility in obese individuals with the PNPLA3148M variant. PLoS ONE 7(6):e39362
Perttila J, Huaman SC, Caron S, Tanhuanpaa K, Staels B, Yki-Järvinen H, Olkkonen MV (2012) PNPLA3 is regulated by glucose in human hepatocytes, and its I148M slows down triglyceride hydrolysis. Am J Physiol Endocrinol Metab 302:1063–1069
Pignitore P, Pirazzi C, Mancina RM, Motta BM, Indiveri C, Pujia A, Montalcini T, Hedfalk K, Romeo S (2013) Recombinant PNPLA3 protein shows triglyceride hydrolase activity and its I148M mutation results in loss of function. Biochim Biophys Acta S1388-1981(13)00279-5
Pirazzi C, Adiels M, Burza AM, Mancina MR, Levin M, Stahlman M, Taskinen MR, Melander MO, Perman J, Pujia A, Andersson L, Maglio C, Montalcini T, Wiklund O, Borén J, Romeo S (2012) Patatin like phospholipase domain containing 3 (PNPLA3) I148M (rs738409) affects hepatic VLDL secretion in humans and in vitro. J Hepatol 57(6):1276–1282
Polson AD, Thompson MP (2003) Adiponutrin mRNA expression in white adipose tissue is rapidly induced by meal-feeding a high sucrose diet. Biochem Biophys Res Commun 301:261–266
Qiao A, Liang J, Ke Y, Li C, Cui Y, Shen L, Zhang H, Cui A, Liu X, Liu C, Chen Y, Zhu Y, Guan Y, Fang F, Chang Y (2010) Mouse Patatin-Like Phospholipase Domain-Containing 3 influences systemic lipid and glucose homeostasis. Hepatology 54:509–521
Riboli E, Elmstahl S, Saracci R, Gullberg B, Lindgarde F (1997) The Malmo Food Study: validity of two dietary assessment methods for measuring nutrient intake. Int J Epidemiol 26:S161–S173
Romeo S, Kozlitina J, Xing C, Pertsemlidis A, Cox D, Pennacchio LA, Boerwinkle E, Cohen JC, Hobbs HH (2008) Genetic variation in PNPLA3 confers susceptibility to non-alcoholic fatty liver disease. Nat Genet 40:1461–1465
Romeo S, Sentinelli F, Cambuli MV, Incani M, Congiu T, Matta V, Pilia S, Doran HI, Cossu E, Loche S, Baroni MG (2010a) The 148M allele of the PNPLA3 gene is associated with indices of liver damage early in life. J Hepatol 53:335–338
Romeo S, Sentinelli F, Dash S, Yeo GSH, Savage DB, Leonetti F, Capoccia D, Incani M, Maglio C, Iacovino M, O’Rahilly S, Baroni MG (2010b) Morbid obesity exposes the association between PNPLA3 I148M (rs738409) and indices of hepatic injury in individuals of European descent. Inter J Obes 34:190–194
Santoro N, Kursawe R, D’Adamo E, Dykas JD, Zhang KC, Bale EA, Cali MA, Narayan D, Shaw MM, Pierpont B, Savoye M, Lartaud D, Eldrich S, Cushman SW, Zhao H, Shulman GI, Caprio S (2010) A common variant in the patatin-like phospholipase gene (PNPLA3) is associated with fatty liver disease in obese children and adolescents. Hepatology 52:1281–1290
Santoro N, Savoya M, Kim G, Marotto K, Shaw MM, Pierpont B, Caprio S (2012) Hepatic fat accumulation is modulated by the interaction between the rs738409 variant in the PNPLA3 gene and the dietary omega6/omega3 PUFA intake. PLoS ONE 7(5):e37827
Santoro N, Caprio S, Giannini C, Kim G, Kursawe R, Pierpont B, Shaw MM, Feldstein AE (2014) Oxidized fatty acids: A potential pathogenic link between fatty liver and type 2 diabetes in obese adolescents? Antioxid Redox Signal 20(2):383–389
Sevastinova K, Santos A, Kotronen A, Hakkarainen A, Makkonen J, Silander K, Peltonen M, Romeo S, Lundbom J, Lundbom N, Olkkonen VM, Gylling H, Fielding BA, Rissanen A, Yki-Järvinen H (2012) Effect of short term carbohydrate overfeeding and long-term weight loss on liver fat in overweight humans. Am J Clin Nutr 96:727–734
Speliotes EK, Yerges-Armstrong LM, Wu J, Hernaez R, Kim JL, Palmer CD, Gudnason V, Eiriksdottir G, Garcia ME, Launer LJ, Nalls MA, Clark JM, Mitchell BD, Shuldiner AR, Butler JL, Tomas M, Hoffmann U, Hwang SJ, Massaro JM, O’Donnell CJ, Sahani DV, Salomaa V, Schadt EE, Schwartz SM, Siscovick DS, Voight BF, Carr JJ, Feitosa MF, Harris TB, Fox CS, Smith AV, Kao WHL, Hirschhorn JN, Borecki IB (2011) Genome wide association analysis identifies variants associated with nonalcoholic fatty liver disease that have distinct effects on metabolic traits. PLoS Genet 7(3):e1001324
Valenti L, Alisi A, Galmozzi E, Bartuli A, Del Menico B, Alterio A, Dongiovanni P, Fargion S, Nobili V (2010) I148M patatin-like phospholipase domain-containing gene variant and severity of pediatric nonalcoholic fatty liver disease. Hepatology 52:1274–1280
Wirfalt E, Mattisson I, Johansson U, Gullberg B, Wallstrom P, Berglund G (2001) A methodological report from the Malmo Diet and Cancer study: development and evaluation of altered routines in dietary data processing. Nutr J 1(1):3
This study is supported by MO-M research grants from Swedish Research Council, Swedish Heart and Lung Foundation, Novo Nordic Foundation, Påhlsson Foundation, Swedish Diabetes Foundation, and by an equipment grant from Knut and Alice Wallenberg foundation, and Linneus foundation for the Lund University Diabetes Center. None of the funders had any role in the study design, data collection and analyses, interpretation of data, decision to publish or preparation of the manuscript. We would like to thank the participants of the MDC-CC study. We are also very grateful to Mrs. Malin Svensson for excellent technical assistance. Author contributions: IS, UE, MR, SR and MO-M designed the project; IS and GR performed genotyping; IS and UE performed the statistical analyses; IS and MO-M wrote the first draft of the paper; UE, GR and SR critically reviewed the manuscript and MO-M had primary responsibility for the final contents of the paper.
Conflict of interest
Ivana A. Stojkovic, Ulrika Ericson, Gull Rukh, Martin Riddestråle, Stefano Romeo and Marju Orho-Melander declare that they have no conflicts of interest.
The ethical committee at Lund University approved MDCS and MPP, and all the participants have given their written informed consent.
Electronic supplementary material
Below is the link to the electronic supplementary material.