We tested genetic variants in GCK, GCKR, and PNPLA3 in a large sample of self-identified Mexican Americans from the BetaGene Study for association with type 2 diabetes-related phenotypes under the hypothesis that they may regulate metabolic balance across the liver and contribute to hepatic steatosis and insulin resistance. We further tested whether interactions with dietary fructose and total sugar contributes to the observed associations. GCK rs1799831 was not associated with any type 2 diabetes-related phenotypes either alone or with any interaction tested. We replicated previous associations reported for GCKR rs780094 and PNPLA3 rs738409. We also show the interaction between GCKR rs780094 and dietary fructose is associated with both glucose effectiveness and glucose effectiveness at zero insulin, measures reflective of hepatic glucose uptake. We further show the interaction between GCKR rs780094 and PNPLA3 rs738409 is associated with type 2 diabetes-related traits, including insulin sensitivity. We conclude variations in GCKR and PNPLA3 and their interactions with each other and dietary fructose are partial determinants of hepatic fat, likely due to alterations in relative contributions of different metabolic pathways in the liver. These findings point to both GCKR and PNPLA3 as important therapeutic targets to mitigate hepatic metabolic dysfunction.

We examined genetic variation in GCKR and PNPLA3 for association with type 2 diabetes-related traits in Mexican Americans without diabetes. We addressed three questions: Are these loci associated with type 2 diabetes-related traits? Is there interaction between these loci? Is there interaction with dietary sugars? These loci individually, and in an interaction, are associated with type 2 diabetes-related traits; these loci also interact with dietary fructose to alter hepatic glucose uptake. Our findings support the hypothesis that these loci, coupled with dietary fructose and total sugars, alter hepatic metabolic balance and may contribute to hepatic fat and insulin resistance.

The PPP1R3B gene, encoding PPP1R3B protein, is critical for liver glycogen synthesis and maintaining blood glucose levels. Genetic variants affecting PPP1R3B expression are associated with several metabolic traits and liver disease, but the precise mechanisms are not fully understood. We studied the effects of both Ppp1r3b overexpression and deletion in mice and cell models and found that both changes in Ppp1r3b expression result in dysregulated metabolism and liver damage, with overexpression increasing liver glycogen stores, while deletion resulted in higher liver lipid accumulation. These patterns were confirmed in humans where variants increasing PPP1R3B expression had lower liver fat and decreased plasma lipids, whereas putative loss-of-function variants were associated with increased liver fat and elevated plasma lipids. These findings support that PPP1R3B is a crucial regulator of hepatic metabolism beyond glycogen synthesis and that genetic variants affecting PPP1R3B expression levels influence if hepatic energy is stored as glycogen or triglycerides.

Non-alcoholic fatty liver disease (NAFLD) represents the most widespread liver disease globally, ranging from non-alcoholic fatty liver (NAFL) and steatohepatitis (NASH) to fibrosis/cirrhosis, with potential progression to hepatocellular carcinoma (HCC). Genome-wide association studies (GWASs) have identified several single nucleotide polymorphisms (SNPs) associated with NAFLD. However, numerous GWAS signals associated with NAFLD locate in non-coding regions, posing a challenge for interpreting their functional annotation.

In this study, we utilized the Activity-by-Contact (ABC) model to construct the enhancer-gene maps of liver by integrating epigenomic data from 15 liver tissues and cell lines. We constructed the most comprehensive genome-wide regulatory maps of the liver, identifying 543,486 enhancer-gene connections, including 267,857 enhancers and 16,872 target genes. Enrichment analyses revealed that the ABC SNPs are significantly enriched in active chromatin regions and active chromatin state. By combining the ABC regulatory maps and NAFLD GWAS data, we systematically identified ABC SNPs associated with NAFLD risk. Through the functional annotations, such as pathway enrichment and drug-gene interaction analyses, we identified 6 genes (GGT1, ACTG1, SPP1, EPHA2, PROZ and SHMT1) as candidate NAFLD genes, with SHMT1 previously reported. Among the SNPs connected to the candidate genes, the ABC SNP rs2017869 (odds ratio [OR] for the C allele = 1.10, 95% CI = 1.04-1.16, P = 5.97 × 10- 4) had the highest ABC score. According to the ABC maps, rs2017869 links to GGT1, and several drugs targeting this gene, such as liothyronine, showed potential benefits to patients with NAFLD. Furthermore, we identified that another novel gene, EPHA2, may play a crucial role in NAFLD by regulating the GGT levels.

Our study provides the most comprehensive ABC regulatory maps of the liver to date. This resource offers a valuable reference for identifying regulatory variants and prioritizing susceptibility genes of liver diseases, such as NAFLD.

Carbohydrate response element-binding protein (ChREBP) is a transcription factor activated by glucose metabolites that orchestrates the expression of genes involved in glycolysis, de novo lipogenesis, and ATP homeostasis. Inadequate ChREBP activity impairs the cellular adaptations to glucose exposure and in humans associates with dyslipidemia, fatty liver disease, and type 2 diabetes. ChREBP activity is regulated by cytosolic-nuclear translocation involving its low-glucose inhibitory domain (LID). Whether this domain is targeted by post-translational lysine acetylation is unknown. Here we report a novel LID acetylation site that controls activity and protein interactions of ChREBP. Mutation of this residue increased glucose-induced activity and target gene expression of ChREBP. Mechanistically, mutant ChREBP protein showed more nuclear localization and enhanced genomic binding to a target promoter. Interactions with proteins that exhibit differential binding upon glucose exposure were attenuated by the mutation, demonstrating the importance of the LID in the formation of the protein interactome. Particularly interactions with 14-3-3 proteins, factors that regulate cytosolic/nuclear trafficking of ChREBP, were reduced, whereas interactions with proteins of the nucleosome remodeling deacetylase complex (NuRD) were increased. These molecular insights may shape new therapeutic strategies to target ChREBP activity and counteract metabolic diseases.

