GBA3 promotes fatty acid oxidation and alleviates non-alcoholic fatty liver by increasing CPT2 transcription

Juyi Li; Yingqun Ni; Yuanyuan Zhang; Huaizhen Liu

doi:10.18632/aging.205616

Research Paper Volume 16, Issue 5 pp 4591—4608

GBA3 promotes fatty acid oxidation and alleviates non-alcoholic fatty liver by increasing CPT2 transcription

Juyi Li^1, , Yingqun Ni^2, , Yuanyuan Zhang^1, , Huaizhen Liu^1, ,

¹ Department of Endocrinology, Geriatrics Center, The First Affiliated Hospital of Anhui University of Chinese Medicine, Hefei 230001, Anhui, China
² Department of Endocrinology, The First Affiliated Hospital of Anhui University of Chinese Medicine, Hefei 230001, Anhui, China

Received: October 19, 2023 Accepted: January 23, 2024 Published: February 29, 2024

https://doi.org/10.18632/aging.205616
How to Cite

Copyright: © 2024 Li et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Background: Excessive lipids accumulation and hepatocytes death are prominent characteristics of non-alcoholic fatty liver disease (NAFLD). Nonetheless, the precise pathophysiological mechanisms are not fully elucidated.

Methods: HepG2 cells stimulated with palmitic acids and rats fed with high-fat diet were used as models for NAFLD. The impact of Glucosylceramidase Beta 3 (GBA3) on fatty acid oxidation (FAO) was assessed using Seahorse metabolic analyzer. Lipid content was measured both in vitro and in vivo. To evaluate NAFLD progression, histological analysis was performed along with measurements of inflammatory factors and liver enzyme levels. Western blot and immunohistochemistry were employed to examine the activity levels of necroptosis. Flow cytometry and reactive oxygen species (ROS) staining were utilized to assess levels of oxidative stress.

Results: GBA3 promoted FAO and enhanced the mitochondrial membrane potential without affecting glycolysis. These reduced the lipid accumulation. Rats supplemented with GBA3 exhibited lower levels of inflammatory factors and liver enzymes, resulting in a slower progression of NAFLD. GBA3 overexpression reduced ROS and the ratio of cell apoptosis. Phosphorylation level was reduced in the essential mediator, MLKL, implicated in necroptosis. Mechanistically, as a transcriptional coactivator, GBA3 promoted the expression of Carnitine Palmitoyltransferase 2 (CPT2), which resulted in enhanced FAO.

Conclusions: Increased FAO resulting from GBA3 reduced oxidative stress and the production of ROS, thereby inhibiting necroptosis and delaying the progression of NAFLD. Our research offers novel insights into the potential therapeutic applications of GBA3 and FAO in the management and treatment of NAFLD.

Introduction

Non-alcoholic fatty liver disease (NAFLD) is a prevalent and widespread chronic liver disorder on a global scale [1, 2]. It is estimated that the prevalence of NAFLD has reached 25%, even higher in certain regions [3]. As a consequence of such a high prevalence, NAFLD is expected to become a major cause of liver transplantation and liver cancer [4]. Although still controversial in clinical practice, there has been a proposal to rename NAFLD as metabolic fatty liver disease [5]. As the main hepatic manifestation of metabolic syndrome, NAFLD exhibits a broad spectrum of histological changes, ranging from mild fatty degeneration to more aggressive inflammation and fibrosis. The accumulation of lipids resulting from an imbalance between lipid intake and disposal serves as the starting point for this cascade of pathological events. Elevated intracellular levels of free fatty acids (FFAs), especially saturated fatty acids, induce lipotoxicity and cell death [6]. Cell death represents the endpoint of hepatocyte fate. Subsequent inflammation and loss of hepatic parenchyma drive disease progression. Increasing evidence suggests that the pattern of hepatocyte death during NAFLD involves necroptosis. NAFLD is one of the few human diseases that activate necroptosis without inhibiting apoptosis [7]. The core molecules mediating necroptosis have been shown to play a role in NAFLD [8]. However, the molecular basis of lipid metabolism driving hepatocyte necroptosis remains unclear.

Glucosylceramidase Beta 3 (GBA3), also known as beta-glucosidase, exhibits broad substrate specificity and hydrolyzes many plant β-glycosides in the human diet [9]. GBA3 may play an important role in the interpretation or biotransformation of glycosides [10, 11]. Due to GBA3’s metabolic activity, its expression in the liver shows relatively high levels [12]. Research has found that decreased expression of GBA3 is associated with poor prognosis in hepatocellular carcinoma [13]. However, the role of GBA3 in NAFLD is still unknown.

In this study, we identified GBA3 as a core gene associated with impaired fatty acid oxidation (FAO) in NAFLD through bioinformatics analysis. Further research revealed that GBA3 enhanced FAO by promoting the transcription of Carnitine Palmitoyltransferase 2 (CPT2). Increased FAO reduced lipid accumulation and levels of reactive oxygen species (ROS) in the cells, protecting hepatocytes from necroptosis and improving NAFLD. These findings demonstrate how lipid homeostasis is regulated through the GBA3/CPT2 signaling pathway during the progression of NAFLD. This may provide a potential strategy for the treatment of NAFLD.

Materials and Methods

Acquisition of FAO-related differentially expressed genes (FRDEGs)

The microarray expression data (GSE160016) were downloaded from the Gene Expression Omnibus. The dataset contained liver transcriptome data from five NAFLD donors and six non-NAFLD donors [14]. When multiple probes were associated with the same gene symbol, the mean value was chosen as the representative expression value for that particular gene.

The “limma” R package with |log fold change| > 2 and screening criteria of P < 0.05 was used to obtain differentially expressed genes (DEGs) [15]. Genes related to FAO were obtained from the GeneCards website [16]. FRDEGs were defined as the genes that overlapped between the set of DEGs and the set of FAO-related genes.

Enrichment analysis was conducted using the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) databases to gain insight into the signaling pathways and biological functions of FRDEGs.

Screening of hub FRDEGs

In GSE160016, the FAO activity score was determined by gene set variation analysis [17]. The preset gene set used was REACTOME_MITOCHONDRIAL_FATTY_ACID_BETA_OXIDATION.v2023.1.Hs.gmt [18]. Pearson correlation coefficients were calculated between the expression levels of FRDEGs and the FAO activity score. A receiver operating characteristic curve was employed to explore the potential diagnostic value of the FRDEGs for NAFLD. The mRNA expression levels of hub genes were measured using qRT-PCR in cell and rat models of NAFLD.

