Breviscapine alleviates NASH by inhibiting TGF-β-activated kinase 1-dependent signaling.

BACKGROUND AND AIMS: NAFLD is a key component of metabolic syndrome, ranging from nonalcoholic fatty liver to NASH, and is now becoming the leading cause of cirrhosis and HCC worldwide. However, due to the complex and unclear pathophysiological mechanism, there are no specific approved agents for treating NASH. Breviscapine, a natural flavonoid prescription drug isolated from the traditional Chinese herb Erigeron breviscapus, exhibits a wide range of pharmacological properties, including effects on metabolism. However, the anti-NASH efficacy and mechanisms of breviscapine have not yet been characterized. APPROACH AND
RESULTS: We evaluated the effects of breviscapine on the development of hepatic steatosis, inflammation, and fibrosis in vivo and in vitro under metabolic stress. Breviscapine treatment significantly reduced lipid accumulation, inflammatory cell infiltration, liver injury, and fibrosis in mice fed a high-fat diet, a high-fat/high-cholesterol diet, or a methionine- and choline-deficient diet. In addition, breviscapine attenuated lipid accumulation, inflammation, and lipotoxicity in hepatocytes undergoing metabolic stress. RNA-sequencing and multiomics analyses further indicated that the key mechanism linking the anti-NASH effects of breviscapine was inhibition of TGF-β-activated kinase 1 (TAK1) phosphorylation and the subsequent mitogen-activated protein kinase signaling cascade. Treatment with the TAK1 inhibitor 5Z-7-oxozeaenol abrogated breviscapine-mediated hepatoprotection under metabolic stress. Molecular docking illustrated that breviscapine directly bound to TAK1.
CONCLUSION: Breviscapine prevents metabolic stress-induced NASH progression through direct inhibition of TAK1 signaling. Breviscapine might be a therapeutic candidate for the treatment of NASH.

acetyl CoA carboxylase alpha

actin beta

alanine aminotransferase

aspartate aminotransferase

C‐C motif chemokine ligand 2/5

collagen type I/III alpha 1

carnitine palmitoyltransferase 1a, liver

connective tissue growth factor

C‐X‐C motif chemokine ligand 2/10

fatty acid binding protein 1

fatty acid synthase

glyceraldehyde 3‐phosphate dehydrogenase

gene set enrichment analysis

gene set variation analysis

hematoxylin and eosin

high‐fat diet

high‐fat/high‐cholesterol

inhibitor of kappa light polypeptide gene enhancer in B cells alpha

inhibitor of NF‐κB kinase subunit beta

c‐Jun N‐terminal kinase

Kyoto Encyclopedia of Genes and Genomes

liver weight

liver weight‐to‐body weight ratio

mitogen‐activated protein kinase

methionine‐ and choline‐deficient

methionine‐ and choline‐sufficient

normal chow

oleic acid

phosphorylated

palmitic acid

principal component analysis

PA + OA

peroxisome proliferator‐activated receptor α/γ

protein–protein interaction

picrosirius red

RNA sequencing

stearoyl‐CoA desaturase 1

TGF‐β‐activated kinase 1

total cholesterol

triglyceride

tissue inhibitor of metalloproteinase 1

INTRODUCTION

NAFLD is a major global health problem that is considered to be a chronic hepatic manifestation of metabolic syndrome, ranging from steatosis to steatohepatitis and liver fibrosis.[ 1, 2 ] NASH is a severe form of NAFLD and a risk factor for cirrhosis and HCC.[ ] Steatosis, hepatocyte injury, inflammation, and various degrees of fibrosis are recognized as important components in the pathogenesis of NASH.[ ] However, there is currently no effective pharmacological therapy approved for NASH, and efforts to control complications arising from the condition are far from satisfactory.[ ] Therefore, approaches for developing promising treatment strategies for this highly prevalent disease are urgently needed.

Breviscapine, as a prescription drug, is a flavonoid extract from the traditional Chinese herb Erigeron breviscapus (Vant.) Hand.‐Mazz, whose main ingredient is scutellarin (≥90%).[ ] At present, more than 10 million patients use breviscapine and related drugs each year in China.[ ] Owing to its multiple pharmacological activities, including anti‐inflammatory, antioxidative, antiapoptotic, vasorelaxant, antiplatelet, anticoagulation, and myocardial protective activities, breviscapine has been widely used to treat cerebrovascular disease, cardiovascular disease, and diabetic complications.[ ] Recent studies have demonstrated that breviscapine protects against CCl4___induced liver injury by reducing proinflammatory cytokine secretion and oxidative stress.[ ] Furthermore, scutellarin, the main component of breviscapine, has been shown to regulate lipid metabolism and may exhibit useful therapeutic effects against NAFLD by reducing oxidative stress.[ 10, 11 ] However, whether breviscapine prevents NASH and its underlying mechanism are not yet clear.

The present study was designed to determine the effects of breviscapine on the pathogenesis of NASH, including its impacts on lipid accumulation, hepatic injury, inflammation, and fibrosis. We found that breviscapine attenuated hepatic lipid accumulation, inflammation, and fibrosis in mice fed a high‐fat diet (HFD), a high‐fat/high‐cholesterol (HFHC) diet, or a methionine‐ and choline‐deficient (MCD) diet, individually. Mechanistically, our data identify breviscapine as a promising candidate for the treatment of NASH whose effects are mediated by direct inhibition of TGF‐β‐activated kinase 1 (TAK1).

MATERIALS AND METHODS

Animals and treatments

