PNPLA3 as a Genetic Determinant of Risk for and Severity of Non-alcoholic Fatty Liver Disease Spectrum

doi:10.14218/JCTH.2016.00009

Original Article

OPEN ACCESS

Download PDF

PNPLA3 as a Genetic Determinant of Risk for and Severity of Non-alcoholic Fatty Liver Disease Spectrum

Habeeb Salameh¹,
Muhannad Al Hanayneh¹,
Maen Masadeh²,
Mohammed Naseemuddin³,
Tasnia Matin³,
Angelika Erwin⁴ and
Ashwani K. Singal^*,5

Journal of Clinical and Translational Hepatology 2016;4(3):175-191

doi: 10.14218/JCTH.2016.00009

Received: March 15, 2016

Revised: May 10, 2016

Accepted: May 23, 2016

Published online: July 8, 2016

Author information

Citation: Salameh H, Hanayneh MA, Masadeh M, Naseemuddin M, Matin T, Erwin A, et al. PNPLA3 as a Genetic Determinant of Risk for and Severity of Non-alcoholic Fatty Liver Disease Spectrum. J Clin Transl Hepatol. 2016;4(3):175-191. doi: 10.14218/JCTH.2016.00009.

Abstract

Background and Aims:Patatin-like phospholipase domain protein 3 (PNPLA3) polymorphisms (rs738409 C>G) are associated with non-alcoholic fatty liver disease (NAFLD). We performed a systematic review and meta-analysis to examine the association of PNPLA3 polymorphisms with the spectrum and severity of this disease. Methods: Studies evaluating the association between the PNPLA3 polymorphism spectrum (fatty liver, steatohepatitis, cirrhosis, and hepatocellular carcinoma) and NAFLD were included. Pooled data are reported as odds ratios (ORs) with 95% confidence intervals. Results: Of 393 potentially relevant studies, 35 on NAFLD were included in the analysis. Compared to healthy controls, the pooled ORs for rs738409 CG and GG compared to CC among patients with non-alcoholic fatty liver (NAFL) were 1.46 (1.16–1.85) and 2.76 (2.30–3.13), and were 1.75 (1.24–2.46) and 4.44 (2.92–6.76) among patients with non-alcoholic steatohepatitis respectively. The respective ORs for CG and GG compared to the CC genotype were 2.35 (0.90–6.13) and 5.05 (1.47–17.29) when comparing non-alcoholic hepatocellular carcinoma to NAFL patients. Among the NAFLD patients, the ORs for G allele frequency when comparing steatosis grade 2–3 to grade 0–1 NAFL, when comparing the NAFLD activity score of ≥ 4 to score ≤ 3, when comparing NASH to NAFLD, when comparing the presence of lobular inflammation to absence, and when comparing the presence of hepatocyte ballooning to absence were 2.33 (1.43–3.80), 1.80 (1.36–2.37), 1.66 (1.42–1.94), 1.58 (1.19–2.10), and 2.63 (1.87–3.69) respectively. Subgroup analysis based on ethnicity showed similar results. Conclusions:PNPLA3 polymorphisms have strong association with the risk for and severity of NAFLDs. PNPLA3 polymorphism plays an evolving role in diagnosis and treatment decisions in patients with NAFLD.

Keywords

PNPLA3, Polymorphisms, Non-alcoholic liver disease, rs738409 C>G

Introduction

Non-alcoholic fatty liver disease (NAFLD) is one of the most common causes of liver cirrhosis and end-stage liver disease in the US, behind hepatitis C virus infection and alcoholic liver disease (ALD).^1,2 Within the past few years, multiple genetic factors have been associated with the severity and predisposition to the spectrum of NAFLDs, including the associated increased risk of hepatocellular carcinoma (HCC).³ Recently, the single nucleotide polymorphism rs738409 C>G (causing an isoleucine to methionine substitution at position 148, I148M) in the patatin-like phospholipase domain containing protein 3 (PNPLA3) gene has been shown to be associated with steatosis among NAFLD patients.³ Since then, other genetic polymorphisms of this gene have been shown to be associated with various other liver diseases, including ALD.⁴In vitro, the human adiponutrin (another name for the PNPLA3 gene product) exhibits activity as a triglyceride hydrolase and transacylase, as well as modest activity as a calcium-independent phospholipase A2. However, in the cellular context, these activities appear minor; similarly, overexpression or knockdown of adiponutrin did not affect lipolysis or the cellular TG level.^5,6

We have shown previously that PNPLA3 is associated with ALD and its severity spectrum (including fatty liver (FL), steatohepatitis, cirrhosis, and HCC).⁴ ALD and NAFLD share several histological features and pathophysiologic mechanisms, including oxidative and endoplasmic reticulum–mediated stress, gut derived endotoxin-mediated cytokine release, and mitochondrial defects.⁷ Over the last 5 years, meta-analyses of pooled data on the associations between the PNPLA3 gene and risk of NAFLD have been reported.^8–10 However, given that these previous meta-analyses were limited and did not examine the association of this genetic polymorphism with the entire spectrum of liver disease—from FL and steatohepatitis to fibrosis, cirrhosis, and HCC—amongst NAFLD patients, we performed the meta-analysis presented herein, with the aim of examining the association of PNPLA3 polymorphisms with the predisposition to and severity of the whole spectrum of liver disease among patients with NAFLDs.

Methods

Identification and retrieval of primary studies

Following the Meta-analysis Of Observational Studies in Epidemiology (MOOSE) guidelines, we searched the literature from the PubMed/Medline, Embase and Cochrane search engines to identify full-length articles written in English that examined PNPLA3 polymorphism in NAFLD patients.¹¹ The initial search terms were “non-alcoholic fatty liver disease” and “adiponutrin, human”. The search was then expanded to include the terms “rs738409 C>G” and “patatin-like phospholipase domain-containing 3 protein”. All databases were searched from their inception through April 2016. Meeting abstracts from major gastroenterology conferences within the past 4 years were also searched to identify studies that were potentially missed by our database searches. Articles were selected for full-text review based on initial reading of the title and abstract. Manual search of the retrieved publications was done by three independent investigators (HS, MM, and MA) to avoid missing articles and ensure as comprehensive search of the literature as possible.

Inclusion and exclusion criteria

Studies reporting association of the PNPLA3 variant (rs738409 C>G) among NAFLD patients and published as full-length articles were eligible for inclusion in this analysis. Studies on patients with other liver diseases but reporting separate data on PNPLA3 polymorphisms for NAFLD were also eligible for inclusion. However, studies without available gene frequency data for analysis and studies including subjects with other liver diseases and without separate data on NAFLD patients were excluded.

Definitions

The disease spectrum across studies on NAFLD patients was defined as non-alcoholic fatty liver (NAFL) according to: steatosis on liver ultrasound in the absence of elevated liver enzymes and/or histologic evidence of steatohepatitis or non-alcoholic steatohepatitis (NASH); steatosis with histologic evidence of steatohepatitis or non-alcoholic cirrhosis (NAC); biopsy or clinical evaluation supported by hematological, biochemical, and/or radiologic imaging findings, and presence of non-alcoholic hepatocellular carcinoma (NAHCC); diagnostic imaging findings on triple0phase magnetic resonance imaging or computed tomography or using histological confirmation from the liver tissue. Healthy controls were defined as subjects without liver disease, for the whole spectrum of NAFLD.

Study selection and data extraction

After determining the eligibility for inclusion, three reviewers (HS, MM, and MA) independently extracted data for (a) study characteristics (author and year of publication, and study design (population-based or not, using controls or not)), (b) study population (liver disease spectrum and the sample size), (c) frequencies of PNPLA3 polymorphism genotypes (rs738409 CC, CG, and GG), and (d) odds ratio (OR; for association of PNPLA3 polymorphism and the spectrum of liver disease and for severity of liver disease). Any discrepancy between the reviewers was resolved by jointly reviewing the study in question. Among studies comparing a diseased population with healthy controls, similar data was also extracted for the healthy controls.

Endpoints and outcomes

Our study endpoints for analysis on the data from included NAFLD studies were risks of (a) entire spectrum of liver disease compared to healthy controls and (b) severity of liver disease spectrum among the NAFLD diseased population.

