Graphical Abstract graphic file with name fx1.jpg [53]Open in a new tab Highlights * • LPO-exposed hPSC-derived hepatocyte-like cells resemble human hepatic steatosis * • Steatosis is associated with decreased mitochondrial maximal respiration * • Purine nucleotide cycle activity increases in steatotic hepatocytes * • Excess fumarate production is associated with induction of macrovesicular steatosis __________________________________________________________________ Endocrine System Physiology; Stem Cells Research; Transcriptomics Introduction Nonalcoholic fatty liver disease (NAFLD) is the most common form of hepatic disease and is strongly associated with obesity and type 2 diabetes ([54]WHO, 2006). In the earliest stage of NAFLD, triglyceride (TG) accumulation in hepatocytes leads to the development of hepatic steatosis ([55]Valenti et al., 2016). This is characterized by the development of macrovesicular steatosis, with TG storage in large lipid droplets ([56]Wang and Yu, 2016). Although steatosis is largely benign, it can progress to nonalcoholic steatohepatitis (NASH), which, in turn, increases the risk of developing cirrhosis and hepatocellular carcinoma ([57]Asrih and Jornayvaz, 2015). However, the mechanism(s) underlying this progression have yet to be determined. At present, there are no specific therapeutics available to reverse or treat NAFLD, and the only effective intervention is through the reduction of obesity through diet and exercise or following bariatric surgery ([58]Laursen et al., 2019). In order to understand the pathophysiology of the NAFLD disease spectrum, studies have been undertaken in rodent models, in vitro culture systems, and human liver biopsy samples. Dietary intervention or genetic manipulation leading to the development of obesity and/or hepatic steatosis in rodents often fails to fully recapitulate the NAFLD phenotype. For example, although the leptin-deficient Ob/Ob mouse develops extreme obesity and metabolic abnormalities, the profound abnormalities in circulating leptin mean that this model is not representative of human disease in the general population ([59]Anstee and Goldin, 2006; [60]Rotundo et al., 2018). Likewise, dietary interventions mimic specific stages of NAFLD, such as steatosis or inflammation, but do not reflect the disease progression that occurs in humans ([61]Oligschlaeger and Shiri-Sverdlov, 2020). A common dietary intervention is the methyl-donor-deficient (MDD) diet, which recapitulates the inflammatory responses associated with NASH but not the other elements of the metabolic syndrome ([62]Lyall et al., 2017; [63]Rinella and Green, 2004). In vitro approaches commonly utilize immortalized cell lines derived from hepatocellular carcinoma; however, malignant transformation is known to induce profound changes in metabolic phenotype ([64]Huang et al., 2013; [65]Pavlova and Thompson, 2016), limiting their utility for modeling NAFLD metabolism. Using whole liver tissue from rodent models or humans to study metabolism in distinct cell populations necessitates disaggregation and sorting, a process that induces oxidative stress, potentially confounding analyses ([66]Llufrio et al., 2018). Analysis of primary cells may also be confounded by the cell sorting process ([67]Llufrio et al., 2018). In terms of mechanisms, studies in animal models and humans have shown that NAFLD is associated with mitochondrial dysfunction including altered electron transport chain (ETC) activity ([68]Koliaki et al., 2015; [69]Rector et al., 2010), altered tricarboxylic acid (TCA) cycle activity and anaplerosis ([70]Satapati et al., 2012; [71]Sunny et al., 2011), induced oxidative metabolism, with a proportional increase in oxidative stress ([72]Satapati et al., 2015) and changes in mitochondrial respiration ([73]Koliaki et al., 2015). In these studies, the presence of multiple cell types in tissue samples may confound the results, and this is further complicated by the changes in cell proportions that may occur in disease states. Here, we present a human-relevant model of hepatic steatosis that can be used for high-resolution analysis of metabolic function and that may provide novel insights into human NAFLD. We hypothesized that nutrient excess in NAFLD leads to increases in ETC activity due to increased availability of respiratory substrates. We aimed to determine the impact of lipid accumulation on mitochondrial respiration specifically in hepatocytes, to avoid the confounding effects associated with bulk tissue analysis. Results Lactate, pyruvate, and octanoate treatment induces macrovesicular steatosis in hepatocyte-like cells We and others have previously reported that the high-energy substrate cocktail LPO promotes intracellular lipid accumulation in pluripotent stem-cell-derived hepatocyte-like cells (HLCs) ([74]Lyall et al., 2018) (schematic of this process is shown in [75]Figure 1A) and human hepatocyte or hepatoblastoma cell lines ([76]Gilchrist et al., 2010; [77]Lockman et al., 2012, [78]2016). However, it was not known whether LPO promotes lipid droplet biogenesis or increased accumulation of lipids within existing droplets. Stem-cell-derived HLCs were differentiated and characterized by measuring mRNA levels of a marker of pluripotency (NANOG) and markers of hepatocytes (Albumin; ALB and hepatocyte nuclear factor-4-alpha; HNF4A) ([79]Wang et al., 2017) ([80]Figure 1B). CYP3A4 activity was present in control HLCs and those challenged with LPO, demonstrating that they retained hepatocyte function ([81]Figure 1C). Analysis of gene expression using RT-qPCR showed increased expression of transcripts that are typically associated with intracellular lipid droplet membranes: perilipin 2 (PLIN2), PLIN4, and PLIN5, suggesting an increase in lipid droplet size ([82]Figures 1D and 1F). This was supported by analysis of lipid accumulation using high content imaging, which demonstrated that 48 h of LPO exposure was not associated with a change in the number of intracellular lipid droplets ([83]Figure 1G) but led to a ~2-fold increase in their size ([84]Figure 1H), consistent with the development of macrovesicular steatosis. Figure 1. [85]Figure 1 [86]Open in a new tab LPO treatment induces development of macrovesicular lipid droplets, similar to benign steatosis in humans (A) Schematic overview showing differentiation process and induction of lipid accumulation. (B) H9 hPSCs differentiate to HLCs, losing pluripotency as shown by diminished expression of the pluripotency marker NANOG during the differentiation process. Concomitantly, these cells begin expressing the hepatocyte markers albumin (ALB) and hepatocyte nuclear factor-4-alpha (HNF4A) (n = 3 biological replicates/group). (C) Differentiated HLCs display robust CYP3A4 activity, which is not diminished following challenge with LPO. (D–F) Perilipin 2 (PLIN2), PLIN4, and PLIN5 expression increases in response to LPO treatment (n = 3 biological replicates per group). (G) Number of lipid droplets do not increase, (H) But intracellular lipid droplets increase in size (n = 32 biological replicates/group). Data were analyzed using two-tailed Student t test and expressed as mean ± SD, ∗p < 0.05, ∗∗∗∗p < 0.0001. LPO treatment induces intracellular lipid accumulation and transcriptomic alterations in key mitochondrial respiratory pathways Induction of steatosis in HLCs was associated with widespread transcriptomic changes ([87]Figure 2A). Differential gene expression analysis identified 853 downregulated and 826 upregulated genes (log[2] fold change cut-off >1.5) in LPO exposed cells compared with controls. A selection of candidate genes was validated by RT-qPCR ([88]Figure S1), based on genes that we previously found to be dysregulated in association with steatosis ([89]Lyall et al., 2018). RNA-seq analysis revealed altered expression of a number of genes previously described as having functional roles in the progression of NAFLD, including PLIN2, PPARGC1A, CYP7A1, and HMGCS2 ([90]Figure 2B). Mapping genes with a log[2] fold change >1.5 to the KEGG pathway database identified a number of enriched pathways ([91]Figures 2C and [92]1D), including those related to steroid hormone biosynthesis (18 genes) and ascorbate and aldarate metabolism (7 genes) ([93]Table S4). There was extensive downregulation of the histone structural units H1, H2A/B, H3, and H4 and UDP-glucuronosyltransferases, which were enriched in multiple pathways. Specifically comparing transcriptomic data with the KEGG terms “TCA cycle” and “oxidative phosphorylation” revealed marked gene expression changes within these pathways ([94]Figures 2E and 2F) including in the majority of genes encoding enzymes that catalyze metabolite interconversion ([95]Table S5). Analysis of the oxidative phosphorylation pathway revealed an overall downregulation of transcription of genes encoding components of respiratory complexes I (ND1, ND2, ND4L, NDUFS2, NDUFV1, NDUFA10, NDUFB2 NDUFA2, NDUFB10) and IV (COX1, COX4I1, COX8A, COX6B1), as well as ATP synthase (ATP6V1A, ATP5MF, ATP6V0D1, ATP6V1E1, ATP6, APT5MC2, ATP5MC1) ([96]Table S6). Figure 2. [97]Figure 2 [98]Open in a new tab Steatosis in hepatocytes is associated with transcriptional rewiring reflecting human NAFLD (A) Heatmap analysis of transcriptional changes associated with macrovesicular steatosis. (B) Expression of key NAFLD-associated genes is disrupted in steatotic HLCs. (C–D) Pathway enrichment analysis reveals a number of disrupted pathways associated with macrovesicular steatosis. (E–F) Analysis of KEGG TCA cycle and oxidative phosphorylation pathways reveals extensive disruption of expression of key genes. Data in (A) were analyzed using a two-tailed Student's t test. Data in (B) were analyzed by two-way ANOVA with Sidak post-hoc testing. n = 3 biological replicates per group. Data are expressed as mean ± SD. Macrovesicular steatosis in hepatocytes is associated with electron transport chain dysfunction It has been suggested that NAFLD-associated macrovesicular steatosis results in impairment of mitochondrial respiration ([99]Koliaki et al., 2015), and we therefore proceeded to analyze this in steatotic HLCs ([100]Figure 3A). Basal oxygen consumption rate (OCR), representing combined mitochondrial and non-mitochondrial oxygen consumption, was unchanged following LPO exposure ([101]Figure 3B). In addition, oligomycin A, a complex V inhibitor decreased OCR equally effectively in control and treatment group. Firstly, this indicated that there were no changes in ATP-linked respiration in response to macrovesicular steatosis ([102]Figure 3C). Secondly, when comparing oligomycin-A-induced alterations in OCR with those following addition of rotenone/antimycin A, we could detect no changes in proton leak between groups ([103]Figure 3D). The addition of the ETC uncoupler FCCP revealed a decrease in maximal respiration in the steatotic HLCs, suggesting mild ETC dysfunction ([104]Figure 3E). Using OCR measurements following FCCP treatment, and comparing with the basal OCR, we calculated that there was a decrease in reserve capacity in the LPO-treated cells, compared with controls ([105]Figure 3F). Subsequently, complex I and III were targeted with rotenone/antimycin A, to completely inhibit oxidative phosphorylation, which reduced OCR to a similar level in both groups, suggesting no difference in non-mitochondrial sources of OCR between the control and steatotic HLCs. Figure 3. [106]Figure 3 [107]Open in a new tab LPO treatment results in decreased mitochondrial maximal respiration and reserve capacity, suggestive of respiratory complex impairment. HLCs were injected sequentially from ports A–C with 4 μM oligomycin, 0.5–1 μM FCCP, and 1 μM antimycin combined with 200 nM rotenone. (A) Raw trace of OCR comparing cells with or without LPO treatment; (B) basal respiration (C) ATP-linked respiration and (D) proton leak were unchanged in LPO-treated HLCs. In contrast, maximal respiration (E) and reserve capacity (F) were both diminished following LPO treatment (n = 18 and 16 biological replicates in the control and treatment group, respectively, from two plates). (G) High-content microscopy revealed no changes in fluorescence of mitochondrial in response to treatment (n = 4 and 6 biological replicates in the control and treatment group, respectively). Representative images of mitochondrial content in control and LPO-treated groups 10x magnification. Blue staining = NucBlue, red staining = MitoTracker Deep Red. Scale bar: 25 μm (H). There were no changes in abundance of the mitochondrial proteins succinate dehydrogenase subunit (SDHA) or pyruvate dehydrogenase α1 and α2 subunits (PDH) (n = 6 biological replicates/group). (I) There were no significant changes in the ratio of mitochondrial to genomic DNA, when comparing mitochondrial region 2 (MT2); mitochondrial region 3 (MT3); or genomic region beta-globin (HBB), indicating no alterations in the quantity of mitochondria in response to treatment (n = 4 biological replicates/group). Levels of citrate synthase activity appears unaltered by LPO treatment (n = 3 biological replicates/group). Data were analyzed using two-tailed Student t test for parametric data, or Mann-Whitney U test (D, I) for non-parametric data, and expressed as mean ± SD. Induction of macrovesicular steatosis is not associated with altered mitochondrial number or mitochondrial integrity We wanted to analyze whether changes in ETC function, as demonstrated by decreased cellular maximal respiration and reserve capacity, were due to altered numbers of mitochondria. Because these were the only respiratory measurements to change, we hypothesized that this was not the case. To test this, we used a high-content microscopy-based approach to measure intracellular mitochondrial fluorescence and observed no changes following exposure to LPO ([108]Figure 3G), indicating no change in quantity. We then measured protein levels of succinate dehydrogenase subunit A (SDHA) and pyruvate dehydrogenase (PDH) α1 and α2 subunits ([109]Figure 3H; representative blot [110]Figure S2), which localize to the mitochondria. Protein levels of SDHA and PDH did not change in response to intracellular lipid accumulation. Furthermore, there were no changes in mitochondrial or nuclear DNA ([111]Figure S3) and no change in the ratio of nuclear to mitochondrial DNA ([112]Figure 3I), suggesting that mitochondrial number is not altered by LPO-induced macrovesicular steatosis. To analyze whether the observed changes in OCR could be due to diminished mitochondrial integrity rather than ETC dysfunction, we measured citrate synthase (CS) activity. This has previously been described as a measure of mitochondrial integrity ([113]Boutagy et al., 2015; [114]Short et al., 2005). Our findings suggest that LPO treatment has no marked impact on CS activity ([115]Figure 3J), suggesting that mitochondria remain intact following intracellular lipid accumulation. Taken together, these data demonstrate that mitochondrial number was not altered in response to macrovesicular steatosis. This suggests that the altered respiratory complex activity may arise due to reactive oxygen species (ROS)-mediated dysfunction, rather than changes in protein content, although this requires further investigation. Induction of an NAFLD-like phenotype in HLCs is associated with increased pyruvate carboxylase activity To investigate TCA cycle dynamics in HLCs, lactate in the LPO cocktail was replaced with ^13C[3]-lactate. Lactate can enter the TCA cycle via mitochondrial pyruvate metabolism, through either pyruvate carboxylase (PC) or pyruvate dehydrogenase (PDH). Alternatively, pyruvate can be transported into the cytosol and used in gluconeogenesis ([116]Figure 4A). LPO exposure was associated with increases in steady state levels of pyruvate, aspartate, and citrate, but there were no changes in metabolites associated with the gluconeogenesis pathway ([117]Figure 4B). Although pyruvate was added as part of the LPO cocktail, increased pyruvate generation also occurred as a result of increased lactate dehydrogenase activity ([118]Figure 4C). Although we were unable to directly measure oxaloacetate levels, we used aspartate as a surrogate and found increased m+3 labeling, demonstrating direct synthesis of aspartate from pyruvate ([119]Figure 4D). Furthermore, isotopomer labeling of citrate indicated increased incorporation from pyruvate-derived oxaloacetate and decreased conversion of labeled acetyl-CoA by PDH ([120]Figure 4E). An alternative pathway for pyruvate is conversion to alanine, and although there were no changes in steady state alanine levels, there were moderate, but significant, changes in the incorporation of ^13C ([121]Figure 4F), indicating overall transamination of pyruvate toward alanine. Despite there being no increase in steady state levels of gluconeogenesis-associated metabolites, the increased m+2 labeling of 3-PG ([122]Figure 4G), serine ([123]Figure 4H), and glycine ([124]Figure 4I) suggests increased flux through the gluconeogenesis pathway. Figure 4. [125]Figure 4 [126]Open in a new tab Macrovesicular steatosis is associated with increased PC activity, leading to preferential anaplerosis of pyruvate into the TCA cycle (A) Schematic outlining conversion of lactate to pyruvate and routes by which this can be converted. Black circles denote carbon atoms. (B) Steady state measurements of pyruvate and metabolites that pyruvate can be converted into. Isotopomer labeling patterns of pyruvate (C), aspartate (D), and citrate (E) show anaplerosis into the TCA cycle. There is minimal conversion of pyruvate to alanine (F–G), but increased flux of the m+2 isotopomer through metabolites related to gluconeogenesis, with sustained conversion to 3-PG, serine, and glycine (H–K). For NMR data, labels are: 1 = ^13C labeling; 0 = no ^13C labeling. All GC-MS data consisted of 10 biological replicates and 2 technical replicates. Isotopomer data were calculated by multiplying MID (multiple ion detection) by normalized total ion count. For NMR data (G and K) n = 4 biological replicates/group. Data were analyzed by two-way ANOVA with Sidak post-hoc testing, or (I–J) two-tailed Student's t test. Data are expressed as mean ± SD. Steatosis in HLCs is associated with increased TCA cycle anaplerosis Isotopic labeling of TCA metabolites can produce a number of different isotopomers, depending on the directionality of metabolite synthesis ([127]Figure 5A). In steatotic HLCs, we observed increases in steady state levels of the TCA cycle metabolites αKG, fumarate, and malate ([128]Figure 5B). In contrast, there were no changes in the steady state levels of succinate and no evidence for increased metabolism of glutamate. To further analyze alterations in TCA cycle dynamics we measured ^13C incorporation into metabolites. Although there was increased ^13C incorporation into PC-derived citrate ([129]Figure 4E), this was not the case for αKG ([130]Figure 5C). The decreased levels of ^13C incorporation into glutamate suggests that this does not result from increased cataplerosis through glutamate ([131]Figure 5D). The increased expression of OGDHL and decreased expression of SUCG1 ([132]Table S5) suggest that steatosis may be associated with metabolic reprogramming within the mitochondrial matrix supported by altered protein synthesis. Although decreased incorporation of ^13C into succinate was observed ([133]Figure 5E), a parallel increase in ^13C incorporation into fumarate, malate, and aspartate was also apparent ([134]Figures 5F and 5G). Because the overall directionality of the malate dehydrogenase and fumarate hydratase reactions can vary depending on the metabolic environment ([135]Dasika et al., 2015; [136]Tyrakis et al., 2017), increased incorporation of ^13C into both malate and fumarate indicate increased overall flux of carbons into these metabolites from the oxaloacetate pool. Taken together, this suggests dysregulated TCA cycle flux in steatotic cells, which may be associated with changes in SDH activity ([137]Figure 5H). The increased incorporation of ^13C into oxaloacetate (using aspartate as a surrogate) demonstrates increased PC activity, driving the conversion of pyruvate to oxaloacetate and suggests possible disruption of associated metabolic cycles, including the malate-aspartate shuttle (MAS) and the purine nucleotide cycle (PNC). Figure 5. [138]Figure 5 [139]Open in a new tab Macrovesicular steatosis in HLCs results in truncation of the TCA cycle, inhibiting conversion of succinate to fumarate. Paradoxically, this is associated with increased accumulation of fumarate. (A) Schematic outlining the TCA cycle and possible ^13C labeled isotopomers from a single cycle. (B) Steady state measurements of TCA cycle-associated metabolites that could be measured by GC-MS. (C) PC-derived αKG is unchanged in the presence of steatosis, and PDH-derived αKG is reduced. (D–E) Labeling patterns of αKG are not due to increased cataplerosis through glutamate, as both PC- and PDH-derived glutamate are reduced in response to steatosis. (F) Furthermore, reduced ^13C labeling of succinate indicates inefficient conversion from αKG. (G–H) Despite reduced labeling of succinate, both fumarate and malate show increased PC-derived label incorporation. (I) Proposed model of preferential anaplerosis of pyruvate into the TCA cycle, with concomitant truncation, preventing conversion of succinate to fumarate. For NMR data, labels are: 1 = ^13C labeling; 0 = no ^13C labeling. All GC-MS data consisted of 10 biological replicates and 2 technical replicates. Isotopomer data were calculated by multiplying MID (multiple ion detection) by normalized total ion count. For NMR data (E) n = 4 biological replicates/group. Data were analyzed by two-way ANOVA with Sidak post-hoc testing and are expressed as mean ± SD. Increased PNC and MAS activity drives NAFLD-associated fumarate accumulation Next, we wanted to determine the source of fumarate accumulation in steatotic HLCs. Transcriptomic analysis identified dysregulation of multiple genes associated with the MAS ([140]Figure 6A) and PNC ([141]Figure 6B), indicating that these pathways may be involved in fumarate accumulation. To investigate this further, HLCs incubated with ^13C[3]-lactate-labeled LPO were co-incubated with either 5-Aminoimidazole-4-carboxamide-1-β-D-ribofuranosyl 5′-monophosphate (AICAR) or O-(Carboxymethyl)hydroxylamine hemihydrochloride (AOA) to inhibit the PNC or MAS, respectively. Addition of AICAR increased only the steady state levels of αKG, and neither AICAR nor AOA had a significant impact on the synthesis of pyruvate, citrate, succinate, or malate ([142]Figure 6C). However, both inhibitors significantly reduced steady state levels of fumarate and aspartate. Although addition of AICAR did not impact on the synthesis of PDH-derived citrate ([143]Figure 6D), it did result in increased accumulation of PDH-derived ^13C in αKG ([144]Figure 6E), suggesting impaired conversion to succinyl-CoA or increased anaplerosis from glutamate. AICAR was also able to partially restore the effects of steatosis on succinate levels ([145]Figure 6F), resulting in small but significant increases in PC- and PDH-derived isotopomers. The most profound effect of AICAR was on fumarate, with each isotopomer reduced to levels below those observed in the control group ([146]Figure 6G), suggesting that in steatotic HLCs, fumarate accumulation is primarily driven by the PNC. Following AOA treatment, we observed a moderate decrease in PC-derived ^13C incorporation into fumarate, demonstrating a small contribution from the MAS. AICAR had a moderate impact on incorporation of PC-derived incorporation of ^13C into malate, whereas AOA impacted on both PC- and PDH-derived ^13C incorporation, suggesting limited contribution of the MAS and PNC to malate accumulation ([147]Figure 6H). Both AICAR and AOA reduced incorporation of ^13C into aspartate ([148]Figure 6I). These data show that the metabolic consequences of LPO treatment, likely driven by increased TCA cycle anaplerosis through OAA, include an increase in both MAS and PNC activity ([149]Figure 6J). Figure 6. [150]Figure 6 [151]Open in a new tab In the presence of TCA cycle truncation, the purine nucleotide cycle and malate-aspartate shuttle fuel fumarate accumulation (A and B) Transcriptomic analysis showed increased expression of malate-aspartate shuttle (MAS) and purine nucleotide cycle (PNC) transcripts, indicating perturbed activity. (C) Inhibition of the PNC and MAS reversed steatosis-induced accumulation of fumarate and aspartate but not malate. (D–G) (D) PNC inhibition moderately impacted PC-derived citrate ^13C labeling but not (E) αKG or (F) succinate. (G) In contrast, PNC inhibition profoundly reduced incorporation of ^13C into fumarate. MAS inhibition also limited PC-derived fumarate accumulation, but to a lesser extent. (H) PNC and MAS expression also inhibited generation of PC-derived malate but to a much smaller extent than fumarate. (I) PNC and MAS inhibition also resulted in a significant reduction in the generation of PC-derived aspartate. (J) Schematic outlining of the proposed pathways by which fumarate accumulation occurs in response to steatosis. (K) Inhibition of fumarate through the PNC reduced LPO-induced macrovesicular steatosis, which was partially restored through addition of exogenous monomethyl fumarate. For NMR data, labels are: 1 = ^13C labeling; 0 = no ^13C labeling. AICAR = PNC inhibitor; AOA = MAS inhibitor. All control and LPO group GC-MS data consisted of 10 biological replicates and 2 technical replicates, as shown in [152]Figure 3. GC-MS LPO + AICAR and LPO + AOA groups consisted of six biological replicates/group. Isotopomer data were calculated by multiplying MID by normalized total ion count. For lipid droplet analysis (K), the LPO group is as shown in [153]Figure 1. For the LPO + AICAR and the LPO + AICAR + Fumarate group, n = 31 and 28 biological replicates/group, respectively. Data were analyzed by two-way ANOVA with Sidak post-hoc testing and are expressed as mean ± SD. Inhibition of fumarate accumulation inhibits lipid droplet hypertrophy in LPO-treated HLCs To determine whether manipulation of fumarate levels affected the development of macrovesicular steatosis, we performed HCA microscopy on HLCs. Cells were treated with LPO in the presence or absence of AICAR. Cells that received AICAR did not develop macrovesicular steatosis ([154]Figure 6K). These results were corroborated by the addition of exogenous fumarate to steatotic HLCs treated with AICAR, which resulted in the development of larger lipid droplets. Fumarate accumulation is not associated with widescale alterations of 5hmC in protein-coding regions in steatotic HLCs One potential consequence of fumarate accumulation in steatotic HLCs is inhibition of αKG-dependent dioxygenase enzyme activity. Given our previous data showing altered 5hmC in a mouse model of NAFLD and the potential importance of the TET enzymes in hepatocellular carcinoma ([155]Lyall et al., 2020; [156]Thomson et al., 2016), we used a DIP-seq approach to measure changes in 5hmC across the genome. Within each group, samples were highly correlated, suggesting that the 5hmC patterns are stable between replicate samples ([157]Figure 7A). We identified 3,294 differentially hydroxymethylated regions (DHRs) between LPO exposed and control cells (>2-fold change, see methods), with the majority of changes located in intragenic or intronic regions ([158]Figure 7B). 5hmC in gene body regions may functionally relate to mRNA transcription ([159]Thomson et al., 2016), and we therefore generated heatmaps representing mean changes across the gene body, with cluster 1 showing subtle increases in 5hmC within the TSS region ([160]Figure 7C). Analysis of genic 5hmC patterns revealed that a number of promoter regions displayed changes in 5hmC levels following LPO exposure ([161]Figure 7D). Integration of DHR and RNA-seq data for promoter regions upstream of the transcriptional start site revealed 12 promoter regions with differential enrichment of 5hmC (>2-fold) and where mRNA expression was altered >0.5-fold. Linear regression of these regions identified a moderate but significant negative relationship between 5hmC enrichment and mRNA expression ([162]Figure S4; [163]Table S7). Of these genes 7 (EPHX3, ERO1B, CSGALNACT1, DOC2A, COL6A1, CASP1 and TMEM88) were previously shown to be dysregulated in the pathogenesis of NAFLD or progression to cirrhosis ([164]Atanasovska et al., 2017; [165]Cazanave et al., 2017; [166]Parafati et al., 2018; [167]Revill et al., 2013; [168]Wilson and Kumar, 2018). Conversely, to our knowledge, there are no described associations between NAFLD and the remaining genes in this list (DDN, OGFRL1, OCIAD2, ABAT, CCDC150). Our findings may indicate a role for 5hmC in the regulation of these genes in NAFLD pathogenesis. Figure 7. [169]Figure 7 [170]Open in a new tab Macrovesicular steatosis in HLCs does not correlate with wide-scale changes in 5hmC enrichment (A) Correlation heatmap of control versus LPO groups, following hmeDIP sequencing, clustered by Euclidean distance. (B) Proportion of DHRs associated with different regions of the genome, showing the majority located in intragenic and intronic regions. (C) Heatmap of DHRs, with an FDR of <0.05 and k-means clustering. (D) Sliding window analysis of the transcriptional start site (TSS), gene body, and transcriptional end site (TES) shows minimal changes between control and LPO groups. For both control and LPO, n = 3 biological replicates/group. Discussion NAFLD is a challenging disease to study in humans. This is due to the difficulty of obtaining tissues that are not compromised by confounding conditions, making high-resolution analysis of transcriptional and metabolic rewiring challenging. In this study we sought to explore the utility of a human-HLC-based model for the study of transcriptomic, epigenomic, and metabolomic dysfunction associated with hepatic steatosis. Using a combination of high content microscopy and gene expression analyses, we showed that treatment of HLCs with LPO induces intracellular lipid accumulation, recapitulating the macrovesicular steatosis observed in NAFLD ([171]Figure 1G). This was in parallel with changes in the expression of genes associated with lipid droplet size ([172]Figures 1D–1F) and is in agreement with findings from human studies ([173]Fujii et al., 2009; [174]Sahini and Borlak, 2016). The disruption of the expression of genes in multiple pathways related to metabolism confirms the findings of a number of previous studies in humans and mice ([175]Collison et al., 2009; [176]Hardwick et al., 2013; [177]Kolwankar et al., 2007; [178]Liu et al., 2016; [179]Nikolaou et al., 2019; [180]Schiöth et al., 2016; [181]Suppli et al., 2019; [182]Yamaguchi and Murata, 2013). Consistent with the altered expression of multiple subunits of respiratory complexes I, IV, and ATPase, respirometry-based analyses revealed decreased maximal mitochondrial respiration and spare capacity ([183]Figures 3E and 3F). Previous data suggest that NALFD associates with inhibition of respiratory activity in complexes I–III in hepatocytes ([184]Koliaki et al., 2015). Our transcriptomic and metabolomic data suggest that this inhibition most likely occurs in complexes I and II. In terms of mechanisms, octanoic acid inhibits complex I–III in rat liver ([185]Scaini et al., 2012) and excess fumarate may also inhibit complex II activity, although manipulation of fumarate did not alter activity of this complex, suggesting that this is not the case. We then questioned whether changes in maximal respiration arose from increased biogenesis or from changes in mitochondrial integrity. Our findings suggest that mitochondrial biogenesis was not affected by LPO treatment ([186]Figure 3), because maximal respiration was the only parameter of respiration that changed, and the ratio of nuclear to mitochondrial DNA, as well as HCA microscopy analysis of mitochondrial mass were not altered. In addition, the finding that CS activity was not adversely affected suggests that mitochondrial integrity was not compromised, although this requires further study. These findings are consistent with human studies of NAFLD that show no changes in mitochondrial mass in steatotic liver tissue ([187]Koliaki et al., 2015). In contrast, liver tissue from patients with NASH contains greater numbers of mitochondria and a decrease in maximal respiration, indicating that a transition in mitochondrial function occurs between these two disease states ([188]Koliaki et al., 2015). The reasons for this are unclear; differences in maximal respiration between NAFLD and NASH could reflect impairment of the ETC and/or changes in baseline mitochondrial respiration. In NASH, the decrease in maximal respiration may be a consequence of increased mitochondrial fragmentation and concurrent decreases in mitophagy ([189]Knott et al., 2008; [190]Koliaki et al., 2015). The differences between studies may also reflect the isolation of mitochondria from the multiple cell types present in whole tissues, which may be altered as a consequence of the disease process. We suggest that the model of hepatic steatosis presented here may mirror the transition between NAFLD and NASH, with respiratory dysfunction in the absence of mitochondrial biogenesis. Future experiments, with longer exposure to LPO and/or the incorporation of non-parenchymal cells into the model, may shed light on this. In agreement with studies in mouse models, and indirect studies in humans, NAFLD was associated with increased hepatic TCA cycle activity ([191]Figure 5B) ([192]Satapati et al., 2012; [193]Sunny et al., 2011). We identified substantial disruption of TCA cycle and oxidative phosphorylation pathways indicating compromised energy metabolism in steatotic HLCs. As TCA cycle flux generates NADH and FADH[2], which transfer electrons to the electron transport chain (ETC), alterations in TCA cycling may have a direct impact on respiration and oxidative stress. In order to delineate detailed information on TCA cycle flux dynamics we utilized stable isotopic tracing with ^13C[3]-lactate, which provides greater metabolite labeling than other substrates ([194]Hui et al., 2017). In the liver, flux of substrates into the TCA cycle predominantly occurs via PC ([195]Lardy et al., 1965), and there is evidence suggesting that PC activity increases in NAFLD ([196]Sunny et al., 2011) and contributes to increased rates of gluconeogenesis through conversion of oxaloacetate to phosphoenolpyruvate ([197]Satapati et al., 2012). Our data support the assertion that PC activity increases in steatotic HLCs and that pyruvate is predominantly utilized to sustain TCA cycle activity. Although we observed sustained synthesis of gluconeogenesis-associated metabolites in steatotic HLCs, there was no increase in gluconeogenesis detected ([198]Figure 4). Increased incorporation of only two ^13C atoms into each of these gluconeogenesis-associated metabolites may suggest that pyruvate-derived acetyl-CoA is being oxidized to OAA via citrate as a means to enter the gluconeogenic pathway. Although increased flux of substrates into the TCA cycle leads to higher levels of TCA cycle activity ([199]Sunny et al., 2011), our data suggest that changes in the TCA cycle are more nuanced than was previously suggested, with enhanced anaplerosis from pyruvate/lactate into the TCA cycle and inhibition of the conversion of succinate to fumarate in steatotic cells. This raised the question of what leads to the diminished incorporation of ^13C into succinate. The decreased cataplerosis through glutamate suggests that this could be due to perturbation of TCA cycle equilibrium, leading to altered interconversion of succinyl-CoA and succinate, which may, in turn, lead to the generation of an increased pool of succinyl-CoA. This is supported by the observation of decreases in PC-derived glutamate and succinate, despite no changes in PC-derived αKG. Because interconversion between succinyl-CoA and succinate is in near equilibrium and readily reversible ([200]Lynn and Guynn, 1978), we propose that this effect could be a result of an alteration in the equilibrium of this reaction, with enhanced succinyl-CoA synthesis, which could be elicited through branched-chain amino acid catabolism. Despite truncation of the TCA cycle, steatotic HLCs rewire their metabolic circuitry to compensate for this, generating increased levels of fumarate through the PNC, and to a lesser extent the MAS. A similar metabolic bypass has been reported in cardiac ischemia ([201]Chouchani et al., 2014), as well as in human and mouse in vitro models of tumorigenesis ([202]Tyrakis et al., 2017). In ischemic reperfusion injury, increased PNC activity results in fumarate overflow, driving reversal of succinate dehydrogenase activity and accumulation of succinate ([203]Chouchani et al., 2014). We did not observe the same phenomenon here, but rather found evidence for inhibition of succinate dehydrogenase activity in steatotic HLCs. Although fumarate can be generated from reversal of fumarate hydratase (FH) activity ([204]Chouchani et al., 2014), we were unable to directly manipulate the activity of this enzyme in order to assess its contribution to the fumarate pool in steatotic HLCs. However, as FH operates at equilibrium ([205]Ajalla Aleixo et al., 2019) it is possible that PNC-fueled fumarate accumulation prevents reverse catalytic activity. Defects in complex I can lead to a reduction in the levels of NAD^+ ([206]Porcelli et al., 2010) and regeneration of fumarate via the PNC may be one mechanism of increasing the NADH pool and maintaining the hydrogen ion gradient of the ETC. Our data suggest that fumarate is important for the development of macrovesicular steatosis; inhibition of the PNC prevented development of macrovesicular lipid droplets, and this was partially restored following the addition of exogenous fumarate ([207]Figure 6K). This correlates with findings in oligodendrocytes and CD8+ T cells, in which exposure to exogenous fumarate resulted in perturbed lipid metabolism ([208]Bhargava et al., 2019; [209]Huang et al., 2015), although the mechanism by which this occurs is unknown. Fumarate accumulation can impact the activity of the αKG-dependent dioxygenase enzyme family, including the ten-eleven translocation (TET) enzymes, which are involved in the conversion of 5-methylcytosine (5mC) into 5-hydroxymethylcytosine (5hmC) and subsequent steps in the DNA demethylation pathway ([210]Tahiliani et al., 2009). We have previously shown that the TET enzymes play an important role in the development of hepatocellular carcinoma and that 5hmC enrichment is reversibly altered at specific loci in response to hepatic fat accumulation ([211]Lyall et al., 2020; [212]Thomson et al., 2016). 5hmC is a stable epigenetic modification and studies suggest that 5hmC enrichment at transcriptional start sites (TSSs) associates with transcriptional repression ([213]Wu et al., 2011). We did not observe wide-scale changes in 5hmC enrichment, and 5hmC changes at specific loci did not, in general, correlate with transcriptional changes in steatotic cells, except at a limited number of promoter regions. This human-cell-based model displays many of the transcriptional and metabolic characteristics associated with human hepatic steatosis, and we propose it as a useful addition to the field. Although there are a number of current in vitro approaches for the study of steatosis, these often use HCC cell lines, which have perturbed metabolic function arising from the process of malignant transformation and immortalization ([214]Huang et al., 2013; [215]Pavlova and Thompson, 2016). This may limit their utility for studying steatosis-associated metabolic changes in non-malignant cells. Furthermore, culturing of HCC cell lines is associated with loss of 5hmC, making them unsuitable for studying the interplay between mitochondrial metabolism and TET-mediated modulation of the epigenome ([216]Chen et al., 2013; [217]Hernandez-Vargas et al., 2010; [218]Liu et al., 2013; [219]Thomson et al., 2016; [220]Xu et al., 2011). As an alternative, primary mouse and human hepatocytes can be used to model hepatic steatosis, but these cells have a limited lifespan, cannot be passaged or expanded in vitro, and human cells are often isolated from fundamentally abnormal, e.g. fatty livers ([221]Hirsova et al., 2016; [222]Lyall et al., 2018). The process of cell sorting may also confound metabolomic, transcriptomic, and epigenomic analyses in cell lines ([223]Llufrio et al., 2018), and the necessity for both tissue disaggregation and cell sorting, which are needed to enable the study downstream effects in individual cell types from whole tissue biopsies, can make the interpretation of results even more challenging ([224]Llufrio et al., 2018). The model system presented here overcomes many of these challenges, as HLCs display a similar metabolic phenotype to human hepatocytes and cell sorting is not required ([225]Wang et al., 2017). Furthermore, these cells retain 5hmC in patterns similar to those observed in humans ([226]Ivanov et al., 2013), making them well suited to studying the interplay between mitochondrial metabolism and the epigenome. The system, allowing the study of transcriptomic, epigenomic, metabolomic, and proteomic effects in a single cell type, can be adapted for the analysis of effects induced by a variety of alternative agents that may induce steatosis. In addition, this model system can be scaled up to model liver injury ([227]Lucendo-Villarin et al., 2020a; [228]Szkolnicka et al., 2014), genetically modified to study gene function ([229]Wang et al., 2019), and moved into 3D culture for longer-term experimentation ([230]Lucendo-Villarin et al., 2020b). These approaches will be valuable to further determine which genes link metabolic disruption and steatosis and for high-throughput drug screening to identify targets for preclinical analysis. In conclusion, we propose that steatotic HLCs represent a new and attractive tool to develop a deeper understanding of steatosis-associated biology and for the development of new therapeutics. Limitations of the study One limitation of this study is the use of AICAR to block PNC activity and limit fumarate accumulation ([231]Chouchani et al., 2014). AICAR is also an AMP-activated protein kinase (AMPK) agonist ([232]Zhao and Saltiel, 2020). AMPK is a regulator of cellular energy homeostasis, and its activation by Metformin, a commonly used treatment for type 2 diabetes, inhibits gluconeogenesis in the liver ([233]Zhou et al., 2001), which may explain why we did not observe increased activity in this pathway. AICAR may have a number of additional effects, notably through downstream pathways, which make the interpretation of its effects problematic, and further studies exploring the effects of overexpression and/or silencing of metabolic enzymes using siRNA on lipid accumulation and the metabolome will be useful to delineate the precise pathways involved. Similarly, there is some evidence that AOA may exert effects upstream of the MAS, leading to inhibition of gluconeogenesis ([234]Subramanian et al., 2007), supporting the need for additional experiments to dissect the precise effect of this inhibitor in the model presented here. Octanoic acid increases oxidative stress through ROS production ([235]Scaini et al., 2012), so future studies should additionally measure oxidative stress and ROS levels in steatotic HLCs. Further studies are required to determine if the mitochondrial dysfunction is due to defects in individual respiratory complexes or changes in a combination of complexes and whether other aspects of mitochondrial respiration may become compromised with longer LPO exposure. Finally, LPO may induce ER stress ([236]Lockman et al., 2016), which may influence mitochondrial oxygen consumption and contribute to the development of NAFLD ([237]Liu and Green, 2019; [238]Wang et al., 2011). Resource availability Lead contact David C Hay: Center for Regenerative Medicine, University of Edinburgh, Institute for Regeneration and Repair, Edinburgh BioQuarter, 5 Little France Crescent, Edinburgh, EH16 4UU. [239]davehay@talktalk.net. Materials availability This study did not generate new unique reagents. Data and code availability The accession number for transcriptomic sequencing data reported in this paper is GEO: [240]GSE138052. The accession number for 5hmC sequencing data reported in this paper is GEO: [241]GSE144955. Original raw GC-MS data have been deposited to Mendeley Data: [242]https://data.mendeley.com/datasets/gd5k5psvf4/1; [243]https://doi.org/10.17632/gd5k5psvf4.1. Methods All methods can be found in the accompanying [244]Transparent methods supplemental file. Acknowledgments