Healthy 8‐ to 10‐week‐old male C57BL/6J mice were housed in cages under standard laboratory conditions (22–24°C temperature, 55%–60% relative humidity, and 12 h light/dark cycle). A mouse model of NASH was established by feeding the mice an HFHC (protein, 14%; fat, 42%; carbohydrates, 44%; cholesterol, 2%; IMA2019001; TrophicDiet, Nantong, China) diet for 16 weeks, and breviscapine (15 and 30 mg/kg) (190620, Shanghai Winherb Medical Science Co., Ltd., Shanghai, China) was administered daily by oral gavage after 8 weeks. The animal experiments were approved by the Institutional Animal Care and Use Committee of Renmin Hospital of Wuhan University according to the Guide for the Care and Use of Laboratory Animals prepared by the National Academy of Sciences and published by the National Institutes of Health. All mice received humane care.

RNA‐sequencing analysis

Total RNA was isolated using TRIzol reagent (T9424; Sigma‐Aldrich, St. Louis, MO), and the integrity of the RNA was assessed with agarose electrophoresis and an RNA 6000 Nano Kit (5067‐1511; Agilent). For RNA sequencing (RNA‐seq), complementary DNA libraries were created using the MGIEasy RNA Library Prep Kit (1000006384; MGI, Shenzhen, China) according to the manufacturer’s instructions. The single‐ended library was sequenced using a BGISEQ‐500 instrument (MGI), and the read length was 50 bp. The reads were mapped to the Ensembl human (hg38/GRCh38) and mouse (mm10/GRCm38) reference genomes with HISAT2 (version 2.1.0) software. Then, SAMtools (version 1.4) was used to convert the files obtained in the above steps into a binary BAM format that could store the comparison information. Next, the fragments per kilobase of exon per million mapped fragments value for each identified gene was calculated using the default parameters of StringTie (version 1.3.3b). Differentially expressed genes (DEGs) were identified with DESeq2 (version 1.2.10) according to the following two criteria: (1) a fold change >1.5 and (2) a corrected p value < 0.05. Finally, the ggplot2 package was used to create a volcano plot showing the fold changes and p values of all genes. The data have been submitted to Sequence ReadArchive (SRA) database with PRJNA number PRJNA781787 and PRJNA782086.

Principal component analysis

Principal component analysis (PCA) was used to compare the differences between the two groups. PCA was implemented using the fast.prcomp function of the R language and then visualized with the ggplot2 package.

Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis

The Kyoto Encyclopedia of Genes and Genomes (KEGG) is a comprehensive database that integrates genomic, chemical, and system functional information. KEGG pathway enrichment analysis was performed using Fisher’s exact test with our in‐house R script, and the KEGG pathway annotations for all genes in the selected genome were downloaded from the KEGG database. Pathways with p < 0.05 were considered significantly enriched.

Gene set variation analysis

Gene set variation analysis (GSVA) was executed using the GSVA R package (version 1.32.0). For each KEGG pathway, a gene set based on the proteins detected in proteomic analysis was generated. Then, the gene expression matrix of the KEGG pathway was used to calculate the single‐sample GSVA scores.

Gene set enrichment analysis

Gene set enrichment analysis (GSEA) was executed on the Java GSEA platform (version 3.0) using the Signal2Noise metric. The genes involved in each KEGG biological pathway were treated as a gene set, and then a ranked list and a gene set permutation of the gene set were generated. Gene sets with p < 0.05 and false discovery rate values < 0.25 were considered statistically significant.

Proteomic analysis

Proteomic analysis was performed using mass spectrometry in data‐dependent acquisition mode. The abundance of peptide fractions in the control group and the treatment group was quantified using unlabeled quantitative methods, and the difference in the abundance of phosphorylated peptide fractions between the control group and the treatment group was tested for significance by t test. The screening criterion for differentially expressed proteins was p < 0.05. The data have been submitted to NCBI PRIDE database with PXD number PXD029867.

Phosphoproteomic analysis

Phosphorylated polypeptides from cell or tissue lysates were enriched with TiO2 magnetic beads and analyzed by liquid chromatography‐tandem mass spectrometry with a Proxeon Easy‐NLC 1000 (Thermo Scientific, Shanghai Applied Protein Technology Co., Ltd., Shanghai, China). The abundance of phosphorylated peptides in the control group and the treatment group was quantified using an unlabeled quantitative method, and the difference in the abundance of phosphorylated peptide fractions between the control group and the treatment group was tested for significance by t test. The screening criterion for differentially expressed phosphorylated proteins was p < 0.05. The data have been submitted to NCBI PRIDE database with PXD number PXD029975.

The detailed materials and methods are provided in the Supporting Information.

RESULTS

Breviscapine attenuates HFD‐induced hepatic injury and metabolic disorders in mice

