Correlation between PNPLA3 rs738409 polymorphism and hepatocellular carcinoma: a meta-analysis of 10,330 subjects

Abstract

Purpose:

The correlation between patatin-like phospholipase domain-containing protein 3 (PNPLA3) rs738409 polymorphism and hepatocellular carcinoma was investigated by several pilot studies, but the results of these studies were controversial. Therefore, we performed this study to better assess the relationship between PNPLA3 rs738409 polymorphism and the likelihood of hepatocellular carcinoma.

Methods:

Eligible studies were searched in PubMed, Medline, EMBASE, and Web of Science. Odds ratios (ORs) and 95% confidence intervals (CIs) were calculated to assess the relationship between PNPLA3 rs738409 polymorphism and hepatocellular carcinoma.

Results:

A total of 17 studies with 10,330 participants were analyzed. A significant association with the likelihood of hepatocellular carcinoma was detected for the PNPLA3 rs738409 polymorphism in dominant (P = 0.0001; OR 0.66; 95% CI 0.53, 0.82), recessive (P < 0.0001; OR 2.32; 95% CI 1.76, 3.06) and allele (P < 0.0001; OR 0.64; 95% CI 0.53, 0.77) comparisons. Further subgroup analyses revealed that the PNPLA3 rs738409 polymorphism was significantly associated with the likelihood of hepatocellular carcinoma in Caucasians (dominant model: P < 0.0001, OR 0.57, 95% CI 0.45, 0.71; recessive model: P < 0.0001, OR 2.74, 95% CI 2.02, 3.71; allele model: P < 0.0001, OR 0.56, 95% CI 0.46, 0.67). However, no positive results were detected in Asians.

Conclusions:

Our findings indicated that the PNPLA3 rs738409 polymorphism may serve as a potential biological marker of hepatocellular carcinoma in Caucasians.

Keywords

gene polymorphism hepatocellular carcinoma (HCC)meta-analysis

Introduction

Hepatocellular carcinoma (HCC) refers to malignancy that originates from hepatocytes. It is the fifth most commonly seen cancer worldwide.¹ Moreover, despite rapid advances in early diagnosis and surgical treatment over the past few decades, HCC still accounts for approximately 9% of cancer-related deaths, making it the second most frequent cause of cancer-related mortality globally.² To date, the exact cause of HCC remains unclear. Although excessive alcohol intake and chronic viral infection have been identified as potential pathogenic factors of HCC,^3-4 the fact that a great inter-individual variability in disease susceptibility existed in those exposed to the abovementioned carcinogenic factors suggests that genetic factors also a play crucial part in the pathogenesis of HCC.

Patatin-like phospholipase domain-containing protein 3 (PNPLA3) is a multi-functional enzyme encoded by the human PNPLA3 gene located on chromosome 22.⁵ Previous experimental studies have shown that PNPLA3 could regulate metabolism of carbohydrates and lipids in the liver,^6-7 and altered expression levels of PNPLA3 have also been proved to be correlated with multiple chronic liver diseases, such as alcoholic liver disease and nonalcoholic fatty liver disease.^8-9 Consequently, functional PNPLA3 polymorphisms were thought to be implicated in the development of end-stage liver diseases like HCC.

The rs738409 polymorphism is located on exon 3 of PNPLA3 gene; the C to G substitution at this locus could lead to reduced enzymatic function, and this variant has been consistently shown to be associated with liver steatosis and fibrosis.^10-11 So far, some pilot studies have been conducted to investigate the possible correlation between PNPLA3 rs738409 polymorphism and the likelihood of HCC. But the results of these studies were controversial and the sample size of individual studies was inadequate to draw a definite conclusion.^12-14 Therefore, we conducted this meta-analysis to better analyze the role of PNPLA3 rs738409 polymorphism in HCC.

Materials and methods

Literature search and inclusion criteria