Assessment of study quality

The quality of included studies was assessed independently by three authors (HS, MM, and MA) using the Newcastle-Ottawa Quality Assessment Scale for case-control studies.¹² This scale has one instrument for assessing case–control studies and another one for cohort studies. Each of these instruments includes measures of quality in selection, comparability, and exposure domains. While one point is granted for each of the areas measured within the selection and exposure domains, a maximum of two points can be assigned within the comparability domain with highest possible total score of nine. A score of seven or greater denotes a high-quality study. Any discrepancies between the three coauthors were addressed by a joint reevaluation of the original article.¹³

Since deviation from Hardy-Weinberg equilibrium (HWE) in healthy controls has been associated with problems in the design and conduct of genetic association studies,^14,15 we examined the studies on the gene allele frequencies of healthy controls for any deviation from HWE. HWE analysis on gene frequencies among healthy controls was examined using the chi-square test, with the value of < 3.84 indicating allele frequency to be in HWE.

Statistical analysis

The strength of association between rs738409 C>G and NAFLD prevalence was expressed by OR and the corresponding 95% confidence interval (CI). Random effects modeling was used for analyzing the pooled data for all the analyses.¹⁶ Heterogeneity was measured using the I² statistic for inter-study variance, with the chi-square test used for statistical analysis.¹⁷ Heterogeneity was defined with I² ≥ 50% or chi-square p < 0.10.¹⁷ At least two studies are needed to examine and report heterogeneity; in our investigations of heterogeneous data and the source of heterogeneity, sensitivity analyses were performed in a stepwise fashion for (a) study quality, (b) excluding studies on the pediatric population for NAFLD analyses, and (c) excluding studies with highest and lowest OR. Publication bias was assessed using Egger’s regression and the Begg and Mazumdar’s rank correlation tests.^18–20 Egger’s test is a regression method that evaluates the association between effect sizes and standard error and uses actual effect size for each study.²⁰ Begg and Mazumdar’s test is a rank correlation test that evaluates the association between effect estimates (taken as a rank and not exact effect size) and sampling variance or standard error.¹⁹ At least three studies are needed for examining and reporting publication bias; in our analyses of publication bias, the analyses were repeated using the Duval and Tweedie trim-and-fill method, a nonparametric (rank-based) data augmentation technique.²¹ Most extreme results on one side of the funnel plot are suppressed and this tool used to estimate the number of studies missing from a meta-analysis. Subsequently, it augments the observed data so that the funnel plot is more symmetric and re-computes the summary estimate based on the complete data.²² This method should be regarded as a way of examining the sensitivity of the results to one particular selection mechanism rather than a way of yielding a more “valid” estimate of the overall effect or outcome. To examine the effect of ethnicity on the association of rs738409 with NAFLD, subgroup analysis was performed by grouping studies with similar populations. All statistical analyses were performed using R software (Foundation for Statistical Computing, Stanford, CA, USA) and its associated metaphor package, or Comprehensive Meta-analysis software (Biostat, Englewood, NJ, USA).

Results

Baseline study characteristics

Of the 393 citations retrieved by the initial search, 462 were reviewed with full-text reading (Fig. 1). Among these, 27 studies were excluded for including patients with ALD (n = 18) and other reasons (n = 9),^3,23–30 including 5 studies due to unavailability of needed data for analysis and 4 studies published as abstracts. The included 35 studies encompassed a total of 13,817 patients (mean age of 52 yrs, 53% male, mean body mass index (BMI) of 29) and reported on PNPLA3 polymorphisms data in NAFLD (Table 1).^31–55 Based on the Newcastle-Ottawa scale, 23 NAFLD studies were of high quality, with study quality score of 7 or more, and the remaining 12 studies were of poor quality (Table 2). Twenty-five studies had control populations.^{31–37,39–43,45–52} Genotype frequency (to perform HWE analysis in healthy controls) was available for only 19 of the studies (Table 3); among these, the control population in only one study deviated from HWE.³¹

Fig. 1. Literature search and study inclusion.

Table 1.

Baseline characteristics of patients from studies included in this analysis

Author (Yr)^Ref	Study design	Controls, n	Cases
			n	Ethnicity	% Males	Mean age	Mean BMI	rs738409 genotype count, CC:CG:GG
			n	Ethnicity	% Males	Mean age	Mean BMI	NAFL	NASH	NAC	NAHCC
Kantartzis (2009)^³¹	P#	50	330	C	39	45	29.51	188:111:31	–	–	–
Sookoian (2009)^³²	P#	94	103	Mixed	28	55	32.5	10:18:12	8:22:33	–	–
Goran (2010)^³⁶Ped	P	117	71	H	46	14	30.3	26:30:15	–	–	–
Hotta (2010)^⁴⁰	P#	578	253	A	48	52	27.7	19:26:19	26:85:78	–	–
Petit (2010)^³³	P	79	139	NR	47	62	34.9	70.4%:29.6%*	–	–	–
Rotman (2010)^³⁵	P#	336	894	Mixed	38	48	34.4	48.2%:51.8%*	50.8%:49.2%*	–	–
Speliotes (2010)^³⁷	P	1405	592	C	37	NR	NR	–	50%:50%*	–	–
Valenti (2010)^³⁴	P#	179	IC: 253	C	70	46	30.5	103:114:36	44:55:23	–	–
Valenti (2010)^³⁴	P#	179	UK&IC: 321	C	62	50	34	142:140:39	65:83:29	–	–
Valenti (2010)^{³⁸ Ped}	P	NA	149	C	62	10	NR	65:61:23	2:45:23	20:29:17	–
Lin (2011)^{⁴¹ Ped}	P#	418	102	A	67	11	26.1	24:52:26	–	–	–
Wagenknecht (2011)^⁵⁵	P#	NA	843	H	37	43	29	324:367:152	–	–	–
Wagenknecht (2011)^⁵⁵	P#	NA	371	AA	38	43	30.1	248:111:12	–	–	–
Wang (2011)^³⁹	P#	723	156	Mixed	60	48	26.9	40:80:36	–	–	–
Kawaguchi (2012)^⁴³	P#	932	527	A	47	55	27.3	49:100:52	39:136:151	–	–
Li (2012)^⁴⁶	P#	202	203	A	51	47	26.6	67:85:51	–	–	–
Peng (2012)^⁴²	P#	553	552	A	72	45	25.4	183:276:93	–	–	–
Petta (2012)^⁴⁴	P#	NA	160	C	68	46	29.7	20:41 **	27:72 **	–	–
Zain (2012)^⁴⁵	P	198	144	A	53	51	28.7	65%:35%*	52%:48%*	–	–
Kitamoto (2013)^⁴⁷	P	1,012	564	A	52	50	27.8	96:241:227	72:183:187	–	–
Verrijken (2013)^⁵⁴	P#	NA	287	C	30	44	37.6	164:106:17	66:69:16	–	–
Guichelaar (2014)^⁴⁹	P	12	132	NR	15	48	46.6	43:13:04	36:28:08	–	–
Islek (2014)^⁴⁸	P#	303	80	C	51	NR	NR	–	23:43:14	–	–
Kanth (2014)^⁵¹	P#	150	156	A	89	37	27	60%:40%*	–	–	–
Lin (2014)^{⁵² Ped}	P#	606	191	A	70	12	NR	58:95:38	–	–	–
Liu (2014)^⁵³	P	NA	375	C	57	56	33.8	125:117:33	–	–	028:043:029
Baclig (2014)^⁷²	P#	36	32	A	50	NA	28.5	14:18 **	–	–	–
Lee (2014)^⁷³	P#	184	155	NA	65.2	44.7	27.5	20:48:31	–	–	–
Park (2015)^⁷⁴	P#	NA	1363	A	100	49	24.4	172:293:137	–	–	–
Wei (2015)^⁷⁵	P#	NA	911	A	41.9	48	NA	359:550 **	–	–	–
Xia (2015)^⁷⁶	P#	2915	1385	A	38.40	62	25.9	486:684:215
Oniki (2015)^⁷⁷	P#	472	119	A	69.70	61.2	25.6	23:69:27	–	–	–
Nishioji (2015)^⁷⁸	P#	NA	824	A	66.5	NA	NA	–	–	–	–
Akuta (2015)^⁵⁶	P#	NA	211	A	57.8	52	25.9	21:60:59	–	–	–
Tai (2015)^⁷⁹	P#	29	152	A	32	30.2	41.8	31:24:05	35:36:21	–	–
Tai (2016)^⁸⁰	P#	NA	177	A	38	29.8	42.7	77:73:27
Honda (2016)^⁸¹ non-obese	P#	782	134	A	38.4	5.51	25.9	23:47:64
Honda (2016)^⁸¹-obese	P#	230	406	A	38.4	5.51	25.9	78:180:148
Summary	–	12,595	13,817	–	51.85	53.2	28.9	–