To identify nongenetic factors influences on DNA methylation (DNAm) variations associated with blood Alanine Aminotransferase (ALT) concentration, this study conducted an epigenome-wide association study (EWAS) on Chinese monozygotic twins.

A total of 61 pairs of Chinese monozygotic twins involved in this study. Whole blood samples were analyzed for DNAm profiling using the Reduced Representation Bisulfite Sequencing (RRBS) technique. We examined the relationship between DNAm levels at each CpG site and serum ALT using a linear mixed-effects model. Enrichment analysis and causal inference analysis was conducted, and differentially methylated regions (DMRs) were further identified. Candidate CpGs were validated in a community sample. Genome-wide significance were calculated by Bonferroni correction (p < 2.14 × 10-7).

We identified 85 CpGs reaching genome-wide significance (p < 2.14 × 10-7), located in 16 genes including FLT4, ADARB2, MRPS31P2, and RELB. Causal inference suggested that DNAm at 61 out of 85 significant CpGs within 14 genes influenced ALT level. 52 DMRs and 1765 pathways such as low voltage-gated calcium channel activity and focal adhesion were identified having influences on ALT levels. Further validation using community population found four CpGs mapped to FLT4 and three to RELB showing hypomethylation and hypermethylation in cases with abnormal ALT (ALT > 40 U/L), respectively.

This study identified several differentially methylated CpG sites associated with serum ALT in the Chinese population, particularly within FLT4 and RELB. These findings provide new insights into the epigenetic modifications underlying liver function.

Metabolic dysfunction-associated steatotic liver disease (MASLD) is characterized by increased hepatic steatosis with cardiometabolic disease and is a leading cause of advanced liver disease. We review here the genetic basis of MASLD. The genetic variants most consistently associated with hepatic steatosis implicate genes involved in lipoprotein input or output, glucose metabolism, adiposity/fat distribution, insulin resistance, or mitochondrial/ER biology. The distinct mechanisms by which these variants promote hepatic steatosis result in distinct effects on cardiometabolic disease that may be best suited to precision medicine. Recent work on gene-environment interactions has shown that genetic risk is not fixed and may be exacerbated or attenuated by modifiable (diet, exercise, alcohol intake) and nonmodifiable environmental risk factors. Some steatosis-associated variants, notably those in patatin-like phospholipase domain-containing 3 (PNPLA3) and transmembrane 6 superfamily member 2 (TM6SF2), are associated with risk of developing adverse liver-related outcomes and provide information beyond clinical risk stratification tools, especially in individuals at intermediate to high risk for disease. Future work to better characterize disease heterogeneity by combining genetics with clinical risk factors to holistically predict risk and develop therapies based on genetic risk is required.

Hypercholesterolemia, characterized by elevated levels of low-density lipoprotein cholesterol (LDL-C), along with inflammation, is a well-known risk factor for developing atherosclerosis and coronary artery disease (CAD). Many patients with hypercholesterolemia may carry inherited genetic variants that are not part of the commonly recognized mutations in the LDLR, APOB, LDLRAP1, and PCSK9 genes. These genetic variants may have cumulative effects that contribute to increased LDL-C levels and CAD development. The polygenic risk score (PRS) may provide an essential tool for evaluating an individual's genetic predisposition to these conditions. This pilot study aimed to investigate the impact of the PRS calculated from specific single nucleotide polymorphisms (SNPs) associated with LDL cholesterol (LDL-C)-namely, CELSR2 rs629301, APOB rs1367117, ABCG8 rs6544713, LDLR rs6511720, APOE rs429358, and rs7412-on LDL-C levels in both healthy individuals with elevated LDL-C levels (>2.6 mmol/L) and those diagnosed with ST-segment elevation myocardial infarction (STEMI). A total of 61 healthy individuals with high LDL-C levels (>2.6 mmol/L) and 93 STEMI patients were selected for the study. The High-Resolution Melting Polymerase Chain Reaction (HRM PCR) method was adopted and sequencing techniques were employed to identify the specific single nucleotide polymorphisms (SNPs) of interest. The patient group exhibited a PRS of 0.824 (with a range of -0.62 to 1.174) compared to 0.674 (range: -0.176 to 0.974) in healthy individuals, indicating a higher genetic predisposition to elevated LDL-C levels (p = 0.001) in patients. Interestingly, patients had lower LDL-C concentrations than healthy individuals. Additionally, a more significant number of patients were past smokers and statin users. The PRS calculations revealed that patients with a higher PRS had increased odds of experiencing an MI, with an odds ratio of 12.044 (95% confidence interval: 1.551-93.517, p = 0.017). Similarly, smokers showed even higher odds, with an odds ratio of 24.962 (95% CI: 7.171-86.890, p < 0.001). Among healthy individuals, those with a higher PRS had increased odds of having an LDL-C concentration greater than 4.9 mmol/L (odds ratio: 20.391, 95% CI: 1.116-358.486, p = 0.039). However, no significant association was found between the PRS and LDL-C levels in the patient group during hospitalization (p = 0.782). This pilot study shows that PRS can be employed to evaluate the risk of MI and to estimate concentrations greater than 4.9 mmol/L LDL-C in healthy individuals.