The current meta-analysis was performed according to the Preferred Reporting Items for Systematic Reviews and Meta-analyses (PRISMA) statement.¹⁵ Potentially related articles were searched in PubMed, Medline, and EMBASE using the combination of following key words: “patatin-like phospholipase domain-containing 3,” “PNPLA3,” “adiponutrin,” “ADPN,” “polymorphism,” “variant,” “mutation,” “genotype,” “allele,” “hepatocellular carcinoma,” and “HCC.” The initial literature search was conducted in March 2018 and the latest update was performed in August 2018. The reference lists of all retrieved publications were also screened to identify other potentially relevant articles.

To test the research hypothesis of this meta-analysis, included studies should meet all the following criteria: (a) case-control study on the correlation between PNPLA3 rs738409 polymorphism and the likelihood of HCC; (b) provide adequate data to calculate odds ratios (ORs) and 95% confidence intervals (CIs); and (c) full text in English is available. Studies were excluded if one of the following criteria was fulfilled: (a) not relevant to PNPLA3 rs738409 polymorphism and HCC; (b) family-based association studies; (c) case reports or case series; and (d) abstracts, reviews, comments, letters, and conference presentations. For duplicate reports, only the study with the largest sample size was included.

Data extraction and quality assessment

The following data were extracted from all included studies: (a) name of first author; (b) year of publication; (c) country and ethnicity of participants; (d) the number of cases and controls; and (e) the genotypic distribution of PNPLA3 rs738409 polymorphism in cases and controls. Additionally, the probability value (P value) of Hardy–Weinberg equilibrium (HWE) was also calculated on the basis of genotypic frequency of PNPLA3 rs738409 polymorphism in the control group.

The Newcastle–Ottawa scale (NOS) was used to assess the quality of eligible studies.¹⁶ The NOS has a score range of zero to nine, and studies with a score of more than seven were thought to be of high quality.

Two reviewers conducted data extraction and quality assessment independently. When necessary, the reviewers wrote to the corresponding authors for extra information or raw data. Any disagreement between two reviewers was solved by discussion until a consensus was reached.

Statistical analysis

All statistical analyses were conducted with Review Manager Version 5.3.3 (The Cochrane Collaboration, Software Update, Oxford, UK). ORs and 95% CIs were calculated to assess the potential association of PNPLA3 rs738409 polymorphism with the likelihood of HCC in the dominant (CC vs. CG + GG), recessive (GG vs. CC + CG), additive (CG vs. CC + GG), and allele (C vs. G) models, and a P value of 0.05 or less was considered to be statistically significant. Between-study heterogeneities were evaluated based on the Q test and the I² statistic. If the P value of the Q test was less than 0.1, or if the I² was greater than 50%, random-effect models (REMs) would be used for analyses on account of obvious between-study heterogeneities. Otherwise, fixed-effect models (FEMs) would be applied for analyses. Subgroup analyses by ethnicity of participants were subsequently conducted to obtain more specific results. Sensitivity analyses were carried out to test the stability of synthetic results. Funnel plots were applied to evaluate possible publication bias.

Results

Characteristics of included studies

The literature search generated 675 results. After exclusion of irrelevant and duplicate articles by reading titles and abstracts, 38 articles were retrieved for further evaluation. Another 21 articles were subsequently excluded after reading the full text. Finally, a total of 17 eligible studies containing 3154 cases and 7176 controls were included for analyses (see Figure 1). All eligible studies were published in English between 2011 and 2018. The NOS score of eligible articles ranged from 7 to 8, which suggested that all included studies were of relatively high quality. The characteristics of included studies are summarized in Table 1.

Figure 1.

Flowchart of study selection for the present study.

Table 1.

The characteristics of included studies for PNPLA3 rs738409 polymorphism and HCC.