Animals

Male Sprague-Dawley rats aged 4 weeks were obtained from Vital River (Beijing, China). The rats were fed a high-fat diet (60 kcal % fat) for 16 weeks to induce NAFLD [19]. Rats were subjected to adenoviral transduction by administering 1 x 10¹¹ viral copies per rat through tail vein injection. We used etomoxir (1 mg/kg, intraperitoneal, every 2 days) to inhibit FAO [20].

Animals were housed in specific pathogen-free environments with controlled lighting (12 hours light and 12 hours dark cycles), temperature (24±2° C), and humidity (50±10%). Animals had unrestricted access to food and water throughout the experiment. Rats were euthanized by administering excessive isoflurane inhalation before liver harvesting.

Cell culture and transfection

The HepG2 cell line was obtained from Procell Biotechnology Co., Ltd. (Wuhan, China). The cells were cultured in MEM supplemented with 10% fetal bovine serum (v/v), penicillin G (100 units/mL), and streptomycin (100 μg/mL). The cell culture was maintained at 37° C in a humidified atmosphere with 5% CO₂. A NAFLD cell model was established using palmitic acid (PA) according to previous research [21].

Lentivirus containing pLVX-GBA3-Puro and pLVX-Puro (vector) was purchased from Genechem (Shanghai, China). Stably transfected cells were selected with puromycin and then validated using qRT-PCR and Western blotting.

SiRNA targeting EP300 was purchased from GenePharma (Suzhou, China). Analysis of knockdown efficiency and other experiments were conducted 48 h post-transfection.

Western blot and co-immunoprecipitation

Cells were prepared with RIPA buffer (Cell Signaling Technology, USA) or Pierce IP lysis buffer (Thermo Fisher Scientific, USA). According to previous study, Western blot and co-immunoprecipitation were performed [22, 23]. In brief, the lysed cells were collected and then centrifuged at 14,000 × g to eliminate cell debris. Protein concentrations were quantified using the BCA Protein Assay Kit (P0009, Beyotime, China). Each sample, containing 10 μg of total protein, was subjected to separation by SDS-PAGE gel and then transferred onto PVDF membranes. The PVDF membrane was incubated with the primary antibody overnight at 4° C. Then, the PVDF membranes were incubated with peroxidase conjugated avidin goat anti-rabbit IgG or goat anti-mouse IgG (1:5000) for one hour at room temperature. The membranes were subsequently scanned, and the protein levels were normalized to β-actin (1:1000) as a control. The Tanon-5200 chemiluminescence imaging system (Shanghai, China) and ImageJ software (NIH, USA) were utilized for recording and quantifying signal intensities. For co-immunoprecipitation, the extract was incubated with the corresponding antibody overnight at 4° C. Protein A&G beads (Bersinbio, Guangzhou, China) were then added and incubated at 4° C for 4 h. The coprecipitated proteins were then washed with SDS loading buffer for 5 min at 95° C. The subsequent results were obtained through Western blotting as previously described. All antibodies were listed in Supplementary Table 1.

Immunofluorescence and PCR

For immunofluorescence, all procedures were performed as previously reported [24]. The results were visualized by Zeiss laser confocal microscope 900.

Cells or tissues were lysed using TRIzol reagent (Life Technologies, USA). Based on previous study, PCR was performed to determine mRNA expression [21]. The cDNA was synthesized using the ReverTra Ace qPCR RT Kit (Toyobo, Osaka, Japan) and reverse transcription was conducted on a MasterCycler (Eppendorf, GRE, Hamburg, Germany). Subsequently, quantitative real-time PCR (qRT-PCR) was carried out on a LightCycler 480 System (Roche, CH, Basel, Switzerland) with the LightCycler 480 software 1.5.1.62 SP3. B-actin was set as control. The relative mRNA expression was determined using the 2^–ΔΔCt method. The primer sequences can be found in Supplementary Table 2.

Cell metabolic analysis

The Seahorse XF Glycolytic Rate Assay Kit was utilized to detect glycolysis in cells, as reported in published literature [25]. To measure FAO through oxygen consumption (OCR) in cells, the Seahorse XF Palmitate Oxidation Stress Test Kit was employed.

The mitochondrial membrane potential (Δψm) was measured using JC-10 (Solarbio, Beijing, China). The relative content of tissue and intracellular triglycerides and cholesterol was quantified using commercial kits (Applygen, Beijing, China). The relative content of FFA was measured using a colorimetric assay kit (Elabscience, Wuhan, China). Intracellular lipid droplets were visualized using BODIPY 493/503 (MedChemExpress, USA) and Zeiss laser confocal microscope 900 [26].

Histological analysis

Livers were sectioned. H&E and Oil Red O stainings were performed using a commercial kit (Solarbio, China). In the IHC process, deparaffinized tissue sections were initially treated with 3% H₂O₂ for 15 minutes, followed by a 10-minute incubation in 5% BSA. Subsequently, the sections were subjected to overnight incubation at 4° C with primary antibodies. This was followed by incubation with biotin-conjugated secondary antibodies and visualization using the streptavidin-biotin staining technique. In the case of inadequate staining due to technical issues, samples were excluded.

The relative levels of tumor necrosis factor-a (TNF-α), interleukin 1β (IL1β), and interleukin 10 (IL10) were measured by enzyme-linked immunosorbent assay (ELISA) (R&D Systems, USA).

The relative activities of aminotransferase (ALT), aspartate aminotransferase (AST) and lactate dehydrogenase (LDH) in the liver were measured using commercial kits (Jiancheng, China).

Cell death and viability assays

A density of 1 × 10³ cells/well was used to seed cells in a 96-well plate for CCK-8 assays. At specific timepoints, 10 μl of CCK-8 reagent was added. The optical density (OD) was measured at 450 nm after incubating at 37° C for one hour. The fraction of DNA-replicating cells, representing cell proliferation status, was assessed using the EdU Detection Kit (RiboBio, Guangzhou, China).