$ Genotype counts were reported as ratios (CC wild genotype: CG heterozygote genotype: GG homozygote genotype) or * as allele frequency percentage (C allele%: G allele%) or ** as ratios (CC wild genotype: CG and GG genotypes).

Abbreviations: A, Asian; AA, African American; BMI, body mass index; C, Caucasian; H, Hispanic; IC, Italian cohort; NA, not applicable; NAC, non-alcoholic cirrhosis; NAHCC, non-alcoholic hepatocellular carcinoma; NAFLD, non-alcoholic liver disease; NASH, non-alcoholic steatohepatitis; NR, not reported; P, prospective; Ped, pediatric studies; R, retrospective; P#, population-based studies; UK&I, United Kingdom and Italian.

Table 2.

Newcastle-Ottawa scale quality score of the included studies

Author (year)^Ref	Selection				Comparability	Exposure			Total
Author (year)^Ref	Independent validation	Case representation	Controls selection	Controls definition	Case and control design/analysis	Ascertainment of exposure	Same method of ascertainment	Same response rate	Total
Kantartzis (2009)³¹	*			*	**	*	*	*	7
Sookoian (2009)³²	*	*	*	*	**	*	*		8
Petit (2010)³³	*	*		*	*	*	*		6
Valenti (2010)³⁴	*	*	*	*	**	*	*		8
Rotman (2010)³⁵	*		*	*	**	*	*		7
Goran (2010)³⁶	*	*	*		**	*	*		7
Speliotes (2010)³⁷	*	*	*		**	*	*		7
Valenti (2010)³⁸	*	*			*	*	*		5
Hotta (2010)⁴⁰	*	*	*	*	**	*	*		8
Wang (2011)³⁹	*	*	*		**	*	*		7
Lin (2011)⁴¹	*	*	*	*	**	*	*		8
Wagenknecht (2011)⁵⁵	*	*			**	*			5
Peng (2012)⁴²	*	*	*	*	**	*	*		8
Kawaguchi (2012)⁴³	*	*	*	*	**	*	*		8
Petta (2012)⁴⁴	*	*			*	*			4
Zain (2012)⁴⁵	*	*	*	*	**	*	*		8
Li (2012)⁴⁶	*		*	*	**	*	*		7
Kitamoto (2013)⁴⁷	*		*	*	**	*	*		7
Verrijken (2013)⁵⁴	*	*			**	*			5
Islek (2014)⁴⁸	*		*	*		*	*		5
Guichelaar (2014)⁴⁹	*	*			**	*			5
Tai (2014/2015)⁵⁰	*		*	*	**	*	*		7
Kanth (2014)⁵¹	*		*	*	**	*	*		7
Lin (2014)⁵²	*	*	*	*	**	*	*		8
Liu (2014)⁵³	*	*	*	*	**	*	*		8
Baclig (2014)⁷²	*	*	*	*	**	*	*		8
Lee (2014)⁷³	*	*	*	*	**	*	*		8
Park(2015)⁷⁴	*	*			**	*	*		6
Wei (2015)⁷⁵	*	*			**	*	*		6
Xia (2015)⁷⁶	*	*	*	*	**		*		7
Oniki (2015)⁷⁷	*	*	*	*	**	*	*		8
Nishioji (2015)⁷⁸	*	*			**	*	*		6
Akuta (2015)⁵⁶	*	*				*	*		4
Tai (2016)⁸⁰	*	*			**	*	*		6
Honda (2016)⁸¹	*	*	*	*	**	*	*	*	9

Table 3.

Observed versus expected PNPLA3 rs738409 genotype frequency in the control populations included in this meta-analysis. Expected genotype frequency is based on the Hardy-Weinberg equilibrium (HWE), which was calculated using the chi-square test. The chi-square value reflects the difference between expected and observed genotype counts. Values < 3.84 indicate that the allele frequency in the population does not deviate from HWE

Study	Observed genotype frequencies CC:CG:GG	Expected genotype frequencies CC:CG:GG	Chi-square value
Kantartzis (2009)³¹	137:70:18	131.48:81.03:12.48	4.17
Valenti (2010)³⁴	118:56:5	119.08:53.83:6.08	0.29
Hotta (2010)⁴⁰	175:296:104	181.44:283.12:110.44	1.19
Goran (2010)³⁶	38:60:19	39.52:56.96:20.52	0.33
Wang (2011)³⁹	269:335:119	263.53:345.94:113.53	0.72
Lin (2011)⁴¹	167:192:59	165.48:195.05:57.48	0.1
Li (2012)⁴⁶	95:91:16	97.72:85.55:18.72	0.81
Kawaguchi (2012)⁴³	247:468:217	248.24:465.52:218.24	0.03
Peng (2012)⁴²	235:259:59	240.25:248.49:64.25	0.99
Kitamoto (2013)⁴⁷	300:513:199	306.02:500.96:205.02	0.58
Islek (2014)⁴⁸	156:129:18	160.46:120.07:22.46	1.67
Guichelaar (2014)⁴⁹	9:3:0	9.19:2.63:0.19	0.25
Tai (2014)⁵⁰	12:15:2	13.11:12.78:3.11	0.87
Lin (2014)⁵²	238:293:75	243.96:281.08:80.96	1.09
Lee (2014)⁷³	30:50:20	30.1:49.5:20.3	0.01
Oniki (2015)⁷⁷	143:249:80	151.6:231.8:88.6	2.60
Park (2015)⁷⁴	280:364:117	280.5:363.0:117.5	0.00
Xia (2015)⁷⁶	1200:1363:352	1214.4:1334.2:366.4	1.36
Honda(2016)⁸¹ non-obese	237:389:156	238.1:386.8:157.1	0.02
Honda (2016)⁸¹ obese	63:124:43	67.9:114.1:47.9	1.72

PNPLA3 polymorphisms in NAFLD compared to healthy controls

NAFL

Among the 26 included studies on 10,236 NAFLD patients, the pooled ORs for rs738409 C> -CG and GG- genotypes compared to the CC genotype in NAFL patients were 1.46 (1.16–1.85) and 2.76 (2.30–3.13) respectively (Fig. 2A, 2B). The data were heterogeneous for the CG analysis but not for the GG analysis, with respective values of I² = 67.4%, p < 0.0001 and I² = 27.9%, p = 0.16. No publication bias was found by Egger’s test (p = 0.14 and 0.28) or Begg and Mazumdar’s test (p = 0.63 and 0.74). The pooled OR for G allele frequency for NAFL compared to healthy controls was 1.91 (1.64–2.21) (Fig. 2C), with heterogeneous data (I² = 72.6%, p < 0.0001) and publication bias (by Egger’s test, p = 0.0023; by Begg and Mazumdar’s test, p = 0.009). Sensitivity analysis using the trim-and-fill method showed homogeneous data, with similar effect size of 1.65 (1.55–1.77).

Forest plots for analysis of non-alcoholic fatty liver studies comparing patients to healthy controls for association of <italic>PNPLA3</italic> polymorphisms with non-alcoholic fatty liver: A) CG vs. CC, B) GG vs. CC, and C) G allele.

Fig. 2. Forest plots for analysis of non-alcoholic fatty liver studies comparing patients to healthy controls for association of PNPLA3 polymorphisms with non-alcoholic fatty liver: A) CG vs. CC, B) GG vs. CC, and C) G allele.

The bottom line in the “statistics for each study” represents the pooled effect size that was analyzed using the random effects model. An OR > 1 denotes risk for the respective outcome or positive association and an OR < 1 indicates a protective effect or negative association. A 95% confidence interval not crossing 1 indicates a significant association. NO, non-obese; O, obese.