First author, year	Country	Ethnicity	Sample size	Genotype distribution		P-value for HWE	NOS score
First author, year	Country	Ethnicity	Sample size	Cases	Controls	P-value for HWE	NOS score
Balasus¹⁷ 2016	Italy	Caucasian	170/304	71/64/35	148/128/28	0.966	8
Casper¹⁸ 2016	Germany	Caucasian	44/174	13/18/13	104/62/8	0.747	8
Ezzikouri et al. ¹³ 2014	Morocco	Caucasian	101/132	43/35/23	66/53/12	0.842	8
Falleti et al.¹² 2016	Italy	Caucasian	150/361	45/64/41	131/170/60	0.700	8
Friedrich¹⁹ 2014	Germany	Caucasian	421/1950	202/177/42	1135/718/97	0.223	8
Gao²⁰ 2017	China	Asian	408/358	166/185/57	143/173/42	0.344	8
Guyot²¹ 2013	France	Caucasian	93/160	54/26/13	86/49/25	<0.001	7
Hai²² 2017	Japan	Asian	131/545	94/35/2	406/130/9	0.702	8
Hassan²³ 2013	USA	Mixed	257/494	130/94/33	274/197/23	0.094	8
Khlaiphuengsin²⁴ 2015	Thailand	Asian	200/388	91/88/21	158/178/52	0.867	8
Liu²⁵ 2014	UK	Caucasian	100/275	28/43/29	125/117/33	0.488	8
Nakaoka et al.¹⁴ 2015	Japan	Asian	48/183	12/22/14	60/96/27	0.253	7
Nischalke²⁶ 2011	Germany	Caucasian	161/190	57/73/31	112/69/9	0.693	8
Seko²⁷ 2017	Japan	Asian	10/120	7/3/0	45/57/18	0.994	8
Stickel²⁸ 2018	Switzerland	Caucasian	751/1165	205/346/200	508/505/152	0.135	8
Ueyama²⁹ 2016	Japan	Asian	59/198	6/23/30	76/78/44	0.007	7
Valenti³⁰ 2011	Italy	Caucasian	50/179	17/21/12	118/56/5	0.590	8

HCC, hepatocellular carcinoma; HWE, Hardy-Weinberg equilibrium; NA, not available; NOS, Newcastle–Ottawa scale; PNPLA3, patatin-like phospholipase domain-containing 3.

Overall and subgroup analysis

A total of 17 studies with 10,330 participants were analyzed. A significant association with the likelihood of HCC was detected for the PNPLA3 rs738409 polymorphism in dominant (P = 0.0001; OR = 0.66; 95% CI 0.53, 0.82), recessive (P < 0.0001; OR = 2.32; 95% CI 1.76, 3.06), and allele (P < 0.0001; OR = 0.64; 95% CI 0.53, 0.77) comparisons. Further subgroup analyses according to ethnicity of participants revealed that the PNPLA3 rs738409 polymorphism was significantly associated with the likelihood of HCC in Caucasians (dominant model: P < 0.0001, OR = 0.57, 95% CI 0.45, 0.71; recessive model: P < 0.0001, OR = 2.74, 95% CI 2.02, 3.71; allele model: P < 0.0001, OR = 0.56, 95% CI 0.46, 0.67). However, no positive results were detected in Asians (see Table 2 and Supplementary file 1).

Table 2.

Results of overall and subgroup analyses for PNPLA3 rs738409 polymorphism and HCC.

Population	Sample size	Dominant comparison			Recessive comparison			Additive comparison			Allele comparison
		P value	OR (95%CI)	I² statistic	P value	OR (95%CI)	I² statistic	P value	OR (95%CI)	I² statistic	P value	OR (95%CI)	I² statistic
Overall	3154/7176	0.0001	0.66(0.53, 0.82)	79%	<0.0001	2.32 (1.76, 3.06)	73%	0.57 1.	03 (0.94, 1.12)	12%	<0.0001	0.64 (0.53, 0.77)	85%
Asian	856/1792	0.50	0.86 (0.56, 1.33)	76%	0.45	1.45 (0.80, 2.62)	73%	0.39 0.	93 (0.78, 1.11)	0%	0.37	0.83 (0.56, 1.24)	86%
Caucasian	2041/4890	<0.0001	0.57(0.45, 0.71)	70%	<0.0001	2.74 (2.02, 3.71)	67%	0.13 1.	09 (0.98, 1.21)	26%	<0.0001	0.56 (0.46, 0.67)	79%

