Genome-wide association study of non-alcoholic fatty liver and steatohepatitis in a histologically characterised cohort☆
Authors: Quentin M. Anstee, Rebecca Darlay, Simon Cockell, Marica Meroni, Olivier Govaere, Dina Tiniakos, Alastair D. Burt, Pierre Bedossa, Jeremy Palmer, Yang-Lin Liu, Guruprasad P. Aithal, Michael Allison, Hannele Yki-Järvinen, Michele Vacca, Jean-Francois Dufour, Pietro Invernizzi, Daniele Prati, Mattias Ekstedt, Stergios Kechagias, Sven Francque, Salvatore Petta, Elisabetta Bugianesi, Karine Clement, Vlad Ratziu, Jörn M. Schattenberg, Luca Valenti, Christopher P. Day, Heather J. Cordell, Ann K. Daly
Journal: Journal of Hepatology (2020)

•Genome-wide association study involved 1,483 biopsied NAFLD cases and 17,781 controls.•Main analysis shows genome-wide significance for PNPLA3, TM6SF2, HSD17B13 and GCKR.•Sub-analyses show significance near LEPR for NASH and near PYGO1 for steatosis.•Except for GCKR, the genome-wide significant signals were replicated. Background & AimsGenetic factors associated with non-alcoholic fatty liver disease (NAFLD) remain incompletely understood. To date, most genome-wide association studies (GWASs) have adopted radiologically assessed hepatic triglyceride content as the reference phenotype and so cannot address steatohepatitis or fibrosis. We describe a GWAS encompassing the full spectrum of histologically characterised NAFLD.MethodsThe GWAS involved 1,483 European NAFLD cases and 17,781 genetically matched controls. A replication cohort of 559 NAFLD cases and 945 controls was genotyped to confirm signals showing genome-wide or close to genome-wide significance.ResultsCase-control analysis identified signals showing p values ≤5 × 10−8 at 4 locations (chromosome [chr] 2 GCKR/C2ORF16; chr4 HSD17B13; chr19 TM6SF2; chr22 PNPLA3) together with 2 other signals with p <1 × 10−7 (chr1 near LEPR and chr8 near IDO2/TC1). Case-only analysis of quantitative traits showed that the PNPLA3 signal (rs738409) had genome-wide significance for steatosis, fibrosis and NAFLD activity score and a new signal (PYGO1 rs62021874) had close to genome-wide significance for steatosis (p = 8.2 × 10−8). Subgroup case-control analysis for NASH confirmed the PNPLA3 signal. The chr1 LEPR single nucleotide polymorphism also showed genome-wide significance for this phenotype. Considering the subgroup with advanced fibrosis (≥F3), the signals on chr2, chr19 and chr22 maintained their genome-wide significance. Except for GCKR/C2ORF16, the genome-wide significance signals were replicated.ConclusionsThis study confirms PNPLA3 as a risk factor for the full histological spectrum of NAFLD at genome-wide significance levels, with important contributions from TM6SF2 and HSD17B13. PYGO1 is a novel steatosis modifier, suggesting that Wnt signalling pathways may be relevant in NAFLD pathogenesis.Lay summaryNon-alcoholic fatty liver disease is a common disease where excessive fat accumulates in the liver and may result in cirrhosis. To understand who is at risk of developing this disease and suffering liver damage, we undertook a genetic study to compare the genetic profiles of people suffering from fatty liver disease with genetic profiles seen in the general population. We found that particular sequences in 4 different areas of the human genome were seen at different frequencies in the fatty liver disease cases. These sequences may help predict an individual's risk of developing advanced disease. Some genes where these sequences are located may also be good targets for future drug treatments. Genetic factors associated with non-alcoholic fatty liver disease (NAFLD) remain incompletely understood. To date, most genome-wide association studies (GWASs) have adopted radiologically assessed hepatic triglyceride content as the reference phenotype and so cannot address steatohepatitis or fibrosis. We describe a GWAS encompassing the full spectrum of histologically characterised NAFLD. The GWAS involved 1,483 European NAFLD cases and 17,781 genetically matched controls. A replication cohort of 559 NAFLD cases and 945 controls was genotyped to confirm signals showing genome-wide or close to genome-wide significance. Case-control analysis identified signals showing p values ≤5 × 10−8 at 4 locations (chromosome [chr] 2 GCKR/C2ORF16; chr4 HSD17B13; chr19 TM6SF2; chr22 PNPLA3) together with 2 other signals with p <1 × 10−7 (chr1 near LEPR and chr8 near IDO2/TC1). Case-only analysis of quantitative traits showed that the PNPLA3 signal (rs738409) had genome-wide significance for steatosis, fibrosis and NAFLD activity score and a new signal (PYGO1 rs62021874) had close to genome-wide significance for steatosis (p = 8.2 × 10−8). Subgroup case-control analysis for NASH confirmed the PNPLA3 signal. The chr1 LEPR single nucleotide polymorphism also showed genome-wide significance for this phenotype. Considering the subgroup with advanced fibrosis (≥F3), the signals on chr2, chr19 and chr22 maintained their genome-wide significance. Except for GCKR/C2ORF16, the genome-wide significance signals were replicated. This study confirms PNPLA3 as a risk factor for the full histological spectrum of NAFLD at genome-wide significance levels, with important contributions from TM6SF2 and HSD17B13. PYGO1 is a novel steatosis modifier, suggesting that Wnt signalling pathways may be relevant in NAFLD pathogenesis.