Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2024 May;56(5):827-837.
doi: 10.1038/s41588-024-01720-y. Epub 2024 Apr 17.

Integrative common and rare variant analyses provide insights into the genetic architecture of liver cirrhosis

Jonas Ghouse  1   2 Gardar Sveinbjörnsson  3 Marijana Vujkovic  4   5   6 Anne-Sofie Seidelin  7 Helene Gellert-Kristensen  7 Gustav Ahlberg  8 Vinicius Tragante  3 Søren A Rand  9   8 Joseph Brancale  10 Silvia Vilarinho  10 Pia Rengtved Lundegaard  8 Erik Sørensen  11 Christian Erikstrup  12 Mie Topholm Bruun  13 Bitten Aagaard Jensen  14 Søren Brunak  15 Karina Banasik  16 Henrik Ullum  17 DBDS Genomic ConsortiumNiek Verweij  18 Luca Lotta  18 Aris Baras  18 Regeneron Genetics CenterTooraj Mirshahi  19 David J Carey  19 Geisinger-Regeneron DiscovEHR CollaborationVA Million Veteran ProgramDavid E Kaplan  4   5 Julie Lynch  20   21 Timothy Morgan  22   23 Tae-Hwi Schwantes-An  22   24 Daniel R Dochtermann  25 Saiju Pyarajan  25   26 Philip S Tsao  27   28 Estonian Biobank Research TeamTriin Laisk  29 Reedik Mägi  29 Julia Kozlitina  30 Anne Tybjærg-Hansen  7 David Jones  31 Kirk U Knowlton  32   33 Lincoln Nadauld  31   34 Egil Ferkingstad  3 Einar S Björnsson  35   36 Magnus O Ulfarsson  3   37 Árni Sturluson  3 Patrick Sulem  3 Ole B Pedersen  38   39 Sisse R Ostrowski  11   38 Daniel F Gudbjartsson  3   40 Kari Stefansson  3 Morten Salling Olesen  9   8 Kyong-Mi Chang  4   5 Hilma Holm  3 Henning Bundgaard  9   38 Stefan Stender  41   42
Collaborators, Affiliations

Integrative common and rare variant analyses provide insights into the genetic architecture of liver cirrhosis

Jonas Ghouse et al. Nat Genet. 2024 May.

Abstract

We report a multi-ancestry genome-wide association study on liver cirrhosis and its associated endophenotypes, alanine aminotransferase (ALT) and γ-glutamyl transferase. Using data from 12 cohorts, including 18,265 cases with cirrhosis, 1,782,047 controls, up to 1 million individuals with liver function tests and a validation cohort of 21,689 cases and 617,729 controls, we identify and validate 14 risk associations for cirrhosis. Many variants are located near genes involved in hepatic lipid metabolism. One of these, PNPLA3 p.Ile148Met, interacts with alcohol intake, obesity and diabetes on the risk of cirrhosis and hepatocellular carcinoma (HCC). We develop a polygenic risk score that associates with the progression from cirrhosis to HCC. By focusing on prioritized genes from common variant analyses, we find that rare coding variants in GPAM associate with lower ALT, supporting GPAM as a potential target for therapeutic inhibition. In conclusion, this study provides insights into the genetic underpinnings of cirrhosis.

PubMed Disclaimer

Conflict of interest statement

The authors who are affiliated with deCODE genetics/Amgen declare competing financial interests as employees. H.B. receives lecture fees from Bristol-Myers Squibb, Merck Sharp and Dohme. J.G. has received lecture fee from Illumina. S.B. is a board member for Proscion A/S and Intomics A/S. N.V., L.L. and A.B. are employees at Regeneron Genetics Center. All other authors have no conflict of interest to declare.

Figures