TM6SF2 -rs58542926-T is associated with increased cirrhosis and modestly decreased coronary artery disease prevalence. However, relative effects of TM6SF2 genotype on major adverse cardiovascular events (MACE) vs liver-related events (LRE) are not known.

We used the UK Biobank, a prospective cohort with genetic and inpatient diagnosis data. The primary predictor was TM6SF2 -rs58542926 genotype, and the primary outcomes were MACE and LRE. Effects were reported as subhazard ratios (sHRs) and 10-year cumulative incidence by Fine-Gray competing risk analyses.

More than 430,000 individuals met inclusion criteria. TM6SF2 -rs58542926-TT genotype (vs CC) was associated with higher incidence of LRE (adjusted sHR 3.16, 95% confidence interval 1.86-5.37) and lower incidence of MACE (adjusted sHR for TT vs CC genotype 0.76, 95% confidence interval 0.63-0.91). In individuals with fibrosis-4 (FIB4) < 1.3, 1.3-2.67, and > 2.67, 10-year LRE incidence in TM6SF2 -rs58542926-TT vs CC individuals was 0.08% vs 0.06% ( P > 0.05), 0.81% vs 0.20% ( P < 0.0001), and 10.5% vs 3.4% ( P = 0.00094), respectively. The corresponding values for MACE were 3.8% vs 5.1% ( P = 0.032), 6.4% vs 8.2% ( P = 0.040), and 17.1% vs 12.4% ( P > 0.05). The absolute decrease in MACE with rs58542926-TT (vs CC) genotype exceeded the absolute increase in LRE in all groups but FIB4 > 2.67. Associations of TM6SF2 genotype with LRE/MACE were significant in men but not women. TM6SF2 -rs58542926-T allele was also associated with increased hepatic steatosis and corrected T1 time by magnetic resonance imaging, with greater effect sizes in men than women.

TM6SF2 genotype has opposite effects on LRE vs MACE incidence, and absolute effects on MACE were greater except in those with highest FIB4 scores. Effects were strongest in men. These findings clarify implications of TM6SF2 genotype based on personalized clinical risk.

The precise pathomechanisms underlying the development of non-alcoholic steatohepatitis (NASH, also known as metabolic dysfunction-associated steatohepatitis [MASH]) remain incompletely understood. In this study, we investigated the potential role of EF-hand domain family member D2 (EFHD2), a novel molecule specific to immune cells, in the pathogenesis of NASH.

Hepatic EFHD2 expression was characterized in patients with NASH and two diet-induced NASH mouse models. Single-cell RNA sequencing (scRNA-seq) and double-immunohistochemistry were employed to explore EFHD2 expression patterns in NASH livers. The effects of global and myeloid-specific EFHD2 deletion on NASH and NASH-related hepatocellular carcinoma were assessed. Molecular mechanisms underlying EFHD2 function were investigated, while chemical and genetic investigations were performed to assess its potential as a therapeutic target.

EFHD2 expression was significantly elevated in hepatic macrophages/monocytes in both patients with NASH and mice. Deletion of EFHD2, either globally or specifically in myeloid cells, improved hepatic steatosis, reduced immune cell infiltration, inhibited lipid peroxidation-induced ferroptosis, and attenuated fibrosis in NASH. Additionally, it hindered the development of NASH-related hepatocellular carcinoma. Specifically, deletion of myeloid EFHD2 prevented the replacement of TIM4+ resident Kupffer cells by infiltrated monocytes and reversed the decreases in patrolling monocytes and CD4+/CD8+ T cell ratio in NASH. Mechanistically, our investigation revealed that EFHD2 in myeloid cells interacts with cytosolic YWHAZ (14-3-3ζ), facilitating the translocation of IFNγR2 (interferon-γ receptor-2) onto the plasma membrane. This interaction mediates interferon-γ signaling, which triggers immune and inflammatory responses in macrophages during NASH. Finally, a novel stapled α-helical peptide targeting EFHD2 was shown to be effective in protecting against NASH pathology in mice.

Our study reveals a pivotal immunomodulatory and inflammatory role of EFHD2 in NASH, underscoring EFHD2 as a promising druggable target for NASH treatment.

Non-alcoholic steatohepatitis (NASH) represents an advanced stage of non-alcoholic fatty liver disease (NAFLD); however, not all patients with NAFLD progress to NASH. A key challenge is identifying the factors that trigger inflammation, which propels the transition from simple fatty liver to NASH. Our research pinpointed EFHD2 as a pivotal driver of NASH, orchestrating the over-activation of interferon-γ signaling within the liver during NASH progression. A stapled peptide designed to target EFHD2 exhibited therapeutic promise in NASH mice. These findings support the potential of EFHD2 as a therapeutic target in NASH.

Metabolic dysfunction-associated steatotic liver disease (MASLD) comprises a spectrum of liver diseases that span simple steatosis, metabolic dysfunction-associated steatohepatitis (MASH) and fibrosis and may progress to cirrhosis and cancer. The pathogenesis of MASLD is multifactorial and is driven by environmental, genetic, metabolic and immune factors. This review will focus on the role of the type 3 cytokines IL-17 and IL-22 in MASLD pathogenesis and progression. IL-17 and IL-22 are produced by similar adaptive and innate immune cells such as Th17 and innate lymphoid cells, respectively. IL-17-related signaling is upregulated during MASLD resulting in increased chemokines and proinflammatory cytokines in the liver microenvironment, enhanced recruitment of myeloid cells and T cells leading to exacerbation of inflammation and liver disease progression. IL-17 may also act directly by activating hepatic stellate cells resulting in increased fibrosis. In contrast, IL-22 is a pleiotropic cytokine with a dominantly protective signature in MASLD and is currently being tested as a therapeutic strategy. IL-22 also exhibits beneficial metabolic effects and abrogates MASH-related inflammation and fibrosis development via inducing the production of anti-oxidants and anti-apoptotic factors. A sex-dependent effect has been attributed to both cytokines, most importantly to IL-22 in MASLD or related conditions. Altogether, IL-17 and IL-22 are key effectors in MASLD pathogenesis and progression. We will review the role of these two cytokines and cells that produce them in the development of MASLD, their interaction with host factors driving MASLD including sexual dimorphism, and their potential therapeutic benefits.