Apoptosis was assessed using the Annexin V-FITC Apoptosis Detection Kit (Vazyme, China) along with flow cytometry. Cells located in the upper right quadrants (positive for PI and Annexin V) were classified as necrotic. Relative activity of caspase 8 was determined using the Caspase 8 Activity Assay Kit (Beyotime, China).

Oxidative stress determination

Tissue or cellular levels of glutathione (GSH) and oxidized glutathione disulfide (GSSG) were quantified by employing commercial kits (Beyotime, China).

Intracellular levels of ROS were measured using a commercial kit (Beyotime, China) and flow cytometry, as well as Zeiss Axio Vert A1 [27]. The hepatic ROS content was quantitatively determined using a commercial kit (Bestbio, China).

Chromatin immunoprecipitation (ChIP)

To investigate the interaction between EP300 and CPT2, Sonication ChIP Kit (RK20258, ABclonal, USA) was used according to the manufacturer’s procedure in our study.

Statistical analysis

All analyses were performed in R language (version 4.2.3). All data represented mean ± SEM from three or more experiments. Group differences were assessed using either a 2-tailed Student’s t-test or a Mann-Whitney U test (both 2-tailed). Statistical significance was set at P < 0.05. Commercial kits used in this study were listed entirely in Supplementary Table 3.

Data sharing statement

The data generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Results

FRDEGs are involved in NAFLD progression

We identified 41 downregulated genes and 17 upregulated genes from our analysis (Figure 1A). Out of these, 29 genes are associated with FAO. The expression patterns of these genes are depicted in Figure 1B.

Figure 1. Screening of FRDEGs from NAFLD. (A) Volcano map of DEGs. (B) Heat plot showed the normalized expression profile of 29 FRDEGs obtained from the overlap of DEGs in NAFLD and fatty acid oxidation related genes from GeneCards. (C, D) Functional enrichment analysis of FRDEGs in KEGG and GO.

Enrichment analysis unveiled the correlation between FRDEGs and diverse metabolic substances, such as linoleic acid and amino acids (Figure 1C). Additionally, FRDEGs play a role in the regulation of inflammation (such as IL17 signaling pathway) and cell death (Figure 1D). This implies that the FAO may impact inflammation and cell death through specific pathways in NAFLD.

GBA3 expression is significantly reduced in NAFLD

First, we calculated the FAO activity score. Consistent with previous studies [28], the activity of FAO was reduced in NAFLD (Figure 2A). This led us to focus on the downregulated FRDEGs that were positively correlated with FAO activity score.

Identification of GBA3 as a hub gene. (A) FAO scores in GSE160016, P-value was calculated by Mann-Whitney U test (*, P B) Pearson correlation coefficient between the expression of FRDEGs and FAO score (*, P C) The receiver operating characteristic curves of six statistically significant genes in GSE160016. (D) Relative mRNA expression of six genes in cells and livers. P-value was calculated by 2-tailed Student’s t test (*, P E) Immunoblotting detection of GBA3 in cells. Protein expression was normalized to actin levels and shown as relative values. P-value was calculated by 2-tailed Student’s t test (*, P F) Immunohistochemical images of GBA3 in livers.

Figure 2. Identification of GBA3 as a hub gene. (A) FAO scores in GSE160016, P-value was calculated by Mann-Whitney U test (*, P < 0.05). (B) Pearson correlation coefficient between the expression of FRDEGs and FAO score (*, P < 0.05). (C) The receiver operating characteristic curves of six statistically significant genes in GSE160016. (D) Relative mRNA expression of six genes in cells and livers. P-value was calculated by 2-tailed Student’s t test (*, P < 0.05; **, P < 0.01; ***, P < 0.001; ns, not significant). (E) Immunoblotting detection of GBA3 in cells. Protein expression was normalized to actin levels and shown as relative values. P-value was calculated by 2-tailed Student’s t test (*, P < 0.05). (F) Immunohistochemical images of GBA3 in livers.

Through Pearson test, we identified six genes that were significant (Figure 2B). Receiver operating characteristic curves demonstrated that these six genes could potentially serve as biomarkers for NAFLD (Figure 2C). To determine the core genes, we used PCR to assess the expression levels of these six genes in the NAFLD model (Figure 2D). The data showed that GBA3 mRNA levels were significantly decreased both in vivo and in vitro. Therefore, further investigations were conducted. Western blotting and immunohistochemistry revealed that the protein levels of GBA3 were also significantly reduced in NAFLD (Figure 2E, 2F).

GBA3 facilitates FAO and decreases lipid accumulation

First, we stably overexpressed GBA3 in cells (Figure 3A). Previous studies have shown that GBA3 can hydrolyze various substrates, including β-glucose linked to hydrophobic groups [29]. This led us to explore its relationship with sugar metabolism before discussing the role of GBA3 in FAO. The data indicated that GBA3 had no effect on extracellular acidification rate, neither in glycolytic capacity nor in reserves (Figure 3B). In the absence of glucose, we found that GBA3 promoted FAO. The mitochondrial oxygen consumption rate significantly increased (Figure 3C). After supplementing with exogenous PA as a substrate, mitochondrial respiration capacity further enhanced (Figure 3D). The improvement in FAO led to an increase in JC10 aggregates (Figure 3E), suggesting a significant elevation in mitochondrial membrane potential, resulting in more energy production.