To examine the pharmacological effects of breviscapine on the hepatic steatosis and injury, we used an HFD‐induced obesity mouse model that destroys the balance between the synthesis and degradation of lipids, which results in lipid accumulation and toxicity in hepatocytes. Mice were fed an HFD for 16 weeks and then administered breviscapine intragastrically (15 and 30 mg/kg/day) for an additional 8 weeks with continuous HFD feeding (Figure 1A). During the whole process of the experiment, the body weight of the HFD‐fed mice (control‐HFD mice) was significantly higher than that of the normal chow (NC)–fed control mice (control‐NC mice), while there were no differences between breviscapine‐treated mice and HFD‐fed mice (Figure S1A). The increased liver weight (LW) and LW‐to‐body weight ratio (LW/BW) in HFD‐fed mice was remarkably decreased by breviscapine (Figure 1B). In addition, hematoxylin and eosin (H&E) staining and oil red O staining showed that lipid accumulation in HFD‐fed mice was significantly reduced by breviscapine (Figure 1C), which was supported by the liver and serum triglyceride (TG) and total cholesterol (TC) levels (Figure 1D). Additionally, reduced hepatocyte ballooning was observed in mice treated with breviscapine (Figure 1C). Quantitative PCR assays showed that hepatic mRNA levels of genes related to fatty acid uptake and synthesis (Cd36, fatty acid binding protein 1 [Fabp1], fatty acid synthase [Fasn], stearoyl CoA desaturase 1 [Scd1], and peroxisome proliferator–activated receptor γ [Pparγ]) were significantly diminished by breviscapine treatment. Moreover, the mRNA levels of lipid β‐oxidation genes (Pparα and carnitine palmitoyltransferase 1a [Cpt1α]) were increased in breviscapine‐treated mice (Figure 1E). Western blot assays also confirmed that the up‐regulation of adipogenic genes (FASN and PPARγ) in HFD‐fed mice was significantly abolished by breviscapine (Figure S1B). The serum levels of alanine aminotransferase (ALT) and aspartate aminotransferase (AST) were lower in breviscapine‐treated mice than in control‐HFD mice (Figure 1F). These results suggest that breviscapine protects mice from hepatic steatosis and injury induced by metabolic stress.

FIGURE 1

Breviscapine alleviates hepatic steatosis and injury in mice fed the HFD. (A) Schematic diagram of the experimental procedure used to examine the protective role of breviscapine in mice fed the NC diet or the HFD for 24 weeks. HFD‐fed mice were intragastrically administered saline or breviscapine (15 or 30 mg/kg) once daily beginning at week 16 for 8 weeks. n = 8 per group. (B) LW and LW/BW values of mice. n = 8 per group. (C) Representative images of H&E and oil red O staining of liver sections from mice. NAFLD activity score and the statistics of oil red O–positive areas are shown. n = 6 per group. Scale bar, 100 µm. (D) Lipid (TG and TC) levels in the liver and serum. n = 8 per group. (E) Quantitative PCR was performed to determine the hepatic mRNA levels of genes related to fatty acid metabolism (Cd36, Fabp1, Fasn, Scd1, Pparγ, Pparα, and Cpt1a) in mice from the indicated groups. Gene expression was normalized to actin beta (Actb) mRNA levels. n = 4–6 per group. (F) Serum levels of ALT and AST were measured in mice after 24 weeks of NC diet feeding or HFD challenge. n = 8 per group. The data are presented as the mean ± SD. Significant difference between the control‐NC group and the control‐HFD group, *p < 0.05, **p < 0.01; significant difference between the control‐HFD group and the Bre‐HFD group, # p < 0.05, ## p < 0.01. The data were analyzed with one‐way ANOVA. Abbreviations: Bre, breviscapine; i.g., intragastric administration; NAS, NAFLD activity score; ns, no significant difference between the control‐HFD group and the Bre‐HFD group

Breviscapine attenuates HFHC diet–induced steatohepatitis

To further assess the potential benefit of breviscapine treatment against the progression of NASH, we next tested the effect of breviscapine in an HFHC diet–induced NASH mouse model. HFHC diet feeding induced more severe hepatic steatosis, inflammation, and fibrosis than HFD feeding, eventually leading to steatohepatitis with features similar to the pathological features of human NASH.[ ] Mice were fed an HFHC diet for 8 weeks and then treated with breviscapine (15 and 30 mg/kg/day) by oral gavage with continuous HFHC diet feeding for another 8 weeks (Figure 2A). The LW and the LW/BW of breviscapine‐treated mice fed the HFHC diet were lower than those of the control HFHC diet–fed mice (Figure 2B). However, there was no difference in body weight after breviscapine treatment (Figure S2A). Histological examination of liver sections showed that breviscapine treatment significantly improved liver histology, resulting in less hepatic steatosis, ballooning degeneration, and inflammatory cell infiltration in treated mice than in control‐HFHC mice (Figure 2C). Consistently, hepatic and serum lipid levels were significantly decreased in breviscapine‐treated mice (Figure 2D). Moreover, breviscapine treatment strikingly decreased the mRNA levels of fatty acid uptake–related and synthesis‐related genes (Cd36, Fabp1, Fasn, Scd1, and Pparγ) (Figure 2E) and alleviated the abnormal protein expression of FASN and PPARγ (Figure S2B) in mice fed the HFHC diet. Additionally, the mRNA levels of β‐oxidation genes (Pparα and Cpt1α) were notably increased in breviscapine‐treated mice compared to control HFHC mice (Figure 2E). Furthermore, the serum levels of ALT and AST in mice fed the HFHC diet were significantly reduced by breviscapine treatment (Figure 2F).

FIGURE 2

Breviscapine represses hepatic steatosis and injury in mice fed an HFHC diet. (A) Schematic diagram of the experimental procedure used to examine the protective effects of breviscapine in mice fed NC or the HFHC diet for 16 weeks. HFHC‐fed mice were intragastrically administered saline or breviscapine (15 or 30 mg/kg) once daily beginning at week 8 for 8 weeks. n = 9 per group. (B) LW and LW/BW values of mice. n = 9 per group. (C) Representative images of H&E and oil red O staining of liver sections from mice. NAFLS activity scores and the statistics of oil red O–positive areas are shown. n = 6 per group. Scale bar, 100 µm. (D) Lipid (TG and TC) levels in the liver and serum. n = 9 per group. (E) Quantitative PCR was performed to determine the hepatic mRNA levels of genes related to fatty acid metabolism (Cd36, Fabp1, Fasn, Scd1, Pparγ, Pparα, and Cpt1a) in mice from the indicated groups. Gene expression was normalized to Actb mRNA levels. n = 6 per group. (F) Serum levels of ALT and AST were measured in mice after 16 weeks of HFHC challenge. n = 9 per group. The data are presented as the mean ± SD. Significant difference between the control‐HFHC group and the Bre‐HFHC group, *p < 0.05, **p < 0.01. The data were analyzed with one‐way ANOVA. Abbreviations: Bre, breviscapine; i.g., intragastric administration; NAS, NAFLD activity score