Gamma-glutamyltransferase (GGT) is known as an oxidative stress marker, induced by alcohol consumption and metabolic disorders, and is reported as a predictor of hepatocellular carcinoma (HCC) development after hepatitis C virus (HCV) elimination. However, it is not clear whether GGT serves simply as a surrogate marker for overlapping metabolic diseases or reflects HCV-specific carcinogenicity. We investigated the association between GGT and hepatocarcinogenesis after achieving a sustained viral response (SVR), accounting for drinking habits or diabetes, and examined predisposing factors associated with GGT levels after SVR.

This is a prospective, multicenter, and observational study using the database of 1001 patients after HCV eradication with direct-acting antiviral agents. The association of GGT at SVR with cumulative HCC development was examined in a multivariate analysis using Cox proportional hazard models after adjustment for covariates including alcohol and diabetes. The association between oxidative stress markers or genetic factors and GGT levels was analyzed.

High GGT levels at SVR were associated with HCC development (HR] 2.38, 95% CI 1.10-5.17). This association was also significant when restricted to patients without alcohol consumption or diabetes (HR 8.38, 95% CI 2.87-24.47). GGT levels were correlated with serum growth differentiation factor 15 levels, a marker of mitochondrial dysfunction. Single-nucleotide polymorphisms of ZNF827 and GDF15 were associated with high GGT levels.

High GGT levels at SVR were associated with HCC development after accounting for alcohol consumption and diabetes. GGT levels are influenced by genetic predisposition and may reflect mitochondrial dysfunction after HCV eradication.

Although the inherited risk factors associated with fatty liver disease are well understood, little is known about the genetic background of metabolic dysfunction-associated steatotic liver disease (MASLD) and its related health impacts. Compared to non-alcoholic fatty liver disease (NAFLD), MASLD presents significantly distinct diagnostic criteria, and epidemiological and clinical features, but the related genetic variants are yet to be investigated. Therefore, we conducted this study to assess the genetic background of MASLD and interactions between MASLD-related genetic variants and metabolism-related outcomes.

Participants from the UK Biobank were grouped into discovery and replication cohorts for an MASLD genome-wide association study (GWAS), and base and target cohorts for polygenic risk score (PRS) analysis. Autosomal genetic variants associated with NAFLD were compared with the MASLD GWAS results. Kaplan-Meier and Cox regression analyses were used to assess associations between MASLD and metabolism-related outcomes.

Sixteen single-nucleotide polymorphisms (SNPs) were identified at genome-wide significance levels for MASLD and duplicated in the replication cohort. Differences were found after comparing these SNPs with the results of NAFLD-related genetic variants. MASLD cases with high PRS had a multivariate-adjusted hazard ratio of 3.15 (95% confidence interval, 2.54-3.90) for severe liver disease (SLD), and 2.81 (2.60-3.03) for type 2 diabetes mellitus. The high PRS amplified the impact of MASLD on SLD and extrahepatic outcomes.

High PRS of MASLD GWAS amplified the impact of MASLD on SLD and metabolism-related outcomes, thereby refining the process of identification of individuals at high risk of MASLD. Supplementation of this process with relevant genetic backgrounds may lead to more effective MASLD prevention and management.

Genetic susceptibility to metabolic associated fatty liver disease (MAFLD) is complex and poorly characterized. Accurate characterization of the genetic background of hepatic fat content would provide insights into disease etiology and causality of risk factors. We performed genome-wide association study (GWAS) on two noninvasive definitions of hepatic fat content: magnetic resonance imaging proton density fat fraction (MRI-PDFF) in 16,050 participants and fatty liver index (FLI) in 388,701 participants from the United Kingdom (UK) Biobank (UKBB). Heritability, genetic overlap, and similarity between hepatic fat content phenotypes were analyzed, and replicated in 10,398 participants from the University Medical Center Groningen (UMCG) Genetics Lifelines Initiative (UGLI). Meta-analysis of GWASs of MRI-PDFF in UKBB revealed five statistically significant loci, including two novel genomic loci harboring CREB3L1 (rs72910057-T, P = 5.40E-09) and GCM1 (rs1491489378-T, P = 3.16E-09), respectively, as well as three previously reported loci: PNPLA3, TM6SF2, and APOE. GWAS of FLI in UKBB identified 196 genome-wide significant loci, of which 49 were replicated in UGLI, with top signals in ZPR1 (P = 3.35E-13) and FTO (P = 2.11E-09). Statistically significant genetic correlation (rg) between MRI-PDFF (UKBB) and FLI (UGLI) GWAS results was found (rg = 0.5276, P = 1.45E-03). Novel MRI-PDFF genetic signals (CREB3L1 and GCM1) were replicated in the FLI GWAS. We identified two novel genes for MRI-PDFF and 49 replicable loci for FLI. Despite a difference in hepatic fat content assessment between MRI-PDFF and FLI, a substantial similar genetic architecture was found. FLI is identified as an easy and reliable approach to study hepatic fat content at the population level.

