- 1Department of Endocrinology and Metabolism, Affiliated Hospital of Weifang Medical University, Weifang, China
- 2Branch of Shandong Provincial Clinical Research Center for Diabetes and Metabolic Diseases, Weifang, China
- 3Clinical Research Center, Affiliated Hospital of Weifang Medical University, Weifang, China
- 4Department of Pathology, Affiliated Hospital of Weifang Medical University, Weifang, China
Empagliflozin is a novel type of sodium-glucose cotransporter two inhibitor with diverse beneficial effects in the treatment of nonalcoholic fatty liver disease (NAFLD). Although empagliflozin impacts NAFLD by regulating lipid metabolism, the underlying mechanism has not been fully elucidated. In this study, we investigated transcriptional regulation pathways affected by empagliflozin in a mouse model of NAFLD. In this study, NAFLD was established in male C57BL/6J mice by administration of a high-fat diet; it was then treated with empagliflozin and whole transcriptome analysis was conducted. Gene expression levels detected by transcriptome analysis were then verified by quantitative real-time polymerase chain reaction, protein levels detected by Western Blot. Differential expression genes screened from RNA-Seq data were enriched in lipid metabolism and synthesis. The Gene Set Enrichment Analysis (GSEA) results showed decreased lipid synthesis and improved lipid metabolism. Empagliflozin improved NAFLD through enhanced triglyceride transfer, triglyceride lipolysis and microsomal mitochondrial β-oxidation. This study provides new insights concerning the mechanisms by which sodium-glucose cotransporter two inhibitors impact NAFLD, particularly in terms of liver lipid metabolism. The lipid metabolism-related genes identified in this experiment provide robust evidence for further analyses of the mechanism by which empagliflozin impacts NAFLD.
Introduction
Nonalcoholic fatty liver disease (NAFLD) is a common chronic liver disease, characterized by a broad spectrum of clinical manifestations (e.g., hepatic steatosis, fibrosis, cirrhosis, and hepatocellular carcinoma) (Rinella, 2015; Manne et al., 2018). Importantly, NAFLD affects one-quarter of adults worldwide (Zhou et al., 2019). Metabolic diseases (e.g., obesity and type 2 diabetes) have been associated with greater NAFLD risk in previous studies (Perumpail et al., 2017). Considering the increasing morbidity of obesity and diabetes, the prevalences of NAFLD and advanced liver disease are expected to increase (Cotter and Rinella, 2020). Lipid metabolism disorders have essential roles in NAFLD occurrence and progression (Mato et al., 2019; Di et al., 2021). For example, hepatic steatosis (a typical symptom of NAFLD) results from lipid acquisition that exceeds lipid disposal (Ipsen et al., 2018). Considering that the liver has a central role in lipid synthesis and metabolism, it is important to study drug prevention and treatment of NAFLD from the perspective of liver lipid metabolism. Recently, SGLT2i has become a focus because of beneficial effect on reduced serum transaminase activities (Seko et al., 2018). Furthermore, several pilot studies showed SGLT2i had a significant reduction in liver transaminases, body weight, and the fatty liver index in NAFLD patients (Komiya et al., 2016; Seko et al., 2017). Therefore, the efficacy of SGLT2 inhibitor on NASH/NAFLD can be expected.
Sodium-glucose cotransporter 2 (SGLT2) inhibitors are effective hypoglycemic drugs that function by inhibiting glucose reabsorption in proximal renal tubules. SGLT2 inhibitors improve glycemic status in patients with type 2 diabetes; they also contribute to weight loss and body fat reduction (Szekeres et al., 2021). Multiple studies have reported that SGLT2 inhibitor treatments reduce serum cholesterol and triglyceride (TG) levels (Hayashi et al., 2017; Calapkulu et al., 2019); they also increase lipolysis, liver fatty acid oxidation, and ketogenesis (Osataphan et al., 2019). Empagliflozin (EMP), a novel anti-hyperglycemic agent (Supplementary Figure S1), is a type of SGLT2 inhibitor (Frampton, 2018). Although EMP was designed as a hypoglycemic drug, it also has an important effect on lipid metabolism and has been shown to reduce liver fat in NAFLD (Kuchay et al., 2018). Studies by Xu et al. revealed that EMP contributes to weight loss through increased fat utilization and browning. Furthermore, reductions of obesity-related inflammation by EMP have been associated with M2 macrophage polarization (Xu et al., 2017; Xu et al., 2019). Additionally, recent studies indicated SGLT2i decreased body weight and lipid profiles, including total cholesterol (TC) and TG (Al-Sharea et al., 2018; Matsubayashi et al., 2020; Liu et al., 2021). However, there is a few researches illustrating how EMP regulated the progress of TG accumulation and associated triglyceride molecules.
Triglyceride molecule is the main form of storage and transport of fatty acids in cells and plasma. Fatty acid activation is catalyzed by long-chain acyl-CoA synthetase (ACSL) family of enzymes and microsomal G3P acyltransferase (GPAT) enzymes catalyzed G3P pathway, which is the main way to synthesize TG (Cohen and Fisher, 2013). The triglyceride transfer protein (MTP) combined with TG to form VLDL particles (Cohen and Fisher, 2013). Furthermore, adipose triglyceride lipase (ATGL, also known as PNPLA2), a rate-limiting step of TG lipolysis in adipocytes was activated by the comparative gene identification-58 (CGI-58) (Lass et al., 2006). Liver specific resection of CGI-58 can lead to NAFLD phenotype in mice, including hepatic steatosis and hepatic fibrosis (Guo et al., 2013). The patatin-like phospholipase domain-containing 3 (PNPLA3) has a genetic association with NALFD, and mice showed accumulation of inactive PNPLA3 in lipid droplets and increased steatosis in high-fat diet feeding (Smagris et al., 2015). Hepatic TG can be oxidized by multiple pathways. Mitochondrial β-oxidation is the main route in hepatocytes (Musso et al., 2009). The expression of genes including mitochondrial β-oxidation is regulated largely by PPARα activity. PGC1α and BAF60a form a complex to regulate activity of PPARα transcriptional in the liver (Li et al., 2008).
Moreover, with the rapid development of RNAseq technology, transcriptome analysis of NAFLD is increasing. Previous investigations suggested lipid metabolism plays a major role in NAFLD (Xu et al., 2021). However, the molecular mechanisms underlying EMP-mediated NAFLD protection of hepatic lipid metabolism have not been fully elucidated. Here, we investigated signaling pathways impacted by EMP in a mouse model of NAFLD and explored the mechanisms and triglyceride molecules underlying the effects of EMP on liver lipid metabolism and synthesis.
Materials and Methods
Animal Model and Drug Administration
All animal experiments in this study were approved by the Animal Ethics Committee of Weifang Medical University and adhered to national guidelines for experiments involving laboratory animals. Six-week-old male C57BL/6J mice (weight, 21 ± 0.8 g) were purchased from Pengyue Co., Ltd. (Jinan, China). All mice were housed under a 12-h day/night cycle at constant temperature (approximately 22°C) and suitable humidity (45–70%). As previous study (Han et al., 2020), the CT group (n = 8) was fed a regular diet (320 kcal per 100 g) (Pengyue Co., Ltd. Jinan, China), and the HFD (n = 8) and EMP groups (n = 8) were fed a high-fat diet (54.05% fat, 15.10% protein, and 30.85% carbohydrate; 529.8 kcal per 100 g) (Fanbo Biotechnology Co., Ltd., Shanghai, China) for 20 weeks. After the 12 weeks of feeding, the EMP group underwent further intragastric administration of EMP (10 mg/kg/day) for 8 weeks, whereas CT and HFD groups received same volume of saline for 8 weeks. After successive drug administration for 8 weeks, the mice were given 2% barbiturate sodium (40 mg/kg, ip) after 12 h of fasting. Body composition of the mice were measured by Body Composition Analysis (Bruker Minispec LF50, German). Blood samples were collected. Livers were also rapidly isolated and cut into several small pieces; two pieces of liver were fixed with formalin, then subjected to hematoxylin-eosin (HE) staining and Oil Red O staining. Plasma triglycerides and liver triglycerides were measured using triglyceride assay kit from Jiancheng (Nanjing, China). The remaining liver tissues were stored at −80°C for RNA extraction.
Total RNA Extraction, cDNA Library Construction, and RNA Sequencing
Total RNA was extracted and purified from liver tissue using TRIzol, in accordance with the manufacturer’s instructions (Invitrogen, United States). Bioanalyzer 2,100 (Agilent, United States) and NanoDrop ND-1000 (NanoDrop, United States) equipment were used to assess total RNA quantity and quality. Poly (A) RNA was purified twice using Dynabeads Oligo (dT) (Thermo Fisher, United States). The resulting RNA was fragmented using a Magnesium RNA Fragmentation Module (NEB, United States) at 94°C for 5–7 min. RNA was used as a template for cDNA synthesis by SuperScript™ II Reverse transcriptase (Invitrogen). Single- or dual-index adapters were ligated to the fragments; AMPureXP beads were used for size selection. The mean insert size in the final cDNA library was 300 ± 50 bp. RNA-Seq was performed by Illumina NovaSeq™ 6,000 (LC Bio Technology Co., Ltd. China); 150-bp paired-end sequences were generated.
Quantification of RNA and Analysis of Differential Expression
FastQC was used to assess sequence quality. After removal of low-quality reads, the remaining 150-bp paired-end sequences were reassembled and mapped to the Mus musculus genome by HISAT2 (Kim et al., 2015). The transcripts were merged from sequences by StringTie (Pertea et al., 2015). To assure quality of RNA, we finally chosen five samples per group for further study. Genes with low expression levels were discarded when their minimum expression thresholds were below one count per million (CPM) in all 15 samples. The expression data were normalized by trimmed mean of M-values (TMM) using the edgR package in Bioconductor (Robinson et al., 2010); counts of sequences were transformed to log2(CPM) values for subsequent estimation of gene expression levels. Differentially expressed genes (DEGs) were identified using the limma package in R (Ritchie et al., 2015) with the following settings: fold change >2 or < −2, and p < 0.01.
Analysis of Enrichment, PPI and GSEA
Gene Ontology (GO) (Young et al., 2010) and Kyoto Encyclopedia of Genes and Genomes (KEGG) (Kanehisa et al., 2008) analyses of DEGs were conducted using the clusterProfiler package in R (Yu et al., 2012). DEGs were integrated with mouse protein-protein interaction (PPI) networks using STRING. PPI networks were analyzed by Cytoscape 3.8.0 (Shannon et al., 2003). Gene Set Enrichment Analysis (GSEA) was used for additional pathway analyses; forward regression was then performed (Subramanian et al., 2005).
Real-Time Quantitative PCR
DEGs related to lipid metabolism and synthesis were subjected to validation by qPCR. Total RNA was extracted as described above for RNA-Seq analysis. cDNA was obtained using the PrimerScript™ RT reagent Kit with gDNA Eraser (TaKaRa, Japan). qPCR was performed using TB Green® Premix Ex Taq™ (TaKaRa). Primers of target genes were synthesized by Sangon Biotech (China); Table 1 shows the sequences of primers used in this study. All relative expression levels of RNA were calculated based on Ct values, then normalized to the levels of β-actin (Liu et al., 2011) as follows: 2 [ΔCt = Ct (target genes) − Ct (β-actin)].
Western Blot
Total protein from liver tissues were lysed in ice-cold Lysis buffer (Whole Cell Lysis Assay, KeyGEN, China), and the total protein concentration was quantified by a BCA Protein Assay Kit (Solarbio, China). Equal amounts of protein were electrophoresed on 10% SDS-PAGE gels and transferred to PVDF membranes. After blocking with 5% BSA, the membranes were incubated with primary antibodies (FDPS, HMGCS1 1:1,000, Proteintech) overnight at 4°C. After washing, membranes were soaked in the secondary antibodies (Anti-rabbit IgG HRP-linked, 1:5,000, Anti-mouse IgG HRP-linked, 1:5,000) for 1 h β-actin was used as an internal control. Immunoblots were then visualized by enhanced chemiluminescence (ECL) using ImageQuant LAS 5000 (GE Healthcare) and analyzed using ImageJ.
Oil Red O Staining and H&E Staining
Mice were anesthetized by intraperitoneal injection of 2% sodium pentobarbital and 4% buffered neutral formalin solution was used to fix their liver tissues for 24 h. Tissue sections were subjected to H&E staining and Oil Red O staining, in accordance with standard protocols. The images were obtained using a microscope (CKX53, Olympus).
Statistics
Statistical analysis of this study was conducted by GraphPad Prism 7.0. Normality test was performed by Shapiro-Wilk normality test. Normally distributed data were presented as the means ± SEM. One-way ANOVA followed by Dunnett’s post-hoc test was performed when data involved in all three groups. Abnormally distributed data are displayed as the median and interquartile intervals, and the distributions between groups were compared with the non-parametric Kruskal-Wallis test. p values <0.05 were considered significant.
Results
Effects of EMP on Mouse Body Weight and Liver Pathology
All mice were acclimated for 1 week, then randomly separated into CT, HFD, and EMP groups. However, the CT + EMP group is not included in our study. Because our previous study had reported EMP does not have an effect on body weight, fat mass, plasma lipids (Sun et al., 2020). Also, a previous study consistent with our study indicated EMP does not alter liver metabolism and biochemical parameters (Petito-da-Silva et al., 2019). Therefore, we did not establish the CT + EMP group. After 20 weeks of high-fat diet treatment, the weight and fat mass were significantly higher in the HFD and EMP groups than in the CT group. After 8 weeks of intragastric administration of EMP, the weight and fat mass were decreased in the EMP group, compared with the HFD group (Figure 1A, B, Table 2). To access the liver metabolism caused by HFD, we measured alanine aminotransferase (ALT), aspartate aminotransferase (AST) and plasma TG levels. The results showed EMP alleviated increased liver injury with reducing AST, ALT and plasma TG levels (Table 2). Liver changes were observed by H&E and Oil Red O staining (Figures 1C,D). Liver tissue structure was normal in the CT group, with distinguishable edges and clear outlines. HFD mice showed prominent diffuse hepatic steatosis with nuclear condensation, cytoplasmic looseness, and increased lipid contents; all of these abnormalities were ameliorated by EMP treatment (Figure 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/793586/fphar-12-793586-HTML/image_m/fphar-12-793586-g001.jpg)
FIGURE 1. Effect of EMP on mice. The impact of EMP on body weight (A) and fat mass (B). Liver histological examination with H&E (200×) (C) and Oil Red O staining (D). Data are presented as means ± SEM, *p < 0.05, Scale bar: 30 μm.
Identification of Differentially Expressed Genes Across Disease and Drug
Here, we performed RNA-Seq using liver tissue from NAFLD model mice that had received EMP to explore the mechanisms underlying the effects of EMP on lipid metabolism and synthesis. We first investigated the liver tissue transcriptome profile in mice with regular and high-fat diets to determine gene expression changes during NAFLD. In total, 404 DEGs were observed in the HFD group, compared with the CT group; these consisted of 251 upregulated genes and 153 downregulated genes (Figures 2A,C). To understand how EMP treatment affects NAFLD, we compared transcriptome profiles between the HFD and EMP groups. Compared with the HFD group, 123 DEGs were observed in the EMP group; these consisted of 50 upregulated genes and 73 downregulated genes (Figures 2A,D). We investigated repeatability among samples using principal components analysis (PCA). The PCA results were displayed in a two-dimensional image format (Figure 2B), showing the degrees of separation among samples and groups. Notably, spatial separation trends were observed between HFD and CT groups, as well as between EMP and HFD groups. The findings indicated significant differences in gene expression patterns among the three groups; sample correlation was high within each group. The union and shared DEGs (492 and 35 genes) were displayed in the Venn diagram (Figure 3A). The gene expression profiles of union and shared DEGs exhibited apparent gene expression variations among the CT, HFD, and EMP groups (Figures 3B,C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/793586/fphar-12-793586-HTML/image_m/fphar-12-793586-g002.jpg)
FIGURE 2. Comparison of transcriptome profiling. (A) The gene number of up-regulation and down-regulation in HFD/CT and EMP/HFD. (B) PCA analysis of three groups. (C) The volcano plots for the distribution of differently expressed genes between HFD and CT group. Blue represents a down-regulation in expression, red represents upregulation and gray represents no significance when compared with control. (D) The volcano plots for the distribution of differently expressed genes between EMP and HFD group.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/793586/fphar-12-793586-HTML/image_m/fphar-12-793586-g003.jpg)
FIGURE 3. The heatmap representing alteration of hepatic union (A) and shared (B) genes expression of HFD/CT and EMP/HFD. (C) The Venn diagram shows an overlap among the HFD/CT and EMP/HFD genes.
GO and KEGG Pathway Analysis
The expression levels of 404 genes were significantly changed in the HFD group compared with the CT group. GO enrichment analysis showed that the DEGs affected by HFD were mainly involved in the “sterol biosynthetic process,” “sterol metabolic process,” “steroid metabolic process,” “cholesterol biosynthetic process,” “cholesterol metabolic process,” “secondary alcohol biosynthetic process,” and “steroid biosynthetic process” categories (Figure 4A). Moreover, KEGG pathway analysis showed that these DEGs were mainly enriched in “steroid biosynthesis” and “terpenoid backbone biosynthesis” pathways (Figure 4B). These results suggested that the HFD treatment primarily affected lipid synthesis and metabolism. Our results are consistent with previous studies which also revealed lipid biosynthetic and lipid metabolic process were enriched in the biological process (Chen et al., 2021; Xu et al., 2021). These results indicate a large proportion of dysregulated protein-coding genes in the transcriptome of a NAFLD liver, associated with glycolipid metabolism. However, there is no study exploring the effect of EMP on HFD-induced NAFLD.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/793586/fphar-12-793586-HTML/image_m/fphar-12-793586-g004.jpg)
FIGURE 4. Main GO terms and pathways enriched by DEGs between HFD and CT group, (A) GO terms; (B) KEGG pathways. Enrichment of GO and pathways based on DEGs between EMP and HFD group, (C) GO terms; (D) KEGG pathways. Gene Number represents the number of the differential expressed genes enriched in each term.
Thus, we established HFD-EMP group to further analysis the potential biological processes. We identified 123 DEGs that were significantly changed in the EMP group, compared with the HFD group. GO enrichment analysis showed that the DEGs were mainly involved in “response to redox state,” “isoprenoid biosynthetic process,” “isoprenoid metabolic process,” “terpenoid biosynthetic process,” “steroid metabolic process,” and “sterol biosynthetic process” categories (Figure 4C). The KEGG pathway analysis indicated that these DEGs were mainly enriched in “terpenoid backbone biosynthesis” and “JAK-STAT signaling” pathways (Figure 4D). These results indicated that EMP mainly affects lipid oxidation and metabolism.
PPI Network and GSEA Analysis
Gene expression data and PPI networks were used in combination to identify the regulated parts of the network in the drug treatment. This process created a PPI network involving DEGs of the EMP and HFD groups, such that the edges connecting co-expressed genes were preserved. To identify the “communities” and “hubs” of the subnetworks, we analyzed the differentially expressed parts of the transcriptome. This analysis revealed multiple communities that contained either up- or downregulated genes associated with EMP treatment (Figure 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/793586/fphar-12-793586-HTML/image_m/fphar-12-793586-g005.jpg)
FIGURE 5. Network of protein-protein interaction (PPI) of DEGs between EMP and HFD group (orange node: upregulated genes, blue node: downregulated genes).
Next, we used GSEA to identify the differentially expressed pathways associated with EMP treatment. Although PPI network analysis identified co-regulated gene communities, GSEA showed additional pathways that were differentially associated with histological severity. In total, 12 pathways were affected by EMP treatment (FDR <0.05, Figure 6A). The full set of pathways identified in this experiment is provided in Table 3. Four pathways were associated with lipid metabolism and synthesis: “ovarian steroidogenesis” (upregulated), “arachidonic acid metabolism” (upregulated), “terpenoid backbone biosynthesis” (downregulated), and “linoleic acid metabolism” (upregulated) (Figures 6B–E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/793586/fphar-12-793586-HTML/image_m/fphar-12-793586-g006.jpg)
FIGURE 6. Metabolic pathways affected by EMP. (A) The top Reactome gene sets that are up- and down-regulated with respect to EMP treatment based on normalized enrichment scores (NES) for EMP/HFD genes. (B–E) Pathways related to lipid metabolism and synthesis which are selected from (A). (orange: upregulated pathways, blue: downregulated pathways, FDR ≤0.01).
EMP Reduced Lipid Metabolism Through Increasing TG Lipolysis and Mitochondrial β-oxidation
Consistent with previous studies, HFD mice showed significantly high TG content both in serum and liver (Figure 7A/B, Table 2), however, there was a decrease tend between HFD mice and EMP mice but no significant difference (p > 0.05). To further explore the underlying triglyceride molecules among groups, we detect the genes expressions of triglyceride metabolism and transport. As showed in Figure Figure7C, EMP had no effect on these genes (ACSL3, ACSL5, GPAT1) involved in the synthesis of TG. Compared with control group, HFD mice had lower expression of MTP, CGI-58, HSL, PNPLA3, meaning decreased triglyceride transfer and TG lipolysis (Figures 7D,E). However, these effects were reversed by EMP treatment. Moreover, EMP also enhanced these decreased-mitochondrial β-oxidation factors (including PPARα, PGC1α), which were overtly suppressed by HFD (Figure7F). To further explore whether RNA-seq data was consistent with protein level, we choose the most relevant genes in lipid metabolism to detect protein level through Western Blot analysis. As expected, the protein level of FDPS and HMGCS1was downregulated in mice after treatment of EMP (Figures 7G,H).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/793586/fphar-12-793586-HTML/image_m/fphar-12-793586-g007.jpg)
FIGURE 7. Empagliflozin reduced triglyceride level both in serum (A) and liver (B) through enhanced triglyceride transfer (D), lipolysis (E) and microsomal mitochondrial β-oxidation (F). Empagliflozin had no effect on the synthesis genes of triglyceride (C). Western-blot analysis showed the protein levels of FDPS and HMGCS1 (G and H), which were the most relevant genes in lipid metabolism. Equal loading of protein was verified by probing β-actin. Data represent means ± SEM. *p < 0.05, **p < 0.01, n = (4–6).
qPCR Assessment of DEGs
To confirm the gene expression levels observed in RNA-Seq analysis, we chose nine genes that were mainly associated with lipid metabolism and synthesis (Smpd3, Fdps, Hmgcs1, Rdh1, Sox9, HR, Sqle, Fam131c, Cyp2g1). qPCR results were generally consistent with the RNA-Seq findings (Figure 8A). As showed in Figure 8B, the mRNA expression of Fam131c, Sox9, Smpd3 was upregulated in HFD groups compared with control group, these effects were reversed after treatment of EMP. Moreover, we found the expression of Fdps, Hmgcs1, HR, Sqle were downregulated in HFD-induced mice, however, these results were not reversed in EMP group.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/793586/fphar-12-793586-HTML/image_m/fphar-12-793586-g008.jpg)
FIGURE 8. Quantitative real-time PCR validation of RNA-seq analysis of DEGs. Graphs represent log2(CPM) values of gene expression levels (A-C) in RNA-seq data and mRNA levels normalized to β-actin, *p < 0.05 vs. CT, #p < 0.05 vs. HFD, n = 5, (D-F) from quantitative real-time PCR data. Data are presented as means ± SEM, *p < 0.05 vs. CT, #p < 0.05 vs. HFD, n = (4–6).
Discussion
This study investigated the effects of EMP on NAFLD. We found that EMP attenuated liver fat content and improved lipid metabolism by regulating genes related to lipid generation, transfer, lipolysis and oxidation. Our results provide a snapshot of the transcriptionally regulated pathways affected by EMP treatment for NAFLD; they may be important candidate pathways for alleviating disease. Additionally, the results provide insights concerning how these pathways interact to inhibit disease.
EMP is a novel type of antihyperglycemic drug that inhibits glucose reabsorption in proximal renal tubules. In addition to its effects on glycemic status in patients with type 2 diabetes, it contributes to weight loss and body fat reduction (Sun et al., 2020; Szekeres et al., 2021). Recently, clinical trials have reported SGLT-2i had beneficial effects on NAFLD, as evidenced by the remarkably reduced ALT, AST, triglycerides, hepatic insulin sensitivity indices (Ranjbar et al., 2019; Kahl et al., 2020). Of note, in EMPA-REG OUTCOME trial and E-LIFT trial also showed that treatment with empagliflozin significantly reduced liver enzymes and liver fat (Sattar et al., 2018; Kim and Lee, 2020). It also increases energy consumption, thermogenesis, and expression of uncoupling protein one in brown fat (Xu et al., 2017). Importantly, EMP attenuates NAFLD by activating autophagy and reducing ER stress and apoptosis (Nasiri-Ansari et al., 2021). Thus, these effects have emerged as important underlying mechanisms in NAFLD development and progression. In the present study, we found that 8-weeks treatment with EMP reduced body weight and fat mass, while alleviating HFD-induced liver injury. These results confirm the protective effects of EMP in a model of NAFLD.
The liver is essential for lipid synthesis and regulation of lipid metabolism. Healthy liver function involves balanced adjustment of lipid metabolism through multiple biological processes. In NAFLD, the compensatory enhancement of fatty acid oxidation is insufficient to normalize the lipid level; it may promote cell and tissue damage, as well as disease progression, by enhancing oxidative stress (Ipsen et al., 2018). In our study, we found that EMP reduced triglyceride level both in serum and liver through enhanced triglyceride transfer, TG lipolysis and microsomal mitochondrial β-oxidation. Previous studies identified DEGs that were associated with NAFLD onset and progression. Moreover, some studies examined gene expression using transcriptome data from liver biopsies of affected patients (Wruck et al., 2015). The findings revealed that many functional pathways and genes exhibited significant associations with NAFLD; many of these pathways and genes were related to lipid metabolism. Here, we used RNA-Seq to investigate the major pathways and genes involved in the effects of EMP on lipid synthesis and metabolism. Subsequent bioinformatics assays identified 123 DEGs with significant changes in expression levels after EMP treatment.
A previous study confirmed that overexpression of sex-determining region Y-box 9 (SOX9) enhanced lipogenesis and expression of PPARγ in sebocytes (Shi et al., 2017). In addition, HR lysine demethylase and nuclear receptor corepressor (HR) may be a adipogenic transcription factor by regulating the expression of PPARγ to generate white adipogenesis (Kumpf et al., 2012). PPARγ is the primary inducer of adipogenesis (Li et al., 2018) and regulates adipocyte differentiation. In our study, SOX9 mRNA levels were significantly reduced in the EMP group. We speculate EMP treatment inhibited SOX9 expression, which led to reduced lipogenesis in the liver. These changes reduced PPARγ expression, thereby inhibiting adipocyte differentiation. Sphingomyelins (SMs) have been reported as potential biomarkers for diagnosis of hepatic steatosis (Li et al., 2014). Additionally, a transcriptome analysis involved in mouse models of NAFLD and liver tissues from patients had reported the expression of SOX9 and SMPD3 (Teufel et al., 2016). In this study, sphingomyelin phosphodiesterase 3 (SMPD3) was inhibited in the EMP group, which may have led to increased SM expression. This finding conflicts with the results of a previous study (Zhang et al., 2019), in which SMPD3 had a lower expression level in tyloxapol-treated mice, then increased after treatment. This inconsistency is presumably because the animal models were established using different methods. The final biosynthetic product of the 3-hydroxy-3-methylglutaryl-CoA synthase 1 (HMGCS1) pathway is cholesterol. HMGCS1 is upregulated by galectin-7 (Gal-7); subsequent HMGCS1 activity increases the accumulation of cellular cholesterol (Fujimoto et al., 2021). Squalene monooxygenase (SQLE) is a critical regulatory point in the cholesterol synthesis pathway (Howe et al., 2017). Farnesyl diphosphate synthase (FDPS) and cytochrome P450 family 2g1 (CYP2G1) are both associated with lipid and cholesterol biosynthesis (Wang et al., 2017). In the present study, the expression levels of HMGCS1, SQLE, and FDPS were highest in the CT group and lowest in the EMP group. These gene expression differences provide important insights into the mechanisms by which EMP affects the liver. Family with sequence similarity 131c (FAM131C) is a group of human N-myristoylated proteins; protein N-myristoylation is required for proper targeting of SAMM50 to mitochondria (Takamitsu et al., 2015). Some studies have shown that overexpression of SAMM50 enhances fatty acid oxidation and reduces intracellular lipid accumulation (Li et al., 2021). In contrast, we found that the Fam131c expression level increased in the HFD group but decreased in the EMP group. Retinol dehydrogenase 1 (RDH1) is one of several enzymes that catalyze the conversion of retinol into all-trans-retinoic acid (atRA). RDH1 can suppress adiposity by promoting brown adipose generation (Krois et al., 2019). In this study, RDH1 expression levels differed between RNA-Seq and qPCR, presumably because of differences in the technical protocols.
Overall, we found that EMP treatment could reverse the activities of NAFLD-associated pathways. These transformations mainly involved reduction of lipid metabolism and synthesis pathways in the liver. Our results indicate that EMP treatment alleviates fat accumulation and slightly affects inflammation in NAFLD. GO and KEGG enrichment were performed to analyze these DEGs, with the aim of clarifying the mechanism by which EMP affects lipid synthesis and metabolism in the liver. Our findings showed that EMP mainly affects the response to redox state, triglyceride transfer, TG lipolysis and microsomal mitochondrial β-oxidation, and JAK-STAT signaling pathways. Importantly, EMP reduced lipid accumulation and alleviated pathological changes involved in NAFLD. These changes may be related to the regulation of lipid oxidation-associated gene expression by EMP. However, we only established the relation between mechanisms of EMP and NAFLD, but which cell type (hepatocytes, macrophages or hepatic stellate cells) in liver for lipid metabolism and synthesis is not determined in this work. Therefore, we will focus this limitation in the next research.
In summary, our results provide snapshot of the transcriptionally regulated pathways affected by EMP treatment for NAFLD; these may be important candidate pathways for treatment to slow the progression of NAFLD. Additionally, our findings provide insights concerning how these pathways interact to inhibit disease. Our research provides new insights concerning the mechanisms by which SGLT inhibitors impact NAFLD, particularly in terms of liver metabolism. The genes identified in this experiment provide robust evidence for further analyses of the mechanism by which EMP impacts 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/, PRJNA770493.
Ethics Statement
The animal study was reviewed and approved by Animal Ethics Committee of Weifang Medical University.
Author Contributions
YM, CK and HQ collected the data, conducted the analysis, and drafted the manuscript. XBC designed the study and conducted the analysis. YL, FH and NH directed the study’s analytic strategy and reviewed the manuscript. XS and JS designed the entire study and revised the manuscript. All authors read and approved the final manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (81870593), Natural Science Foundation of Shandong Province of China (ZR2020MH106, ZR202102240146), Shandong Province Higher Educational Science and Technology Program for Youth Innovation (2020KJL004) and Quality Improvement of Postgraduate Education in Shandong Province (SDYAL19156).
Conflict of Interest
The 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/fphar.2021.793586/full#supplementary-material
References
Al-Sharea, A., Murphy, A. J., Huggins, L. A., Hu, Y., Goldberg, I. J., and Nagareddy, P. R. (2018). SGLT2 Inhibition Reduces Atherosclerosis by Enhancing Lipoprotein Clearance in Ldlr-/- Type 1 Diabetic Mice. Atherosclerosis 271, 166–176. doi:10.1016/j.atherosclerosis.2018.02.028
Calapkulu, M., Cander, S., Gul, O. O., and Ersoy, C. (2019). Lipid Profile in Type 2 Diabetic Patients with New Dapagliflozin Treatment; Actual Clinical Experience Data of Six Months Retrospective Lipid Profile from Single center. Diabetes Metab. Syndr. 13, 1031–1034. doi:10.1016/j.dsx.2019.01.016
Chen, X., Tang, Y., Chen, S., Ling, W., and Wang, Q. (2021). IGFBP-2 as a Biomarker in NAFLD Improves Hepatic Steatosis: an Integrated Bioinformatics and Experimental Study. Endocr. Connect. 10, 1315–1325. doi:10.1530/EC-21-0353
Cohen, D. E., and Fisher, E. A. (2013). Lipoprotein Metabolism, Dyslipidemia, and Nonalcoholic Fatty Liver Disease. Semin. Liver Dis. 33, 380–388. doi:10.1055/s-0033-1358519
Cotter, T. G., and Rinella, M. (2020). Nonalcoholic Fatty Liver Disease 2020: The State of the Disease. Gastroenterology 158, 1851–1864. doi:10.1053/j.gastro.2020.01.052
Di Mauro, S., Salomone, F., Scamporrino, A., Filippello, A., Morisco, F., Guido, M., et al. (2021). Coffee Restores Expression of lncRNAs Involved in Steatosis and Fibrosis in a Mouse Model of NAFLD. Nutrients 13, 2952. doi:10.3390/nu13092952
Frampton, J. E. (2018). Empagliflozin: A Review in Type 2 Diabetes. Drugs 78, 1037–1048. doi:10.1007/s40265-018-0937-z
Fujimoto, N., Akiyama, M., Satoh, Y., and Tajima, S. (2021). Interaction of Galectin-7 with HMGCS1 In Vitro May Facilitate Cholesterol Deposition in Cultured Keratinocytes. J. Invest. Dermatol. S0022-202X, 02085–02086. doi:10.1016/j.jid.2021.04.038
Guo, F., Ma, Y., Kadegowda, A. K., Betters, J. L., Xie, P., Liu, G., et al. (2013). Deficiency of Liver Comparative Gene Identification-58 Causes Steatohepatitis and Fibrosis in Mice. J. Lipid Res. 54, 2109–2120. doi:10.1194/jlr.M035519
Han, X., Ding, C., Zhang, G., Pan, R., Liu, Y., Huang, N., et al. (2020). Liraglutide Ameliorates Obesity-Related Nonalcoholic Fatty Liver Disease by Regulating Sestrin2-Mediated Nrf2/HO-1 Pathway. Biochem. Biophys. Res. Commun. 525, 895–901. doi:10.1016/j.bbrc.2020.03.032
Hayashi, T., Fukui, T., Nakanishi, N., Yamamoto, S., Tomoyasu, M., Osamura, A., et al. (2017). Dapagliflozin Decreases Small Dense Low-Density Lipoprotein-Cholesterol and Increases High-Density Lipoprotein 2-cholesterol in Patients with Type 2 Diabetes: Comparison with Sitagliptin. Cardiovasc. Diabetol. 16, 8. doi:10.1186/s12933-016-0491-5
Howe, V., Sharpe, L. J., Prabhu, A. V., and Brown, A. J. (2017). New Insights into Cellular Cholesterol Acquisition: Promoter Analysis of Human HMGCR and SQLE, Two Key Control Enzymes in Cholesterol Synthesis. Biochim. Biophys. Acta Mol. Cel Biol Lipids 1862, 647–657. doi:10.1016/j.bbalip.2017.03.009
Ipsen, D. H., Lykkesfeldt, J., and Tveden-Nyborg, P. (2018). Molecular Mechanisms of Hepatic Lipid Accumulation in Non-alcoholic Fatty Liver Disease. Cell. Mol. Life Sci. 75, 3313–3327. doi:10.1007/s00018-018-2860-6
Kahl, S., Gancheva, S., Straßburger, K., Herder, C., Machann, J., Katsuyama, H., et al. (2020). Empagliflozin Effectively Lowers Liver Fat Content in Well-Controlled Type 2 Diabetes: A Randomized, Double-Blind, Phase 4, Placebo-Controlled Trial. Diabetes Care 43, 298–305. doi:10.2337/dc19-0641
Kanehisa, M., Araki, M., Goto, S., Hattori, M., Hirakawa, M., Itoh, M., et al. (2008). KEGG for Linking Genomes to Life and the Environment. Nucleic Acids Res. 36, D480–D484. doi:10.1093/nar/gkm882
Kim, D., Langmead, B., and Salzberg, S. L. (2015). HISAT: a Fast Spliced Aligner with Low Memory Requirements. Nat. Methods 12, 357–360. doi:10.1038/nmeth.3317
Kim, K. S., and Lee, B. W. (2020). Beneficial Effect of Anti-diabetic Drugs for Nonalcoholic Fatty Liver Disease. Clin. Mol. Hepatol. 26, 430–443. doi:10.3350/cmh.2020.0137
Komiya, C., Tsuchiya, K., Shiba, K., Miyachi, Y., Furuke, S., Shimazu, N., et al. (2016). Ipragliflozin Improves Hepatic Steatosis in Obese Mice and Liver Dysfunction in Type 2 Diabetic Patients Irrespective of Body Weight Reduction. PLoS ONE 11, e0151511. doi:10.1371/journal.pone.0151511
Krois, C. R., Vuckovic, M. G., Huang, P., Zaversnik, C., Liu, C. S., Gibson, C. E., et al. (2019). RDH1 Suppresses Adiposity by Promoting Brown Adipose Adaptation to Fasting and Re-feeding. Cel. Mol. Life Sci. 76, 2425–2447. doi:10.1007/s00018-019-03046-z
Kuchay, M. S., Krishan, S., Mishra, S. K., Farooqui, K. J., Singh, M. K., Wasir, J. S., et al. (2018). Effect of Empagliflozin on Liver Fat in Patients with Type 2 Diabetes and Nonalcoholic Fatty Liver Disease: A Randomized Controlled Trial (E-LIFT Trial). Diabetes Care 41, 1801–1808. doi:10.2337/dc18-0165
Kumpf, S., Mihlan, M., Goginashvili, A., Grandl, G., Gehart, H., Godel, A., et al. (2012). Hairless Promotes PPARγ Expression and Is Required for white Adipogenesis. EMBO Rep. 13, 1012–1020. doi:10.1038/embor.2012.133
Lass, A., Zimmermann, R., Haemmerle, G., Riederer, M., Schoiswohl, G., Schweiger, M., et al. (2006). Adipose Triglyceride Lipase-Mediated Lipolysis of Cellular Fat Stores Is Activated by CGI-58 and Defective in Chanarin-Dorfman Syndrome. Cell Metab 3, 309–319. doi:10.1016/j.cmet.2006.03.005
Li, J. F., Qu, F., Zheng, S. J., Wu, H. L., Liu, M., Liu, S., et al. (2014). Elevated Plasma Sphingomyelin (D18:1/22:0) Is Closely Related to Hepatic Steatosis in Patients with Chronic Hepatitis C Virus Infection. Eur. J. Clin. Microbiol. Infect. Dis. 33, 1725–1732. doi:10.1007/s10096-014-2123-x
Li, S., Liu, C., Li, N., Hao, T., Han, T., Hill, D. E., et al. (2008). Genome-wide Coactivation Analysis of PGC-1alpha Identifies BAF60a as a Regulator of Hepatic Lipid Metabolism. Cel Metab 8, 105–117. doi:10.1016/j.cmet.2008.06.013
Li, Y., Jin, D., Xie, W., Wen, L., Chen, W., Xu, J., et al. (2018). PPAR-γ and Wnt Regulate the Differentiation of MSCs into Adipocytes and Osteoblasts Respectively. Curr. Stem Cel Res Ther 13, 185–192. doi:10.2174/1574888X12666171012141908
Li, Z., Shen, W., Wu, G., Qin, C., Zhang, Y., Wang, Y., et al. (2021). The Role of SAMM50 in Non-alcoholic Fatty Liver Disease: from Genetics to Mechanisms. FEBS Open Bio 11, 1893–1906. doi:10.1002/2211-5463.13146
Liu, W., Xu, G., Ma, J., Jia, W., Li, J., Chen, K., et al. (2011). Osteopontin as a Key Mediator for Vasculogenic Mimicry in Hepatocellular Carcinoma. Tohoku J. Exp. Med. 224, 29–39. doi:10.1620/tjem.224.29
Liu, Y., Xu, J., Wu, M., Xu, B., and Kang, L. (2021). Empagliflozin Protects against Atherosclerosis Progression by Modulating Lipid Profiles and Sympathetic Activity. Lipids Health Dis. 20, 5. doi:10.1186/s12944-021-01430-y
Manne, V., Handa, P., and Kowdley, K. V. (2018). Pathophysiology of Nonalcoholic Fatty Liver Disease/Nonalcoholic Steatohepatitis. Clin. Liver Dis. 22, 23–37. doi:10.1016/j.cld.2017.08.007
Mato, J. M., Alonso, C., Noureddin, M., and Lu, S. C. (2019). Biomarkers and Subtypes of Deranged Lipid Metabolism in Non-alcoholic Fatty Liver Disease. World J. Gastroenterol. 25, 3009–3020. doi:10.3748/wjg.v25.i24.3009
Matsubayashi, Y., Yoshida, A., Suganami, H., Osawa, T., Furukawa, K., Suzuki, H., et al. (2020). Association of Increased Hepatic Insulin Clearance and Change in Serum Triglycerides or β-hydroxybutyrate Concentration via the Sodium/glucose-Cotransporter 2 Inhibitor Tofogliflozin. Diabetes Obes. Metab. 22, 947–956. doi:10.1111/dom.13980
Musso, G., Gambino, R., and Cassader, M. (2009). Recent Insights into Hepatic Lipid Metabolism in Non-alcoholic Fatty Liver Disease (NAFLD). Prog. Lipid Res. 48, 1–26. doi:10.1016/j.plipres.2008.08.001
Nasiri-Ansari, N., Nikolopoulou, C., Papoutsi, K., Kyrou, I., Mantzoros, C. S., Kyriakopoulos, G., et al. (2021). Empagliflozin Attenuates Non-alcoholic Fatty Liver Disease (NAFLD) in High Fat Diet Fed ApoE(-/-) Mice by Activating Autophagy and Reducing ER Stress and Apoptosis. Int. J. Mol. Sci. 22, 818. doi:10.3390/ijms22020818
Osataphan, S., Macchi, C., Singhal, G., Chimene-Weiss, J., Sales, V., Kozuka, C., et al. (2019). SGLT2 Inhibition Reprograms Systemic Metabolism via FGF21-dependent and -independent Mechanisms. JCI Insight 4, e123130. doi:10.1172/jci.insight.123130
Pertea, M., Pertea, G. M., Antonescu, C. M., Chang, T. C., Mendell, J. T., and Salzberg, S. L. (2015). StringTie Enables Improved Reconstruction of a Transcriptome from RNA-Seq Reads. Nat. Biotechnol. 33, 290–295. doi:10.1038/nbt.3122
Perumpail, B. J., Khan, M. A., Yoo, E. R., Cholankeril, G., Kim, D., and Ahmed, A. (2017). Clinical Epidemiology and Disease burden of Nonalcoholic Fatty Liver Disease. World J. Gastroenterol. 23, 8263–8276. doi:10.3748/wjg.v23.i47.8263
Petito-da-Silva, T. I., Souza-Mello, V., and Barbosa-da-Silva, S. (2019). Empaglifozin Mitigates NAFLD in High-Fat-Fed Mice by Alleviating Insulin Resistance, Lipogenesis and ER Stress. Mol. Cel. Endocrinol. 498, 110539. doi:10.1016/j.mce.2019.110539
Ranjbar, G., Mikhailidis, D. P., and Sahebkar, A. (2019). Effects of Newer Antidiabetic Drugs on Nonalcoholic Fatty Liver and Steatohepatitis: Think Out of the Box!. Metabolism 101, 154001. doi:10.1016/j.metabol.2019.154001
Rinella, M. E. (2015). Nonalcoholic Fatty Liver Disease: a Systematic Review. JAMA 313, 2263–2273. doi:10.1001/jama.2015.5370
Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). Limma powers Differential Expression Analyses for RNA-Sequencing and Microarray Studies. Nucleic Acids Res. 43, e47. doi:10.1093/nar/gkv007
Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edgeR: a Bioconductor Package for Differential Expression Analysis of Digital Gene Expression Data. Bioinformatics 26, 139–140. doi:10.1093/bioinformatics/btp616
Sattar, N., Fitchett, D., Hantel, S., George, J. T., and Zinman, B. (2018). Empagliflozin Is Associated with Improvements in Liver Enzymes Potentially Consistent with Reductions in Liver Fat: Results from Randomised Trials Including the EMPA-REG OUTCOME® Trial. Diabetologia 61, 2155–2163. doi:10.1007/s00125-018-4702-3
Seko, Y., Sumida, Y., Sasaki, K., Itoh, Y., Iijima, H., Hashimoto, T., et al. (2018). Effects of Canagliflozin, an SGLT2 Inhibitor, on Hepatic Function in Japanese Patients with Type 2 Diabetes Mellitus: Pooled and Subgroup Analyses of Clinical Trials. J. Gastroenterol. 53, 140–151. doi:10.1007/s00535-017-1364-8
Seko, Y., Sumida, Y., Tanaka, S., Mori, K., Taketani, H., Ishiba, H., et al. (2017). Effect of Sodium Glucose Cotransporter 2 Inhibitor on Liver Function Tests in Japanese Patients with Non-alcoholic Fatty Liver Disease and Type 2 Diabetes Mellitus. Hepatol. Res. 47, 1072–1078. doi:10.1111/hepr.12834
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 13, 2498–2504. doi:10.1101/gr.1239303
Shi, G., Wang, T. T., Quan, J. H., Li, S. J., Zhang, M. F., Liao, P. Y., et al. (2017). Sox9 Facilitates Proliferation, Differentiation and Lipogenesis in Primary Cultured Human Sebocytes. J. Dermatol. Sci. 85, 44–50. doi:10.1016/j.jdermsci.2016.10.005
Smagris, E., BasuRay, S., Li, J., Huang, Y., Lai, K. M., Gromada, J., et al. (2015). Pnpla3I148M Knockin Mice Accumulate PNPLA3 on Lipid Droplets and Develop Hepatic Steatosis. Hepatology 61, 108–118. doi:10.1002/hep.27242
Subramanian, A., Tamayo, P., Mootha, V. K., Mukherjee, S., Ebert, B. L., Gillette, M. A., et al. (2005). Gene Set Enrichment Analysis: a Knowledge-Based Approach for Interpreting Genome-wide Expression Profiles. Proc. Natl. Acad. Sci. U S A. 102, 15545–15550. doi:10.1073/pnas.0506580102
Sun, X., Han, F., Lu, Q., Li, X., Ren, D., Zhang, J., et al. (2020). Empagliflozin Ameliorates Obesity-Related Cardiac Dysfunction by Regulating Sestrin2-Mediated AMPK-mTOR Signaling and Redox Homeostasis in High-Fat Diet-Induced Obese Mice. Diabetes 69, 1292–1305. doi:10.2337/db19-0991
Szekeres, Z., Toth, K., and Szabados, E. (2021). The Effects of SGLT2 Inhibitors on Lipid Metabolism. Metabolites 11, 87. doi:10.3390/metabo11020087
Takamitsu, E., Otsuka, M., Haebara, T., Yano, M., Matsuzaki, K., Kobuchi, H., et al. (2015). Identification of Human N-Myristoylated Proteins from Human Complementary DNA Resources by Cell-free and Cellular Metabolic Labeling Analyses. PLoS ONE 10, e0136360. doi:10.1371/journal.pone.0136360
Teufel, A., Itzel, T., Erhart, W., Brosch, M., Wang, X. Y., Kim, Y. O., et al. (2016). Comparison of Gene Expression Patterns between Mouse Models of Nonalcoholic Fatty Liver Disease and Liver Tissues from Patients. Gastroenterology 151, 513–e0. doi:10.1053/j.gastro.2016.05.051
Wang, C. C., Yen, J. H., Cheng, Y. C., Lin, C. Y., Hsieh, C. T., Gau, R. J., et al. (2017). Polygala Tenuifolia Extract Inhibits Lipid Accumulation in 3T3-L1 Adipocytes and High-Fat Diet-Induced Obese Mouse Model and Affects Hepatic Transcriptome and Gut Microbiota Profiles. Food Nutr. Res. 61, 1379861. doi:10.1080/16546628.2017.1379861
Wruck, W., Kashofer, K., Rehman, S., Daskalaki, A., Berg, D., Gralka, E., et al. (2015). Multi-omic Profiles of Human Non-alcoholic Fatty Liver Disease Tissue Highlight Heterogenic Phenotypes. Sci. Data 2, 150068. doi:10.1038/sdata.2015.68
Xu, L., Nagata, N., Chen, G., Nagashimada, M., Zhuge, F., Ni, Y., et al. (2019). Empagliflozin Reverses Obesity and Insulin Resistance through Fat browning and Alternative Macrophage Activation in Mice Fed a High-Fat Diet. BMJ Open Diabetes Res. Care 7, e000783. doi:10.1136/bmjdrc-2019-000783
Xu, L., Nagata, N., Nagashimada, M., Zhuge, F., Ni, Y., Chen, G., et al. (2017). SGLT2 Inhibition by Empagliflozin Promotes Fat Utilization and Browning and Attenuates Inflammation and Insulin Resistance by Polarizing M2 Macrophages in Diet-Induced Obese Mice. EBioMedicine 20, 137–149. doi:10.1016/j.ebiom.2017.05.028
Xu, L., Yin, L., Qi, Y., Tan, X., Gao, M., and Peng, J. (2021). 3D Disorganization and Rearrangement of Genome Provide Insights into Pathogenesis of NAFLD by Integrated Hi-C, Nanopore, and RNA Sequencing. Acta Pharm. Sin B 11, 3150–3164. doi:10.1016/j.apsb.2021.03.022
Young, M. D., Wakefield, M. J., Smyth, G. K., and Oshlack, A. (2010). Gene Ontology Analysis for RNA-Seq: Accounting for Selection Bias. Genome Biol. 11, R14. doi:10.1186/gb-2010-11-2-r14
Yu, G., Wang, L. G., Han, Y., and He, Q. Y. (2012). clusterProfiler: an R Package for Comparing Biological Themes Among Gene Clusters. OMICS 16, 284–287. doi:10.1089/omi.2011.0118
Zhang, T., Zhao, Q., Xiao, X., Yang, R., Hu, D., Zhu, X., et al. (2019). Modulation of Lipid Metabolism by Celastrol. J. Proteome Res. 18, 1133–1144. doi:10.1021/acs.jproteome.8b00797
Keywords: empaglifiozin, NAFLD, obesity, lipid, transcriptome
Citation: Ma Y, Kan C, Qiu H, Liu Y, Hou N, Han F, Shi J and Sun X (2021) Transcriptomic Analysis Reveals the Protective Effects of Empagliflozin on Lipid Metabolism in Nonalcoholic Fatty Liver Disease. Front. Pharmacol. 12:793586. doi: 10.3389/fphar.2021.793586
Received: 12 October 2021; Accepted: 06 December 2021;
Published: 21 December 2021.
Edited by:
Francisco Javier Cubero, Complutense University of Madrid, SpainReviewed by:
Stefania Di Mauro, University of Catania, ItalyCarlos Sanz Garcia, Complutense University of Madrid, Spain
Copyright © 2021 Ma, Kan, Qiu, Liu, Hou, Han, Shi and Sun. 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: Junfeng Shi, jfshi@wfmc.edu.cn; Xiaodong Sun, xiaodong.sun@wfmc.edu.cn
†These authors have contributed equally to this work and share first authorship