Elisabeth Rohbeck1,2,3
Corinna Niersmann1,2,3
Karl Köhrer4Thorsten Wachtmeister4
Michael Roden5,1,2Jürgen Eckel1,3
Tania Romacho1,6*- 1German Diabetes Center, Institute for Clinical Diabetology, Leibniz Center for Diabetes Research at Heinrich-Heine-University Düsseldorf, Düsseldorf, Germany
- 2German Center for Diabetes Research (DZD), Partner Düsseldorf, München-Neuherberg, Germany
- 3CureDiab Metabolic Research GmbH, Düsseldorf, Germany
- 4Biological and Medical Research Centre (BMFZ), Medical Faculty, Heinrich-Heine-University Düsseldorf, Düsseldorf, Germany
- 5Department of Endocrinology and Diabetology, Medical Faculty and University Hospital Düsseldorf, Heinrich-Heine-University Düsseldorf, Düsseldorf, Germany
- 6Chronic Complications of Diabetes Lab (ChroCoDiL), Department of Nursing Sciences, Physiotherapy and Medicine, Faculty of Health Sciences, University of Almería, Almería, Spain
Introduction: We have previously shown that the novel positive allosteric modulator of the GABAA receptor, HK4, exerts hepatoprotective effects against lipotoxicity-induced apoptosis, DNA damage, inflammation and ER stress in vitro. This might be mediated by downregulated phosphorylation of the transcription factors NF-κB and STAT3. The current study aimed to investigate the effect of HK4 on lipotoxicity-induced hepatocyte injury at the transcriptional level.
Methods: HepG2 cells were treated with palmitate (200 μM) in the presence or absence of HK4 (10 μM) for 7 h. Total RNA was isolated and the expression profiles of mRNAs were assessed. Differentially expressed genes were identified and subjected to the DAVID database and Ingenuity Pathway Analysis software for functional and pathway analysis, all under appropriate statistical testing.
Results: Transcriptomic analysis showed substantial modifications in gene expression in response to palmitate as lipotoxic stimulus with 1,457 differentially expressed genes affecting lipid metabolism, oxidative phosphorylation, apoptosis, oxidative and ER stress among others. HK4 preincubation resulted in the prevention of palmitate-induced dysregulation by restoring initial gene expression pattern of untreated hepatocytes comprising 456 genes. Out of the 456 genes, 342 genes were upregulated and 114 downregulated by HK4. Enriched pathways analysis of those genes by Ingenuity Pathway Analysis, pointed towards oxidative phosphorylation, mitochondrial dysregulation, protein ubiquitination, apoptosis, and cell cycle regulation as affected pathways. These pathways are regulated by the key upstream regulators TP53, KDM5B, DDX5, CAB39 L and SYVN1, which orchestrate the metabolic and oxidative stress responses including modulation of DNA repair and degradation of ER stress-induced misfolded proteins in the presence or absence of HK4.
Discussion: We conclude that HK4 specifically targets mitochondrial respiration, protein ubiquitination, apoptosis and cell cycle. This not only helps to counteract lipotoxic hepatocellular injury through modification of gene expression, but - by targeting transcription factors responsible for DNA repair, cell cycle progression and ER stress - might even prevent lipotoxic mechanisms. These findings suggest that HK4 has a great potential for the treatment of non-alcoholic fatty liver disease (NAFLD).
1 Introduction
No specific pharmacological therapy has been approved to treat non-alcoholic fatty liver disease (NAFLD), a complex disease with many environmental, but also genetic factors contributing to its origin and progression (Byrne and Targher, 2015; Targher et al., 2021). Several previous studies targeting the GABAergic system have shown promising hepatoprotective effects against liver toxicity in vitro and in vivo (Norikura et al., 2007; Shilpa et al., 2012; Hata et al., 2019). It has been reported that GABA improved mitochondrial function in parallel to an attenuation of apoptotic cell death in mice with severe acute liver injury (Hata et al., 2019). Therefore, we aimed to determine if the positive allosteric modulator (PAM) of the GABAA receptor, HK4, could prevent the deleterious effects of lipotoxicity in hepatocytes. In our recent publication (Rohbeck et al., 2022) we were able to confirm the presence of the GABAA receptor in HepG2, as described by Minuk and colleagues (Minuk et al., 1987; Rohbeck et al., 2022). Moreover, we showed that the GABAA receptor PAM, HK4, could reduce palmitate (PA)-induced DNA fragmentation, inflammation, cell death and especially apoptosis, in vitro (Rohbeck et al., 2022). These effects were mediated by activation of the two transcription factors nuclear factor kappa-light-chain-enhancer of activated B Cells (NF-κB) and signal transducers and activators of transcription 3 (STAT3) (Rohbeck et al., 2022). Those transcription factors control the expression of a large number of downstream genes related to cell proliferation, survival, stress responses and inflammation (He and Karin, 2011). Therefore we hypothesize that besides these two transcription factors, HK4 might also counteract gene expression patterns induced by lipotoxic stimuli such as PA exposure.
The development of NAFLD in humans strongly correlates with the dysregulation of transcriptional regulators that affect lipid metabolism (CAR, ChREBP, C/EBPα, FXR, LXR, PPARα/γ/δ, SREBP1c, STAT5), inflammation (c-Jun, C/EBPβ, IRF1/3, NF-κB, RELA, SHP, STAT1/3), metabolic stress (ATF4/6, CYP2E1, eIF2α, IRE1α, Nrf2, Xbp1), and fibrosis (AEBP1, RUNX2, Smad/TGFβ, YAP) (Steensels et al., 2020). Besides the aforementioned transcription factors, there are plenty of genes whose expression correlates positively or negatively with aberrant pathways involved in NAFLD (Greco et al., 2008; Jonas and Schürmann, 2021). Therefore we aimed to identify PA-regulated genes and to further understand the underlying signalling pathways regulated by HK4 as a protection against PA-induced lipotoxicity. In fact, lipid intake enriched in PA, has frequently been associated with obesity, type 2 diabetes mellitus and NAFLD (Hernández et al., 2017). In detail, ingestion of saturated fat rapidly increases energy metabolism, hepatic lipid storage, and insulin resistance (Nowotny et al., 2013). It has been postulated that this metabolic change is accompanied by regulation of hepatic gene expression and signalling that may contribute to the development of NAFLD (Hernández et al., 2017). Thus, several studies have identified candidates that may target transcriptional networks as a response to palmitate (Das et al., 2010; Piccolis et al., 2019; Pérez-Schindler et al., 2022), while other studies identified a typical gene expression signature in NAFLD progression in humans mostly related to lipid metabolism, oxidative phosphorylation, inflammation, endoplasmic reticulum (ER) and oxidative stress and apoptosis pathways (Greco et al., 2008; Ryaboshapkina and Hammar, 2017). RNA-sequencing can further identify and unravel specific gene patterns and involved pathway regulations.
Therefore, the global aim of the present study was to further explore the effects of HK4 on palmitate-induced lipotoxicity in hepatocytes at transcriptional level. In detail, we have examined if the previously described preventive effect of HK4 over palmitate can be translated to the transcriptional level by taking well-known transcriptional features of NAFLD into account.
2 Materials and methods
2.1 Cell culture
HepG2 cells were obtained from Merck KGaA (Darmstadt, Germany; ECACC certified) and cultured according to the corresponding manufacturer’s instructions. During experiments, cells were cultured in serum-free medium. Cells were treated with 200 μM bovine serum albumin-conjugated palmitate (Cayman chemical, Ann Arbor, MI) alone or in combination with 10 μM of HK4 (Taros Chemicals, Dortmund, Germany) up to 30 min prior to palmitate exposure for 7 h.
2.2 RNA isolation
Total RNA of 5 independent experiments, each one comprising 3 treatments (untreated, PA, PA + HK4), from HepG2 cells were isolated with an automated RNA isolation machine Innupure C16 using an RNA Isolation Kit (Analytic Jena, Jena, Germany) following the manufacturer’s instructions.
2.3 3′-RNA-seq analysis
The total RNA samples used for 3′-mRNA Seq analyses were quantified by fluorometric measurement using the Qubit device and a RNA High Sensitivity assay (Thermo Fisher Scientific Inc. Massachusetts, USA) and quality measured by capillary electrophoresis using the Fragment Analyzer and the “Total RNA Standard Sensitivity Assay” (Agilent Technologies Inc. Santa Clara, USA). All samples in this study showed high quality in RNA Quality Numbers (RQN 10). 100 ng total RNA per sample were used for library preparation performed according to the manufacturer’s protocol using the QuantSeq 3′-mRNA-Seq Library Prep Kit FWD (Lexogen®, Vienna, Austria). Bead purified libraries were normalized and finally sequenced on the NextSeq550 system (Illumina Inc. San Diego, USA), using single-end sequencing with a read length of 76bp. The bcl2fastq2 tool was used to convert the bcl files to fastq files as well as for adapter trimming and demultiplexing.
2.4 Bioinformatic analysis
Expression analyses were conducted with CLC Genomics Workbench (version 22.0.1, QIAGEN, Venlo. NL). The reads of all samples were adapter trimmed and quality trimmed (using the default parameters: bases below Q13 were trimmed from the end of the reads, ambiguous nucleotides maximal 2). Mapping was done against the Homo sapiens (GRCh38) genome sequence. After grouping of samples according to their respective experimental condition, multi-group comparisons were made and statistically determined using the in-build algorithm. The resulting p-values were corrected for multiple testing by FDR and Bonferroni-correction. A p-value of ≤0.05 was considered significant.
Ingenuity Pathway Analysis (IPA, Qiagen) was performed to understand interaction networks within differentially expressed genes (DEGs). Absolute values of fold change (|FC|) ≥ 1.5, p-values ≤0.05 and the Fisher’s exact test for p-values have been considered for the identification of canonical pathways. The activation z-score helps to infer the activation states of implicated biological functions, while values ≥2 indicates an activation and a z-score ≤2 points towards the inhibition of potential upstream regulators or canonical pathways. If meaningful, a trend towards activation or inhibition (|z-score| ≥ 1.5) have been considered. Either networks have been designed by selecting corresponding DEGs of mitochondrial respiration or protein ubiquitination pathway or by adapting a proposal of IPA for cell death and cell cycle.
DAVID (Database for Annotation, Visualization and Integrated Discovery) Enrichment Analysis was used for the functional annotation of DEGs (2021 online version, https://david.ncifcrf.gov/) (Da Huang et al., 2009; Sherman et al., 2022). To ascertain the gene-enriched pathways and the potential Gene Ontology (GO) classification, terms comprising biological process, molecular functions, and signalling pathways concerning specifically from Kyoto Encyclopedia of Genes and Genomes (KEGG) database were used. The modified Fisher exact p-value (EASE score) ≤ 0.05 are considered strongly enriched.
2.5 Statistical analyses
Statistical analyses were performed using the GraphPad Prism software (Version 8.0.1, San Diego, USA). Two-way ANOVA (post hoc test: Dunnett’s) were used to assess statistical significant differences. Transcript levels (expressed as reads per kilobase of transcript per million mapped reads (RPKM)) were represented as mean values with standard error of the mean (SEM). p-values ≤0.05 were considered as statistically significant.
3 Results
3.1 Functional annotation of the gene expression profile of hepatocytes under lipotoxicity
Hepatocyte lipid metabolism, oxidative phosphorylation, inflammation, ER and oxidative stress as well as apoptosis pathways are altered during NAFLD progression (Byrne and Targher, 2015; Dewidar et al., 2020). In parallel, modular enrichment analysis of our lipotoxicity model showed that PA altered pathways related to lipid metabolism, cell cycle, NAFLD pathways in general, oxidative phosphorylation, apoptosis, oxidative stress response, peroxisome proliferator-activated receptors (PPAR) signalling, tumor protein p53 (TP53) and unfolded protein response (Table 1). Between 77 and 10 DEG are involved in the aforementioned pathways addressed by modular enrichment analysis with an increasing fold enrichment value.
TABLE 1. Modular enrichment analysis of genes differentially expressed by palmitate exposure.
3.2 Treatment with palmitate and HK4 results in differentially regulated gene profiles
3′-mRNA sequencing and expression analysis identified differentially regulated protein-coding genes among untreated HepG2 cells, in PA-treated cells and PA + HK4-treated cells (Figure 1). In a differential expression analysis, genes have been filtered by a p-value ≤0.05 and |FC| ≥ 1.5. 1,457 genes were differentially expressed when comparing untreated and PA-treated hepatocytes, while only half of the number of genes (723) remained differentially expressed between PA alone compared to PA in combination with HK4. Interestingly, a much lower number of genes (271) are differently expressed between untreated and PA + HK4. From the 784 upregulated genes in untreated cells compared to PA-treated and the 472 upregulated genes in PA + HK4 compared to PA, 342 genes are identic. Thus, these 342 genes are downregulated when treated with PA and reversibly upregulated in the presence of HK4. In contrast, from 673 downregulated genes in untreated cells compared to PA and 251 downregulated genes in PA + HK4 compared to PA, 114 genes are identic. Thus, 114 genes are upregulated when treated with PA and reversibly downregulated in the presence of HK4. We consider those 456 genes (342 upregulated, 114 downregulated in PA + HK4) as key regulated genes by HK4 to restore initial physiological conditions from lipotoxicity.
FIGURE 1. Workflow and overview of the number of genes in untreated, palmitate (PA) or PA + HK4-treated HepG2 cells. After 7 h treatment with 200 μM PA and 10 μM HK4, RNA has been isolated from five independent HepG2 cell lysates per treatment group and 3′-mRNA sequencing has been performed with Illumina NextSeq550 system. Thick arrows indicate up or downregulated genes in the first treatment group of the comparison. Differential expression (DE) of genes have been filtered by p-value ≤0.05 and an absolute value of fold change (|FC|) ≥ 1.5.
3.3 Determining the transcriptome of palmitate-induced lipotoxicity and restoring of initial gene expression pattern as in untreated hepatocytes by HK4
Hierarchical clustering showed distinctive gene expression profiles of untreated hepatocytes, hepatocytes exposed to PA alone or PA combined with HK4 (Figure 2). A heat map with differentially regulated genes after 7 h treatment, showed an opposite expression pattern in the majority of genes in the PA-treated group compared to untreated and PA + HK4 group. Thus, gene expression pattern of PA + HK4-treated cells is more similar to the untreated group than the PA-treated group, except for one PA + HK4 sample. In total, the heat map considered 1,570 transcripts (Supplementary Figure S1).
FIGURE 2. Overview of hierarchical clustering and heat map of transcription profile of HepG2 cells in response to palmitate (PA, pink) and HK4 (blue) or untreated cells (green). Euclidean metric for distance measurement was used. Changes in the abundance of 1,570 individual genes are shown, while each rectangle represents a gene. The intensity of the red and blue colors correlates with the degree of up- and downregulation, respectively. Differentially expressed genes were defined by ANOVA with a p-value ≤0.05 and an absolute fold change ≥1.5.
3.4 HK4-mediated modification in gene expression related to mitochondrial respiration, protein ubiquitination, cell cycle and apoptosis are crucial for protecting cells from palmitate-induced lipotoxicity
Besides the clustering of genes, we used bioinformatic approaches (IPA) to evaluate the overrepresentation of a group of genes mapping to a specific pathway compared to the total reference set of genes. Thus, we obtained an overview of the biological pathways modified by HK4 in lipotoxic conditions (Figure 3). The 456 DEGs with an up- or downregulation in PA and a reversible regulation in the presence of HK4 were analysed.
FIGURE 3. Enriched pathways of 456 differentially expressed genes in untreated or PA + HK4-treated HepG2 cells compared to PA-treated cells. The length of the bars is proportional to the significance of the association between the set of genes and the pathway, expressed by the negative logarithm of the p-value. Only pathways with p ≤ 0.05 (dotted threshold line) are shown. Upregulated pathways in untreated or PA + HK4-treated cells associated with a positive z-score are colored in orange, downregulated pathway with a negative z-score are shown in blue and pathways with no change or an unknown pattern are marked in white or grey. Selected canonical pathways which are subsequently considered in more detail are highlighted in red.
Regarding affected canonical pathways in IPA with a |z-score| ≥ 2, several DEGs are targeting and activating oxidative phosphorylation (z-score = 3, p = 0.0004, 9 genes), cell cycle control (z-score = 2.24, p = 0.006, 5 genes) or on the contrary inhibiting the caspase-independent cell death granzyme A signalling pathway (z-score = −2.24, p = 0.016, 5 genes) in PA + HK4 or untreated hepatocytes compared to PA-treated cells. While several other pathways related to processes such as senescence (z-score = 1.67, p = 0.046, 11 genes), ferroptosis (z-score = 1.13, p = 0.005, 8 genes), apoptosis (z-score = 0.81, p = 0.02, 6 genes), cell cycle regulation (z-score = −1, p = 0.02, 4 genes) and sirtuin signalling pathway (z-score = −0.7, p = 0.007, 13 genes) were significantly addressed without passing threshold |z-score| of 2. Since the underlying network also includes findings without associated directional attributes, several pathways with a considerable number of DEGs show up with unknown, unchanged or controversial pattern. Among those, mitochondrial dysfunction (p = 0.0007, 11 genes) protein ubiquitination (p = 0.002, 14 genes), iron homeostasis (p = 0.007, 8 genes), proline biosynthesis (p = 0.003, 2 genes), glucose degradation (p = 0.043, 2 genes) and granzyme B (p = 0.043, 2 genes) signalling were found.
To continue assessing the influence of HK4 on the gene expression profile, associated genes of some of the aforementioned canonical pathways were identified by IPA. Selected canonical pathways were mitochondrial respiration, protein ubiquitination, apoptosis and cell cycle related pathways (highlighted in red in Figure 3). We have taken a closer look at oxidative phosphorylation and mitochondrial dysfunction (clustered together as mitochondrial respiration), as well as protein ubiquitination, since they are the most relevant pathways according to their p-values and number of included genes. Additionally, we have focused on apoptosis, because we previously demonstrated the protective effect of HK4 on PA-induced apoptosis (Rohbeck et al., 2022). Since cell cycle plays a pivotal role in our lipotoxicity model with the highest number of related DEGs (Table 1) and it is closely connected to apoptosis (Caldez et al., 2020), we also took cell cycle into account. We listed genes related to the selected pathways (mitochondrial respiration, protein degradation, apoptosis, and cell cycle) suggested by IPA and DAVID together with the fold change of their expression in Supplementary Table S1. Further identified pathways were ER stress, inflammation and lipid metabolism which display few related DEGs, such as CREB3L3, CHOP, NFKB2, NFKBIL1, IL17RC, ABHD3, PLA2G2A (Supplementary Table S1). Out of the 456 DEGs, 14 genes (1 down-, 13 upregulated) can strongly be linked to mitochondrial respiration, 22 genes (3 down-, 19 upregulated) to protein ubiquitination pathway and 14 genes both to apoptosis (5 down-, 9 upregulated) and cell cycle (2 down-, 12 upregulated). The genes ANAPC10, BIRC3, CDC26, FZR1, SIAH1 and SKP2 are matched to several pathways of IPA (protein ubiquitination and apoptosis or cell cycle) at the same time, since the pathways strongly interact with each other. Therefore, we only depict them once for apoptosis or cell cycle, when plotting RPKM values of all associated genes in Figures 4A–D. Remarkably, in all signalling pathways considerable fewer genes are downregulated in PA + HK4-treated hepatocytes.
FIGURE 4. Reads per kilobase of transcript per million mapped reads (RPKM) expression values of genes related to mitochondrial respiration (A), protein ubiquitination (B), apoptosis (C) and cell cycle (D). Genes were selected of the 456 DEGs among untreated (white), palmitate (PA, black) and PA + HK4 (black striped)-treated HepG2 cells and sorted by descending expression values of PA. p-values were calculated by two-way ANOVA; ***p ≤ 0.001, **p ≤ 0.01, *p ≤ 0.05 vs. untreated, #p ≤ 0.05 vs. PA + HK4-treated cells. Abbreviations of genes are explained in Supplementary Table S1.
Most of the genes associated with mitochondrial respiration (Figure 4A) are coding components of the respiratory chain complexes like the NADH dehydrogenase (NDUFA4, NDUFB3, NDUFB5), the succinate dehydrogenase (ETFA, SDHB, SDHC), the cytochrome bc1 complex (CYB5A), the cytochrome c oxidase (COX7B, COX17) or the ATP synthase (ATP5F1C, ATP6V0A1). Additionally, the pyrophosphatase-coding gene PPA2 and peroxiredoxin 3-coding gene PRDX3 are listed as associated genes of oxidative phosphorylation and mitochondrial dysfunction. While PPA2 regulates energy metabolism of the cell, PRDX3 exerts an antioxidative function. From all 14 genes, only ATP6V0A1 is downregulated in PA + HK4-treated hepatocytes compared to PA-treated cells alone.
Genes associated with protein ubiquitination (Figure 4B) depict mainly components of the E3 ubiquitin-protein ligase (COP1, MARCHF5, PELI3, SYVN1, UHRF2), the E2 ubiquitin-conjugated enzyme (UBE2T), the E1 ubiquitin-activating enzyme (ATG7), the proteasome (PSMA2, PSMB8) or as heat shock proteins (DNAJC4, HSPA1B, HSPB11, HSPE1). Furthermore, TAP2 is coding a member of ATP transporter and MED20 and PAN2 are involved in the deubiquitination. In total, only three genes (PELI3, ATG7 and DNAJC4) are downregulated in PA + HK4-treated hepatocytes compared to PA alone.
Regarding the genes associated with apoptosis (Figure 4C), 5 genes are downregulated in the presence of HK4 combined with PA compared to PA-treated cells. The most downregulated gene in this pathway, referred to its fold change (Supplementary Table S1), is the anti-apoptotic gene BIRC3, followed by the pro-apoptotic genes HRAS, BAX, TP73 and FASTK. 9 genes (TNFRSF10B, STEAP3, SIAH1, CTSC, BCL2L11, LMNB2, NRAS, DDIAS, DFFB) show an upregulation when treated with PA in combination with HK4.
Regarding genes related to cell cycle (Figure 4D), two genes (CDKN1C, FZR1) are downregulated in untreated or PA + HK4-treated cells compared to PA-treated hepatocytes, while 11 genes are upregulated in the presence of PA + HK4. Among those, there are the cell division cycle–related genes (CDCA7, CDC26 and CDC45), CDK7, CHEK2, MCM2 and ORC5, which are mainly regulating the progression steps of the interphase, while ANAPC10 and ESPL1 play a role in the anaphase of the mitosis. PPM1D is a negative regulator of cellular stress responses.
All genes of mitochondrial respiration, protein degradation, apoptosis and cell cycle have been included in networks designed with IPA with the purpose to understand how groups of data set molecules might interact (Figures 5A–D). B-cell receptor (BCR) complex, cytochrome-c oxidase, Interferon-α, IKK complex, IL 6 and CCND1, as well as the upstream regulators TP53, TP73, NF-κB, STAT5A, MAPK9, RB1 and DDX5 have been listed as additional nodes to connect the genes of the pathways.
FIGURE 5. Gene interaction network map generated with IPA for mitochondrial respiration (A), protein ubiquitination (B), apoptosis (C) and cell cycle (D) pathways. Genes associated with one of the four pathways are represented by nodes with their shape representing the type of molecule/functional class. The causal network results of the 456 DEG identified the master regulator SYVN1 targeting the highest number of molecules in our dataset through the intermediate regulators HTT and TP53 (E). Network includes up and down regulation of downstream dataset molecules. Nodes in red are upregulated in untreated/PA + HK4-treated HepG2 cells, orange nodes are predicted to be activated from Ingenuity knowledge base, green colored nodes show downregulation and blue ones are predicted to be inhibited. Orange lines between the nodes indicate an activating relationship, blue lines an inhibition. Pointed arrowheads indicate that the downstream node is expected to be activated if the upstream node connected to it is activated, while blunt arrowheads indicate that the downstream node is expected to be inhibited if the upstream node that connects to it is activated. Predicted legend is shown on the right side.
3.5 Upstream regulators of PA and HK4-modified genes affect DNA repair and cell cycle
Upstream regulators cover a spectrum of molecule types (e.g., cytokines, kinases, microRNA, receptors and transcription factors) that affect the expression, transcription, or phosphorylation of other (downstream) molecules. To unravel critical mediators of lipotoxicity upstream of the 456 genes modified by PA and HK4, we identified upstream regulators with two different analyses in the IPA software. The first approach of the upstream regulator analysis (URA) determines likely upstream regulators (Table 2) that are connected to dataset genes through a set of direct or indirect relationships. URA revealed 4 potential upstream regulators of the gene expression profile by PA and HK4 with passing |z-score| threshold of 2 (Table 2). Among those potential regulators, 2 are predicted to be activated (CAB39L (z-score = 2), DDX5 (z-score = 2.65)) and 2 to be inhibited (KDM5B (z-score = −2.82), TP53 (z-score = −2,11)) with HK4 pretreatment. The most relevant upstream regulator is the tumor suppressor TP53, affecting 30 targets. The encoded protein responds to diverse ways of cellular stress to regulate expression of target genes, thereby inducing cell cycle arrest, apoptosis, senescence, DNA repair or changes in metabolism (Vogelstein et al., 2000; Marcel et al., 2010). The other transcription regulator, KDM5B, is affecting 8 target molecules from our dataset. Interestingly, KDM5B plays a role in DNA repair and the transcriptional repression of certain tumor suppressor genes (Klein et al., 2014; Li et al., 2014). DDX5, showing the highest activation z-score (2.646), is an established co-activator of TP53 and therefore having a pivotal role in orchestrating the cellular response to DNA damage and repair (Nicol et al., 2013). CAB39L, the last upstream regulator of Table 2, exerts a tumor suppressive effect by inducing apoptosis and cell cycle arrest (Li et al., 2018).
TABLE 2. Potential upstream regulators of palmitate and HK4-induced gene expression profile.
The second approach is the causal network analysis (CNA) connecting upstream regulators to dataset molecules with the advantage of taking paths into account that involve more than one link (i.e., through intermediate regulators). Therefore, the CNA can been seen as a generalization of URA, which are used to generate a more complete picture of possible root causes for the observed expression changes (Krämer et al., 2014). In total, 19 master regulators (data not shown) are passing |z-score| threshold of 2. We plotted causal network analysis results for the “root” regulator SYVN1 which affects the most regulators (16) and target molecules (55) in our dataset (Figure 5E). SYVN1 is affecting many biological processes like response to unfolded protein, protein ubiquitination and intrinsic apoptotic signalling pathway in response to ER stress through the two intermediate regulators HTT and TP53. HTT (huntingtin) is known to be involved in Huntington’s disease signalling, but also required for normal development, including vesicle transport, protein trafficking, and transcriptional regulation (Rodriguez-Lebron et al., 2005). The upstream regulator TP53 not only affects 30 target molecules of our dataset in the URA, but it is also predicted to activate or inhibit 12 other upstream regulators (PRKCE, JNK, MEK, NOTCH1, NF-κB (complex), PRKAA1, GLI1, TSC2, TP73, RB1, NPM1, MAPK9). Some important activated upstream regulators by TP53, namely NF-κB, MAPK9, TP73 and RB1 are shown to interact together with STAT5A and DDX5 with selected downstream genes of mitochondrial respiration, protein ubiquitination, apoptosis and cell cycle (Figure 5E).
4 Discussion
PA is one of the most abundant saturated fatty acids in diet (Murru et al., 2022). Analysing the transcriptome in hepatocytes treated with PA to induce lipotoxicity, shows that short term palmitate exposure induces a specific gene expression profile. Although the effect of lipotoxicity on the induction of transcriptome responses has been analysed in hepatocytes, the impact of the novel GABAA receptor PAM, HK4, on lipotoxicity has not been previously studied. Therefore, the purpose of this study was to describe gene expression changes occurring under conditions resembling lipotoxicity and to analyse the potential protection by HK4. We found that mitochondrial respiration, protein ubiquitination, apoptosis and cell cycle pathways are major targets of HK4, mainly addressed by the upstream regulators SYVN1 and TP53.
4.1 Differential gene expression profile of hepatocytes under PA-induced lipotoxicity and its reversion by HK4
Under conditions mimicking hepatic lipotoxicity with PA treatment, 1,457 genes were differentially expressed. In line, in HepG2 cells it has been reported that 776 genes were affected by 1 mM PA after 6 h exposure (Das et al., 2010), while others detected changes in only 11 genes after 24 h exposure of 50 μM PA (Vock et al., 2007). These differences in differential gene expression induced by palmitate are likely due concentration used, treatment duration and/or hepatocyte cell type.
A total of 456 genes were differentially expressed in PA treatment compared to untreated cells and reverted by HK4 pre-treatment. Putting these DEGs into pathophysiological context, the canonical pathways in IPA pointed towards most relevant DEGs related to mitochondrial respiration, protein ubiquitination, apoptosis and cell cycle pathways. These pathways are regulated in hepatocytes undergoing an adaptive response similarly to other studies resembling pathophysiology in NAFLD (Das et al., 2010; Mota et al., 2016; Piccolis et al., 2019).
Both oxidative phosphorylation and mitochondrial dysfunction, are listed as the top two canonical pathways of our dataset according to their p-value comparing PA alone and in combination with HK4.
During NAFLD electron transport chain complex expression and respiratory control ratio are decreased, while β-oxidation and TCA cycle activity are increased (Koliaki and Roden, 2016). This is in line with the assumption that hepatic mitochondria might upregulate their oxidative capacity at the expense of decreased coupling efficiency when transiently adapting to lipid overload (Koliaki et al., 2015; Jelenik et al., 2017). Later on, the loss of mitochondrial adaptation will favor lipid deposition and insulin resistance and in turn accelerate oxidative stress, resulting finally in non-alcoholic steatohepatitis (NASH) with impaired mitochondrial biogenesis (Fromenty and Roden, 2022). We observed that mitochondrial respiration-related genes as nicastrin (NCSTN), NADH dehydrogenase (NDUFA4) and succinate dehydrogenase (SDHB) were downregulated by acute palmitate exposure, exhausted from mitochondrial adaption. In mice fed with a Western style diet, mimicking NAFLD, a reduced succinate-activated respiration has been described, due to reduced SDHB gene expression (Staňková et al., 2021). Reduced half-life of oxidative phosphorylation subunits contributed to mitochondrial impairment in mice with NAFLD (Lee et al., 2018). It has also been described, that impaired activity of the NADH dehydrogenase, due to pathogenic mtDNA mutation, occurs in people with type 2 diabetes mellitus. (Sharma et al., 2009). However description about the gene expression of specific NADH dehydrogenase subunits (NDUFA4, NDUFB3 or NDUFB5) in the context of NAFLD is lacking. In line with our sequencing data, NCSTN protein is overexpressed in liver in HCC leading to an enhanced cell growth and migration through Notch1 and Akt signalling pathways, however little is known about its protective effect in NAFLD (Li et al., 2020). Therefore, it can be postulated that an upregulation of the novel identified genes NCSTN and NDUFA4, as well as SDHB might be a protective mechanism of HK4 to restore mitochondrial capacity and boost energy supply.
The protein ubiquitination pathway is affected in our lipotoxicity model, since dysregulation of ER-protein folding and ER stress response is one of the predominant hallmarks of NAFLD progression. Specifically impaired proliferation and apoptosis in hepatocellular injury is correlated with loss of the proteasome or the inhibition of the ubiquitin-proteasome pathway (Wójcik, 2002). Thus, in the liver of people with NAFLD, inactivation of components of the ubiquitin-proteasome pathway, promotes apoptotic cell death (Joshi-Barve et al., 2003). It has also been suggested, that ATP deficiency due to reduced mitochondrial respiration contributed to inhibition of ubiquitin-proteasome and proper protein degradation and therefore activate mitophagy (Lee et al., 2018). This explains the lower expression of the protein ubiquitination-regulated genes, like HSPA1B and TAP after PA exposure and the protective effect of HK4 by elevating their expression. In line, expression of HSP70 members (like HSPA1B) has been proposed to take over an anti-inflammatory role in tissue-resident macrophages (like Kupffer cells) upon metabolic challenge (Brykczynska et al., 2020). There is no literature on the gene expression of TAP2 in NAFLD. However, TAP2 downregulation has been confirmed in hepatocyte-like cells from obese people, leading to abnormally metabolism and liver regeneration (Li, Y. et al., 2021).
Ji and colleagues have shown that palmitic acid could induce mitochondrial-mediated apoptosis in HepG2 by upregulating Bax and downregulating Bcl-2 (Ji et al., 2005; Panasiuk et al., 2006). Our data are in line to what was previously described, but also show that HK4 can prevent the upregulation of Bax and downregulation of Bcl-2 induced by PA. Additionally, the lower expression of BIRC3 in HK4-treated cells might also play a protective role in our cell model, since inhibition of BIRC3 reduced hepatocellular carcinoma and progression of metastases (Fu et al., 2019; Frazzi, 2021). These findings underline the results of our recent publication where HK4 exerted an anti-apoptotic effect in HepG2 cells after 24 h of PA exposure (Rohbeck et al., 2022).
Cell division is essential for organismal growth and tissue homeostasis, especially when tissue is damaged in NAFLD (Caldez et al., 2020). Aberrations in cell cycle proteins or their regulators might even lead to HCC (Bisteau et al., 2014). In line, gene expression of CDKN1C in our dataset reflects protein downregulation of CDKN1C in HCC and cirrhotic liver samples (Fornari et al., 2008). Furthermore, decreased levels of Mcm2 as a sign of dysregulated cell cycle progression has been shown in hepatocytes in NAFLD (Dabravolski et al., 2021). Inactivity or deletion of CDC26 might result in a reduced activity of anaphase-promoting complex-cyclosome, thus impairing a variety of cellular processes such as cell division, differentiation, genome stability, energy metabolism, cell death and autophagy (Zhou et al., 2016). Although less is known about the specific expression of cell cycle-related genes in patients with NAFLD, impairment (by gene mutation) and absence of cell cycle-related genes in general repress cell proliferation, which inhibits liver regeneration after acute injury or chronic damage (Caldez et al., 2020). Therefore, we can postulate that cell cycle control also plays a pivotal role in the protective effect of HK4.
4.2 The upstream regulators SYVN1 and TP53 control and integrate HK4 protective effects against lipotoxicity
The causal network analysis suggests that the E3 ubiquitin-protein ligase synoviolin (SYVN1, aka HDR1) is an upstream regulator. The protein encoded by this gene is involved in the ubiquitin-dependent degradation of misfolded proteins which are accumulated during ER stress and proposed as a liver metabolic regulator and a potential drug target for fatty liver disease and progressed hepatocellular carcinoma (Kikkert et al., 2004; Ji et al., 2021). In cirrhotic liver tissues upregulated SYVN1 has been implicated in the downregulation of NRF2, a transcription factor that combats oxidative stress (Wu et al., 2014; Bathish et al., 2022). Overexpression of SYVN1 has been proposed to ameliorate hepatic steatosis and enhanced insulin sensitivity in db/db mice (Li, K. et al., 2021). Our data suggest that upregulation of SYVN1 by HK4 pretreatment might be one mediator of the protective effect of HK4. On the contrary, other works points towards SYVN1 deletion specifically in the liver as protection against HFD-induced obesity and liver steatosis and insulin resistance in mice (Wei et al., 2018). Genome-wide mRNA sequencing revealed that SYVN1 deficiency reprograms liver metabolic gene expression profiles, including suppressing genes involved in glycogenesis and lipogenesis and upregulating genes involved in glycolysis and fatty acid oxidation (Wei et al., 2018).
Our results point towards a negative regulation of SYVN1 over HTT (Yang et al., 2007) and TP53 (Yamasaki et al., 2007). Mutant htt is predicted to inhibit the 26S ubiquitin proteasome system, composed of the 20S catalytic core complex and the 19S regulatory units, resulting in protein aggregation and ER stress (Valera et al., 2005; Qu et al., 2021). Thus an upregulation of the components of the proteasome complex PSMA2 and PSMB8, in the presence of HK4 might lead to proper protein degradation and avoid accumulation of misfolded proteins, which could otherwise cause hepatic steatosis (Valera et al., 2005; Rutkowski et al., 2008).
Regarding TP53, it has been described that SYVN1 sequestrates and metabolizes TP53 in the cytoplasm and negatively regulates its cellular level and biological functions, including transcription, cell cycle regulation, senescence, DNA repair and apoptosis (Yamasaki et al., 2007). As a double-edged sword TP53 takes over a dual function in both the aggravation and amelioration of NAFLD (Yan et al., 2018). Evidence also suggests that induction of TP53 improves pathophysiological conditions associated to NAFLD (Gross et al., 2017; Humpton et al., 2022). Thus, TP53 is a relevant upstream regulator in our setting orchestrating cellular signalling, like ER stress response by impaired mitochondrial function through counteraction of RB1 and NF-κB (Johnson et al., 2011). Cells lacking Rb1 exhibit defective mitochondria and decreased oxygen consumption (Váraljai et al., 2015). Most likely, also the upstream regulator DDX5, interacting with TP53, directly regulates transcriptional activity of mitochondrial-related genes by binding at or near their promoter (Xing et al., 2020; Xu et al., 2022). Similar to our results, DDX5 expression is reported to be downregulated in palmitate-stimulated hepatocytes but also in patients with NASH (Zhang et al., 2022). Additionally, modulation of TP53-dependent pathways during prolonged metabolic stress has been also described as linked to the protein ubiquitination pathway (French and Bardag-Gorce, 2005; Lee et al., 2012). In this regard, we can speculate that through inhibition of NF-κB activation (Wadgaonkar et al., 1999), TP53 controls ubiquitination via the upregulation of PSMA2, PSMB8 and TAP2 in HK4-treated cells. Furthermore, we can postulate that TP53 might also control apoptosis either via interaction with TP73 and JNK pathway (MAPK9) or directly interfering with BCL2L11 (Brockhaus and Brüne, 1999; Dhanasekaran and Reddy, 2008; Kunst et al., 2016). TP53, one of the key upstream regulators, is a well-known member of the p53 tumor suppressor family, that can inhibit cell cycle-related genes and arrest cell cycle (Kunst et al., 2016; Engeland, 2022). This might explain, why all cell cycle regulating genes except for its negative regulating inhibitor of cyclin-dependent kinases, CDKN1C, and FZR1 are downregulated in lipotoxic conditions and again upregulated in the presence of HK4 in our setting (Fornari et al., 2008).
We have previously shown that the protective effects of HK4 can be mediated by downregulated phosphorylation of NF-κB and STAT3 (Rohbeck et al., 2022). Substantial evidence supports that activation of the transcription factor NF-κB and downstream inflammatory signalling pathways are involved in hepatic insulin resistance (Tilg et al., 2017). In line, upregulated NF-κB gene expression in PA-treated HepG2 cells compared to untreated (p = 0.004) or PA and HK4-co-treated cells (p = 0.089) fits with protein expression analysis of our previous study (Rohbeck et al., 2022). Also in line with our previous study, the current data indicate an involvement of STAT signaling, since STAT5A is differentially expressed in our dataset. STAT5A, likewise our recently recognized mediator of HK4 protective effect, STAT3, is known to control cell survival, differentiation, proliferation, and metabolism in response to extracellular stimuli (Wingelhofer et al., 2018). A very recent study suggests that non-canonical NF-κB activation can attenuate the hepatoprotective JAK2/STAT5 signalling (Vesting et al., 2022). We conclude from our study that mitochondrial respiration, protein ubiquitination, apoptosis and cell cycle are major targets of HK4 to minimize hepatocellular injury under a lipotoxic stimulus as PA. Our findings strongly suggest that transcription factors responsible for DNA repair and ER stress such as TP53 play a central regulatory role in hepatocyte lipotoxicity with this novel drug. Targeting those transcription factors and interfere in gene expression pattern holds great therapeutic potential for HK4 in the treatment of NAFLD.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, PRJNA902305.
Author contributions
ER, CN, JE, MR, and TR contributed to conception and design of the study. ER performed experiments. KK and TW performed sequencing. ER, KK, and TW performed data analysis. ER, CN, JE, TR interpreted data. ER wrote the first draft of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.
Funding
TR and JE are supported by KomIT (Center of Competence for Innovative Diabetes Therapy) funded by the European Funds for Regional Development EFRE-0400192. TR is the recipient of a Next-Generation EU Maria Zambrano fellowship. This work was supported by grants from the German Federal Ministry of Health (BMG) and the Ministry of Culture and Science of the State North Rhine-Westphalia (MKW NRW) to the German Diabetes Center (DDZ), the German Federal Ministry of Education and Research (BMBF) to German Center for Diabetes Research (DZD e.V.).
Acknowledgments
We are grateful to Alexandra Stein for excellent technical support and thank Algiax Pharmaceuticals GmbH for providing HK4. Computational infrastructure and support were provided by the Center for Information and Media Technology at Heinrich-Heine-University Düsseldorf.
Conflict of interest
Authors ER, CN, and JE are employed by CMR CureDiab Metabolic Research GmbH.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2023.1106075/full#supplementary-material
References
Bathish B., Robertson H., Dillon J. F., Dinkova-Kostova A. T., Hayes J. D. (2022). Nonalcoholic steatohepatitis and mechanisms by which it is ameliorated by activation of the CNC-bZIP transcription factor Nrf2. Free Radic. Biol. Med. 188, 221–261. doi:10.1016/j.freeradbiomed.2022.06.226
Bisteau X., Caldez M. J., Kaldis P. (2014). The complex relationship between liver cancer and the cell cycle: A story of multiple regulations. Cancers 6 (1), 79–111. doi:10.3390/cancers6010079
Brockhaus F., Brüne B. (1999). p53 accumulation in apoptotic macrophages is an energy demanding process that precedes cytochrome c release in response to nitric oxide. Oncogene 18 (47), 6403–6410. doi:10.1038/sj.onc.1203058
Brykczynska U., Geigges M., Wiedemann S. J., Dror E., Böni-Schnetzler M., Hess C., et al. (2020). Distinct transcriptional responses across tissue-resident macrophages to short-term and long-term metabolic challenge. Cell. Rep. 30 (5), 1627–1643. doi:10.1016/j.celrep.2020.01.005
Byrne C. D., Targher G. (2015). Nafld: A multisystem disease. J. hepatology 62 (1), S47–S64. doi:10.1016/j.jhep.2014.12.012
Caldez M. J., Bjorklund M., Kaldis P. (2020). Cell cycle regulation in NAFLD: When imbalanced metabolism limits cell division. Hepatol. Int. 14 (4), 463–474. doi:10.1007/s12072-020-10066-6
Da Huang W., Sherman B. T., Lempicki R. A. (2009). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 4 (1), 44–57. doi:10.1038/nprot.2008.211
Dabravolski S. A., Bezsonov E. E., Orekhov A. N. (2021). The role of mitochondria dysfunction and hepatic senescence in NAFLD development and progression. Biomed. Pharmacother. = Biomedecine Pharmacother. 142, 112041. doi:10.1016/j.biopha.2021.112041
Das S. K., Mondal A. K., Elbein S. C. (2010). Distinct gene expression profiles characterize cellular responses to palmitate and oleate. J. lipid Res. 51 (8), 2121–2131. doi:10.1194/jlr.M004275
Dewidar B., Kahl S., Pafili K., Roden M. (2020). Metabolic liver disease in diabetes - from mechanisms to clinical trials. Metabolism 111S, 154299. clinical and experimental. doi:10.1016/j.metabol.2020.154299
Dhanasekaran D. N., Reddy E. P. (2008). JNK signaling in apoptosis. Oncogene 27 (48), 6245–6251. doi:10.1038/onc.2008.301
Engeland K. (2022). Cell cycle regulation: p53-p21-RB signaling. Cell. death Differ. 29 (5), 946–960. doi:10.1038/s41418-022-00988-z
Fornari F., Gramantieri L., Ferracin M., Veronese A., Sabbioni S., Calin G. A., et al. (2008). MiR-221 controls CDKN1C/p57 and CDKN1B/p27 expression in human hepatocellular carcinoma. Oncogene 27 (43), 5651–5661. doi:10.1038/onc.2008.178
Frazzi R. (2021). BIRC3 and BIRC5: Multi-faceted inhibitors in cancer. Cell. & Biosci. 11 (1), 8. doi:10.1186/s13578-020-00521-0
French S., Bardag-Gorce F. (2005). “Ubiquitin-proteasome pathway in the pathogenesis of liver disease,” in Signaling pathways in liver diseases: With 15 tables. Editor J.-F. Dufour (Berlin, Heidelberg: Springer), 377–389.
Fromenty B., Roden M. (2022). Mitochondrial alterations in fatty liver diseases. J. hepatology 78, 415–429. doi:10.1016/j.jhep.2022.09.020
Fu P.-Y., Hu B., Ma X.-L., Yang Z.-F., Yu M.-C., Sun H.-X., et al. (2019). New insight into BIRC3: A novel prognostic indicator and a potential therapeutic target for liver cancer. J. Cell. Biochem. 120 (4), 6035–6045. doi:10.1002/jcb.27890
Greco D., Kotronen A., Westerbacka J., Puig O., Arkkila P., Kiviluoto T., et al. (2008). Gene expression in human NAFLD. Am. J. physiology. Gastrointest. liver physiology 294 (5), G1281–G1287. doi:10.1152/ajpgi.00074.2008
Gross B., Pawlak M., Lefebvre P., Staels B. (2017). PPARs in obesity-induced T2DM, dyslipidaemia and NAFLD. Nat. Rev. Endocrinol. 13 (1), 36–49. doi:10.1038/nrendo.2016.135
Hata T., Rehman F., Hori T., Nguyen J. H. (2019). GABA, γ-aminobutyric acid, protects against severe liver injury. J. Surg. Res. 236, 172–183. doi:10.1016/j.jss.2018.11.047
He G., Karin M. (2011). NF-κB and STAT3 - key players in liver inflammation and cancer. Cell. Res. 21 (1), 159–168. doi:10.1038/cr.2010.183
Hernández E. Á., Kahl S., Seelig A., Begovatz P., Irmler M., Kupriyanova Y., et al. (2017). Acute dietary fat intake initiates alterations in energy metabolism and insulin resistance. J. Clin. investigation 127 (2), 695–708. doi:10.1172/JCI89444
Humpton T. J., Hall H., Kiourtis C., Nixon C., Clark W., Hedley A., et al. (2022). p53-mediated redox control promotes liver regeneration and maintains liver function in response to CCl4. Cell. death Differ. 29 (3), 514–526. doi:10.1038/s41418-021-00871-3
Jelenik T., Kaul K., Séquaris G., Flögel U., Phielix E., Kotzka J., et al. (2017). Mechanisms of insulin resistance in primary and secondary nonalcoholic fatty liver. Diabetes 66 (8), 2241–2253. doi:10.2337/db16-1147
Ji J., Zhang L., Wang P., Mu Y.-M., Zhu X.-Y., Wu Y.-Y., et al. (2005). Saturated free fatty acid, palmitic acid, induces apoptosis in fetal hepatocytes in culture. Exp. Toxicol. pathology official J. Gesellschaft fur Toxikologische Pathologie 56 (6), 369–376. doi:10.1016/j.etp.2005.02.003
Ji F., Zhou M., Sun Z., Jiang Z., Zhu H., Xie Z., et al. (2021). Integrative proteomics reveals the role of E3 ubiquitin ligase SYVN1 in hepatocellular carcinoma metastasis. Cancer Commun. Lond. Engl. 41 (10), 1007–1023. doi:10.1002/cac2.12192
Johnson R. F., Witzel I.-I., Perkins N. D. (2011). p53-dependent regulation of mitochondrial energy production by the RelA subunit of NF-κB. Cancer Res. 71 (16), 5588–5597. doi:10.1158/0008-5472.CAN-10-4252
Jonas W., Schürmann A. (2021). Genetic and epigenetic factors determining NAFLD risk. Mol. Metab. 50, 101111. doi:10.1016/j.molmet.2020.101111
Joshi-Barve S., Barve S. S., Butt W., Klein J., McClain C. J. (2003). Inhibition of proteasome function leads to NF-kappaB-independent IL-8 expression in human hepatocytes. Hepatol. Baltim. Md 38 (5), 1178–1187. doi:10.1053/jhep.2003.50470
Kikkert M., Doolman R., Dai M., Avner R., Hassink G., van Voorden S., et al. (2004). Human HRD1 is an E3 ubiquitin ligase involved in degradation of proteins from the endoplasmic reticulum. J. Biol. Chem. 279 (5), 3525–3534. doi:10.1074/jbc.M307453200
Klein B. J., Piao L., Xi Y., Rincon-Arano H., Rothbart S. B., Peng D., et al. (2014). The histone-H3K4-specific demethylase KDM5B binds to its substrate and product through distinct PHD fingers. Cell. Rep. 6 (2), 325–335. doi:10.1016/j.celrep.2013.12.021
Koliaki C., Roden M. (2016). Alterations of mitochondrial function and insulin sensitivity in human obesity and diabetes mellitus. Annu. Rev. Nutr. 36, 337–367. doi:10.1146/annurev-nutr-071715-050656
Koliaki C., Szendroedi J., Kaul K., Jelenik T., Nowotny P., Jankowiak F., et al. (2015). Adaptation of hepatic mitochondrial function in humans with non-alcoholic fatty liver is lost in steatohepatitis. Cell. metab. 21 (5), 739–746. doi:10.1016/j.cmet.2015.04.004
Krämer A., Green J., Pollard J., Tugendreich S. (2014). Causal analysis approaches in ingenuity pathway analysis. Bioinforma. Oxf. Engl. 30 (4), 523–530. doi:10.1093/bioinformatics/btt703
Kunst C., Haderer M., Heckel S., Schlosser S., Müller M. (2016). The p53 family in hepatocellular carcinoma. Transl. Cancer Res. 5 (6), 632–638. doi:10.21037/tcr.2016.11.79
Lee I. H., Kawai Y., Fergusson M. M., Rovira I. I., Bishop A. J. R., Motoyama N., et al. (2012). Atg7 modulates p53 activity to regulate cell cycle and survival during metabolic stress. Sci. (New York, N.Y.) 336 (6078), 225–228. doi:10.1126/science.1218395
Lee K., Haddad A., Osme A., Kim C., Borzou A., Ilchenko S., et al. (2018). Hepatic mitochondrial defects in a nonalcoholic fatty liver disease mouse model are associated with increased degradation of oxidative phosphorylation subunits. Mol. Cell. proteomics MCP 17 (12), 2371–2386. doi:10.1074/mcp.RA118.000961
Li X., Liu L., Yang S., Song N., Zhou X., Gao J., et al. (2014). Histone demethylase KDM5B is a key regulator of genome stability. Proc. Natl. Acad. Sci. U. S. A. 111 (19), 7096–7101. doi:10.1073/pnas.1324036111
Li W., Wong C. C., Zhang X., Kang W., Nakatsu G., Zhao Q., et al. (2018). CAB39L elicited an anti-Warburg effect via a LKB1-AMPK-PGC1α axis to inhibit gastric tumorigenesis. Oncogene 37 (50), 6383–6398. doi:10.1038/s41388-018-0402-1
Li H., Lan T., Xu L., Liu H., Wang J., Li J., et al. (2020). NCSTN promotes hepatocellular carcinoma cell growth and metastasis via β-catenin activation in a Notch1/AKT dependent manner. J. Exp. Clin. cancer Res. CR 39 (1), 128. doi:10.1186/s13046-020-01638-3
Li K., Zhang K., Wang H., Wu Y., Chen N., Chen J., et al. (2021). Hrd1-mediated ACLY ubiquitination alleviate NAFLD in db/db mice. Metabolism Clin. Exp. 114, 154349. doi:10.1016/j.metabol.2020.154349
Li Y., Lin Y., Han X., Li W., Yan W., Ma Y., et al. (2021). GSK3 inhibitor ameliorates steatosis through the modulation of mitochondrial dysfunction in hepatocytes of obese patients. iScience 24 (3), 102149. doi:10.1016/j.isci.2021.102149
Marcel V., Perrier S., Aoubala M., Ageorges S., Groves M. J., Diot A., et al. (2010). Δ160p53 is a novel N-terminal p53 isoform encoded by Δ133p53 transcript. FEBS Lett. 584 (21), 4463–4468. doi:10.1016/j.febslet.2010.10.005
Minuk G. Y., Bear C. E., Sarjeant E. J. (1987). Sodium-independent, bicuculline-sensitive 3HGABA binding to isolated rat hepatocytes. Am. J. physiology 252 (5), G642–G647. doi:10.1152/ajpgi.1987.252.5.G642
Mota M., Banini B. A., Cazanave S. C., Sanyal A. J. (2016). Molecular mechanisms of lipotoxicity and glucotoxicity in nonalcoholic fatty liver disease. Metabolism Clin. Exp. 65 (8), 1049–1061. doi:10.1016/j.metabol.2016.02.014
Murru E., Manca C., Carta G., Banni S. (2022). Impact of dietary palmitic acid on lipid metabolism. Front. Nutr. 9, 861664. doi:10.3389/fnut.2022.861664
Nicol S. M., Bray S. E., Black H. D., Lorimore S. A., Wright E. G., Lane D. P., et al. (2013). The RNA helicase p68 (DDX5) is selectively required for the induction of p53-dependent p21 expression and cell-cycle arrest after DNA damage. Oncogene 32 (29), 3461–3469. doi:10.1038/onc.2012.426
Norikura T., Kojima-Yuasa A., Opare Kennedy D., Matsui-Yuasa I. (2007). Protective effect of gamma-aminobutyric acid (GABA) against cytotoxicity of ethanol in isolated rat hepatocytes involves modulations in cellular polyamine levels. Amino acids 32 (3), 419–423. doi:10.1007/s00726-006-0381-3
Nowotny B., Zahiragic L., Krog D., Nowotny P. J., Herder C., Carstensen M., et al. (2013). Mechanisms underlying the onset of oral lipid-induced skeletal muscle insulin resistance in humans. Diabetes 62 (7), 2240–2248. doi:10.2337/db12-1179
Panasiuk A., Dzieciol J., Panasiuk B., Prokopowicz D. (2006). Expression of p53, Bax and Bcl-2 proteins in hepatocytes in non-alcoholic fatty liver disease. World J. gastroenterology 12 (38), 6198–6202. doi:10.3748/wjg.v12.i38.6198
Pérez-Schindler J., Vargas-Fernández E., Karrer-Cardel B., Ritz D., Schmidt A., Handschin C. (2022). Characterization of regulatory transcriptional mechanisms in hepatocyte lipotoxicity. Sci. Rep. 12 (1), 11477. doi:10.1038/s41598-022-15731-4
Piccolis M., Bond L. M., Kampmann M., Pulimeno P., Chitraju C., Jayson C. B. K., et al. (2019). Probing the global cellular responses to lipotoxicity caused by saturated fatty acids. Mol. Cell. 74 (1), 32–44. doi:10.1016/j.molcel.2019.01.036
Qu J., Zou T., Lin Z. (2021). The roles of the ubiquitin-proteasome system in the endoplasmic reticulum stress pathway. Int. J. Mol. Sci. 22 (4), 1526. doi:10.3390/ijms22041526
Rodriguez-Lebron E., Denovan-Wright E. M., Nash K., Lewin A. S., Mandel R. J. (2005). Intrastriatal rAAV-mediated delivery of anti-huntingtin shRNAs induces partial reversal of disease progression in R6/1 Huntington's disease transgenic mice. Mol. Ther. J. Am. Soc. Gene Ther. 12 (4), 618–633. doi:10.1016/j.ymthe.2005.05.006
Rohbeck E., Hasse B., Koopmans G., Romero A., Belgardt B.-F., Roden M., et al. (2022). Positive allosteric γ-aminobutyric acid type A receptor modulation prevents lipotoxicity-induced injury in hepatocytes in vitro. Diabetes, Obes. metabolism 24 (8), 1498–1508. doi:10.1111/dom.14719
Rutkowski D. T., Wu J., Back S.-H., Callaghan M. U., Ferris S. P., Iqbal J., et al. (2008). UPR pathways combine to prevent hepatic steatosis caused by ER stress-mediated suppression of transcriptional master regulators. Dev. Cell. 15 (6), 829–840. doi:10.1016/j.devcel.2008.10.015
Ryaboshapkina M., Hammar M. (2017). Human hepatic gene expression signature of non-alcoholic fatty liver disease progression, a meta-analysis. Sci. Rep. 7 (1), 12361. doi:10.1038/s41598-017-10930-w
Sharma L. K., Lu J., Bai Y. (2009). Mitochondrial respiratory complex I: Structure, function and implication in human diseases. Curr. Med. Chem. 16 (10), 1266–1277. doi:10.2174/092986709787846578
Sherman B. T., Hao M., Qiu J., Jiao X., Baseler M. W., Lane H. C., et al. (2022). David: A web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic acids Res. 50, W216–W221. doi:10.1093/nar/gkac194
Shilpa J., Roshni B. T., Chinthu R., Paulose C. S. (2012). Role of GABA and serotonin coupled chitosan nanoparticles in enhanced hepatocyte proliferation. J. Mater. Sci. Mater. Med. 23 (12), 2913–2921. doi:10.1007/s10856-012-4754-8
Staňková P., Kučera O., Peterová E., Elkalaf M., Rychtrmoc D., Melek J., et al. (2021). Western diet decreases the liver mitochondrial oxidative flux of succinate: Insight from a murine NAFLD model. Int. J. Mol. Sci. 22, 6908. doi:10.3390/ijms22136908
Steensels S., Qiao J., Ersoy B. A. (2020). Transcriptional regulation in non-alcoholic fatty liver disease. Metabolites 10, 283. doi:10.3390/metabo10070283
Targher G., Corey K. E., Byrne C. D., Roden M. (2021). The complex link between NAFLD and type 2 diabetes mellitus - mechanisms and treatments. Nat. Rev. Gastroenterology hepatology 18 (9), 599–612. doi:10.1038/s41575-021-00448-y
Tilg H., Moschen A. R., Roden M. (2017). NAFLD and diabetes mellitus. Nat. Rev. Gastroenterology hepatology 14 (1), 32–42. doi:10.1038/nrgastro.2016.147
Valera A. G., Díaz-Hernández M., Hernández F., Ortega Z., Lucas J. J. (2005). The ubiquitin-proteasome system in Huntington's disease. Neurosci. a Rev. J. bringing Neurobiol. neurology psychiatry 11 (6), 583–594. doi:10.1177/1073858405280639
Váraljai R., Islam A. B. M. M. K., Beshiri M. L., Rehman J., Lopez-Bigas N., Benevolenskaya E. V. (2015). Increased mitochondrial function downstream from KDM5A histone demethylase rescues differentiation in pRB-deficient cells. Genes. & Dev. 29 (17), 1817–1834. doi:10.1101/gad.264036.115
Vesting A. J., Jais A., Klemm P., Steuernagel L., Wienand P., Fog-Tonnesen M., et al. (2022). NIK/MAP3K14 in hepatocytes orchestrates NASH to hepatocellular carcinoma progression via JAK2/STAT5 inhibition. Mol. Metab. 66, 101626. doi:10.1016/j.molmet.2022.101626
Vock C., Gleissner M., Klapper M., Döring F. (2007). Identification of palmitate-regulated genes in HepG2 cells by applying microarray analysis. Biochimica biophysica acta 1770 (9), 1283–1288. doi:10.1016/j.bbagen.2007.07.001
Vogelstein B., Lane D., Levine A. J. (2000). Surfing the p53 network. Nature 408 (6810), 307–310. doi:10.1038/35042675
Wadgaonkar R., Phelps K. M., Haque Z., Williams A. J., Silverman E. S., Collins T. (1999). CREB-binding protein is a nuclear integrator of nuclear factor-kappaB and p53 signaling. J. Biol. Chem. 274 (4), 1879–1882. doi:10.1074/jbc.274.4.1879
Wei J., Yuan Y., Chen L., Xu Y., Zhang Y., Wang Y., et al. (2018). ER-associated ubiquitin ligase HRD1 programs liver metabolism by targeting multiple metabolic enzymes. Nat. Commun. 9 (1), 3659. doi:10.1038/s41467-018-06091-7
Wingelhofer B., Neubauer H. A., Valent P., Han X., Constantinescu S. N., Gunning P. T., et al. (2018). Implications of STAT3 and STAT5 signaling on gene regulation and chromatin remodeling in hematopoietic cancer. Leukemia 32 (8), 1713–1726. doi:10.1038/s41375-018-0117-x
Wójcik C. (2002). Regulation of apoptosis by the ubiquitin and proteasome pathway. J. Cell. Mol. Med. 6 (1), 25–48. doi:10.1111/j.1582-4934.2002.tb00309.x
Wu T., Zhao F., Gao B., Tan C., Yagishita N., Nakajima T., et al. (2014). Hrd1 suppresses Nrf2-mediated cellular protection during liver cirrhosis. Genes. & Dev. 28 (7), 708–722. doi:10.1101/gad.238246.114
Xing Z., Russon M. P., Utturkar S. M., Tran E. J. (2020). The RNA helicase DDX5 supports mitochondrial function in small cell lung cancer. J. Biol. Chem. 295 (27), 8988–8998. doi:10.1074/jbc.RA120.012600
Xu K., Sun S., Yan M., Cui J., Yang Y., Li W., et al. (2022). DDX5 and DDX17-multifaceted proteins in the regulation of tumorigenesis and tumor progression. Front. Oncol. 12, 943032. doi:10.3389/fonc.2022.943032
Yamasaki S., Yagishita N., Sasaki T., Nakazawa M., Kato Y., Yamadera T., et al. (2007). Cytoplasmic destruction of p53 by the endoplasmic reticulum-resident ubiquitin ligase 'Synoviolin. EMBO J. 26 (1), 113–122. doi:10.1038/sj.emboj.7601490
Yan Z., Miao X., Zhang B., Xie J. (2018). p53 as a double-edged sword in the progression of non-alcoholic fatty liver disease. Life Sci. 215, 64–72. doi:10.1016/j.lfs.2018.10.051
Yang H., Zhong X., Ballar P., Luo S., Shen Y., Rubinsztein D. C., et al. (2007). Ubiquitin ligase Hrd1 enhances the degradation and suppresses the toxicity of polyglutamine-expanded huntingtin. Exp. Cell. Res. 313 (3), 538–550. doi:10.1016/j.yexcr.2006.10.031
Zhang Y., Ye S., Lu W., Zhong J., Leng Y., Yang T., et al. (2022). RNA helicase DEAD-box protein 5 alleviates nonalcoholic steatohepatitis progression via tethering TSC complex and suppressing mTORC1 signaling. Hepatol. Online ahead of print. doi:10.1002/hep.32651
Keywords: hepatocytes, lipotoxicity, GABAA receptor, NAFLD, RNA sequencing (RNAseq)
Citation: Rohbeck E, Niersmann C, Köhrer K, Wachtmeister T, Roden M, Eckel J and Romacho T (2023) Positive allosteric GABAA receptor modulation counteracts lipotoxicity-induced gene expression changes in hepatocytes in vitro. Front. Physiol. 14:1106075. doi: 10.3389/fphys.2023.1106075
Received: 23 November 2022; Accepted: 31 January 2023;
Published: 13 February 2023.
Edited by:
Sonia Michael Najjar, Ohio University, United StatesReviewed by:
Komal Ramani, Cedars Sinai Medical Center, United StatesSivarajan Kumarasamy, Ohio University, United States
Xavier Vidal-Gomez, INSERM U1046 Physiologie et Médecine Expérimentale du Coeur et des Muscles, France
Copyright © 2023 Rohbeck, Niersmann, Köhrer, Wachtmeister, Roden, Eckel and Romacho. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Tania Romacho, tromacho@ual.es