Among studies of Asian populations, the pooled ORs for the CG (8 studies with 2,610 patients) and GG (10 studies with 2,748 patients) genotypes as compared to the CC genotype were 1.50 (1.13–2.00) and 2.67 (2.16–3.31) respectively. The data were heterogeneous for the CG analysis (I² = 70%, p = 0.001) and homogeneous for the GG analysis (I² = 32%, p = 0.15). No publication bias was found by Egger’s test (p = 0.17 and 0.59) or Begg and Mazumdar’s test (p = 0.27 and 1.0). Similarly, among studies of Caucasian populations, the pooled ORs for the CG (2 studies with 904 patients) and GG (2 studies with 904 patients) genotypes as compared to the CC genotype were 1.29 (0.29–3.49) and 3.51 (2.00–6.14) respectively. The data were heterogeneous for both analyses (I² = 84%, p = 0.01 and I² = 89%, p < 0.0001); however, publication bias could not be assessed, as there were only two studies in each analysis. The pooled OR for the G allele frequency in 12 Asian studies with 5,319 patients was 1.78 (1.52–2.08), with heterogeneous data (I² = 70, p < 0.0001) and publication bias as assessed by Egger’s test (p = 0.04) but not by Begg and Mazumdar’s test (p = 0.12). Other ethnicities (Caucasians, Blacks, and Hispanics) were reported in one study each (Table 4). Combined analysis on the non-Asian population for the G allele frequency was 3.29 (1.89–5.54), with homogenous data (I² = 51, p = 0.13) and no publication bias (Egger’s test, p = 0.37; Begg and Mazumdar’s test, p = 0.5).

Table 4.

Ethnicity-based analysis on the association between PNPLA3 and risk for and severity of disease among NAFLD studies

Comparison	Ethnicity	No. of studies	Genotypes	OR	95% CI
NAFL vs. controls	Caucasians	2	CG vs CC	1.29	0.29–3.49
	Asians	8	CG vs CC	1.50	1.13–2.00
	Caucasians	2	GG vs CC	3.51	2.00–6.14
	Hispanics	1		4.70	2.39–9.25
	Asians	10		2.67	2.16–3.31
	Caucasians	1	G allele	2.38	1.36–4.17
	Blacks	1		9.71	2.78–33.84
	Hispanics	1		2.97	1.90–4.80
	Asians	12		1.78	1.52–2.08
NASH vs. controls	Caucasian	2	G allele	3.20	2.72–3.75
NASH vs. controls	Asian	3	G allele	2.25	1.94–2.61
NAFL 2–3 vs. NAFL 0–1	Asian	1	G allele	3.32	1.56–7.01
NAFL 2–3 vs. NAFL 0–1	Caucasian	4	G allele	2.69	1.27–5.69
NASH vs. NAFL	Caucasian	3	G allele	1.42	1.01–2.00
NASH vs. NAFL	Asian	5	G allele	1.84	1.53–2.21
NAFLD activity score	Mixed	1	G allele	1.56	1.12–2.17
NAFLD activity score	Asian	3	G allele	2.05	1.35–3.12
Lobular inflammation	Mixed	1	G allele	1.84	1.33–2.55
	Caucasian	2		1.47	1.19–1.83
	Asian	2		1.72	0.52–5.69
Hepatocyte ballooning	Caucasian	2	G allele	3.17	2.68–3.76
Hepatocyte ballooning	Asian		G allele	1.81	1.07–3.03
Fibrosis/Cirrhosis	Caucasian	2	G allele	1.41	1.13–1.74
Fibrosis/Cirrhosis	Asian	2	G allele	2.16	1.90–2.58

Abbreviations: NAFL, non-alcoholic fatty liver; NAFLD, non-alcoholic fatty liver disease; NASH, non-alcoholic steatohepatitis.

NASH

For the two studies on 353 NAFLD patients, the pooled ORs for the rs738409 C>G -CG and GG- genotypes compared to the CC genotype for presence of NASH were 1.75 (1.24–2.46) and 4.44 (2.92–6.76) respectively (Fig. 3A, 3B). The data were very robust for both analyses, with respective values of I² = 0%, p = 0.55 and I² = 0%, p = 0.38. The Pooled OR for the G allele frequency from five studies with 1,596 NASH patients was 2.54 (2.03–3.16) (Fig. 3C). The data were heterogeneous (I² = 67%, p = 0.019) and there was no publication bias (Egger’s test, p = 0.56; Begg and Mazumdar’s test, p = 0.5). Ethnicity-based analysis among the 3 Asian studies encompassing 924 patients showed similar effect size of 2.25 (1.94–2.61), with homogeneous data (I² = 0%, p = 0.73). In additional, two Caucasian studies encompassing 672 patients showed an effect size was 3.20 (2.72–3.75), with homogeneous data (I² = 0%, p = 0.45) (Table 4).

Fig. 3. Forest plots for analysis of non-alcoholic fatty liver studies comparing patients to healthy controls for association of PNPLA3 polymorphisms with non-alcoholic steatohepatitis: A) CG vs. CC, B) GG vs. CC, and C) G allele frequency.

The effect size is reported as odds ratio with 95% confidence interval. The bottom line in the “statistics for each study” represents the pooled effect size that was analyzed using the random effects model. An OR > 1 denotes risk for the respective outcome or positive association and an OR < 1 indicates a protective effect or negative association. A 95% confidence interval not crossing 1 indicates a significant association.

Association of PNPLA3 polymorphisms with the severity of NAFLD

Steatosis grade 2–3 vs. 0–1

The ORs of the rs738409 C>G -CG and GG- genotypes compared to the CC genotype for NAFL grade 2–3 vs. grade 0–1 were 1.17 (0.35–3.92) and 2.15 (0.96–4.84) respectively (Fig. 4A, 4B). The data for the GG analysis were heterogeneous (I² = 67.8%, p < 0.045) and there was no publication bias (Egger’s test, p = 0.61; Begg and Mazumdar’s test, p = 1.0). After excluding the study with extremely high OR and wide CI effect size, the results remained similar, with OR of 1.64 (1.10–2.46) and homogenous data (I² = 8.9%, p = 0.30). The OR for the G allele frequency pooled from six studies encompassing 2,521 patients was 2.33 (1.43–3.80) (Fig. 4C), with heterogeneous data (I² = 86%, p < 0.0001) and evidence of publication bias by Egger’s test (p = 0.025) but not by Begg and Mazumdar’s test (p = 0.13). Ethnicity based analysis for four Caucasian studies encompassing 1,475 patients showed the OR for the G allele frequency to be 2.69 (1.27–5.69) (Table 4), with homogeneous data. However, after excluding one study on a pediatric population,³⁸ the pooled OR for the G allele frequency became 1.54 (1.19–1.98), with homogeneous data (I² = 48.9%, p = 0.1) and evidence of publication bias by the Egger’s test (p = 0.019). Using the trim-and-fill method, the effect size remained unchanged at 2.33 (1.43–3.80), without any publication bias.

Forest plots for analysis of the risk of spectrum of non-alcoholic fatty liver disease on the association of <italic>PNPLA3</italic> polymorphism with steatosis grades 2–3 vs. 0–1: A) CG vs. CC, B) GG vs. CC, and C) G allele frequency.

Fig. 4. Forest plots for analysis of the risk of spectrum of non-alcoholic fatty liver disease on the association of PNPLA3 polymorphism with steatosis grades 2–3 vs. 0–1: A) CG vs. CC, B) GG vs. CC, and C) G allele frequency.

NASH in NAFLD