CI, confidence interval; HCC, hepatocellular carcinoma; NA, not available; OR, odds ratio; PNPLA3, patatin-like phospholipase domain-containing 3.

The values in bold represent statistically significant differences between cases and controls.

Sensitivity analysis

Sensitivity analyses were carried out to examine the stability of synthetic results by eliminating studies that deviated from HWE. No changes of results were observed in any comparisons, which indicated that our findings were statistically reliable.

Publication bias

Potential publication bias in the current study was evaluated with funnel plots. No obvious asymmetry of funnel plots was observed in any comparisons, which suggested that our findings were unlikely to be influenced by a severe publication bias.

Discussion

To the best of our knowledge, so far this is the largest meta-analysis on the correlation between the PNPLA3 rs738409 polymorphism and HCC. Our overall and subgroup analyses revealed that the PNPLA3 rs738409 polymorphism was significantly associated with the likelihood of HCC in Caucasians. The stability of synthetic results was subsequently evaluated in sensitivity analyses, and no changes of results were observed in any comparisons, which indicated that our findings were quite stable and reliable. As for evaluation of heterogeneities, obvious between-study heterogeneities were detected in dominant, recessive, and allele comparisons, and thus analyses were performed with REMs. In further stratified analyses, a reduction tendency of heterogeneity was found in the Caucasian subgroup, which suggested that differences in ethnic background could partially explain observed heterogeneities between studies.

There are several points that need to be addressed regarding this meta-analysis. First, previous experimental studies have shown that the mutant allele of rs738409 polymorphism was correlated with reduced enzymatic activity and hepatic steatosis, and this may partially explain positive findings in the current study.^6-9 Second, the pathogenic mechanism of HCC is highly complex, and hence it is unlikely that a single gene polymorphism can significantly contribute to its development. Therefore, to better illustrate the potential correlations of certain gene polymorphisms with HCC, we strongly recommend further studies to perform haplotype analyses and to explore potential gene–gene interactions.

As with all meta-analyses, this study certainly has some limitations. First, our results were derived from unadjusted analyses due to lack of raw data, and failure to conduct further stratified analyses according to age, gender, and co-morbidity conditions may impact the reliability of our findings.³¹ Second, obvious heterogeneities were detected in certain subgroup comparisons, which indicated that the inconsistent results of included studies could not be fully explained by the differences in ethnic background. Other unmeasured characteristics of participants may also contribute to between-study heterogeneities.³² Third, the association between the PNPLA3 rs738409 polymorphism and HCC also may be influenced by gene–gene and gene–environmental interactions. However, the majority of studies did not consider these potential interactions, which accordingly hindered the performance of relevant analyses.³³ Considering these limitations, the results of the current study should be interpreted with caution.

Overall, our meta-analysis suggests that the PNPLA3 rs738409 polymorphism may serve as a potential biological marker of HCC in Caucasians. However, further well-designed studies are warranted to confirm our findings. Additionally, future investigations are needed to explore the potential roles of other PNPLA3 gene polymorphisms in the development of HCC.

Supplemental Material

Supplementary_file_1-Forest_plots_for_PNPLA3_rs738409_polymorphism – Supplemental material for Correlation between PNPLA3 rs738409 polymorphism and hepatocellular carcinoma: a meta-analysis of 10,330 subjects

Supplemental material, Supplementary_file_1-Forest_plots_for_PNPLA3_rs738409_polymorphism for Correlation between PNPLA3 rs738409 polymorphism and hepatocellular carcinoma: a meta-analysis of 10,330 subjects by Zongsheng Huang, Xianwen Guo, Guo Zhang, Liexin Liang and Bing Nong in The International Journal of Biological Markers