Fig. 1
Fig. 1. Study design and main results.
a, Overall study design, with stage 1 representing the European-specific GWAS meta-analysis and stage 2 representing the cross-ancestry GWAS meta-analysis. b, Miami plot of cirrhosis GWAS. The x axis is the chromosomal location of SNPs and the y axis is the strength of association −log10(P). Note that the y axis includes a break at 50. The lead SNPs and SNPs located within ±1 Mb are highlighted, and the nearest genes are annotated. The top plot shows results from the European-specific GWAS meta-analysis, whereas the bottom plot displays results from the cross-ancestry GWAS meta-analysis. The dashed red line represents the threshold for genome-wide significance (P < 5 × 10−8). P values were two-sided and based on an IVW fixed-effects meta-analysis, and not adjusted for multiple testing.
Fig. 2
Fig. 2. Cross-trait associations between cirrhosis variants and selected metabolic and hepatobiliary traits relevant for cirrhosis.
The heatmap shows associations of variants identified through cirrhosis GWAS meta-analysis or endophenotype-driven analyses with 41 binary and quantitative traits sampled from meta-analysis of data from CHB-CID/DBDS, deCODE, Intermountain Healthcare, FinnGen, UKB and external sources, where available. The number of cases for binary traits and sample size for quantitative traits are shown in parenthesis following each trait. Shown are variants and phenotypes with significant associations after correcting for multiple testing using an FDR of <0.05. P values (two-sided) were derived from linear and logistic regression models. Hierarchical clustering was performed on a variant level using the complete linkage method based on Euclidian distance. Coloring represents z scores for each respective trait or disease, oriented toward the cirrhosis risk-increasing allele. Red indicates an increase in the trait or disease risk, while blue indicates a decrease in the trait or disease risk. SHBG, sex hormone-binding globulin; IGF-1, insulin growth factor 1; ApoA, apolipoprotein A; ApoB, apolipoprotein B; COPD, chronic obstructive pulmonary disease; WHRadjBMI, waist-to-hip-ratio adjusted for BMI; LDL-C, low-density lipoprotein cholesterol; T1D, type 1 diabetes.
Fig. 3
Fig. 3. Comparison of genetic associations with NAFLD, ALD and cirrhosis.
a, Shown are the effects of 18 previously reported NAFLD variants and 36 cirrhosis variants identified in this study, totaling 38 distinct signals. The NAFLD effects were derived from meta-analysis of data from deCODE, UKB, CHB, Intermountain Healthcare and FinnGen (n = 22,944). Cirrhosis effects were derived from the cross-ancestry meta-analysis. Variants that had stronger effects (PHet < 0.05/38, that is corrected for 38 tests) on NAFLD compared with cirrhosis are colored blue, while variants with stronger effects on cirrhosis compared with NAFLD are colored red. b, Shown are the effects of the 36 cirrhosis variants on NAFLD (n = 22,944) and ALD (n = 2,931). Variants with stronger effects (PHet < 0.05/36, that is corrected for 36 tests) on ALD compared with NAFLD are colored red. For both a and b, points refer to effect estimates (log(OR), measure of center), error bars represent 95% CI and the solid line represents the line of best fit. The dashed identity line (y = x) is shown for reference. PHet were two-sided and obtained using a likelihood ratio test (Cochran’s Q).
Fig. 4
Fig. 4. Interaction between PNPLA3 p.Ile148Met (rs738409), environment and risk of liver-related outcomes in the UKB.
Shown are the association between rs738409 carrier status and risk of five liver outcomes according to BMI categories (<30 versus ≥30 kg m2), weekly alcohol intake (≤14 versus >14 weekly drinks) and T2D (no versus yes). Points refer to OR (measure of center), error bars represent 95% CIs and P represents P value for interaction. We used logistic regression models, adjusted for age, sex and ten principal components (PCs). Interactions between the variant and environmental factors were evaluated using likelihood ratio tests, comparing a main-effect model (variant + environmental factor) with a model including an interaction term (variant × environmental factor). The number of exposed individuals and the number of outcomes within each subcategory are listed in Supplementary Table 15.
Fig. 5
Fig. 5. PRSs and hepatobiliary outcomes.
Associations between six different PRS (PRSALT, PRS5-SNP, PRS15-SNP, PRS36-SNP, PRSEUR and PRSCA) and hepatobiliary outcomes in the UKB. Points refer to ORs (measure of center) per s.d. increase in PRS; error bars represent 95% CIs. Total number of cases is provided following each outcome. Logistic regression models were used, adjusted for age, sex and ten PCs.
Fig. 6
Fig. 6. PRSs and disease progression.
a,b, Risk of HCC in individuals with liver cirrhosis according to polygenic risk percentile. c,d, Risk of cirrhosis in individuals with NAFLD diagnosis, according to polygenic risk percentile. We used the PRS15-SNP to assign individuals to different risk groups. Cumulative incidence (solid line, measure of center) was estimated using Fine-Gray regression, which takes the competing risk of death into account. The lighter shades represent the respective 95% CIs. Number of individuals at risk according to each exposure group and events are given below each plot.

Similar articles

References

    1. Ginès P, et al. Liver cirrhosis. Lancet Lond. Engl. 2021;398:1359–1376. doi: 10.1016/S0140-6736(21)01374-X. - DOI - PubMed
    1. Huang DQ, et al. Global epidemiology of cirrhosis—aetiology, trends and predictions. Nat. Rev. Gastroenterol. Hepatol. 2023;20:388–398. doi: 10.1038/s41575-023-00759-2. - DOI - PMC - PubMed
    1. Loomba R, et al. Heritability of hepatic fibrosis and steatosis based on a prospective twin study. Gastroenterology. 2015;149:1784–1793. doi: 10.1053/j.gastro.2015.08.011. - DOI - PMC - PubMed
    1. Emdin CA, et al. Association of genetic variation with cirrhosis: a multi-trait genome-wide association and gene–environment interaction study. Gastroenterology. 2021;160:1620–1633. doi: 10.1053/j.gastro.2020.12.011. - DOI - PMC - PubMed
    1. Gellert-Kristensen H, et al. Combined effect of PNPLA3, TM6SF2, and HSD17B13 variants on risk of cirrhosis and hepatocellular carcinoma in the general population. Hepatology. 2020;72:845–856. doi: 10.1002/hep.31238. - DOI - PubMed