The pooled ORs for presence of NASH among NAFLD patients from four studies encompassing 1,146 patients for the CG and GG genotypes compared to the CC genotype were 1.92 (1.43–2.57) and 3.53 (2.02–6.15) respectively (Fig. 5A, 5B); the data were homogeneous, with respective values of I² = 0%, p = 0.51 and I² = 0%, p = 0.8. There was no publication bias for the CG analysis (Egger’s test, p = 0.7; Begg and Mazumdar’s test, p = 1.0). There was also no publication bias for either the CG or the GG analysis (Egger’s test, p = 0.43 and 0.2; Begg and Mazumdar’s test, p = 1.0 and 0.29). Ethnicity-based analysis was not performed, as all studies were carried out with Asian populations. The OR for the G allele frequency pooled from 10 studies encompassing 2,814 patients was 1.66 (1.24–1.94) (Fig. 5C). The data were homogenous (I² = 21.4%, p = 0.24) and there was no publication bias (Egger’s test, p = 0.096; Begg and Mazumdar’s test, p = 0.15). Ethnicity-based analysis showed that the OR for the G allele frequency in three Caucasian studies encompassing 1,326 patients was 1.42 (1.01–2.00) and in five Asian studies was 1.84 (1.53–2.21); the data were homogenous for both analyses, with respective values of I² = 54%, p = 0.11 and I² = 0%, p = 0.29. In addition, there was no publication bias in either analysis (Egger’s test, p = 0.39 and 0.29; Begg and Mazumdar’s test, p = 0.30 and 0.46). One study reported data in a population of mixed ethnicities, without separate ethnicity-based data.³²

Fig. 5. Forest plots for analysis of the risk of spectrum of non-alcoholic fatty liver disease on the association of PNPLA3 polymorphism with non-alcoholic steatohepatitis vs. non-alcoholic fatty liver: A) CG vs. CC, B) GG vs. CC, and C) G allele frequency.

NAFLD activity score and histopathologic severity

NAS

The pooled OR for the G allele frequency in the presence of NAS > 4 vs. < 3 among NAFLD patients from four studies encompassing 1,586 patients was 1.80 (1.36–2.37) (Fig. 6A); the data were homogenous (I² = 19.8%, p = 0.12) and there was no publication bias (Egger’s test, p = 0.12; Begg and Mazumdar’s test, p = 0.31). However, after excluding a study that used a population of mixed ethnicities,³⁵ the OR from the three Asian studies, now encompassing 692 patients, was 2.05 (1.35–3.12), with no publication bias (Egger’s test, p = 0.32; Begg and Mazumdar’s test, p = 1.0).

Fig. 6. Forest plots for analysis of the risk of spectrum of non-alcoholic fatty liver disease on the association of PNPLA3 polymorphism with G allele frequency in NAFLD activity score in panel (A), in lobular inflammation in panel (B) and in hepatocyte ballooning in panel (C).

Lobular Inflammation

The pooled OR for the G allele frequency for the presence of lobular inflammation in five studies encompassing 2,313 patients was 1.58 (1.19–2.10) (Fig. 6 B); the data were heterogeneous (I² = 56.9%, p = 0.054) and there was no publication bias (Egger’s test, p = 0.46; Begg and Mazumdar’s test, p = 0.46). Ethnicity-based analysis (Table 4) showed that the OR for the G allele frequency in two studies of Caucasian populations encompassing 879 patients was 1.47 (1.19–1.83) and in two studies of Asian populations encompassing 540 patients was 1.72 (0.52–5.69). The data were homogenous in the Caucasian analysis (I² = 0%, p = 0.40) but not in the Asian analysis (I² = 84%, p = 0.012). No analysis of publication bias was performed, as only two studies were included in each analysis.

Hepatocyte ballooning

The pooled OR for the G allele frequency for hepatocyte ballooning in four studies encompassing 1,419 patients was 2.63 (1.87–3.69) (Fig. 6C); the data were homogenous (I² = 43.3%, p = 0.15) and there was no publication bias (Egger’s test, p = 0.24; Begg and Mazumdar’s test, p = 0.31). Ethnicity-based analysis (Table 4) showed that the OR for the G allele frequency in two studies of Caucasian populations encompassing 879 patients was 3.17 (2.68–3.76) and in two studies of Asian populations encompassing 540 patients was 1.81 (1.07–3.03). The data were homogeneous in both analyses (I² = 0%, p = 0.66 and I² = 0%, p = 0.32 respectively)

Portal inflammation

Only one study encompassing 894 patients reported the OR of the G allele frequency for portal inflammation, and the value was 1.57 (1.24–1.99).

Fibrosis and cirrhosis

Among the eight studies of advanced cirrhosis and fibrosis encompassing 3,571 patients, the pooled ORs for the rs738409 C>G CG and GG genotypes compared to the CC genotype were 18.80 (2.12–166.75) (from one study⁵⁶) and 2.72 (0.84–8.83) respectively (Fig. 7A). The data were heterogeneous for the GG analysis (I² = 76.8%, p = 0.013) and there was no publication bias (Egger’s test, p = 0.22; Begg and Mazumdar’s test, p = 0.29). When the pediatric study³⁸ was excluded, the pooled OR showed a similar effect size of 1.43 (0.94–2.17). The data were homogenous (I² = 0%, p = 0.61). The pooled OR for the G allele frequency was 1.67 (1.37–2.03) (Fig. 7B), with heterogeneous data (I² = 55.9%, p = 0.045) and no publication bias (Egger’s test, p = 0.61; Begg and Mazumdar’s test, p = 1.0). The pooled OR for the G allele frequency after excluding the pediatric study³⁸ was 1.64 (1.32–2.04) (data not shown in figure form), with heterogeneous data (I² = 64%, p = 0.03) and no publication bias (Egger’s test, p = 0.49; Begg and Mazumdar’s test, p = 1.0). Ethnicity-based analysis showed that among two studies on Caucasian populations encompassing 1,166 patients the pooled OR for the G allele frequency was 1.41 (1.13–1.74), with homogeneous data (I² = 0%, p = 0.61), and among the two studies on Asian populations encompassing 682 patients the pooled OR for the G allele frequency was 2.16 (1.90–2.58), with homogeneous data (I² = 0%, p = 0.60) (Table 4).

Fig. 7. Forest plots for analysis of the risk of spectrum of non-alcoholic fatty liver disease on the association of PNPLA3 polymorphism with fibrosis/cirrhosis vs. NAFLD: A) GG vs CC, and B) G allele frequency.

The bottom line in the “statistics for each study” represents the pooled effect size that was analyzed using the random effects model. An OR > 1 denotes risk for the respective outcome or positive association and an OR < 1 indicates protective effect or negative association. A 95% confidence interval not crossing 1 indicates a significant association. NO, non-obese; O, obese. Valenti 10 refers to the pediatric study.³⁸

HCC in NAFLD

The ORs of the rs738409 C>G -CG and GG- genotypes as compared to the CC genotype for the risk of NAHCC among NAFL patients without cirrhosis were 2.35 (0.90–6.13) and 5.05 (1.47–17.29) respectively. The pooled OR of the rs738409 C>G CG as compared to the CC genotype for NAHCC among patients with NAC was 2.06 (1.07–3.94). All these analyses were reported in one study each, and the forest plots are not shown for these analyses.⁵³

Discussion

This meta-analysis includes data for both adult and pediatric NAFLD, by which we demonstrate that PNPLA3 polymorphisms are associated with a) the risk of NAFL and NASH compared to healthy controls and b) more severe disease on the spectrum of NAFLD, where the polymorphism correlated with the grade of steatosis, NAS, histopathologic findings, cirrhosis and when NASH or NAHCC were compared to NAFL.

The PNPLA3 gene, a transmembrane protein encoding for adiponutrin in human, is highly expressed in the liver and adipose tissues. Its expression in subcutaneous and intra-abdominal adipose tissues correlates with obesity.⁵⁷ The biochemical function of adiponutrin is uncertain, but it is considered to have lipogenic transacetylase activity, likely facilitating energy mobilization and lipid storage in adipose and liver tissues.^57,58 Recently, it has been reported that the 148M adiponutrin allele is a loss-of function variant that predisposes an individual to steatosis by decreasing triglyceride hydrolysis in hepatocytes.⁵⁹ However, this genetic variant does not appear to correlate with visceral or subcutaneous fat content, insulin sensitivity, or peripheral blood lipid levels.^37,38 Consistent with the phenotype of NAFLD patients, overexpression of PNPLA3 G allele increases the size of lipid droplets, as compared with the overexpression of PNPLA3 C allele.⁶⁰ The results of this meta-analysis are compatible with the mechanism of action of PNPLA3 and its presumed effect on liver toxicity. The exact mechanism by which PNPLA3 polymorphism increases the risk of HCC is yet to be determined. PNPLA3 was recently linked to increasing circulating levels of the intercellular adhesion molecule 1 (ICAM-1),⁶¹ which is a proinflammatory marker and characterized as a determinant of the malignant nature of tumors;^62,63 in addition, it has been shown to be associated with decreased levels of adiponectin, a protein with oncosuppressive properties.⁶⁴