Footnotes

Authors’ contributions

Zongsheng Huang and Bing Nong conceived the study and participated in its design. Zongsheng Huang and Xianwen Guo conducted the systematic literature review. Guo Zhang and Liexin Liang performed the data analyses. Zongsheng Huang and Bing Nong drafted the manuscript. All authors have read and approved the final manuscript.

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Ethical approval

This article does not contain any studies with human participants or animals performed by any of the authors.

Funding

This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors.

Informed consent

For this type of study formal consent is not required.

Supplemental material

Supplemental material for this article is available online.

References

Siegel

Zou

et al . Cancer statistics, 2014. CA Cancer J Clin 2014; 64: 9–29.

Ferlay

Soerjomataram

Dikshit

et al . Cancer incidence and mortality worldwide: Sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer 2015; 136: E359–E386.

Forner

Llovet

Bruix

Hepatocellular carcinoma. Lancet 2012; 379: 1245–1255.

El-Serag

HB.

Hepatocellular carcinoma. N Engl J Med 2011; 365: 1118–1127.

Collins

Goward

Cole

et al . Reevaluating human gene annotation: a second-generation analysis of chromosome 22. Genome Res 2003; 13: 27–36.

Wilson

Gardner

Lambie

et al . Characterization of the human patatin-like phospholipase family. J Lipid Res 2006; 47: 1940–1949.

Dubuquoy

Robichon

Lasnier

et al . Distinct regulation of adiponutrin/PNPLA3 gene expression by the transcription factors ChREBP and SREBP1c in mouse and human hepatocytes. J Hepatol 2011; 55: 145–153.

Trépo

Romeo

Zucman-Rossi

et al . PNPLA3 gene in liver diseases. J Hepatol 2016; 65: 399–412.

Trépo

Contribution of PNPLA3 gene to the natural history of liver diseases. Acta Gastroenterol Belg 2017; 80: 43–51.

10.

Chen

Xin

Geng

et al . PNPLA3 I148M variant in nonalcoholic fatty liver disease: demographic and ethnic characteristics and the role of the variant in nonalcoholic fatty liver fibrosis. World J Gastroenterol 2015; 21: 794–802.

11.

Dongiovanni

Donati

Fares

et al . PNPLA3 I148M polymorphism and progressive liver disease. World J Gastroenterol 2013; 19: 6969–6978.

12.

Falleti

Cussigh

Cmet

et al . PNPLA3 rs738409 and TM6SF2 rs58542926 variants increase the risk of hepatocellular carcinoma in alcoholic cirrhosis. Dig Liver Dis 2016; 48: 69–75.

13.

Ezzikouri

Alaoui

Tazi

et al . The adiponutrin I148M variant is a risk factor for HCV-associated liver cancer in North-African patients. Infect Genet Evol 2014; 21: 179–183.

14.

Nakaoka

Hashimoto

Kawabe

et al . PNPLA3 I148M associations with liver carcinogenesis in Japanese chronic hepatitis C patients. Springerplus 2015; 4: 83.

15.

Moher

Liberati

Tetzlaff

et al . Preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement. Ann Intern Med 2009; 151: 264–269.

16.

Stang

Critical evaluation of the Newcastle–Ottawa scale for the assessment of the quality of nonrandomized studies in meta-analyses. Eur J Epidemiol 2010; 25: 603–605.

17.

Balasus

Way

Fusilli

et al . The association of variants in PNPLA3 and GRP78 and the risk of developing hepatocellular carcinoma in an Italian population[J]. Oncotarget 2016; 7(52): 86791–86802.

18.

Casper

Krawczyk

Behrmann

et al . Variant PNPLA3 increases the HCC risk: Prospective study in patients treated at the Saarland University Medical Center[J]. Z Gastroenterol 2016; 54(6): 585–586.

19.

Friedrich