Liver inflammation and fibrosis exacerbate NASH progression.[ ] Next, we examined the effects of breviscapine on NASH‐associated inflammation and fibrosis. Immunofluorescence staining of the macrophage marker CD11b showed markedly less macrophage infiltration in liver sections from breviscapine‐treated mice than in those from control‐HFHC mice (Figure 3A). Quantitative PCR assays showed that breviscapine significantly reduced the mRNA levels of proinflammatory genes (C‐C motif chemokine ligands 2 and 5 [Ccl2, Ccl5], C‐X‐C motif chemokine ligands 2 and 10 [Cxcl2, Cxcl10], and Tnfα) (Figure 3B). Furthermore, breviscapine treatment significantly down‐regulated the hepatic expression of phosphorylated inhibitor of NF‐κB kinase subunit beta (p‐IKKβ) and nuclear p‐p65 in mice fed the HFHC diet, whereas it up‐regulated the expression of inhibitor of kappa light polypeptide gene enhancer in B cells alpha (IKBα; Figure 3C). We next measured the effect of breviscapine on NASH‐associated fibrosis in mice fed the HFHC diet. Picrosirius red (PSR) staining showed less collagen deposition in the liver sections of breviscapine‐treated mice than in those of control HFHC mice (Figure 3D). Additionally, quantitative PCR assays confirmed that breviscapine treatment significantly decreased the hepatic mRNA levels of profibrotic genes (collagen type I/III alpha 1 [Col1a1, Col3a1], connective tissue growth factor (Ctgf), and tissue inhibitor of metalloproteinase 1 [Timp1]) in mice fed the HFHC diet (Figure 3E). These results indicate that breviscapine ameliorates HFHC diet–induced steatohepatitis in mice.

FIGURE 3

Breviscapine attenuates hepatic inflammation and fibrosis in mice fed an HFHC diet. (A) Immunofluorescence staining of CD11b (in red) in the livers of HFHC‐fed mice. Nuclei were labeled with DAPI (in blue). The ratios of the numbers of CD11b‐positive cells are shown. n = 3–4 per group. Scale bar, 50 µm. (B) Quantitative PCR was performed to determine hepatic mRNA levels of genes related to inflammation (Ccl2, Cxcl2, Ccl5, Cxcl10, and Tnfa) in mice from the indicated groups. Gene expression was normalized to Actb mRNA levels. n = 6 per group. (C) Expression of p‐IKKβ, IKKβ, IKBα, p‐p65, p65, and glyceraldehyde 3‐phosphate dehydrogenase (GAPDH) was analyzed by western blotting. GAPDH served as a loading control. n = 3 per group. (D) Representative histological images of liver sections were examined by PSR staining. n = 6 per group. Scale bar, 100 µm. (E) Quantitative PCR was performed to determine hepatic mRNA levels of genes related to fibrosis indicators (Col1a1, Col3a1, Ctgf, and Timp1) in mice from the indicated groups. Gene expression was normalized to Actb mRNA levels. n = 6 per group. The data are presented as the mean ± SD. *Significant difference between the control‐HFHC group and the breviscapine‐HFHC group, *p < 0.05, **p < 0.01. The data were analyzed with one‐way ANOVA. Abbreviations: Bre, breviscapine; HPF, high‐power field

Breviscapine attenuates MCD diet–induced steatohepatitis

To further verify the effect of breviscapine on steatohepatitis, another mouse steatohepatitis model induced by MCD diet feeding for 4 weeks was employed.[ ] MCD diet–fed mice were administered breviscapine (30 mg/kg) daily for the entire experimental duration (Figure 4A). Consistent with the results obtained with the HFHC mouse model, serum levels of ALT and AST were significantly reduced by breviscapine treatment (Figure 4B). Hepatic TG content was lower in breviscapine‐treated mice than in control‐MCD mice (Figure 4C). Moreover, H&E staining and oil red O staining showed that breviscapine‐treated mice had less severe steatosis and inflammatory cell infiltration than control‐MCD mice (Figure 4D). Additionally, CD11b and PSR staining showed that MCD‐induced hepatic inflammation and collagen deposition were abolished by breviscapine treatment (Figure 4E). Quantitative PCR assays also confirmed that breviscapine decreased the mRNA levels of genes related to inflammation (Ccl2, Cxcl2, and Tnfα) and fibrosis (Col1a1, Col3a1, and Ctgf) (Figure 4F). These data further confirm that breviscapine protects mice against experimental steatohepatitis.

FIGURE 4

