Genome-wide interaction study with BMI identifies CYP7A1 and GIPR as genetic modulators of MASLD

Metabolic dysfunction-associated steatotic liver disease (MASLD) may progress to liver inflammation, fibrosis, cirrhosis and hepatocellular carcinoma. So far, genome-wide association studies explain a small fraction of MASLD heritability. We sought to identify novel genetic determinants of MASLD by exploring interactions between genetic variants and body mass index (BMI). First, we examined genome-wide interactions with BMI for circulating alanine aminotransferase (ALT) levels using UK Biobank data. For identified loci, we next examined associations with hepatic proton density fat fraction (PDFF) in 35,146 independent UK Biobank participants. Associations with PDFF were replicated in four independent European cohorts, followed by a phenome-wide association study. Finally, we used human liver epigenomic maps and CRISPR/Cas9 experiments in vitro and in vivo to functionally characterize the CYP7A1 locus. Thirteen loci interact with BMI for ALT (P<5E-8), including eight well-known genetic modulators of MASLD. Two loci - UBXN2B/CYP7A1 and GIPR - are additionally associated with PDFF. For the intronic rs34783010 in GIPR, the minor T allele is associated with lower BMI and higher HbA1c and liver triglyceride content in humans. The UBXN2B/CYP7A1 locus is associated with PDFF in four additional European cohorts. Epigenomic data and in vitro experiments in human liver cells prioritise rs10504255 and CYP7A1 as the functional effectors in this locus. Perturbation of CYP7A1 orthologues using CRISPR/Cas9 results in less liver fat in 10-day-old, metabolically challenged zebrafish larvae. A genome-wide SNPxBMI design fuelled identification of two MASLD genes: CYP7A1 and GIPR.