A splice variant in HSD17B13 gene is demonstrated to protect against nonalcoholic fatty liver disease (NAFLD), and mitigate the effect of Patatin-like phospholipase domain-containing 3 (PNPLA3-I148M). It is being explored as a putative drug target and in polygenic risk scores. Based on whole exome sequencing (WES) in our cohort of biopsy proven NAFLD and limited data on the variant in our ethnicity, we sought to explore its role.

This is a cross-sectional study that recruited 1,020 individuals with ultrasound/biopsy-confirmed NAFLD and matched controls. Liver enzymes and lipid profiles were estimated (Beckman coulter LX750/DXH800); WES was performed in NAFLD patients and controls (Illumina; HiSeqX); HSD17B13-A-INS/I148M-PNPLA3 variants were genotyped (sequencing/qR T-PCR); HSD17B13 protein expression was estimated (immunohistochemistry); the Student's t-test/Mann-Whitney U/Chi-square test and odds ratio (95% confidence interval) were used.

There was no significant difference (Odds ratio = 0.76; 95% CI -0.57 to 1.03; P = 0.76) in the frequency of the rs72613567-A-INS between controls and patients (17.8% vs. 14.4%). No difference in the ALT (Alanine transaminase; 72.24 ± 65.13 vs. 73.70 ± 60.06; P = 0.51) and AST levels (Aspartate aminotransferase; 60.72 ± 55.59 vs. 61.63 ± 60.33; P = 0.91) between HSD17B13-wild and variant carriers were noted. Significantly elevated liver enzymes were seen in PNPLA3-148-variant/HSD17B13-wild compared with PNPLA3-148-variant/HSD17B13-variant (90.44 ± 59.0 vs. 112.32 ± 61.78; P = 0.02). No difference in steatosis (P = 0.51) between HSD17B13-wild and variant carriers was noted. No other variants in the intron-exon boundaries were identified. Although, protein expression differences were noted between wild and variant carriers, no difference in the extent of steatosis was seen.

Our study reports lack of association of the splice variant with reduced risk of NAFLD, and mitigating the effect of PNPLA3 variant. Ethnicity-based validation must be carried out before including it in assessing protection against NAFLD.

Mitochondrial (MT) dysfunction is a hallmark of liver diseases. However, the effects of functional variants such as protein truncating variants (PTVs) in MT-related genes on the risk of liver diseases have not been extensively explored.

We extracted 60,928 PTVs across 2466 MT-related nucleus genes using whole-exome sequencing data obtained from 442,603 participants in the UK Biobank. We examined their associations with liver dysfunction that represented by the liver-related biomarkers and the risks of chronic liver diseases and liver-related mortality.

96.10% of the total participants carried at least one PTV. We identified 866 PTVs that were positively associated with liver dysfunction at the threshold of P value < 8.21e - 07. The coding genes of these PTVs were mainly enriched in pathways related to lipid, fatty acid, amino acid, and carbohydrate metabolisms. The 866 PTVs were presented in 1.07% (4721) of participants. Compared with participants who did not carry any of the PTVs, the carriers had a 5.33-fold (95% CI 4.15-6.85), 2.82-fold (1.69-4.72), and 4.41-fold (3.04-6.41) increased risk for fibrosis and cirrhosis of liver, liver cancer, and liver disease-related mortality, respectively. These adverse effects were consistent across subgroups based on age, sex, body mass index, smoking status, and presence of hypertension, diabetes, dyslipidemia, and metabolic syndrome.

Our findings revealed a significant impact of PTVs in MT-related genes on liver disease risk, highlighting the importance of these variants in identifying populations at risk of liver diseases and facilitating early clinical interventions.

Pseudokinases are catalytically inactive proteins in the human genome that lack the ability to transfer phosphate from ATP to their substrates. The Tribbles family of pseudokinases contains three members: Tribbles 1, 2, and 3. Tribbles 1 has recently gained importance because of its involvement in various diseases, including cancer. It acts as a scaffolding protein that brings about the degradation of its substrate proteins, such as C/EBPα/β, MLXIPL, and RAR/RXRα, among others, via the ubiquitin proteasome system. It also serves as an adapter protein, which sequesters different protein molecules and activates their downstream signaling, leading to processes, such as cell survival, cell proliferation, and lipid metabolism. It has been implicated in cancers such as AML, prostate cancer, breast cancer, CRC, HCC, and glioma, where it activates oncogenic signaling pathways such as PI3K-AKT and MAPK and inhibits the anti-tumor function of p53. TRIB1 also causes treatment resistance in cancers such as NSCLC, breast cancer, glioma, and promyelocytic leukemia. All these effects make TRIB1 a potential drug target. However, the lack of a catalytic domain renders TRIB1 "undruggable", but knowledge about its structure, conformational changes during substrate binding, and substrate binding sites provides an opportunity to design small-molecule inhibitors against specific TRIB1 interactions.