Wannhoff

Kattner

et al

PNPLA3 in end-stage liver disease: Alcohol consumption, hepatocellular carcinoma development, and transplantation-free survival.[J]. J Gastroenterol Hepatol, 2014; 29(7): 1477–1484.

20.

Gao

Liu

Yang

et al

Association between PNPLA3 gene polymorphisms and risk of hepatitis B virus-related hepatocellular carcinoma in Han population in China: A case-control study[J]. Scand J Gastroenterol 2017: 1–8.

21.

Guyot

Sutton

Rufat

et al . PNPLA3 rs738409, hepatocellular carcinoma occurrence and risk model prediction in patients with cirrhosis[J]. J Hepatol 2013; 58(2): 312–318.

22.

Hai

Tamori

Thuy

LTT

et al . Polymorphisms in MICA, but not in DEPDC5, HCP5 or PNPLA3, are associated with chronic hepatitis C-related hepatocellular carcinoma[J]. Sci Rep 2017; 7(1): 11912.

23.

Hassan

Kaseb

Etzel

et al . Genetic variation in the\r, PNPLA3\r, gene and hepatocellular carcinoma in USA: Risk and prognosis prediction[J]. Mol Carcinog 2013; 52(S1): 139–147.

24.

Khlaiphuengsin

Kiatbumrung

Payungporn

et al

Association of PNPLA3 Polymorphism with Hepatocellular Carcinoma Development and Prognosis in Viral and Non-Viral Chronic Liver Diseases.[J]. Asian Pac J Cancer Prev 2015; 16(18): 8377–8382.

25.

Liu

Patman

Leathart

JBS

et al . Carriage of the PNPLA3 rs738409 C >G polymorphism confers an increased risk of non-alcoholic fatty liver disease associated hepatocellular carcinoma[J]. J Hepatol 2014; 61(1): 75–81.

26.

Nischalke

Berger

Luda

et al . The PNPLA3 rs738409 148M/M Genotype Is a Risk Factor for Liver Cancer in Alcoholic Cirrhosis but Shows No or Weak Association in Hepatitis C Cirrhosis[J]. PLOS ONE 2011; 6(11): e27087.

27.

Seko

Sumida

Tanaka

et al . Development of hepatocellular carcinoma in Japanese patients with biopsy-proven nonalcoholic fatty liver disease: In relation to PNPLA3 genotype and hepatocarcinogenesis/fibrosis progression[J]. Hepatol Res 2017.

28.

Stickel

Buch

Nischalke

et al . Genetic variants in PNPLA3 and TM6SF2 predispose to the development of hepatocellular carcinoma in individuals with alcohol-related cirrhosis[J]. Am J Gastroenterol 2018.

29.

Ueyama

Nishida

Korenaga

et al . The impact of PNPLA3 and JAZF1 on hepatocellular carcinoma in non-viral hepatitis patients with type 2 diabetes mellitus[J]. J Gastroenterol 2016; 51(4): 370–379.

30.

Valenti

Rumi

Galmozzi

et al . Patatin-Like phospholipase domain-containing 3 I148M polymorphism, steatosis, and liver damage in chronic hepatitis C[J]. Hepatology 2011; 53(3): 791–799.

31.

Xie

Shi

Liu

The roles of TLR gene polymorphisms in atherosclerosis: a systematic review and meta-analysis of 35,317 subjects. Scand J Immunol 2017; 86: 50–58.

32.

Shi

Xie

Jia

et al . Associations of insulin receptor and insulin receptor substrates genetic polymorphisms with polycystic ovary syndrome: A systematic review and meta-analysis. J Obstet Gynaecol Res 2016; 42: 844–854.

33.

Xie

Shi

Xun

et al . Endothelial nitric oxide synthase gene single nucleotide polymorphisms and the risk of hypertension: A meta-analysis involving 63,258 subjects. Clin Exp Hypertens 2017; 39: 175–182.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

0.17 MB