GBA3 enhanced FAO and reduced the intracellular accumulation of lipids. (A) Constructing the cell line overexpressing GBA3. P-value was calculated by 2-tailed Student’s t test (*, P B) Extracellular acidification rate (ECAR) was analyzed, P-value was calculated by 2-tailed Student’s t test (ns, not significant). (C, D) Respiratory rates expressed as oxygen consumption rate (OCR) were measured using Seahorse metabolic analyzer. P-value was calculated by 2-tailed Student’s t test (***, P E) Representative images of JC-10 staining (Scale bars, 50 μm). (F) Lipid droplet was visualized by Bodipy 493/503 (green fluorescence), nuclei were counterstained blue with DAPI (blue fluorescence) (Scale bars, 20 μm). (G) Relative levels of free fatty acid, triglyceride and cholesterol were quantified in cells, P-value was calculated by 2-tailed Student’s t test (***, P

Figure 3. GBA3 enhanced FAO and reduced the intracellular accumulation of lipids. (A) Constructing the cell line overexpressing GBA3. P-value was calculated by 2-tailed Student’s t test (*, P < 0.01). (B) Extracellular acidification rate (ECAR) was analyzed, P-value was calculated by 2-tailed Student’s t test (ns, not significant). (C, D) Respiratory rates expressed as oxygen consumption rate (OCR) were measured using Seahorse metabolic analyzer. P-value was calculated by 2-tailed Student’s t test (***, P < 0.001). (E) Representative images of JC-10 staining (Scale bars, 50 μm). (F) Lipid droplet was visualized by Bodipy 493/503 (green fluorescence), nuclei were counterstained blue with DAPI (blue fluorescence) (Scale bars, 20 μm). (G) Relative levels of free fatty acid, triglyceride and cholesterol were quantified in cells, P-value was calculated by 2-tailed Student’s t test (***, P < 0.001).

Accumulation of lipids is one of the initiating factors of NAFLD. Images from confocal microscopy showed that GBA3 reduced lipids in cells (Figure 3F). As shown in Figure 3G, GBA3 led to a significant decrease in cellular content of FFA, triglycerides, and cholesterol. However, the effect of GBA3 was blocked by etomoxir. This indicates that the decrease in lipid content in cells is caused by the increased demand for substrates due to the higher FAO induced by GBA3.

GBA3 protects rats from NAFLD

So far, we have identified the crucial role of GBA3 in counteracting lipid overload by stimulating FAO in vitro. To determine the physiological function of GBA3, we transferred it to a NAFLD rat model. Before being fed a high-fat diet, the rats received adenovirus expressing GBA3 (Figure 4A).

GBA3 protected rats from NAFLD. (A) Successful supplementation of GBA3 in the liver of mice was confirmed through PCR and immunohistochemistry. P-value was calculated by 2-tailed Student’s t test (***, P B) H&E staining of livers. (C) Oil Red O staining was used to assess lipid content. (D) Relative levels of free fatty acid, triglyceride, cholesterol, inflammatory factors, ALT, AST and LDH were detected in livers, P-value was calculated by 2-tailed Student’s t test (*, P

Figure 4. GBA3 protected rats from NAFLD. (A) Successful supplementation of GBA3 in the liver of mice was confirmed through PCR and immunohistochemistry. P-value was calculated by 2-tailed Student’s t test (***, P < 0.001). (B) H&E staining of livers. (C) Oil Red O staining was used to assess lipid content. (D) Relative levels of free fatty acid, triglyceride, cholesterol, inflammatory factors, ALT, AST and LDH were detected in livers, P-value was calculated by 2-tailed Student’s t test (*, P < 0.05; **, P < 0.01; ***, P < 0.001).

Histological analysis showed that GBA3 alleviated hepatic steatosis induced by a high-fat diet (Figure 4B, 4C). Compared to the control group, GBA3 reduced lipid accumulation. Additionally, levels of inflammatory factors and transaminases in the liver were significantly decreased. Etomoxir partially reversed the effect of GBA3 (Figure 4D). These data support the physiological protective role of GBA3 in high-fat diet-induced NAFLD.

GBA3 protects cells from necroptosis by reducing ROS levels

Apoptosis of hepatocytes caused by lipid overload often promotes the development of NAFLD [30]. Therefore, we investigated whether the pro-FAO effect of GBA3 is beneficial for cell survival. CCK8 and EdU staining showed that GBA3 had no effect on cell proliferation (Figure 5A, 5B). We then turned our attention to apoptosis. Data from flow cytometry showed that GBA3 overexpression only reduced the percentage of cells in the upper right quadrant (PI positive and Annexin V positive) (Figure 5C). This indicates that the protective effect of GBA3 mainly prevents cells from undergoing late-stage apoptosis.

GBA3 reduced the occurrence of necroptosis. (A) EdU staining (green fluorescence) on cell proliferation, nuclei were counterstained blue with DAPI (Scale bars, 50 μm). (B) Cell proliferation was detected by measuring the optical density at 450 nm, P-value was calculated by 2-tailed Student’s t test (ns, not significant). (C) Percentage of apoptotic cell was assessed via flow cytometry, P-value was calculated by 2-tailed Student’s t test (ns, not significant; **, P D) Immunoblotting detection of MLKL and pMLKL (Ser358) in cells. Protein expression was normalized to actin levels and shown as relative values. P-value was calculated by 2-tailed Student’s t test (*, P E) The activity of caspase 8 was assessed by measuring the optical density at 405 nm. Incremental activity of caspase 8 was expressed as a percentage of control, P-value was calculated by 2-tailed Student’s t test (*, P F) Phosphorylation level of MLKL in livers detected by immunohistochemistry. (G) The infiltrated level of immune cells, P-value was calculated by Mann-Whitney U test (*, P

Figure 5. GBA3 reduced the occurrence of necroptosis. (A) EdU staining (green fluorescence) on cell proliferation, nuclei were counterstained blue with DAPI (Scale bars, 50 μm). (B) Cell proliferation was detected by measuring the optical density at 450 nm, P-value was calculated by 2-tailed Student’s t test (ns, not significant). (C) Percentage of apoptotic cell was assessed via flow cytometry, P-value was calculated by 2-tailed Student’s t test (ns, not significant; **, P < 0.01). (D) Immunoblotting detection of MLKL and pMLKL (Ser358) in cells. Protein expression was normalized to actin levels and shown as relative values. P-value was calculated by 2-tailed Student’s t test (*, P < 0.05). (E) The activity of caspase 8 was assessed by measuring the optical density at 405 nm. Incremental activity of caspase 8 was expressed as a percentage of control, P-value was calculated by 2-tailed Student’s t test (*, P < 0.05; ***, P < 0.001). (F) Phosphorylation level of MLKL in livers detected by immunohistochemistry. (G) The infiltrated level of immune cells, P-value was calculated by Mann-Whitney U test (*, P < 0.05; **, P < 0.01).

Necroptosis is a late-stage apoptosis that induces inflammation and infiltration of immune cells [31]. MLKL is the core mediator of necroptosis. Its phosphorylation is a specific and irreversible molecular event that drives necroptosis [32]. Consequently, we performed related experiments. As expected, GBA3 overexpression reduced the phosphorylation level of MLKL and the activity of Caspase 8 (Figure 5D, 5E). Similar trends were observed in vivo (Figure 5F). We divided all NAFLD samples from GSE160016 into high-expression and low-expression groups based on GBA3 expression level. Then, we compared the levels of immune cell infiltration between the two groups. The data showed that when GBA3 was highly expressed, the infiltration levels of activated B cells, natural killer cells, eosinophils, and T follicular helper cells decreased (Figure 5G). Most of the other changes in immune cells were not significant. These findings suggest that GBA3 can inhibit the transition of cells to necroptosis and prevent the occurrence of inflammation.

However, it is still unclear how increased FAO protects cells from necroptosis. Some studies have drawn our attention to this. Robert O Watson et al. found that excess ROS in mitochondria can trigger necrotic apoptosis [33]. Enhanced FAO promoted the generation of reducing substances and reduced ROS production [34]. In this study, we found that GBA3 increased the GSH/GSSG ratio (Figure 6A, 6B) and reduced ROS levels (Figure 6C–6E), both in vitro and in vivo. In conclusion, the above data suggest that GBA3 protects cells from necroptosis by increasing FAO and reducing ROS levels.

GBA3 maintains mitochondrial homeostasis by reducing ROS. (A, B) Levels of GSH and GSSG were quantified (**, P C) Representative image of ROS detection in cells (Scale bars, 50 μm). (D) ROS was quantified via flow cytometry, P-value was calculated by 2-tailed Student’s t test (*, P E) Hepatic ROS content was quantitatively determined, P-value was calculated by 2-tailed Student’s t test (**, P

Figure 6. GBA3 maintains mitochondrial homeostasis by reducing ROS. (A, B) Levels of GSH and GSSG were quantified (**, P < 0.01; ***, P < 0.001). (C) Representative image of ROS detection in cells (Scale bars, 50 μm). (D) ROS was quantified via flow cytometry, P-value was calculated by 2-tailed Student’s t test (*, P < 0.05). (E) Hepatic ROS content was quantitatively determined, P-value was calculated by 2-tailed Student’s t test (**, P < 0.01; ***, P < 0.001).

GBA3 upregulates the expression of CPT2 to promote FAO

Finally, we focus on the specific mechanism of GBA3 in regulating FAO. FAO is regulated by a series of enzymes. PCR was used to preliminarily detect the expression levels of these enzymes. As shown in Figure 7A, the transcription level of CPT2 increased. WB showed an increase in the expression of CPT2 (Figure 7B). This indicates that CPT2 is involved in GBA3-mediated metabolic reprogramming. We obtained the transcription factors of CPT2 from the Cistrome Data Browser. Their expression in GSE160016 is shown in Figure 7C. As the expression of GBA3 is downregulated in NAFLD, we paid more attention to the transcription factors whose expression levels also decreased. PCR showed that only the expression of EP300 was decreased in the NAFLD cell model (Figure 7D). The results of WB further confirmed this finding (Figure 7E). As shown in Supplementary Figure 1, we demonstrated EP300 was transcription factor of CPT2 via ChIP. We infer the existence of interaction between GBA3 and EP300. Immunoprecipitation (Figure 7F) and immunofluorescence (Figure 7G) demonstrated our hypothesis.

Figure 7. GBA3 interacted with EP300 to promote CPT2 transcription. (A) Relative mRNA expression, P-value was calculated by 2-tailed Student’s t test (***, P < 0.001; ns, not significant). (B) Immunoblotting detection of CPT2. P-value was calculated by 2-tailed Student’s t test (*, P < 0.05). (C) Heat plot showed the normalized expression profile of transcription factors of CPT2 in GSE160016. (D) Relative mRNA expression of 8 transcription factors, P-value was calculated by 2-tailed Student’s t test (*, P < 0.05; **, P < 0.01; ns, not significant). (E) Immunoblotting detection of EP300. P-value was calculated by 2-tailed Student’s t test (*, P < 0.05). (F) Coimmunoprecipitation of GBA3 with endogenous EP300. (G) Representative images of immunofluorescence staining of GBA3 and EP300. Nuclei were counterstained blue with DAPI (Scale bars, 5 μm). (H) Knockdown of EP300 in cells overexpressing GBA3. (I) CPT2 expression was confirmed through PCR and Western blot. P-value was calculated by 2-tailed Student’s t test (***, P < 0.001). (J) Respiratory rates expressed as oxygen consumption rate (OCR) were measured using Seahorse metabolic analyzer. P-value was calculated by 2-tailed Student’s t test (***, P < 0.001).

We used siRNA to interfere with the expression of EP300 to investigate whether EP300 is involved in GBA3-mediated regulation of FAO (Figure 7H). The data showed that when the expression of EP300 was knocked down, the expression of CPT2 decreased (Figure 7I). The promoting effect of GBA3 on FAO was also blocked (Figure 7J). All these data indicate that GBA3 enhances FAO by interacting with EP300 and promoting the transcription of CPT2.

Discussion

Our study provides the first evidence of the crucial contribution of GBA3 in the pathology of NAFLD in liver cells. The data suggest that GBA3 promotes the transcription of CPT2, thus alleviating NAFLD by promoting FAO. Our study provides evidence supporting the toxic effects of fatty acid overload on liver cells [35]. Furthermore, it further indicates that the regulation of lipid metabolism is a potential strategy for managing NAFLD.

Imbalance in FFA metabolism is a key event leading to the occurrence of NAFLD [36]. In the highly reversible stage of simple steatosis, various cellular stress responses are activated to accelerate FFA processing, thus counteracting the progression of non-alcoholic fatty liver [37, 38]. Enhancing FAO in liver cells can greatly help the liver process excess lipids [39]. Consistent with these studies, our work suggests that GBA3 can enhance cytoplasmic FAO, thereby counteracting the pathological changes caused by NAFLD. If the development of NAFLD exceeds the functional adaptability of mitochondria, FAO may collapse, leading to the excessive production of lipotoxic metabolites and reactive oxygen species [40]. This will activate signaling pathways of inflammation, leading to hepatocyte apoptosis [41]. In fact, oxidative DNA damage has been detected in the livers of late-stage NAFLD patients [42].

In addition to FFA overload, hepatocyte necroptosis is also a key event in the process of NAFLD [43]. Necrosis and necroinflammation are typical histological manifestations of NAFLD [44]. Necroptosis simultaneously exhibits morphological and biochemical characteristics of both apoptosis and necrosis. Molecules that promote apoptosis, such as RIPK1, can also mediate necroptosis. Sanjoy Roychowdhury et al. found that the occurrence of necroptosis in hepatocytes may not require the activity of RIPK1. The role of RIPK3, another key molecule in necroptosis, in NAFLD remains controversial [45, 46]. The work of Stefan Wirtz et al. suggests that MLKL can mediate necroptosis independently of RIPK3. Furthermore, there is evidence that the loss of Caspase 8 in hepatocytes contributes to the initiation of necroptosis. Therefore, in our work, the combination test of MLKL and Caspase 8 was used to determine the activation of necroptosis. Given the limited knowledge of the association between lipid metabolism and necroptosis in NAFLD, we used cell and animal models to focus on the regulation of necroptosis by FAO. Our results indicate that FAO reduces the level of ROS, maintains mitochondrial homeostasis, and prevents the triggering of necroptosis.

Currently, research on the functions of GBA3 is still in the exploratory stage. In this study, we have demonstrated through in vitro experiments that there is an interaction between GBA3 and EP300. This contributes to the promotion of CPT2 by EP300. The non-enzymatic function of GBA3 has been previously discovered. The work of Tadashi Suzuki and others has demonstrated a physical interaction between GBA3 and NEU2 [47]. These phenomena suggest that the function of GBA3 may not be limited to enzymes. Future research may further explore the physiological role of GBA3 from a protein perspective.

Conclusions

In summary, this study provides direct evidence for the important role of GBA3 in the pathogenesis of NAFLD. Mouse and in vitro experiments demonstrate that GBA3 interacts with EP300 to promote the transcription of CPT2, thereby enhancing FAO. The resulting reduced ROS levels prevent necrotic apoptosis from being triggered. Our research supports the potential therapeutic value of GBA3 in NAFLD.

Supplementary Materials

Supplementary Figure 1

Supplementary Tables

Author Contributions

Huaizhen Liu planned the study and supervised the analyses. Juyi Li planned the study, performed the experiments, and drafted the manuscript. Yingqun Ni and Yuanyuan Zhang collected data and completed the pictures and tables. All authors contributed to the interpretation of results and revision of the paper, and approved the final manuscript.

Conflicts of Interest

The authors state that they have no conflicts of interest.

Ethical Statement

All experiments were approved by Laboratory Animal Ethics Committee of Anhui University of Traditional Chinese Medicine (No. 2022101) and adhered to the animal care and use guidelines outlined in the National Institutes of Health Publication No. 8023.

Funding

This research was supported by the National Natural Science Foundation of China (82274468) and the Scientific Research Program of Higher Education Institutions in Anhui Province (2023AH040103).

References

1. Pafili K, Roden M. Nonalcoholic fatty liver disease (NAFLD) from pathogenesis to treatment concepts in humans. Mol Metab. 2021; 50:101122. https://doi.org/10.1016/j.molmet.2020.101122 [PubMed]
2. Younossi ZM, Stepanova M, Afendy M, Fang Y, Younossi Y, Mir H, Srishord M. Changes in the prevalence of the most common causes of chronic liver diseases in the United States from 1988 to 2008. Clin Gastroenterol Hepatol. 2011; 9:524–30.e1. https://doi.org/10.1016/j.cgh.2011.03.020 [PubMed]
3. Younossi ZM, Koenig AB, Abdelatif D, Fazel Y, Henry L, Wymer M. Global epidemiology of nonalcoholic fatty liver disease-Meta-analytic assessment of prevalence, incidence, and outcomes. Hepatology. 2016; 64:73–84. https://doi.org/10.1002/hep.28431 [PubMed]
4. Calzadilla Bertot L, Adams LA. The Natural Course of Non-Alcoholic Fatty Liver Disease. Int J Mol Sci. 2016; 17:774. https://doi.org/10.3390/ijms17050774 [PubMed]
5. Eslam M, Newsome PN, Sarin SK, Anstee QM, Targher G, Romero-Gomez M, Zelber-Sagi S, Wai-Sun Wong V, Dufour JF, Schattenberg JM, Kawaguchi T, Arrese M, Valenti L, et al. A new definition for metabolic dysfunction-associated fatty liver disease: An international expert consensus statement. J Hepatol. 2020; 73:202–9. https://doi.org/10.1016/j.jhep.2020.03.039 [PubMed]
6. Palomer X, Pizarro-Delgado J, Barroso E, Vázquez-Carrera M. Palmitic and Oleic Acid: The Yin and Yang of Fatty Acids in Type 2 Diabetes Mellitus. Trends Endocrinol Metab. 2018; 29:178–90. https://doi.org/10.1016/j.tem.2017.11.009 [PubMed]
7. Schwabe RF, Luedde T. Apoptosis and necroptosis in the liver: a matter of life and death. Nat Rev Gastroenterol Hepatol. 2018; 15:738–52. https://doi.org/10.1038/s41575-018-0065-y [PubMed]
8. Tao L, Yi Y, Chen Y, Zhang H, Orning P, Lien E, Jie J, Zhang W, Xu Q, Li Y, Ding Z, Wu C, Ding Q, et al. RIP1 kinase activity promotes steatohepatitis through mediating cell death and inflammation in macrophages. Cell Death Differ. 2021; 28:1418–33. https://doi.org/10.1038/s41418-020-00668-w [PubMed]
9. Lopes-Marques M, Serrano C, Cardoso AR, Salazar R, Seixas S, Amorim A, Azevedo L, Prata MJ. GBA3: a polymorphic pseudogene in humans that experienced repeated gene loss during mammalian evolution. Sci Rep. 2020; 10:11565. https://doi.org/10.1038/s41598-020-68106-y [PubMed]
10. Berrin JG, McLauchlan WR, Needs P, Williamson G, Puigserver A, Kroon PA, Juge N. Functional expression of human liver cytosolic beta-glucosidase in Pichia pastoris. Insights into its role in the metabolism of dietary glucosides. Eur J Biochem. 2002; 269:249–58. https://doi.org/10.1046/j.0014-2956.2001.02641.x [PubMed]
11. Gopalan V, Pastuszyn A, Galey WR Jr, Glew RH. Exolytic hydrolysis of toxic plant glucosides by guinea pig liver cytosolic beta-glucosidase. J Biol Chem. 1992; 267:14027–32. [PubMed]
12. de Graaf M, van Veen IC, van der Meulen-Muileman IH, Gerritsen WR, Pinedo HM, Haisma HJ. Cloning and characterization of human liver cytosolic beta-glycosidase. Biochem J. 2001; 356:907–10. https://doi.org/10.1042/0264-6021:3560907 [PubMed]
13. Ying JF, Zhang YN, Song SS, Hu ZM, He XL, Pan HY, Zhang CW, Wang HJ, Li WF, Mou XZ. Decreased expression of GBA3 correlates with a poor prognosis in hepatocellular carcinoma patients. Neoplasma. 2020; 67:1139–45. https://doi.org/10.4149/neo_2020_190928N980 [PubMed]
14. Hou J, Zhang J, Cui P, Zhou Y, Liu C, Wu X, Ji Y, Wang S, Cheng B, Ye H, Shu L, Zhang K, Wang D, et al. TREM2 sustains macrophage-hepatocyte metabolic coordination in nonalcoholic fatty liver disease and sepsis. J Clin Invest. 2021; 131:e135197. https://doi.org/10.1172/JCI135197 [PubMed]
15. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015; 43:e47. https://doi.org/10.1093/nar/gkv007 [PubMed]
16. Stelzer G, Rosen N, Plaschkes I, Zimmerman S, Twik M, Fishilevich S, Stein TI, Nudel R, Lieder I, Mazor Y, Kaplan S, Dahary D, Warshawsky D, et al. The GeneCards Suite: From Gene Data Mining to Disease Genome Sequence Analyses. Curr Protoc Bioinformatics. 2016; 54:1.30.1-1.30.33. https://doi.org/10.1002/cpbi.5 [PubMed]
17. Hänzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics. 2013; 14:7. https://doi.org/10.1186/1471-2105-14-7 [PubMed]
18. Jassal B, Matthews L, Viteri G, Gong C, Lorente P, Fabregat A, Sidiropoulos K, Cook J, Gillespie M, Haw R, Loney F, May B, Milacic M, et al. The reactome pathway knowledgebase. Nucleic Acids Res. 2020; 48:D498–503. https://doi.org/10.1093/nar/gkz1031 [PubMed]
19. Jiao N, Baker SS, Chapa-Rodriguez A, Liu W, Nugent CA, Tsompana M, Mastrandrea L, Buck MJ, Baker RD, Genco RJ, Zhu R, Zhu L. Suppressed hepatic bile acid signalling despite elevated production of primary and secondary bile acids in NAFLD. Gut. 2018; 67:1881–91. https://doi.org/10.1136/gutjnl-2017-314307 [PubMed]
20. Mørkholt AS, Oklinski MK, Larsen A, Bockermann R, Issazadeh-Navikas S, Nieland JG, Kwon TH, Corthals A, Nielsen S, Nieland JD. Pharmacological inhibition of carnitine palmitoyl transferase 1 inhibits and reverses experimental autoimmune encephalitis in rodents. PLoS One. 2020; 15:e0234493. https://doi.org/10.1371/journal.pone.0234493 [PubMed]
21. Li J, Kou C, Sun T, Liu J, Zhang H. Identification and Validation of Hub Immune-Related Genes in Non-Alcoholic Fatty Liver Disease. Int J Gen Med. 2023; 16:2609–21. https://doi.org/10.2147/IJGM.S413545 [PubMed]
22. Li X, Ma J, Sun Z, Li N, Jiao G, Zhang T, Cao H. Development and validation of a N6 methylation regulator-related gene signature for prognostic and immune response prediction in non-small cell lung cancer. Am J Cancer Res. 2023; 13:2984–97. [PubMed]
23. Zhang L, Jiao G, You Y, Li X, Liu J, Sun Z, Li Q, Dai Z, Ma J, Zhou H, Li G, Meng C, Chen Y. Arginine methylation of PPP1CA by CARM1 regulates glucose metabolism and affects osteogenic differentiation and osteoclastic differentiation. Clin Transl Med. 2023; 13:e1369. https://doi.org/10.1002/ctm2.1369 [PubMed]
24. Yin S, Liu L, Brobbey C, Palanisamy V, Ball LE, Olsen SK, Ostrowski MC, Gan W. PRMT5-mediated arginine methylation activates AKT kinase to govern tumorigenesis. Nat Commun. 2021; 12:3444. https://doi.org/10.1038/s41467-021-23833-2 [PubMed]
25. Mistry JJ, Hellmich C, Moore JA, Jibril A, Macaulay I, Moreno-Gonzalez M, Di Palma F, Beraza N, Bowles KM, Rushworth SA. Free fatty-acid transport via CD36 drives β-oxidation-mediated hematopoietic stem cell response to infection. Nat Commun. 2021; 12:7130. https://doi.org/10.1038/s41467-021-27460-9 [PubMed]
26. Qiu B, Simon MC. BODIPY 493/503 Staining of Neutral Lipid Droplets for Microscopy and Quantification by Flow Cytometry. Bio Protoc. 2016; 6:e1912. https://doi.org/10.21769/BioProtoc.1912 [PubMed]
27. Guo C, Yang M, Jing L, Wang J, Yu Y, Li Y, Duan J, Zhou X, Li Y, Sun Z. Amorphous silica nanoparticles trigger vascular endothelial cell injury through apoptosis and autophagy via reactive oxygen species-mediated MAPK/Bcl-2 and PI3K/Akt/mTOR signaling. Int J Nanomedicine. 2016; 11:5257–76. https://doi.org/10.2147/IJN.S112030 [PubMed]
28. Hou T, Tian Y, Cao Z, Zhang J, Feng T, Tao W, Sun H, Wen H, Lu X, Zhu Q, Li M, Lu X, Liu B, et al. Cytoplasmic SIRT6-mediated ACSL5 deacetylation impedes nonalcoholic fatty liver disease by facilitating hepatic fatty acid oxidation. Mol Cell. 2022; 82:4099–115.e9. https://doi.org/10.1016/j.molcel.2022.09.018 [PubMed]
29. Daniels LB, Coyle PJ, Chiao YB, Glew RH, Labow RS. Purification and characterization of a cytosolic broad specificity beta-glucosidase from human liver. J Biol Chem. 1981; 256:13004–13. [PubMed]
30. Kanda T, Matsuoka S, Yamazaki M, Shibata T, Nirei K, Takahashi H, Kaneko T, Fujisawa M, Higuchi T, Nakamura H, Matsumoto N, Yamagami H, Ogawa M, et al. Apoptosis and non-alcoholic fatty liver diseases. World J Gastroenterol. 2018; 24:2661–72. https://doi.org/10.3748/wjg.v24.i25.2661 [PubMed]
31. Weinlich R, Oberst A, Beere HM, Green DR. Necroptosis in development, inflammation and disease. Nat Rev Mol Cell Biol. 2017; 18:127–36. https://doi.org/10.1038/nrm.2016.149 [PubMed]
32. Sun L, Wang H, Wang Z, He S, Chen S, Liao D, Wang L, Yan J, Liu W, Lei X, Wang X. Mixed lineage kinase domain-like protein mediates necrosis signaling downstream of RIP3 kinase. Cell. 2012; 148:213–27. https://doi.org/10.1016/j.cell.2011.11.031 [PubMed]
33. Weindel CG, Martinez EL, Zhao X, Mabry CJ, Bell SL, Vail KJ, Coleman AK, VanPortfliet JJ, Zhao B, Wagner AR, Azam S, Scott HM, Li P, et al. Mitochondrial ROS promotes susceptibility to infection via gasdermin D-mediated necroptosis. Cell. 2022; 185:3214–31.e23. https://doi.org/10.1016/j.cell.2022.06.038 [PubMed]
34. Choi HJ, Jhe YL, Kim J, Lim JY, Lee JE, Shin MK, Cheong JH. FoxM1-dependent and fatty acid oxidation-mediated ROS modulation is a cell-intrinsic drug resistance mechanism in cancer stem-like cells. Redox Biol. 2020; 36:101589. https://doi.org/10.1016/j.redox.2020.101589 [PubMed]
35. Bessone F, Razori MV, Roma MG. Molecular pathways of nonalcoholic fatty liver disease development and progression. Cell Mol Life Sci. 2019; 76:99–128. https://doi.org/10.1007/s00018-018-2947-0 [PubMed]
36. Lei K. Long non-coding RNAs regulate fatty acid and cholesterol metabolism. Genome Instability and Disease. 2022; 3:70–82. https://doi.org/10.1007/s42764-022-00070-5
37. Sunny NE, Bril F, Cusi K. Mitochondrial Adaptation in Nonalcoholic Fatty Liver Disease: Novel Mechanisms and Treatment Strategies. Trends Endocrinol Metab. 2017; 28:250–60. https://doi.org/10.1016/j.tem.2016.11.006 [PubMed]
38. Begriche K, Massart J, Robin MA, Bonnet F, Fromenty B. Mitochondrial adaptations and dysfunctions in nonalcoholic fatty liver disease. Hepatology. 2013; 58:1497–507. https://doi.org/10.1002/hep.26226 [PubMed]
39. Gariani K, Ryu D, Menzies KJ, Yi HS, Stein S, Zhang H, Perino A, Lemos V, Katsyuba E, Jha P, Vijgen S, Rubbia-Brandt L, Kim YK, et al. Inhibiting poly ADP-ribosylation increases fatty acid oxidation and protects against fatty liver disease. J Hepatol. 2017; 66:132–41. https://doi.org/10.1016/j.jhep.2016.08.024 [PubMed]
40. Guo Q, Li F, Duan Y, Wen C, Wang W, Zhang L, Huang R, Yin Y. Oxidative stress, nutritional antioxidants and beyond. Sci China Life Sci. 2020; 63:866–74. https://doi.org/10.1007/s11427-019-9591-5 [PubMed]
41. Anstee QM, Reeves HL, Kotsiliti E, Govaere O, Heikenwalder M. From NASH to HCC: current concepts and future challenges. Nat Rev Gastroenterol Hepatol. 2019; 16:411–28. https://doi.org/10.1038/s41575-019-0145-7 [PubMed]
42. Koliaki C, Szendroedi J, Kaul K, Jelenik T, Nowotny P, Jankowiak F, Herder C, Carstensen M, Krausch M, Knoefel WT, Schlensak M, Roden M. Adaptation of hepatic mitochondrial function in humans with non-alcoholic fatty liver is lost in steatohepatitis. Cell Metab. 2015; 21:739–46. https://doi.org/10.1016/j.cmet.2015.04.004 [PubMed]
43. Afonso MB, Rodrigues PM, Carvalho T, Caridade M, Borralho P, Cortez-Pinto H, Castro RE, Rodrigues CM. Necroptosis is a key pathogenic event in human and experimental murine models of non-alcoholic steatohepatitis. Clin Sci (Lond). 2015; 129:721–39. https://doi.org/10.1042/CS20140732 [PubMed]
44. Diehl AM, Day C. Cause, Pathogenesis, and Treatment of Nonalcoholic Steatohepatitis. N Engl J Med. 2017; 377:2063–72. https://doi.org/10.1056/NEJMra1503519 [PubMed]
45. Roychowdhury S, McCullough RL, Sanz-Garcia C, Saikia P, Alkhouri N, Matloob A, Pollard KA, McMullen MR, Croniger CM, Nagy LE. Receptor interacting protein 3 protects mice from high-fat diet-induced liver injury. Hepatology. 2016; 64:1518–33. https://doi.org/10.1002/hep.28676 [PubMed]
46. Gautheron J, Vucur M, Reisinger F, Cardenas DV, Roderburg C, Koppe C, Kreggenwinkel K, Schneider AT, Bartneck M, Neumann UP, Canbay A, Reeves HL, Luedde M, et al. A positive feedback loop between RIP3 and JNK controls non-alcoholic steatohepatitis. EMBO Mol Med. 2014; 6:1062–74. https://doi.org/10.15252/emmm.201403856 [PubMed]
47. Wang L, Seino J, Tomotake H, Funakoshi Y, Hirayama H, Suzuki T. Co-Expression of NEU2 and GBA3 Causes a Drastic Reduction in Cytosolic Sialyl Free N-glycans in Human MKN45 Stomach Cancer Cells-Evidence for the Physical Interaction of NEU2 and GBA3. Biomolecules. 2015; 5:1499–514. https://doi.org/10.3390/biom5031499 [PubMed]

Research Paper Volume 16, Issue 5 pp 4591—4608