Fatty acid desaturase 1 (FADS1) is a rate-limiting enzyme in long-chain polyunsaturated fatty acid (LCPUFA) synthesis. Reduced activity of FADS1 was observed in metabolic dysfunction-associated steatotic liver disease (MASLD). The aim of this study was to determine whether adeno-associated virus serotype 8 (AAV8) mediated hepatocyte-specific overexpression of Fads1 (AAV8-Fads1) attenuates western diet-induced metabolic phenotypes in a rat model. Male weanling Sprague-Dawley rats were fed with a chow diet, or low-fat high-fructose (LFHFr) or high-fat high-fructose diet (HFHFr) ad libitum for 8 weeks. Metabolic phenotypes were evaluated at the endpoint. AAV8-Fads1 injection restored hepatic FADS1 protein levels in both LFHFr and HFHFr-fed rats. While AAV8-Fads1 injection led to improved glucose tolerance and insulin signaling in LFHFr-fed rats, it significantly reduced plasma triglyceride (by ~50%) and hepatic cholesterol levels (by ~25%) in HFHFr-fed rats. Hepatic lipidomics analysis showed that FADS1 activity was rescued by AAV8-FADS1 in HFHFr-fed rats, as shown by the restored arachidonic acid (AA)/dihomo-γ-linolenic acid (DGLA) ratio, and that was associated with reduced monounsaturated fatty acid (MUFA). Our data suggest that the beneficial role of AAV8-Fads1 is likely mediated by the inhibition of fatty acid re-esterification. FADS1 is a promising therapeutic target for MASLD in a diet-dependent manner.

Genome-wide association studies (GWASs) have identified hundreds of risk loci for liver disease and lipid-related metabolic traits, although identifying their target genes and molecular mechanisms remains challenging. We predicted target genes at GWAS signals by integrating them with molecular quantitative trait loci for liver gene expression (eQTL) and liver chromatin accessibility QTL (caQTL). We predicted specific regulatory caQTL variants at four GWAS signals located near EFHD1, LITAF, ZNF329, and GPR180. Using transcriptional reporter assays, we determined that caQTL variants rs13395911, rs11644920, rs34003091, and rs9556404 exhibit allelic differences in regulatory activity. We also performed a protein binding assay for rs13395911 and found that FOXA2 differentially interacts with the alleles of rs13395911. For variants rs13395911 and rs11644920 in putative enhancer regulatory elements, we used CRISPRi to demonstrate that repression of the enhancers altered the expression of the predicted target and/or nearby genes. Repression of the element at rs13395911 reduced the expression of EFHD1, and repression of the element at rs11644920 reduced the expression of LITAF, SNN, and TXNDC11. Finally, we showed that EFHD1 is a metabolically active gene in HepG2 cells. Together, these results provide key steps to connect genetic variants with cellular mechanisms and help elucidate the causes of liver disease.

Genome-wide association analyses using high-throughput metabolomics platforms have led to novel insights into the biology of human metabolism1-7. This detailed knowledge of the genetic determinants of systemic metabolism has been pivotal for uncovering how genetic pathways influence biological mechanisms and complex diseases8-11. Here we present a genome-wide association study for 233 circulating metabolic traits quantified by nuclear magnetic resonance spectroscopy in up to 136,016 participants from 33 cohorts. We identify more than 400 independent loci and assign probable causal genes at two-thirds of these using manual curation of plausible biological candidates. We highlight the importance of sample and participant characteristics that can have significant effects on genetic associations. We use detailed metabolic profiling of lipoprotein- and lipid-associated variants to better characterize how known lipid loci and novel loci affect lipoprotein metabolism at a granular level. We demonstrate the translational utility of comprehensively phenotyped molecular data, characterizing the metabolic associations of intrahepatic cholestasis of pregnancy. Finally, we observe substantial genetic pleiotropy for multiple metabolic pathways and illustrate the importance of careful instrument selection in Mendelian randomization analysis, revealing a putative causal relationship between acetone and hypertension. Our publicly available results provide a foundational resource for the community to examine the role of metabolism across diverse diseases.

General obesity is a well-established risk factor for gallstone disease (GSD), but whether central obesity contributes additional independent risk remains controversial. We aimed to comprehensively clarify the effect of body fat distribution on GSD.

We first investigated the observational association of central adiposity, characterized by waist-to-hip ratio (WHR), with GSD risk using data from UK Biobank (N=472,050). We then explored the genetic relationship using summary statistics from the largest genome-wide association study of GSD (ncase=43,639, ncontrol=506,798) as well as WHR, with and without adjusting for body mass index (BMI) (WHR: n=697,734; WHRadjBMI: n=694,649).

Observational analysis demonstrated an increased risk of GSD with one unit increase in WHR (HR=1.18, 95%CI=1.14-1.21). A positive WHR-GSD genetic correlation (r g =0.41, P=1.42×10-52) was observed, driven by yet independent of BMI (WHRadjBMI: r g =0.19, P=6.89×10-16). Cross-trait meta-analysis identified four novel pleiotropic loci underlying WHR and GSD with biological mechanisms outside of BMI. Mendelian randomization confirmed a robust WHR-GSD causal relationship (OR=1.50, 95%CI=1.35-1.65) which attenuated yet remained significant after adjusting for BMI (OR=1.17, 95%CI=1.09-1.26). Furthermore, observational analysis confirmed a positive association between general obesity and GSD, corroborated by a shared genetic basis (r g =0.40, P=2.16×10-43), multiple novel pleiotropic loci (N=11) and a causal relationship (OR=1.67, 95%CI=1.56-1.78).

Both observational and genetic analyses consistently provide evidence on an association of central obesity with an increased risk of GSD, independent of general obesity. Our work highlights the need of considering both general and central obesity in the clinical management of GSD.

Nonalcoholic Fatty Liver Disease (NAFLD) is a complex human disease. Common genetic variation in the patatin-like phospholipase domain containing 3 (PNPLA3) and transmembrane 6 superfamily member 2 (TM6SF2) genes have been associated with an increased risk of developing NAFLD, nonalcoholic steatohepatitis (NASH), and fibrosis in adults. The role of rare genetic variants in the development and progression of NAFLD in children is not well known. We aimed to explore the role of rare genetic variants in pediatric patients with advanced fibrosis.