Breviscapine suppresses hepatic steatosis, inflammation, and fibrosis in mice fed an MCD diet. (A) Schematic diagram of the experimental procedure used to examine the protective effects of breviscapine in mice fed a methionine‐ and choline‐sufficient (MCS) or MCD diet for 4 weeks. MCD diet–fed mice were intragastrically administered saline or breviscapine (15 or 30 mg/kg) once daily for 4 weeks. n = 12–15 per group. (B) Serum levels of ALT and AST measured in mice after 4 weeks of MCD diet challenge. n = 12–15 per group. (C) Liver lipid (TG) levels. n = 10–15 per group. (D) Representative images of hepatic steatosis and lipid accumulation stained with H&E and oil red O. NAFLD activity scores and the ratios of oil red O–positive areas are shown. n = 6 per group. Scale bar, 100 µm. (E) Immunofluorescence staining of CD11b (in red) and representative histological images stained with PSR in liver sections from the indicated mice after MCD diet challenge. Nuclei were labeled with DAPI (in blue). The ratios of the numbers of CD11b‐positive cells are shown. n = 4 per group. Scale bar, 50 µm. n = 6 mice per group. Scale bar, 50 µm. (F) Quantitative PCR was performed to determine the hepatic mRNA levels of genes related to inflammation (Ccl2, Cxcl2, and Tnfa) and fibrosis indicators (Col1a1, Col3a1, and Ctgf) in the liver samples of control‐MCS, control‐MCD, and breviscapine (30 mg/kg)‐MCD mice. Gene expression was normalized to Actb mRNA levels. n = 5 per group. The data are presented as the mean ± SD. Significant difference between the control‐MCS group and the control‐MCD group, *p < 0.05, **p < 0.01; significant difference between the control‐MCD group and the breviscapine‐MCD group, # p < 0.05, ## p < 0.01. The data were analyzed with one‐way ANOVA or two‐tailed Student t test. Abbreviations: Bre, breviscapine; HPF, high‐power field; i.g., intragastric administration; NAS, NAFLD activity score

Breviscapine attenuates lipid accumulation and inflammation in hepatocytes

Primary hepatocytes are considered to be the primary cell type of the liver architecture and the major cells responsible for metabolic disorders. Next, we investigated the effects of breviscapine (50 and 100 μM) on lipid accumulation and inflammation in primary hepatocytes and L02 hepatocytes treated with palmitic acid (PA) + oleic acid (OA; = PO). Lipid accumulation (assayed by oil red O staining) (Figure 5A,B) and hepatic TG (Figure 5C) and TC (Figure S3) levels were significantly increased in primary hepatocytes and L02 hepatocytes treated with PO. However, the induction of lipid accumulation in primary hepatocytes and L02 hepatocytes was blunted by breviscapine treatment (Figure 5A–C). Furthermore, breviscapine treatment reduced the mRNA levels of the lipogenic genes Fasn, Scd1, and acetyl CoA carboxylase alpha (Accα) in PO‐treated primary hepatocytes (Figure 5D). In addition, the expression of proinflammatory genes (Cxcl10, Tnfα, and Ccl2) was significantly decreased by breviscapine in primary hepatocytes cultured in PO medium (Figure 5E). Similar results from the quantitative PCR assay were observed in L02 hepatocytes cultured with PO in the absence or presence of breviscapine (Figure 5D,E). Together, these data indicate that breviscapine attenuates steatohepatitis changes in hepatocytes.

FIGURE 5

Breviscapine reduces lipid accumulation and inflammation in primary hepatocytes and L02 cells treated with PO. Primary hepatocytes were cultured in PO medium containing 0.2 mM PA and 0.4 mM OA for 12 h. L02 cells were cultured in PO medium containing 0.5 mM PA and 1.0 mM OA for 12 h. (A) Representative oil red O staining images of primary hepatocytes and L02 cells stimulated with PO and treated with breviscapine (50 and 100 μM) for 12 h. Scale bar, 100 μm. n = 3 independent experiments per group. (B) Quantitative analysis of oil red O staining in hepatocytes. (C) Levels of TGs and TC in primary hepatocytes and L02 cells in the indicated groups. The data were obtained from three independent experiments per group. (D) mRNA levels of genes related to lipid metabolism (Fasn, Scd1, and Accα) and (E) inflammation (Cxcl10, Tnfα, and Ccl2) in primary hepatocytes and L02 cells in the indicated groups. Gene expression was normalized to Actb mRNA levels. The data were obtained from three independent experiments per group. The data are presented as the mean ± SD. Significant difference between the DMSO‐bovine serum albumin group and the DMSO‐PO group, *p < 0.05, **p < 0.01; significant difference between the DMSO‐PO group and the breviscapine‐PO group, # p < 0.05, ## p < 0.01. The data were analyzed with one‐way ANOVA. Abbreviations: Bre, breviscapine; BSA, bovine serum albumin

Breviscapine prevents NASH by systematically blunting the pathways involved in lipid metabolism, inflammation, fibrosis, and apoptosis

To systemically elucidate how breviscapine prevents NASH, we performed RNA‐seq on the livers of HFHC diet–fed mice treated with breviscapine. Unsupervised PCA and hierarchical clustering clearly separated the samples from the control mice fed the HFHC diet and the breviscapine‐treated mice fed the HFHC diet into two clusters (Figure 6A). GSEA revealed that genes related to lipid metabolism, inflammation, fibrosis, and apoptosis were enriched (Figure 6B) and significantly down‐regulated by breviscapine treatment (Figure 4A). Moreover, a heatmap based on GSEA indicated that hepatic genes involved in lipid metabolism, inflammation, fibrosis, and apoptosis pathways were significantly down‐regulated by breviscapine treatment in mice fed the HFHC diet for 16 weeks in vivo (Figure 6C).

FIGURE 6

Transcriptomic analysis revealed the key differential targets in breviscapine‐treated, HFHC diet–fed mice. (A) PCA of RNA‐seq data from mice fed a control or HFHC diet for 16 weeks. n = 5 per group. (B) GSEA of pathways related to lipid metabolism, inflammation, fibrosis, and apoptosis. n = 5 per group. (C) Heatmaps of gene expression profiles related to lipid metabolism, inflammation, fibrosis, and apoptosis based on the RNA‐seq data set. n = 5 per group. Abbreviations: Bre, breviscapine; Ecm, extracellular matrix; NES, normalized enrichment score; PC, principal component