This is the first comprehensive meta-analysis examining the association of PNPLA3 polymorphisms with the risk for and severity of the spectrum of liver disease among NAFLD patients. The previous meta-analysis by Xu et al⁹ included 23 studies and looked at the association between NAFL/NASH and PNPLA3 polymorphism, but did not include NAHCC, NAC, steatosis grade, NAS, etc. The other most recent meta-analysis published in the literature investigated the association between NAFLD and PNPLA3 in studies of Asian populations only.¹⁰ Our current analysis included a total of 25 NAFLD studies and confirmed these findings, expanding our understanding of the associations with different liver disease phenotypes. In addition, we confirmed the HWE calculations when possible, looked at the study quality assessment based on the Newcastle-Ottawa scale, and performed subgroup analysis based on age and ethnicity to confirm our findings. All of these features of our study make the current analysis more comprehensive than the previous ones for gaining a more detailed understanding of the NAFLD and PNPLA3 genetic association.

Herein we present results from our meta-analysis of ten ALD studies, which also included individual participant data from five studies, showing that the PNPLA3 polymorphism was associated with predisposition to and severity of the ALD spectrum.⁴ When we compared the current meta-analysis with the previous one we noted certain differences between NAFLD and the ALD disease spectrum and the association with the PNPLA3 polymorphism. First of all, there is a stronger association amongst NAFLD, as compared to the ALD studies. This could be due to more studies published and analyzed with larger sample size for the current analyses. Secondly, the association was lacking with alcoholic FL but it was more prominent in the current study for association with NAFL. A similar prevalence of the polymorphisms in drinkers with FL and non-drinkers aligns well with the fact that the FL in alcoholics is a universal phenomenon in alcoholics and is reversible upon abstinence.^65–67

Our study is potentially limited by the possibility of publication bias. In order to minimize this potential limitation, however, and to subsequently overestimate the true effect size due to negative studies’ identification failure,⁶⁸ we combined searches from the PubMed/Medline, Embase and Cochrane databases with manual searches. Unfortunately, we could not examine the association of PNPLA3 polymorphism with NAC and NAHCC compared to controls due to the lack of studies examining this association. Unavailability of individual patient data also precluded our ability to examine the impact of confounders such as age, sex, comorbidities, and BMI on the association with the PNPLA3 polymorphisms.

Another limitation to our study relates to the differences in the definition of NASH used across the studies included in this analysis. Bedossa et al⁶⁹ described various features of hepatic injury in NASH, with the histological finding of lobular inflammation and hepatocyte ballooning being key findings. Different scoring systems were proposed, such as the Steatosis, Activity and Fibrosis system⁷⁰ and the NASH Clinical Research Network system.⁷¹ To resolve this limitation, we grouped studies with almost similar definitions for the sake of homogeneity of the analysis.

In summary, the present study provided unequivocal evidence of rs738409 as a strong modifier of the natural history of NAFLD in different ethnic and age populations. As such, the PNPLA3 gene may be a potential target for therapy in NAFLD. Prospective data are now needed to further understand the association of PNPLA3 polymorphisms, particularly related to a) a response to control of risk factors of NAFLD and b) for prediction of the natural history of the disease.

Abbreviations

ALD:: alcoholic liver disease

BMI:: body mass index

CRN:: clinical research network

FL:: fatty liver

HCC:: hepatocellular carcinoma

HWE:: Hardy-Weinberg equilibrium

ICAM-1:: intercellular adhesion molecule 1

MOOSE:: Meta-analysis Of Observational Studies in Epidemiology

NAC:: non-alcoholic cirrhosis

NAHCC:: non-alcoholic hepatocellular carcinoma

NAFL:: non-alcoholic fatty liver

NAFLD:: non-alcoholic fatty liver disease

NAS:: NAFLD activity score

NASH:: non-alcoholic steatohepatitis

OR:: odds ratio

PNPLA3:: patatin-like phospholipase domain protein 3

Declarations

Conflict of interest

None

Authors’ contributions

Study design (HS, AKS), data collection (HS, MAH, MM, MN, TM, AE), drafting the manuscript (HS, MAH, MM, MN, TM, AE, AKS), data analysis and interpretation (HS, AKS), critical revision and funding (AKS), statistical analysis (AKS, AE).

References

1	Singal AK, Guturu P, Hmoud B, Kuo YF, Salameh H, Wiesner RH. Evolving frequency and outcomes of liver transplantation based on etiology of liver disease. Transplantation 2013;95:755-760 View Article

2	Charlton MR, Burns JM, Pedersen RA, Watt KD, Heimbach JK, Dierkhising RA. Frequency and outcomes of liver transplantation for nonalcoholic steatohepatitis in the United States. Gastroenterology 2011;141:1249-1253 View Article

3	Kotronen A, Johansson LE, Johansson LM, Roos C, Westerbacka J, Hamsten A. A common variant in PNPLA3, which encodes adiponutrin, is associated with liver fat content in humans. Diabetologia 2009;52:1056-1060 View Article

4	Salameh H, Raff E, Erwin A, Seth D, Nischalke HD, Falleti E. PNPLA3 gene polymorphism is associated with predisposition to and severity of alcoholic liver disease. Am J Gastroenterol 2015;110:846-856 View Article

5	Kienesberger PC, Oberer M, Lass A, Zechner R. Mammalian patatin domain containing proteins: a family with diverse lipolytic activities involved in multiple biological functions. J Lipid Res 2009;50:S63-S68 View Article

6	Kershaw EE, Hamm JK, Verhagen LA, Peroni O, Katic M, Flier JS. Adipose triglyceride lipase: function, regulation by insulin, and comparison with adiponutrin. Diabetes 2006;55:148-157 View Article

7	Day CP. From fat to inflammation. Gastroenterology 2006;130:207-210 View Article

8	Sookoian S, Pirola CJ. Meta-analysis of the influence of I148M variant of patatin-like phospholipase domain containing 3 gene (PNPLA3) on the susceptibility and histological severity of nonalcoholic fatty liver disease. Hepatology 2011;53:1883-1894 View Article

9	Xu R, Tao A, Zhang S, Deng Y, Chen G. Association between patatin-like phospholipase domain containing 3 gene (PNPLA3) polymorphisms and nonalcoholic fatty liver disease: a HuGE review and meta-analysis. Sci Rep 2015;5:9284 View Article

10	Zhang L, You W, Zhang H, Peng R, Zhu Q, Yao A. PNPLA3 polymorphisms (rs738409) and non-alcoholic fatty liver disease risk and related phenotypes: a meta-analysis. J Gastroenterol Hepatol 2015;30:821-829 View Article

11	Stroup DF, Berlin JA, Morton SC, Olkin I, Williamson GD, Rennie D. Meta-analysis of observational studies in epidemiology: a proposal for reporting. Meta-analysis Of Observational Studies in Epidemiology (MOOSE) group. JAMA 2000;283:2008-2012 View Article

12	Available at: http://www.ohri.ca/programs/clinical_epidemiology/oxford.asp

13	Ungaro R, Bernstein CN, Gearry R, Hviid A, Kolho KL, Kronman MP. Antibiotics associated with increased risk of new-onset Crohn’s disease but not ulcerative colitis: a meta-analysis. Am J Gastroenterol 2014;109:1728-1738 View Article

14	Minelli C, Thompson JR, Abrams KR, Thakkinstian A, Attia J. How should we use information about HWE in the meta-analyses of genetic association studies?. Int J Epidemiol 2008;37:136-146 View Article

15	Rodriguez S, Gaunt TR, Day IN. Hardy-Weinberg equilibrium testing of biological ascertainment for Mendelian randomization studies. Am J Epidemiol 2009;169:505-514 View Article

16	DerSimonian R, Laird N. Meta-analysis in clinical trials. Control Clin Trials 1986;7:177-188 View Article