Whole exome sequencing data was generated for 229 pediatric patients diagnosed with NAFLD recruited from the NASH Clinical Research Network (NASH CRN). Case-control single variant and gene-based collapsing analyses were used to test for rare variants that were enriched or depleted within the pediatric NAFLD cohort specifically for advanced fibrosis (cases) versus those without fibrosis (controls) or six other histologic characteristics. Exome data from non-NAFLD population controls were also used for additional analyses. All results were adjusted for multiple testing using a Bonferroni correction.

No genome-wide significant associations were found between rare variation and presence of advanced fibrosis or NASH, nor the severity of steatosis, inflammation, or hepatocellular ballooning. Significantly, no enrichment of rare variants in PNPLA3 or TM6SF2 was observed across phenotypes.

In a cohort of children with histologically proven NAFLD, no genome-wide significant associations were found between rare genetic variation and advanced fibrosis or six other histologic features. Of particular interest was the lack of association with genes of interest in adults: PNPLA3 and TM6SF2, though limitations in sample size may reduce the ability to detect associations, particularly with rare variation.

The investigation of environmental effects on clinical measurements using individual samples is challenging because their genetic and environmental factors are different. However, using monozygotic twins (MZ) makes it possible to investigate the influence of environmental factors as they have the same genetic factors within pairs because the difference in the clinical traits within the MZ mostly reflect the influence of environmental factors. We hypothesized that the within-pair differences in the traits that are strongly affected by genetic factors become larger after genetic risk score (GRS) correction. Using 278 Japanese MZ pairs, we compared the change in within-pair differences in each of the 45 normalized clinical measurements before and after GRS correction, and we also attempted to correct for the effects of genetic factors to identify Cytosine-phosphodiester-Guanine (CpG) sites in DNA sequences with epigenetic effects that are regulated by genetic factors. Five traits were classified into the high heritability group, which was strongly affected by genetic factors. CpG sites could be classified into three groups: regulated only by environmental factors, regulated by environmental factors masked by genetic factors, and regulated only by genetic factors. Our method has the potential to identify trait-related methylation sites that have not yet been discovered.

Sorafenib induces ferroptosis, making it a useful treatment against advanced liver hepatocellular carcinoma (LIHC). However, sorafenib resistance is extremely common among LIHC patients. Here, we used a comprehensive approach to investigate the effects of ABHD12, which regulates tumorigenesis and sorafenib resistance in LIHC. We validated ABHD12 expression was upregulated in LIHC tissue, which correlated with worse overall survival and related to tumor size or stage. ABHD12 facilitated a pro-tumorigenic phenotype involving increased cell proliferation, migration, and clonogenicity as well as sorafenib resistance. Knockout of ABHD12 sensitized liver cancer cells to sorafenib-induced ferroptosis. Co-delivery of sorafenib and ABHD12 inhibitor into a nude mouse model enhanced therapeutic efficacy for LIHC. Our study demonstrates that ABHD12 contributes to tumor growth and sorafenib resistance in liver cancer, which indicate the promising potential of ABHD12 in diagnosis and prognosis as well as highlight the potential therapeutic applications for co-delivery of sorafenib and ABHD12 inhibitor.

Genome-wide association studies (GWASs) have identified tens of thousands of genetic loci associated with human complex traits. However, the majority of GWASs were conducted in individuals of European ancestries. Failure to capture global genetic diversity has limited genomic discovery and has impeded equitable delivery of genomic knowledge to diverse populations. Here we report findings from 102,900 individuals across 36 human quantitative traits in the Taiwan Biobank (TWB), a major biobank effort that broadens the population diversity of genetic studies in East Asia. We identified 968 novel genetic loci, pinpointed novel causal variants through statistical fine-mapping, compared the genetic architecture across TWB, Biobank Japan, and UK Biobank, and evaluated the utility of cross-phenotype, cross-population polygenic risk scores in disease risk prediction. These results demonstrated the potential to advance discovery through diversifying GWAS populations and provided insights into the common genetic basis of human complex traits in East Asia.

Elevated Gamma-glutamyl transferase (GGT) levels are often suggestive of cholelithiasis, and previous studies have indicated that GGT is highly expressed in the urinary system. Therefore, we hypothesized that there may be an association between GGT levels and calculus of kidney (CK) incidence. To investigate this potential causal relationship, we employed Mendelian randomization (MR) analysis. Additionally, we analyzed the levels of other liver enzymes, including alanine transaminase (ALT) and alkaline phosphatase (ALP). The relationship between GGT levels and CK incidence was analyzed using two-sample Mendelian randomization. Summary Genome-Wide Association Studies data were utilized for this analysis. 33 single nucleotide polymorphisms known to be associated with GGT levels were employed as instrumental variables. We employed several MR methods including IVW (inverse variance weighting), MR-Egger, weighted median, weighted mode, and MR-PRESSO (Mendelian Randomization Pleiotropy RESidual Sum and Outlier). Furthermore, we conducted tests for horizontal multivariate validity, heterogeneity, and performed leave-one-out analysis to ensure the stability of the results. Overall, several MR methods yielded statistically significant results with a p-value < 0.05. The results from the IVW analysis yielded an odds ratio (OR) of 1.0062 with a 95% confidence interval (CI) of 1.0016-1.0109 (p = 0.0077). Additional MR methods provided supplementary results: MR-Egger (OR 1.0167, 95% CI 1.0070-1.0266, p = 0.0040); weighted median (OR 1.0058, 95% CI 1.0002-1.0115, p = 0.0423); and weighted mode (OR 1.0083, 95% CI 1.0020-1.0146, p- = 0.0188). Sensitivity analyses did not reveal heterogeneity or outliers. Although potential horizontal pleiotropy emerged, we speculate that this could be attributed to inadequate test efficacy. However, subsequent use of MR-PRESSO did not provide evidence of pleiotropy. Our analysis suggests a positive association between elevated GGT levels and CK incidence, indicating an increased risk of CK development. However, no causal relationship was observed between levels of ALP or ALT and CK incidence.