Furthermore, we performed RNA‐seq analysis on PO‐treated L02 cells treated with breviscapine in vitro. PCA and hierarchical clustering clearly separated the samples from the breviscapine‐treated group and the DMSO‐control group (Figure S5A). The GSEA results showed that the lipid metabolism–related, inflammation‐related, and apoptosis‐related pathways were mainly enriched (Figure S5B,C) and that genes were significantly down‐regulated by breviscapine treatment compared with DMSO treatment (Figure S6). Collectively, these results demonstrate that breviscapine treatment protects against NASH by suppressing hepatic damage, steatosis, inflammation, and fibrosis under conditions of metabolic stress.

Breviscapine inhibits the mitogen‐activated protein kinase pathway in hepatocytes and liver tissue under metabolic stress

To investigate the molecular mechanism(s) underlying the anti‐NASH effect of breviscapine, the transcriptomes of breviscapine‐treated hepatocytes treated with PO and mice fed the HFHC diet were combined and analyzed by KEGG analysis and GSEA (Figure 7A). Subsequently, we analyzed the data in combination with the transcriptomes of cells and mice and found that a total of nine signaling pathways were enriched at the cellular level, while five signaling pathways were enriched at the tissue level. Through Venn analysis, we found that breviscapine treatment affected two identical signaling pathways in vitro and in vivo: the mitogen‐activated protein kinase (MAPK) signaling pathway and the TGF‐β signaling pathway (Figure 7B). Moreover, the gene numbers of the MAPK signaling pathway were greater than those of the TGF‐β signaling pathway (Figure 7B). Subsequently, proteomic analysis showed that differential proteins related to lipid metabolism and inflammation were highly enriched in mice fed the HFHC diet but significantly down‐regulated by breviscapine treatment (Figure 7C,D).

FIGURE 7

Multiomics analysis revealed that the MAPK pathway is the downstream target of breviscapine. (A) GSEA and KEGG pathway enrichment analysis of the transcriptomes of cells and HFHC diet–fed mouse liver samples. n = 3 independent experiments per group for cells, n = 5 per group for tissues. (B) Venn diagram of GSEA and KEGG pathway intersections based on transcriptome data both in vivo and in vitro. The histogram of the two key pathways is based on the corresponding scores. n = 3 independent experiments per group for cells, n = 5 per group for tissues. (C) Scatterplot showing the fold changes in all proteins in the proteomic data set, with differentially expressed proteins shown with corresponding font colors (inflammation in blue, lipid metabolism in yellow, fibrosis in red, apoptosis in purple). n = 3 per group. (D) Heatmap showing the degrees of influence of pathways related to lipid metabolism and inflammation based on a proteomic data set. n = 3 per group. (E) PPI network diagram of multiomics joint analysis. n = 3 per group. (F) Bar chart of phosphorylated proteins ranked according to PPI degree. n = 3 per group. (G) Western blot analysis of proteins involved in the MAPK signaling cascade in cells stimulated with PA and mice subjected to HFD or HFHC feeding. GAPDH served as a loading control. n = 3 per group. The data are presented as the mean ± SD. Significant difference between the DMSO/control group and the breviscapine group, *p < 0.05, **p < 0.01. The data were analyzed with one‐way ANOVA. Abbreviations: Bre, breviscapine; FDR, false discovery rate; Hspa8, heat shock protein a8; RIG‐I, retinoic acid–inducible gene I

To uncover the key effector(s) that contributes to breviscapine‐induced protection from NASH, we integrated the results from proteomic analysis, phosphoproteomic analysis, and transcriptomic characterization. We found a strong correlation between the DEGs and the differentially phosphorylated MAPK signaling molecules under breviscapine treatment, as illustrated by the STRING protein interaction network database (Figure 7E). Furthermore, MAPK14, heat shock protein A8, and MAP3K7 were the top 3 phosphoproteins according to the protein–protein interaction (PPI) degree (Figure 7F). MAPK14 is a member of the MAPK family. MAPK signaling is integrally involved in various physiological processes, such as inflammation, apoptosis, oxidative stress, and fibrotic responses.[ ] TAK1 is a MAP3K protein and has been defined as a critical upstream molecule of MAPKs.[ ] Thus, to further investigate whether breviscapine protects against NASH by regulating the MAPK signaling pathway, western blotting was used to analyze the activation of MAPK signaling in vitro and in vivo. The phosphorylation of TAK1, c‐Jun N‐terminal kinase (JNK), and p38 was significantly suppressed by breviscapine treatment after PA, HFD, or HFHC diet treatment (Figure 7G). These data suggest that breviscapine inhibits TAK1 and its downstream signaling pathway in hepatocytes or livers subjected to metabolic stress.

TAK1 inhibition is required for the breviscapine‐induced suppression of lipid accumulation and the proinflammatory response in hepatocytes

