Genotype-First Approach Identifies an Association between rs28374544/FOG2S657G and Liver Disease through Alterations in mTORC1 Signaling

Metabolic dysfunction-associated Fatty Liver Disease (MAFLD) has emerged as one of the leading cardiometabolic diseases. Friend of GATA2 (FOG2) is a transcriptional co-regulator that has been shown to regulate hepatic lipid metabolism and accumulation. Using meta-analysis from several different biobank datasets, we identified a coding variant of FOG2 (rs28374544, A1969G, S657G) predominantly found in individuals of African ancestry (minor allele frequency~20%), which is associated with liver failure/cirrhosis phenotype and liver injury. To gain insight into potential pathways associated with this variant, we interrogated a previously published genomics dataset of 38 human induced pluripotent stem cell (iPSCs) lines differentiated into hepatocytes (iHeps). Using Differential Gene Expression Analysis and Gene Set Enrichment Analysis, we identified the mTORC1 pathway as differentially regulated between iHeps from individuals with and without the variant. Transient lipid-based transfections were performed on the human hepatoma cell line (Huh7) using wild-type FOG2 and FOG2S657G and demonstrated that FOG2S657G increased mTORC1 signaling, de novo lipogenesis, and cellular triglyceride synthesis and mass. In addition, we observed a significant downregulation of oxidative phosphorylation in FOG2S657G cells in fatty acid-loaded cells but not untreated cells, suggesting that FOG2S657G may also reduce fatty acid to promote lipid accumulation. Taken together, our multi-pronged approach suggests a model whereby the FOG2S657G may promote MAFLD through mTORC1 activation, increased de novo lipogenesis, and lipid accumulation. Our results provide insights into the molecular mechanisms by which FOG2S657G may affect the complex molecular landscape underlying MAFLD.


Introduction
Metabolic dysfunction-associated fatty liver disease (MAFLD) is a leading cardiometabolic disease characterized by dysregulation of hepatic lipid metabolism and accumulation of triglyceride (TG) in the liver.If left untreated, it can progress from steatosis to metabolic dysfunction-associated steatohepatitis (MASH), cirrhosis, and, increasingly, hepatocellular carcinoma (HCC).Both globally and in the United States, MAFLD is causing increased disease burden and mortality and will soon be the leading cause of liver transplantation [1].
Friend of GATA 2 (FOG2) is a transcriptional co-regulator encoded by the ZFPM2 gene that is crucial for both development and organ homeostasis in adulthood [2,3].We have recently identified a coding variant of FOG2 (rs28374544, A1969G, Ser657Gly, hereafter referred to as FOG2 S657G ), which is highly present in individuals of African ancestry (MAF~0.20)and associated with coronary microvascular disease [4].In addition to cardiac expression, FOG2 is also present in the developing and adult liver.It has been shown to affect insulin resistance and hepatic fat accumulation through the regulation of Ppara in the liver and to regulate AKT [5].Single nucleotide RNAseq of whole livers in humans found that FOG2 expression was increased in hepatocytes from livers of individuals with MASH as compared to those with MAFLD [6].Higher levels of FOG2 expression were also associated with a more favorable prognosis in a study of hepatocellular carcinoma in humans [7].
Based on these data, we hypothesized that FOG2 S657G may play a role in liver homeostasis.Here, we report human genetic evidence linking FOG2 S657G with liver disease and injury and hepatic lipid metabolism.We used gene expression data from iPSCsdifferentiated to hepatocytes to identify pathways through which FOG2 S657G may regulate hepatocyte biology.We show that FOG2 S657G increases mTORC1 signaling and, in turn, promotes lipid accumulation and affects cellular metabolism.