17	Higgins JP, Thompson SG, Deeks JJ, Altman DG. Measuring inconsistency in meta-analyses. BMJ 2003;327:557-560 View Article

18	Deeks JJ, Altman DG, Bradburn MJ. Systematic reviews in health care: meta-analysis in context. 2nd edition. London: BMJ Books; 2005, p.285-312

19	Begg CB, Mazumdar M. Operating characteristics of a rank correlation test for publication bias. Biometrics 1994;50:1088-1101 View Article

20	Egger M, Davey Smith G, Schneider M, Minder C. Bias in meta-analysis detected by a simple, graphical test. BMJ 1997;315:629-634 View Article

21	Duval S, Tweedie R. Trim and fill: A simple funnel-plot-based method of testing and adjusting for publication bias in meta-analysis. Biometrics 2000;56:455-463 View Article

22	Available at: http://finzi.psych.upenn.edu/library/metafor/html/trimfill.html

23	Hamza S, Petit JM, Masson D, Jooste V, Binquet C, Sgro C. PNPLA3 rs738409 GG homozygote status is associated with increased risk of hepatocellular carcinoma in cirrhotic patients. J Hepatol 2012;56:S281-S282 View Article

24	Way M, McQuillin A, Gurling HMD, Morgan MY. The PNPLA3 I148M mutation significantly increases the risk of developing alcohol-related cirrhosis in alcohol-dependent individuals. J Hepatol 2013;58:S563-S564 View Article

25	Dutta AK. Genetic factors affecting susceptibility to alcoholic liver disease in an Indian population. Ann Hepatol 2013;12:901-907

26	Nguyen-Khac E, Houchi H, Dreher ML, Herpe YE, Naassila M. Is PNPLA3 polymorphism involved in severe acute alcoholic hepatitis. Hepatology 2011;54:976A

27	Bhatt SP, Nigam P, Misra A, Guleria R, Pandey RM, Pasha MA. Genetic variation in the patatin-like phospholipase domain-containing protein-3 (PNPLA-3) gene in Asian Indians with nonalcoholic fatty liver disease. Metab Syndr Relat Disord 2013;11:329-335 View Article

28	Burza MA, Pirazzi C, Maglio C, Sjöholm K, Mancina RM, Svensson PA. PNPLA3 I148M (rs738409) genetic variant is associated with hepatocellular carcinoma in obese individuals. Dig Liver Dis 2012;44:1037-1041 View Article

29	Romeo S, Kozlitina J, Xing C, Pertsemlidis A, Cox D, Pennacchio LA. Genetic variation in PNPLA3 confers susceptibility to nonalcoholic fatty liver disease. Nat Genet 2008;40:1461-1465 View Article

30	Santoro N, Kursawe R, D’Adamo E, Dykas DJ, Zhang CK, Bale AE. A common variant in the patatin-like phospholipase 3 gene (PNPLA3) is associated with fatty liver disease in obese children and adolescents. Hepatology 2010;52:1281-1290 View Article

31	Kantartzis K, Peter A, Machicao F, Machann J, Wagner S, Königsrainer I. Dissociation between fatty liver and insulin resistance in humans carrying a variant of the patatin-like phospholipase 3 gene. Diabetes 2009;58:2616-2623 View Article

32	Sookoian S, Castaño GO, Burgueño AL, Gianotti TF, Rosselli MS, Pirola CJ. A nonsynonymous gene variant in the adiponutrin gene is associated with nonalcoholic fatty liver disease severity. J Lipid Res 2009;50:2111-2116 View Article

33	Petit JM, Guiu B, Masson D, Duvillard L, Jooste V, Buffier P. Specifically PNPLA3-mediated accumulation of liver fat in obese patients with type 2 diabetes. J Clin Endocrinol Metab 2010;95:E430-E436 View Article

34	Valenti L, Al-Serri A, Daly AK, Galmozzi E, Rametta R, Dongiovanni P. Homozygosity for the patatin-like phospholipase-3/adiponutrin I148M polymorphism influences liver fibrosis in patients with nonalcoholic fatty liver disease. Hepatology 2010;51:1209-1217 View Article

35	Rotman Y, Koh C, Zmuda JM, Kleiner DE, Liang TJ, NASH CRN. The association of genetic variability in patatin-like phospholipase domain-containing protein 3 (PNPLA3) with histological severity of nonalcoholic fatty liver disease. Hepatology 2010;52:894-903 View Article

36	Goran MI, Walker R, Le KA, Mahurkar S, Vikman S, Davis JN. Effects of PNPLA3 on liver fat and metabolic profile in Hispanic children and adolescents. Diabetes 2010;59:3127-3130 View Article

37	Speliotes EK, Butler JL, Palmer CD, Voight BF. PNPLA3 variants specifically confer increased risk for histologic nonalcoholic fatty liver disease but not metabolic disease. Hepatology 2010;52:904-912 View Article

38	Valenti L, Alisi A, Galmozzi E, Bartuli A, Del Menico B, Alterio A. I148M patatin-like phospholipase domain-containing 3 gene variant and severity of pediatric nonalcoholic fatty liver disease. Hepatology 2010;52:1274-1280 View Article

39	Wang CW, Lin HY, Shin SJ, Yu ML, Lin ZY, Dai CY. The PNPLA3 I148M polymorphism is associated with insulin resistance and nonalcoholic fatty liver disease in a normoglycaemic population. Liver Int 2011;31:1326-1331 View Article

40	Hotta K, Yoneda M, Hyogo H, Ochi H, Mizusawa S, Ueno T. Association of the rs738409 polymorphism in PNPLA3 with liver damage and the development of nonalcoholic fatty liver disease. BMC Med Genet 2010;11:172 View Article

41	Lin YC, Chang PF, Hu FC, Yang WS, Chang MH, Ni YH. A common variant in the PNPLA3 gene is a risk factor for non-alcoholic fatty liver disease in obese Taiwanese children. J Pediatr 2011;158:740-744 View Article

42	Peng XE, Wu YL, Lin SW, Lu QQ, Hu ZJ, Lin X. Genetic variants in PNPLA3 and risk of non-alcoholic fatty liver disease in a Han Chinese population. PLoS One 2012;7:e50256 View Article

43	Kawaguchi T, Sumida Y, Umemura A, Matsuo K, Takahashi M, Takamura T. Genetic polymorphisms of the human PNPLA3 gene are strongly associated with severity of non-alcoholic fatty liver disease in Japanese. PLoS One 2012;7:e38322 View Article

44	Petta S, Grimaudo S, Cammà C, Cabibi D, Di Marco V, Licata G. IL28B and PNPLA3 polymorphisms affect histological liver damage in patients with non-alcoholic fatty liver disease. J Hepatol 2012;56:1356-1362 View Article

45	Zain SM, Mohamed R, Mahadeva S, Cheah PL, Rampal S, Basu RC. A multi-ethnic study of a PNPLA3 gene variant and its association with disease severity in non-alcoholic fatty liver disease. Hum Genet 2012;131:1145-1152 View Article

46	Li Y, Xing C, Cohen JC, Hobbs HH. Genetic variant in PNPLA3 is associated with nonalcoholic fatty liver disease in China. Hepatology 2012;55:327-328 View Article

47	Kitamoto T, Kitamoto A, Yoneda M, Hyogo H, Ochi H, Nakamura T. Genome-wide scan revealed that polymorphisms in the PNPLA3, SAMM50, and PARVB genes are associated with development and progression of nonalcoholic fatty liver disease in Japan. Hum Genet 2013;132:783-792 View Article

48	Islek EE, Sazci A, Ozel MD, Aygun C. Genetic variants in the PNPLA3 gene are associated with nonalcoholic steatohepatitis. Genet Test Mol Biomarkers 2014;18:489-496 View Article

49	Guichelaar MM, Gawrieh S, Olivier M, Viker K, Krishnan A, Sanderson S. Interactions of allelic variance of PNPLA3 with nongenetic factors in predicting nonalcoholic steatohepatitis and nonhepatic complications of severe obesity. Obesity (Silver Spring) 2013;21:1935-1941 View Article