To verify whether inhibition of TAK1‐dependent signaling is required for the anti‐NASH effects of breviscapine, we used 5Z‐7‐oxozeaenol, which is a potent irreversible selective inhibitor of TAK1. 5Z‐7‐Oxozeaenol prevents inflammation by inhibiting the catalytic activity of TAK1 MAPK kinase kinase.[ ] Notably, oil red O staining showed that 5Z‐7‐oxozeaenol pretreatment almost completely abrogated the reduction in lipid accumulation induced by breviscapine treatment in PO‐treated hepatocytes (Figure 8A,B). In addition, breviscapine significantly decreased the intracellular TG and TC levels, while these effects were abolished by pretreatment with 5Z‐7‐oxozeaenol (Figure 8C). Furthermore, pretreatment with 5Z‐7‐oxozeaenol abolished the reductions in the expression of fatty acid synthesis genes (Fasn, Scd1, and Accα) and proinflammatory genes (Tnfα, Ccl2, and Cxcl10) induced by breviscapine treatment (Figure 8D). Moreover, pretreatment with 5Z‐7‐oxozeaenol almost completely abrogated the breviscapine‐induced phosphorylation of TAK1 and the downstream phosphorylation of JNK and p38 without changing the expression levels of total TAK1, JNK, and p38 in hepatocytes cultured with PA for 12 h (Figure 8E; Figure S7). Finally, we investigated the possibility of a direct interaction between breviscapine and TAK1. Because scutellarin, a prescription drug, is the main ingredient (≥90%) of breviscapine, we performed molecular docking of scutellarin with the TAK1 protein. Scutellarin exhibited suitable steric complementarity with the binding site of TAK1. Hydrogen bond interactions occurred between TAK1 and scutellarin. The oxygen atoms of scutellarin, regarded as hydrogen bond acceptors, formed hydrogen bonds with the backbone nitrogen atoms of Ala107 and Asp175. The nitrogen atom of scutellarin, regarded as a hydrogen bond donor, formed a hydrogen bond with the backbone oxygen atom of Ala107 (Figure 8F). These data suggest that inhibition of TAK1 is required for breviscapine to protect hepatocytes from PO‐induced lipid accumulation and inflammation.

FIGURE 8

TAK1 inhibition is responsible for the effects of breviscapine on lipid accumulation and inflammation in hepatocytes. (A) Representative images of oil red O staining showing the degree of lipid accumulation in primary hepatocytes treated with DMSO or breviscapine (100 μM) after PO (0.2 mM PA and 0.4 mM OA) stimulation for 12 h in the presence or absence of 5Z‐7‐oxozeaenol (2 μM) for 12 h. Scale bar, 200 µm. The data were obtained from three independent experiments per group. (B) Quantitative analysis of the degree of lipid accumulation. (C) TG and TC levels in primary hepatocytes in the indicated groups stimulated with bovine serum albumin or PO (0.5 mM PA and 1.0 mM OA) for 12 h. Four independent experiments per group were performed. (D) Relative mRNA levels of inflammatory factors (Tnfα, IL‐1β, and Ccl10) and lipogenic genes (Fasn, Scd1, and Accα) in primary hepatocytes treated with PO. The relative mRNA expression was normalized to that of the β‐actin‐encoding gene Actb. The data were obtained from three independent experiments per group. (E) Western blot analysis of total and phosphorylated protein levels of apoptosis signal–regulating kinase, JNK, and p38 protein in primary hepatocytes stimulated with PA for 12 h. Actin served as the loading control. The data were obtained from three independent experiments per group. (F) Molecular docking of breviscapine with the protein TAK1. The data are presented as the mean ± SD. Significant difference between the control‐DMSO‐PO group and the control‐breviscapine‐PO group, *p < 0.05, **p < 0.01. The data were analyzed with one‐way ANOVA. Abbreviations: Bre, breviscapine; ns, no significant difference between the 5Z‐7‐oxozeaenol‐DMSO‐PO group and the 5Z‐7‐oxozeaenol‐breviscapine‐PO group; 5Z‐7‐Ox, 5Z‐7‐oxozeaenol

DISCUSSION

Exercise training and dietary interventions remain the primary recommendations for patients with NAFLD and its progressive stage, NASH.[ ] However, for multiple societal, psychological, physical, genetic, and epigenetic reasons, it is challenging for patients to adhere to such lifestyle modifications; thus, pharmacotherapy for NAFLD, especially incurable NASH, is essential. Nevertheless, few medications are currently available to meet the increasing disease burdens of NAFLD and NASH.[ ] Traditional Chinese medicine (TCM) has played a very large role in the prevention and treatment of NAFLD in China due to its unique theories on the etiology, pathogenesis, diagnosis, and treatment of fatty liver, such as the theories of Tiao Gan Qi Shu Hua Zhuo and the liver‐based regulatory system for metabolic homeostasis.[ ] In recent years, accumulating evidence has demonstrated the beneficial role of Chinese herbal medicine in the treatment of NAFLD.[ 12, 21, 22 ] Breviscapine has a wide range of pharmacological effects, such as antioxidative, anti‐inflammatory, antifibrotic, and lipid‐lowering effects (similar to the Hua Zhuo effect of TCM) and is clinically used for the treatment of diabetes, NAFLD, and other metabolic diseases.[ 23, 24 ] These uses show that modern medicine and traditional medicine have the same understanding of these diseases. In this study, we discovered that breviscapine ameliorates NAFLD and NASH through suppression of the TAK1 signaling pathway.

Breviscapine, a mixture of flavonoid glycosides, has strong biological activity and is also a common monomeric compound used in clinical practice. It has been formulated into a variety of dosage forms for the treatment of cardiovascular and cerebrovascular diseases, such as atherosclerosis, coronary heart disease, and insufficient cerebral blood supply.[ ] A clinical study has shown that breviscapine reduces hyperlipidemia in patients.[ ] Additionally, breviscapine reduces hepatic lipid accumulation in HFD‐fed rats[ ] and mice.[ ] Our results revealed that breviscapine significantly attenuated lipid accumulation in response to metabolic stress in both murine models and in vitro hepatocytes. Furthermore, our results confirmed that breviscapine abolished hepatic inflammation and fibrosis in a variety of animal models under metabolic stress. Aberrant lipid metabolism, steatosis, hepatic inflammation, apoptosis, and fibrosis are considered the major pathological factors promoting the progression of NAFLD.[ ] Using an unbiased systemic analysis, we found that breviscapine largely decreased the expression of genes associated with lipid metabolism, inflammation, apoptosis, and fibrosis in mice fed an HFHC diet.