FOG2 S657G Is Associated with Liver Injury and Disease in Multiple Biobank Populations
In order to determine whether FOGS657G was associated with human liver phenotypes, we first used Genebass, which encompasses 4529 phenotypes with gene-based and single-variant testing across 394,841 individuals with exome sequence data from the UK Biobank (UKBB) (https://app.genebass.orgaccessed on 29 March 2024) [8].There was an association between rs28374544 and liver failure/cirrhosis phenotype (p = 0.0053375) in this population.Given the relatively low representation of African ancestry individuals in the UKBB and the prevalence of FOG2 S657G in individuals of African ancestry, it is possible that the signal that we observed was due to ancestry and not the variant.Therefore, we performed a meta-analysis for rs28374544 in the PennMedicine Biobank (PMBB) [9], All of Us Biobank [10], and the Million Veterans Program (MVP) [11].Datasets include only individuals of African ancestry [12].Here, we observed a significant increase in "cirrhosis of the liver without mention of alcohol" (p = 0.01239, β = +0.068)(Table 1).We next examined the levels of serum proteins in 389 PMBB participants of African ancestry who had undergone genotyping and proteomic profiling with SOMAscan (Soma-Logic) for markers of liver disease and liver injury [13].The mean age of participants was 61 ± 13.6, and 58% were male.We found that there was a 30% increase in plasma alanine transaminase (ALT) with at least one copy of FOG2 S657G (Table 1).Taken together, these data support a positive association between FOG2 S657G and liver disease.

Genomic Analysis from a Cohort of Induced Pluripotent Stem Cell Lines Differentiated into Hepatocytes Idenfies an Association between FOG2 S657G and mTORC1 Pathway
To gain insight into mechanisms by which FOG2 S657G may alter hepatocyte gene expression, we examined the effects of rs28374544 on gene expression in a cohort of human Genes 2024, 15, 1098 3 of 12 induced pluripotent stem cell (iPSC) lines differentiated into hepatocytes (iHeps) previously reported [14].This dataset includes RNAseq data from 38 iHeps lines isolated from selfreported African American individuals (25 females, 13 males), 23 of which carry the reference allele (AA) and 15 of which are heterozygous for the variant (AG) (Figure 1a).We performed differential gene expression analysis and found that 152 genes were upregulated and 131 were downregulated in lines from iHeps with FOG2 S657G at a nominal p-value of 0.5 and log2 fold change of 1.2 (Figure 1b).The top differentially expressed gene was IGFBP2 (fold change = 1.38, p = 2.1 × 10 −4 ), which encodes the Insulin growth factor binding protein 2 (IGFBP2) and has been reported to be associated with MAFLD and progression to MASH [15,16] (Table S1) Other top differentially expressed and abundant genes included ENO1, ALDOA, PKM, and PGK1, key regulators of various metabolic pathways, including notably anaerobic respiration, hypoxia, glycolysis, and oxidative phosphorylation (Table S1).
the reference allele (AA) and 15 of which are heterozygous for the variant (AG) (Figure 1a).We performed differential gene expression analysis and found that 152 genes were upregulated and 131 were downregulated in lines from iHeps with FOG2 S657G at a nominal p-value of 0.5 and log2 fold change of 1.2 (Figure 1b).The top differentially expressed gene was IGFBP2 (fold change = 1.38, p = 2.1 × 10 −4 ), which encodes the Insulin growth factor binding protein 2 (IGFBP2) and has been reported to be associated with MAFLD and progression to MASH [15,16] (Table S1) Other top differentially expressed and abundant genes included ENO1, ALDOA, PKM, and PGK1, key regulators of various metabolic pathways, including notably anaerobic respiration, hypoxia, glycolysis, and oxidative phosphorylation (Table S1).
We then performed Gene Set Enrichment Analysis (GSEA) using the differentially expressed genes and Hallmark gene sets, which represent canonical pathways.The most differentially regulated pathway was mTORC1 signaling (Enrichment score 19.8, p = 2.6 × 10 −9 , Figure 1c).Related pathways, including hypoxia, glycolysis, and oxidative phosphorylation, were also differentially regulated.Closer examination of the genes differentially regulated in the mTORC1 pathway was consistent with an increase in mTORC1 signaling (Figure 1d), which has been shown to affect hepatic lipid metabolism in several ways, including driving de novo lipogenesis (DNL).We then performed Gene Set Enrichment Analysis (GSEA) using the differentially expressed genes and Hallmark gene sets, which represent canonical pathways.The most differentially regulated pathway was mTORC1 signaling (Enrichment score 19.8, p = 2.6 × 10 −9 , Figure 1c).Related pathways, including hypoxia, glycolysis, and oxidative phosphorylation, were also differentially regulated.Closer examination of the genes differentially regulated in the mTORC1 pathway was consistent with an increase in mTORC1 signaling (Figure 1d), which has been shown to affect hepatic lipid metabolism in several ways, including driving de novo lipogenesis (DNL).

Functional In Vitro Studies of FOG2 S657G Overexpression in Huh7 Cells Confirm Regulation of mTORC1, DNL, and TG Synthesis
To gain mechanistic insights into the associations identified via genomic analyses, we turned to in vitro functional studies.To determine the role of FOG2 S657G on mTORC1 signaling, we transfected Huh7 cells with vectors expressing empty vector, wild type (WT) FOG2, or FOG2 S657G under EF1alpha promoter.These resulted in similar levels of gene and protein over-expression (Figure 2a,b).Gene expression of IGFBP2 was significantly increased (1.9-fold, p < 0.001) in cells overexpressing FOG2 S657G , similar in direction and magnitude to what was observed in FOG2 S657G expressing iHeps (Figure S1).Using western blot analysis of whole cell lysates, we found that there was significantly higher phosphorylation of S6 and 4EBP in FOG2 S657G -expressing cells as compared to WT (Figure 2c,d).Increased phosphorylation of these mTORC1 targets is consistent with increased mTORC1 signaling in the cells and increased changes in metabolic pathways that are important for cell survival and proliferation and are often altered in disease states [17].
and female (square) individuals of African American Ancestry in Penn Medicine Biobank.(b) Volcano plot of differentially expressed genes as determined by DeSeq2 analysis.Blue dots indicate significantly downregulated and red dots significantly upregulated (nom p < 0.05 and log2 (Fold change) > 1.2).(c) Gene set enrichment analysis of differentially expressed genes using Hallmark pathways.(d) Heat map of differentially expressed genes in the mTORC1 hallmark pathway consistent with an increase in mTORC1 signaling in individuals with the variant.

Functional In Vitro Studies of FOG2 S657G Overexpression in Huh7 Cells Confirm Regulation of mTORC1, DNL, and TG Synthesis
To gain mechanistic insights into the associations identified via genomic analyses, we turned to in vitro functional studies.To determine the role of FOG2 S657G on mTORC1 signaling, we transfected Huh7 cells with vectors expressing empty vector, wild type (WT) FOG2, or FOG2 S657G under EF1alpha promoter.These resulted in similar levels of gene and protein over-expression (Figure 2a,b).Gene expression of IGFBP2 was significantly increased (1.9-fold, p < 0.001) in cells overexpressing FOG2 S657G , similar in direction and magnitude to what was observed in FOG2 S657G expressing iHeps (Figure S1).Using western blot analysis of whole cell lysates, we found that there was significantly higher phosphorylation of S6 and 4EBP in FOG2 S657G -expressing cells as compared to WT (Figure 2c,d).Increased phosphorylation of these mTORC1 targets is consistent with increased mTORC1 signaling in the cells and increased changes in metabolic pathways that are important for cell survival and proliferation and are often altered in disease states [17].In addition to regulating protein and nucleotide translation, nutrient sensing, cell growth, and autophagy, mTORC1 is also known to promote the cleavage and nuclear localization of SREBP1c to drive the transcription of genes related to DNL [18].To determine if the increase in mTORC1 signaling affected lipid metabolism, we measured SREBP1 protein and observed a 35% increase (p < 0.05) in the active cleaved form of SREBP1 in cells overexpressing FOG2 S657G relative to WT with no difference in total SREBP1 protein (Figure 3a,b).Since the nuclear form of SREBP1 can drive the transcription of genes related to DNL, we measured the relative gene expression of key genes involved in DNL, including FASN, ACACA, and GPAM (Figure 3c-e).This was associated with an increase in cell TG mass (Figure 3f).To determine whether this was due to an increase in TG synthesis in the cells, we labeled cells with 3 H glycerol and measured the incorporation of 3 H glycerol into both TG and phospholipid (PL).There was increased intracellular 3 H-TG in FOG2 S657G -expressing cells, indicating an increase in TG synthesis (Figure 3g).There was also increased 3 H-TG in the media which is consistent with an increase in TG secretion due to more cellular TG availability (Figure 3h).There was no difference in cellular or secreted 3 H-PL, demonstrating that there is a specific effect on TG synthesis and not a result of a difference in label uptake (Figure 3i,j).
SREBP1 protein and observed a 35% increase (p < 0.05) in the active cleaved form of SREBP1 in cells overexpressing FOG2 S657G relative to WT with no difference in total SREBP1 protein (Figure 3a,b).Since the nuclear form of SREBP1 can drive the transcription of genes related to DNL, we measured the relative gene expression of key genes involved in DNL, including FASN, ACACA, and GPAM (Figure 3c-e).This was associated with an increase in cell TG mass (Figure 3f).To determine whether this was due to an increase in TG synthesis in the cells, we labeled cells with 3 H glycerol and measured the incorporation of 3 H glycerol into both TG and phospholipid (PL).There was increased intracellular 3 H-TG in FOG2 S657G -expressing cells, indicating an increase in TG synthesis (Figure 3g).There was also increased 3 H-TG in the media which is consistent with an increase in TG secretion due to more cellular TG availability (Figure 3h).There was no difference in cellular or secreted 3 H-PL, demonstrating that there is a specific effect on TG synthesis and not a result of a difference in label uptake (Figure 3i,j).

FOG2 S657G Cells Have Reduced Oxidative Phosphorylation in iHeps and Huh7 Cells
Another differentially regulated pathway in iHep lines from individuals with the FOG2 S657G variant was oxidative phosphorylation (Figure 1c), and the directionality of the differentially expressed genes of this pathway was suggestive of decreased oxidative phosphorylation (Figure 4a).A reduction in oxidative phosphorylation can also contribute to lipid accumulation.We, therefore, measured the oxygen consumption rate (OCR) in transfected cells using the Seahorse bioanalyzer.We found that both FOG2 and FOG2 S657G transfected cells had lower OCR as compared to EV control cells under basal conditions, and there was no significant difference between FOG2 S657G -expressing cells and WT

FOG2 S657G Cells Have Reduced Oxidative Phosphorylation in iHeps and Huh7 Cells
Another differentially regulated pathway in iHep lines from individuals with the FOG2 S657G variant was oxidative phosphorylation (Figure 1c), and the directionality of the differentially expressed genes of this pathway was suggestive of decreased oxidative phosphorylation (Figure 4a).A reduction in oxidative phosphorylation can also contribute to lipid accumulation.We, therefore, measured the oxygen consumption rate (OCR) in transfected cells using the Seahorse bioanalyzer.We found that both FOG2 and FOG2 S657G transfected cells had lower OCR as compared to EV control cells under basal conditions, and there was no significant difference between FOG2 S657G -expressing cells and WT control (Figure 4b,c).However, when cells were lipid-loaded with 0.2 mM oleic acid (OA) for 16 h prior to the assay, we noted a decrease in both basal and maximal respiration in FOG2 S657Gexpressing cells relative to control (Figure 4d,e).This decrease in oxidative phosphorylation in the presence of fatty acid loading suggests that there may be a decrease in fatty acid oxidation, which may contribute to the lipid accumulation phenotype observed in the cells.Prior studies have shown that Fog2 may regulate Ppara expression and affect fatty acid oxidation in the mouse liver.However, we found no difference in the expression of Ppara or its target CPT1 (Figure 4f,g).

Discussion
We have identified an association between FOG2 S657G and liver disease/injury and increased plasma ALTs.Genomic analysis from multiple human iPSC lines differentiated to iHeps suggested an association between FOG2 S657G and metabolic pathways, including mTORC1 signaling, hypoxia, glycolysis, and oxidative phosphorylation.In vitro functional studies showed that FOG2 S657G can activate hepatocyte mTORC1 to increase DNL, increase TG synthesis, and decrease oxidative phosphorylation.Together, these act to increase hepatic TG accumulation.
In addition to regulating hepatic lipid metabolism [5], FOG2 has also been shown to affect metabolic pathways in other cell types.A mutation in Fog2's LXCX retinoblastoma protein binding motif has been shown to regulate metabolic processes as adipogenesis, weight gain, and small intestine homeostasis in mice [19,20].Guo et al. demonstrated that overexpression of Fog2 in the liver protected mice from the development of fatty liver but, interestingly, was associated with an increase in insulin resistance [5].The reciprocal results of increased fatty liver and protection from insulin resistance were observed with Fog2 siRNA knockdown experiments and these effects were mediated through altered expression of Ppara in the liver.Interestingly, this group also observed that Fog2 could affect the phosphorylation of AKT (an upstream regulator of mTORC1) and thus alter the insulin signaling pathway [5].This is consistent with previous results that showed that the drosophila homolog of Fog2, U-Shaped, and human FOG2 inhibits phosphoinositide-3 kinase (PI3K) and results in decreased phosphorylation of AKT and reduced cell growth in drosophila and human Hep3B cells [21].Additionally, they showed that the middle region of FOG2 from amino acid 507-789 is the domain that binds directly to the P85 subunit of PI3K and inhibits its activity.The FOG2 S657G mutation is in this region and suggests that the effect of the amino acid substitution might be to alter the interaction of FOG2 with PI3K in the cytosol.If FOG2 S657G prevented the interaction with PI3K, then this might result in increased activity, increased AKT phosphorylation, and ultimately increased activation of the downstream mTORC1 pathways.
Based on the results of Guo et al., FOG2 S657G -based lipid accumulation would result in a loss of function.FOG2 S657G has previously been reported to serve as loss-of-function in failing to repress GATA4 in the heart [22].However, activation of the mTORC1 pathway in FOG2 S657G hepatocytes is more consistent with the development of insulin resistance, which was seen with the FOG2 overexpression [5].These discrepant results suggest more complex interactions of FOG2 in hepatocytes.For example, FOG2 S657G could mediate changes via transcriptional regulation in the nucleus while at the same time exerting different effects in the cytoplasm through regulation of AKT and downstream pathways.Therefore, it is possible that hepatocyte FOG2 S657G results in a change of function rather than a simple loss or gain of function in regard to hepatocyte metabolism.
Prior human GWAS studies for hepatic steatosis have not identified this variant [23,24].One potential reason is that FOG2 S657G is primarily present in individuals of African ancestry.There are known disparities in the diagnosis of MAFLD, with lower reported incidence of MAFLD in African American relative to Hispanic and Caucasian individuals [25,26].African American individuals have also been traditionally underrepresented in population-based studies of MAFLD [27].Most MAFLD studies are relatively small due to the diagnostic difficulty and include an even smaller number of patients of African ancestry, limiting the ability to detect this specific variant.For example, only 1-2% of individuals in the UK Biobank have African ancestry, limiting the discovery of relevant new variants in this population [28].More recent studies in a larger multi-ancestry cohort also did not identify this variant [29].Additionally, analyses using biobank data are subject to participation bias, particularly when individuals are from a healthcare population, such as PMBB or MVP.These populations may have a higher incidence of risk factors for liver disease, which can confound the signal for liver disease.However, we performed a meta-analysis using three different biobanks to help mediate the effects of bias.Furthermore, our in vitro studies demonstrate that there is biological plausibility for this variant altering hepatic lipid metabolism and contributing to liver disease.
Our analysis identified an association between liver disease phenotypes and ALT, which are more consistent with MASH.Our data from human biobanks and cell lines suggest mild TG accumulation with FOG2 S657G .This phenotype is subtle, consistent with a point mutation resulting in a single amino acid change that is predicted to generate a slight conformational change [4].It is clear that this signal is not driven by steatosis alone.This may be because the variant is exerting only minor effects on lipid metabolism and steatosis, or alternately, steatosis may not be driving liver injury.FOG2 S657G may be driving metabolic changes in other tissues and cell types that may contribute to the development of insulin resistance in the liver.For example, our analysis identified IGFBP2 as differentially regulated by FOG2 S657G in both iHeps and Huh7 cells.IGFBP2 has recently been identified to be the second-most differentially regulated protein in the progression from MAFLD to cirrhosis in an Icelandic population [30].The mechanism by which IGFBP2 may affect this progression is not known; however, recent data suggest that zone 2 hepatocytes that express IGFBP2 play key roles in liver regeneration [31].Additional studies are needed to understand the interactions between FOG2 S657G , IGFBP2, and mechanisms of progression from MAFLD to cirrhosis.
Future studies in pre-clinical models of FOG2 S657G are needed to understand the mechanisms by which FOG2 S5657G drives hepatocyte lipid accumulation and may promote progression to liver injury.This would allow for the examination of other liver cell types and extra-hepatic tissues and their effect on liver injury.Currently, there are limited treatment options for MAFLD and MASH, and more investigation into disease pathogenesis is needed [32].Understanding the mechanisms by which steatosis can develop, and progress is critical to developing new treatments.We have shown that FOG2 S657G is associated with mild TG accumulation and downstream sequelae such as elevated ALT and liver disease phenotypes.Elucidating the role of FOG2 and FOG2 S657G , specifically in the liver, may reveal new pathways to target and novel treatment avenues.

Meta-Analysis of Biobanks
PheWAS for rs28354744 was performed in PMBB, All of Us, and MVP in individuals of African ancestry using SAIGE [12].Each analysis included 10,995, 50,457, and 120,039 participants and 959, 1472, 945 phecodes, respectively.Fixed-effect meta-analysis was performed using the metagen function in the R package meta (version 6.5-0) [33].Results from the PheWAS meta-analysis were combined with their respective phecodes using the R package PheWAS (version 0.99.5-4).

Plasma Protein Measurements
Plasma was collected from 389 African American individuals enrolled in the PMBB and stored at −80 • C before shipping on dry ice to SomaLogic, Inc. (Boulder, CO, USA) for measurement of the relative concentration of 1305 proteins using the multiplexed, aptamer-based SOMAscan assay, as previously described [13].Data were normalized and log-transformed, and plasma ALT levels in the 102 carriers of rs28354744 were compared to non-carriers using linear regression analysis.

Analysis of Published iPSC-Hepatocyte Data
Raw count data were downloaded from GEO and annotated using biomaRt and EnsDb.Hsapiens.v86packages in R 4.2.2.Partek Flow (11.0.23.0918) was used for all genomic analyses, including principal component, differential gene expression, pathway, and clustering analyses.DEseq2 was used to normalize data, generate normalized counts, and perform differential gene expression analysis.Gene Set Enrichment Analysis was performed for genes with a nominal p-value < 0.05 and log2(Fold change) > 1.2, using Hallmark pathways which represent canonical pathways.We report GSEA enrichment scores and all pathways with nominal p value of less than 0.05.Heatmaps represent hierarchical clustering using average linkage Euclidean distance.

Huh7 Cell Culture and Transfection
HuH-7 cells were obtained from the JCRB cell bank and maintained in DMEM with 10% fetal bovine serum (FBS) and 1% antibiotics.Prior to the assay, cells were plated collagen coated 6-well plates.Cells were transfected with empty vector (EV), WT FOG2, or FOG2 S657G (Twist Biosciences, San Francisco, CA, USA) using commercially available reagents (Lipofectamine 2000, OptiMEM, Invitrogen Carlsbad, CA, USA) according to manufacturer instructions.Experiments were performed 72 h after transfection.

Gene Expression
Gene expression was measured by quantitative real-time PCR (qRT-PCR).In brief, RNA was extracted from cell samples with a Qiagen RNA Easy Protect Kit, and 1 µg of RNA was reversed transcribed to cDNA using a High-Capacity cDNA Reverse Transcription Kit (Thermo, Carlsbad, CA, USA).The cDNA was diluted 30 times and used for qRT-PCR analysis using PowerUp SYBR Green Master Mix (Applied Biosystems, Carlsbad, CA, USA) and indicated primers.All samples analyzed were normalized to the average of 3 housekeeping using the delta Ct method and normalized to the empty vector control sample.

Western Blot
Cells were lysed in 1 mL of lysis buffer containing 62.5 mM sucrose, 0.5% sodium deoxycholate, 0.5% Triton X-100, 50 mM Tris-HCl, pH 7.4, 150 mM NaCl, 50 µg/mL leupeptin, 50 µg/mL pepstatin A and 30 µL/mL of a protease inhibitor.Protein concentration was measured by BCA protein reagent.Samples were subject to electrophoresis on 4-12% Bis-Tris gels (Thermo Fisher, Carlsbad, CA, USA) and transferred to nitrocellulose membranes.Membranes were incubated with the following antibodies at manufacturerrecommended dilutions: Protein bands were visualized using the LICOR Odyssey Fc imager (LI-COR).Bands were quantified using Image J software.Phospho-S6, Phospho S6K, and Phospho-4EBP were normalized to their respective total proteins.SREBP1c was normalized to actin.

TG Mass Assay
For detection of TGs by colorimetric enzymatic assay, HuH-7 cells were washed with PBS and incubated with 3:2 Hexane: Isopropanol for 2 h.The extracted lipid was dried under N 2 and reconstituted in 300 µL of 15% Triton X-100 in chloroform and dried under N 2 and reconstituted in 300 µL of water.Then, 20 µL of extracted lipid was used for colorimetric TG measurement by enzymatic assay (Thermo Fisher, Carlsbad, CA, USA) using commercial lipid standard.Total protein was recovered by solubilizing the cells in 0.1 N NaOH and measured using bicinchoninic acid (BCA) assay (Pierce, Thermo Fisher Scientific, Carlsbad, CA, USA).

Figure 1 .
Figure 1.RNAseq analysis of iHeps from individuals with and without FOG2 S657G identifies mTORC1 signaling as the top differentially regulated pathway.(a) Principal component analysis of the normalized RNAseq data in AA (black) and AG (crimson) iHeps derived from both male (round) and female (square) individuals of African American Ancestry in Penn Medicine Biobank.(b) Volcano plot of differentially expressed genes as determined by DeSeq2 analysis.Blue dots indicate significantly downregulated and red dots significantly upregulated (nom p < 0.05 and log2 (Fold change) > 1.2).(c) Gene set enrichment analysis of differentially expressed genes using Hallmark pathways.(d) Heat map of differentially expressed genes in the mTORC1 hallmark pathway consistent with an increase in mTORC1 signaling in individuals with the variant.

Figure 2 .
Figure 2. FOG2 S657G increases mTORC1 signaling in Huh7 cells.Huh7 cells were transiently transfected with vectors expressing control empty vector (EV), FOG2 (WT), or FOG2 S657G (S657G).(a,b) Gene and protein expression measured 72 h post-transfection, using qRT-PCR and Western blot, respectively, confirms increased and similar expression of FOG2 WT and S657G relative to EV. ns= not significant (c) Representative Western Blot of phosphorylated and total fractions of mTORC1 targets in transfected Huh7 cells.(d) Quantification of the ratio of phosphorylated to total protein using Image J analysis, n = 4/group * p < 0.05, ** p < 0.01 by Student's t-test.

Figure 2 .
Figure 2. FOG2 S657G increases mTORC1 signaling in Huh7 cells.Huh7 cells were transiently transfected with vectors expressing control empty vector (EV), FOG2 (WT), or FOG2 S657G (S657G).(a,b) Gene and protein expression measured 72 h post-transfection, using qRT-PCR and Western blot, respectively, confirms increased and similar expression of FOG2 WT and S657G relative to EV. ns = not significant (c) Representative Western Blot of phosphorylated and total fractions of mTORC1 targets in transfected Huh7 cells.(d) Quantification of the ratio of phosphorylated to total protein using Image J analysis, n = 4/group * p < 0.05, ** p < 0.01 by Student's t-test.

Figure 3 .
Figure 3. FOG2 S657G expression in Huh7 cells results in increased de novo lipogenesis and TG synthesis.(a) Representative Western blot of total cell lysate for SREBP1c levels.(b) Quantification of the amount of cleaved SREBP1c as normalized to actin levels using ImageJ, n = 4/group.(c-e) RNA expression of DNL genes in WT and S657G transfected Huh7 cells show significantly increased expression of key regulators of DNL.n = 6/group (f) Total cell TG mass was normalized to total protein as measured by BCA and expressed as triglyceride per mg protein.n= 6/group (g-j) Transfected cells were labeled overnight with 3H glycerol, and the total amount of 3H TG in the cell (g) and the media (h) and 3Hphospholipid in cell (i) and media (j) was measured and normalized to total cell protein.n = 4/group, * p < 0.05 by Student's t-test.

Figure 3 .
Figure 3. FOG2 S657G expression in Huh7 cells results in increased de novo lipogenesis and TG synthesis.(a) Representative Western blot of total cell lysate for SREBP1c levels.(b) Quantification of the amount of cleaved SREBP1c as normalized to actin levels using ImageJ, n = 4/group.(c-e) RNA expression of DNL genes in WT and S657G transfected Huh7 cells show significantly increased expression of key regulators of DNL.n = 6/group (f) Total cell TG mass was normalized to total protein as measured by BCA and expressed as triglyceride per mg protein.n= 6/group (g-j) Transfected cells were labeled overnight with 3H glycerol, and the total amount of 3H TG in the cell (g) and the media (h) and 3Hphospholipid in cell (i) and media (j) was measured and normalized to total cell protein.n = 4/group, * p < 0.05 by Student's t-test.

Figure 4 .
Figure 4. Decreased oxidative phosphorylation with FOG2 S657G expression.(a) Heat map of differentially expressed genes in the iHeps in the Hallmark oxidative phosphorylation pathway is consistent with decreased activity in the pathway.(b-e) Seahorse Bioanalyzer data in Huh7 cells transfected with FOG2 or FOG2 S657G .Oxygen consumption rate (OCR) is measured in cells over 75 min in non-treated cells (b,c) other transfected cells were incubated with 0.2 mM OA overnight (d,e) prior to measurement.All conditions represent 10-12 wells * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001, ns = not significant (f,g) RNA expression of PPARA and CPT1 genes in WT and S657G transfected Huh7 cells.

Table 1 .
Summary of liver-associated phenotypes associated with rs28374544.