50	Tai CM, Huang CK, Tu HP, Hwang JC, Chang CY, Yu ML. PNPLA3 genotype increases susceptibility of nonalcoholic steatohepatitis among obese patients with nonalcoholic fatty liver disease. Surg Obes Relat Dis 2015;11:888-894 View Article

51	Kanth VV, Sasikala M, Rao PN, Steffie Avanthi U, Rao KR, Nageshwar Reddy D. Pooled genetic analysis in ultrasound measured non-alcoholic fatty liver disease in Indian subjects: A pilot study. World J Hepatol 2014;6:435-442 View Article

52	Lin YC, Chang PF, Chang MH, Ni YH. Genetic variants in GCKR and PNPLA3 confer susceptibility to nonalcoholic fatty liver disease in obese individuals. Am J Clin Nutr 2014;99:869-874 View Article

53	Liu YL, Patman GL, Leathart JB, Piguet AC, Burt AD, Dufour JF. Carriage of the PNPLA3 rs738409 C>G polymorphism confers an increased risk of non-alcoholic fatty liver disease associated hepatocellular carcinoma. J Hepatol 2014;61:75-81 View Article

54	Verrijken A, Beckers S, Francque S, Hilden H, Caron S, Zegers D. A gene variant of PNPLA3, but not of APOC3, is associated with histological parameters of NAFLD in an obese population. Obesity (Silver Spring) 2013;21:2138-2145 View Article

55	Wagenknecht LE, Palmer ND, Bowden DW, Rotter JI, Norris JM, Ziegler J. Association of PNPLA3 with non-alcoholic fatty liver disease in a minority cohort: the Insulin Resistance Atherosclerosis Family Study. Liver Int 2011;31:412-416 View Article

56	Akuta N, Kawamura Y, Arase Y, Suzuki F, Sezaki H, Hosaka T. Relationships between genetic variations of PNPLA3, TM6SF2 and histological features of nonalcoholic fatty liver disease in Japan. Gut Liver 2016;10:437-445 View Article

57	Baulande S, Lasnier F, Lucas M, Pairault J. Adiponutrin, a transmembrane protein corresponding to a novel dietary- and obesity-linked mRNA specifically expressed in the adipose lineage. J Biol Chem 2001;276:33336-33344 View Article

58	Wilson PA, Gardner SD, Lambie NM, Commans SA, Crowther DJ. Characterization of the human patatin-like phospholipase family. J Lipid Res 2006;47:1940-1949 View Article

59	He S, McPhaul C, Li JZ, Garuti R, Kinch L, Grishin NV. A sequence variation (I148M) in PNPLA3 associated with nonalcoholic fatty liver disease disrupts triglyceride hydrolysis. J Biol Chem 2010;285:6706-6715 View Article

60	Chamoun Z, Vacca F, Parton RG, Gruenberg J. PNPLA3/adiponutrin functions in lipid droplet formation. Biol Cell 2013;105:219-233 View Article

61	Paré G, Ridker PM, Rose L, Barbalic M, Dupuis J, Dehghan A. Genome-wide association analysis of soluble ICAM-1 concentration reveals novel associations at the NFKBIK, PNPLA3, RELA, and SH2B3 loci. PLoS Genet 2011;7:e1001374 View Article

62	Qu Z, Wu M, Xie T. Expression of ICAM-1 mRNA in HCC using in situ hybridization technique. Zhonghua Bing Li Xue Za Zhi 1997;26:82-84

63	Roland CL, Harken AH, Sarr MG, Barnett CC. ICAM-1 expression determines malignant potential of cancer. Surgery 2007;141:705-707 View Article

64	Trépo E, Nahon P, Bontempi G, Valenti L, Falleti E, Nischalke HD. Association between the PNPLA3 (rs738409 C>G) variant and hepatocellular carcinoma: Evidence from a meta-analysis of individual participant data. Hepatology 2014;59:2170-2177 View Article

65	Edmondson HA, Peters RL, Frankel HH, Borowsky S. The early stage of liver injury in the alcoholic. Medicine (Baltimore) 1967;46:119-129 View Article

66	Naveau S, Giraud V, Borotto E, Aubert A, Capron F, Chaput JC. Excess weight risk factor for alcoholic liver disease. Hepatology 1997;25:108-111 View Article

67	Crabb DW. Pathogenesis of alcoholic liver disease: newer mechanisms of injury. Keio J Med 1999;48:184-188 View Article

68	Thornton A, Lee P. Publication bias in meta-analysis: its causes and consequences. J Clin Epidemiol 2000;53:207-216 View Article

69	Bedossa P. Histological assessment of NAFLD. Dig Dis Sci 2016;61:1348-1355 View Article

70	Bedossa P, Poitou C, Veyrie N, Bouillot JL, Basdevant A, Paradis V. Histopathological algorithm and scoring system for evaluation of liver lesions in morbidly obese patients. Hepatology 2012;56:1751-1759 View Article

71	Kleiner DE, Brunt EM, Van Natta M, Behling C, Contos MJ, Cummings OW. Design and validation of a histological scoring system for nonalcoholic fatty liver disease. Hepatology 2005;41:1313-1321 View Article

72	Baclig MO, Lozano-Kühne JP, Mapua CA, Gopez-Cervantes J, Natividad FF. Genetic variation I148M in patatin-like phospholipase 3 gene and risk of non-alcoholic fatty liver disease among Filipinos. Int J Clin Exp Med 2014;7:2129-2136

73	Lee SS, Byoun YS, Jeong SH, Woo BH, Jang ES, Kim JW. Role of the PNPLA3 I148M polymorphism in nonalcoholic fatty liver disease and fibrosis in Korea. Dig Dis Sci 2014;59:2967-2974 View Article

74	Park JH, Cho B, Kwon H, Prilutsky D, Yun JM, Choi HC. I148M variant in PNPLA3 reduces central adiposity and metabolic disease risks while increasing nonalcoholic fatty liver disease. Liver Int 2015;35:2537-2546 View Article

75	Wei JL, Leung JC, Loong TC, Wong GL, Yeung DK, Chan RS. Prevalence and severity of nonalcoholic fatty liver disease in non-obese patients: a population study using proton-magnetic resonance spectroscopy. Am J Gastroenterol 2015;110:1306-1314

76	Xia MF, Ling Y, Bian H, Lin HD, Yan HM, Chang XX. I148M variant of PNPLA3 increases the susceptibility to non-alcoholic fatty liver disease caused by obesity and metabolic disorders. Aliment Pharmacol Ther 2016;43:631-642 View Article

77	Oniki K, Saruwatari J, Izuka T, Kajiwara A, Morita K, Sakata M. Influence of the PNPLA3 rs738409 polymorphism on non-alcoholic fatty liver disease and renal function among normal weight subjects. PLoS One 2015;10:e0132640 View Article

78	Nishioji K, Mochizuki N, Kobayashi M, Kamaguchi M, Sumida Y, Nishimura T. The Impact of PNPLA3 rs738409 genetic polymorphism and weight gain ≥ 10 kg after age 20 on non-alcoholic fatty liver disease in non-obese Japanese individuals. PLoS One 2015;10:e0140427 View Article

79	Tai CM, Huang CK, Tu HP, Hwang JC, Chang CY, Yu ML. PNPLA3 genotype increases susceptibility of nonalcoholic steatohepatitis among obese patients with nonalcoholic fatty liver disease. Surg Obes Relat Dis 2015;11:888-894 View Article

80	Tai CM, Huang CK, Tu HP, Hwang JC, Yeh ML, Huang CF. Interactions of a PPARGC1A Variant and a PNPLA3 variant affect nonalcoholic steatohepatitis in severely obese Taiwanese patients. Medicine (Baltimore) 2016;95:e3120 View Article

81	Honda Y, Yoneda M, Kessoku T, Ogawa Y, Tomeno W, Imajo K. The characteristics of non-obese NAFLD: Effect of genetic and environmental factors. Hepatol Res 2016 View Article

Copyright © 2016 The Second Affiliated Hospital of Chongqing Medical University. Published by XIA & HE Publishing Inc. All rights reserved. This is an Open Access article distributed under the terms of the Creative Commons Attribution-Noncommercial 4.0 License (CC BY-NC 4.0), permitting all non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.