NAFLD is triggered by an excessive supply of nutrients. Most ingested fat is diverted to adipose tissue or working muscle for storage or oxidation, but fat stored in white adipose tissue undergoes lipolysis to release fatty acids.[ ] Excess fatty acids contribute to accumulation of very large amounts of TGs in hepatocytes and frequently lead to steatosis in the liver.[ ] In this study, breviscapine reduced the hepatic lipid content in PO‐treated hepatocytes and the livers of diet‐induced mice. Moreover, breviscapine inhibited Fasn, Accα, Scd1, Cd36, Fabp1, and Pparγ mRNA expression and promoted Pparα and Cpt1α mRNA expression in vivo and in vitro, suggesting that breviscapine may decrease hepatic lipid accumulation and steatosis by restraining de novo synthesis and uptake of fatty acids while accelerating fatty acid β‐oxidation.

A characteristic of NASH is the presence of hepatic inflammation, and continued inflammation in the liver is thought to drive the development of fibrosis.[ ] Previous studies have demonstrated that breviscapine inhibits lipopolysaccharide‐induced production of proinflammatory mediators.[ ] Recent studies have shown that the anti‐inflammatory effects of breviscapine may be due to suppression of NF‐κB and NLR family pyrin domain containing 3 inflammasome activation.[ 33, 34 ] In the current study, breviscapine inhibited NF‐κB signaling pathway activation and attenuated the secretion of proinflammatory cytokines under metabolic stress. Studies have shown that activation of NF‐κBs, especially the NF‐κB p65 subunit, is closely associated with collagen deposition in liver fibrosis.[ ] In our study, genes associated with collagen production, such as Col1a1, Col3a1, Ctgf, and Timp1, were down‐regulated by breviscapine in HFHC/MCD diet–induced murine models. In addition, analysis of PSR staining representing collagen deposition showed that breviscapine attenuated NASH‐associated fibrosis.

To elucidate the possible mechanism by which breviscapine prevents NASH, a systematic analysis based on RNA sequence data was performed, and our data implicated MAPK signaling in the anti‐NASH effects of breviscapine. The MAPK signaling pathway is involved in the regulation of inflammation and fatty acid metabolism.[ ] TAK1 is an upstream kinase that activates NF‐κB and MAPK signaling during NAFLD progression.[ ] Due to the pivotal role of TAK1 in many physiological processes, such as inflammation, cell differentiation, and apoptosis, inhibition of TAK1 activation induced by metabolic stress may generate a potent hepatoprotective effect.[ ] It has been demonstrated that inhibition of the long noncoding RNA HULC attenuates hepatic fibrosis and hepatocyte apoptosis by inhibiting the MAPK signaling pathway in rats with NAFLD.[ ] Recent studies have pointed out the important roles of TAK1 in regulating hepatocyte lipid metabolism and reducing hepatocyte inflammation in NAFLD.[ 40, 41 ] Restraining TAK1 activation attenuates hepatic steatosis, lipid deposition, and the inflammatory response in animals with HFD‐induced NAFLD.[ 42, 43, 44 ] The results of the multiomics and western blot analyses further suggest that breviscapine inhibits TAK1 phosphorylation and downstream JNK/p38 activation, thereby postponing the pathogenesis of NASH. Our study demonstrates that breviscapine ameliorates NASH under metabolic stress conditions by regulating lipid metabolism, suppressing inflammation, and reducing liver fibrosis by specifically inhibiting TAK1 phosphorylation.

Collectively, our current findings demonstrate that breviscapine, identified as a specific inhibitor of TAK1, ameliorates the pathogenesis of NASH through suppression of hepatic lipid accumulation, inflammation, and fibrogenesis. Therefore, our findings will provide insights supporting the development of therapeutic approaches for NASH.

CONFLICT OF INTEREST

Nothing to report.

AUTHOR CONTRIBUTIONS

Conceptualization: Jiao Guo; Project administration: Jiao Guo; Funding acquisition: Jiao Guo, Tian Lan; Visualization: Tian Lan, Shuo Jiang, Jing Zhang, Qiqing Weng, Yang Yu, and Haonan Li; Writing – original draft: Tian Lan, Shuo Jiang, Jing Zhang, and Qiqing Weng; Writing – review & editing: Tian Lan, Shuo Jiang, Jing Zhang, Qiqing Weng, and Sha Hu; Investigation: Shuo Jiang, Jing Zhang, Qiqing Weng, Yang Yu, Haonan Li, and Song Tian; Data curation: Shuo Jiang, Jing Zhang, Qiqing Weng, Xin Ding, Yiqi Yang, Weixuan Wang, Lexun Wang, Duosheng Luo, Xue Xiao, Shenghua Piao, Qing Zhu, and Xianglu Rong; Supervision: Jiao Guo

Fig S1

Click here for additional data file.

Fig S2

Click here for additional data file.

Fig S3

Click here for additional data file.

Fig S4

Click here for additional data file.

Fig S5

Click here for additional data file.

Fig S6

Click here for additional data file.

Fig S7

Click here for additional data file.

Supplementary Material

Click here for additional data file.