Both psychiatric vulnerability and cortical structure are shaped by the cumulative effect of common genetic variants across the genome. However, the shared genetic underpinnings between psychiatric disorders and brain structural phenotypes, such as thickness and surface area of the cerebral cortex, remains elusive. In this study, we employed pleiotropy-informed conjunctional false discovery rate analysis to investigate shared loci across genome-wide association scans of regional cortical thickness, surface area, and seven psychiatric disorders in approximately 700,000 individuals of European ancestry. Aggregating regional measures, we identified 50 genetic loci shared between psychiatric disorders and surface area, as well as 26 genetic loci shared with cortical thickness. Risk alleles exhibited bidirectional effects on both cortical thickness and surface area, such that some risk alleles for each disorder increased regional brain size while other risk alleles decreased regional brain size. Due to bidirectional effects, in many cases we observed extensive pleiotropy between an imaging phenotype and a psychiatric disorder even in the absence of a significant genetic correlation between them. The impact of genetic risk for psychiatric disorders on regional brain structure did exhibit a consistent pattern across highly comorbid psychiatric disorders, with 80% of the genetic loci shared across multiple disorders displaying consistent directions of effect. Cortical patterning of genetic overlap revealed a hierarchical genetic architecture, with the association cortex and sensorimotor cortex representing two extremes of shared genetic influence on psychiatric disorders and brain structural variation. Integrating multi-scale functional annotations and transcriptomic profiles, we observed that shared genetic loci were enriched in active genomic regions, converged on neurobiological and metabolic pathways, and showed differential expression in postmortem brain tissue from individuals with psychiatric disorders. Cumulatively, these findings provide a significant advance in our understanding of the overlapping polygenic architecture between psychopathology and cortical brain structure.

Genetic factors influence the risk of fatty liver disease (FLD) in adults. The aim of this study was to test if, and when, genetic risk factors known to affect FLD in adults begin to exert their deleterious effects during childhood, adolescence and early adulthood.

We included up to 4018 British children and adolescents from the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort. Three genetic variants known to associate robustly with FLD in adults (PNPLA3 rs738409, TM6SF2 rs58542926 and HSD17B13 rs72613567) were tested for association with plasma levels of alanine transaminase (ALT) and aspartate transaminase (AST) during childhood (mean age: 9.9 years), early adolescence (15.5 years), late adolescence (17.8 years), and early adulthood (24.5 years). We also tested the associations of a 17-variant score and whole-genome polygenic risk scores (PRS) derived from associations in adults with plasma ALT and AST at the same four time points. Associations with elastography-derived liver steatosis and fibrosis were tested in early adulthood.

Genetic risk factors for FLD (individually, combined into a 3-variant score, a 17-variant score and as a genome-wide PRS), were associated with higher liver enzymes, beginning in childhood and throughout adolescence and early adulthood. The ALT-increasing effects of the genetic risk variants became larger with increasing age. The ALT-PRS was associated with liver steatosis in early adulthood. No genetic associations with fibrosis were observed.

Genetic factors that promote FLD in adults associate with elevated liver enzymes already during childhood, and their effects get amplified with increasing age.

Nonalcoholic fatty liver disease (NAFLD) is common and partially heritable and has no effective treatments. We carried out a genome-wide association study (GWAS) meta-analysis of imaging (n = 66,814) and diagnostic code (3,584 cases versus 621,081 controls) measured NAFLD across diverse ancestries. We identified NAFLD-associated variants at torsin family 1 member B (TOR1B), fat mass and obesity associated (FTO), cordon-bleu WH2 repeat protein like 1 (COBLL1)/growth factor receptor-bound protein 14 (GRB14), insulin receptor (INSR), sterol regulatory element-binding transcription factor 1 (SREBF1) and patatin-like phospholipase domain-containing protein 2 (PNPLA2), as well as validated NAFLD-associated variants at patatin-like phospholipase domain-containing protein 3 (PNPLA3), transmembrane 6 superfamily 2 (TM6SF2), apolipoprotein E (APOE), glucokinase regulator (GCKR), tribbles homolog 1 (TRIB1), glycerol-3-phosphate acyltransferase (GPAM), mitochondrial amidoxime-reducing component 1 (MARC1), microsomal triglyceride transfer protein large subunit (MTTP), alcohol dehydrogenase 1B (ADH1B), transmembrane channel like 4 (TMC4)/membrane-bound O-acyltransferase domain containing 7 (MBOAT7) and receptor-type tyrosine-protein phosphatase δ (PTPRD). Implicated genes highlight mitochondrial, cholesterol and de novo lipogenesis as causally contributing to NAFLD predisposition. Phenome-wide association study (PheWAS) analyses suggest at least seven subtypes of NAFLD. Individuals in the top 10% and 1% of genetic risk have a 2.5-fold to 6-fold increased risk of NAFLD, cirrhosis and hepatocellular carcinoma. These genetic variants identify subtypes of NAFLD, improve estimates of disease risk and can guide the development of targeted therapeutics.