- Research
- Open access
- Published:
Revealing potential anti-fibrotic mechanism of Ganxianfang formula based on RNA sequence
Chinese Medicine volume 17, Article number: 23 (2022)
Abstract
Background
Ganxianfang (GXF) formula as a traditional Chinese medicine (TCM) is used for liver fibrosis in clinical practice while its mechanism is unclear. The aim of this study is to explore the molecular mechanism of GXF against CCl4-induced liver fibrosis rats.
Methods
Detected the main compounds of GXF by UPLC-MS/MS. Evaluated the efficacy of GXF (1.58, 3.15, 4.73 g/kg/day) and Fuzheng Huayu (FZHY, positive control, 0.47 g/kg/day) through serum alanine aminotransferase (ALT), aspartate aminotransferase (AST) levels and histopathological changes. Explored the underlying mechanisms by integrating our total liver RNA sequencing (RNA-seq) data with recent liver single-cell sequencing (scRNA-seq) studies. Verified potential pharmacodynamic substances of GXF by hepatic stellate cell (HSC)-T6 line.
Results
Main compounds were identified in GXF by UPLC-MS/MS, including baicalin, wogonoside and matrine etc. With GXF-high dose treatment, the elevation of ALT and AST induced by CCl4 were significantly reduced, and the protective effect of GXF-high dose treatment was better than FZHY. Liver histopathological changes were alleviated by GXF-high dose treatment, the ISHAK scoring showed the incidence of liver cirrhosis (F5/F6) decreased from 76.5 to 55.6%. The results of liver hydroxyproline content were consistent with the histopathological changes. RNA-seq analysis revealed the differential genes (DEGs) were mainly enriched in ECM-receptor interaction and chemokine signaling pathway. GXF effectively inhibited collagen deposition and significantly downregulated CCL2 to inhibit the recruitment of macrophages in liver tissue. Integrating scRNA-seq data revealed that GXF effectively inhibited the expansion of scar-associated Trem2+CD9+ macrophages subpopulation and PDGFRα+PDGFRβ+ scar-producing myofibroblasts in the damaged liver, and remodeled the fibrotic niche via regulation of ligand-receptor interactions including TGFβ/EGFR, PDGFB/PDGFRα, and TNFSF12/TNFRSF12a signaling. In vitro experiments demonstrated that baicalin, matrine and hesperidin in GXF inhibited the activation of hepatic stellate cells.
Conclusions
This study clarified the potential anti-fibrotic effects and molecular mechanism of GXF in CCl4-induced liver fibrosis rats, which deserves further promotion and application.
Background
Liver fibrosis is a histological change caused by acute or chronic liver diseases, such as hepatitis virus infection, autoimmune liver diseases, and alcoholic or non-alcoholic fatty liver disease [1]. The main pathological feature of liver fibrosis is the excessive deposition of extracellular matrix (ECM) component, which disrupts normal liver structure and function, resulting in the decompensation of liver cirrhosis and the occurrence of liver cancer [2, 3]. Therefore, it is important to effectively reverse liver fibrosis. At present, there is no specific medicine for the treatment of liver fibrosis except for the etiological treatment. Due to the complicated pathological mechanism of liver fibrosis, drugs developed for a single target are difficult to be effective in clinical practice [4]. Studies in recent decades have shown that traditional Chinese medicine (TCM) treatments have shown curative advantages in the field of prevention and treatment of liver fibrosis [5,6,7]. In clinical practice in China, several traditional Chinese medicines such as Fuzheng Huayu, Biejia Ruangan tablet and Anluo Huaxian pill are widely used in anti-fibrosis treatment, and Fuzheng Huayu has been included in the guideline for the diagnosis and treatment of liver fibrosis with integrated traditional Chinese and western medicine [8,9,10]. Therefore, traditional Chinese medicines have a broad prospect for the treatment of liver fibrosis.
Ganxianfang (GXF) formula is obtained by summarizing the traditional Chinese medicine prescriptions from "Jin Gui Yao Lue". Combined with the clinical experience of Professor Huabao Liu from Chongqing Traditional Chinese Medicine Hospital, GXF has been further optimized. This prescription contains a total of five herbs: Huang Qin (Scutellariae Radix), Ku Shen (Sophorae Flavescentis Radix), Fang Feng (Saposhnikoviae Radix), Chen Pi (Citrus Reticulata) and Xi Yang Shen (Panacis Quinquefolii Radix). Modern pharmacological studies indicated that Scutellariae Radix and Saposhnikoviae Radix in GXF have been reported to attenuate inflammation [11, 12]. In addition, Sophorae Flavescentis Radix has also been reported to have anti-fibrotic effects [13]. GXF has been commonly used in Chongqing region and exerts remarkable clinical effects on hepatic protection. After a large number of clinical observations, we found that 70.3% of patients with liver cirrhosis, caused by hepatitis B and alcoholic liver disease, and Child–Pugh class A or B can be reversed after GXF treatment. For patients with Child–Pugh class C, GXF treatment can stabilize liver function and reduce portal hypertension. A large cohort of clinical study is currently being established to further scientifically confirm its anti-fibrosis efficacy. Although GXF has a great anti-fibrotic potential, the complex chemical compositions have blocked the systematic understanding of the material basis of classic TCM formulas. More importantly, it is difficult to progress its use on a wider scale due to its unclear molecular mechanism.
The point of this research is to verify the efficacy of GXF and explore the potential molecular mechanism on the premise of identifying the material basis of GXF through UPLC-MS/MS. The rat model involved an artificial liver injury was induced by CCl4, a selective hepatotoxic drug, which can trigger inflammation and fibrosis. This model simulated clinical liver fibrosis in morphology and pathophysiology [14,15,16]. The clinical dosage of GXF is 30 g/person/day, and used this as the standard to set the middle dose (3.15 g/kg) for rats and set high (4.73 g/kg) and low dose (1.58 g/kg) in this study. We explored the potential mechanism of GXF with our RNA sequencing (RNA-seq) data combined with recent scRNA-seq studies [17,18,19]. With the supplementary demonstration of basic experiments, the clinical application of GXF could be better promoted. Figure 1 shows the flow chart of the whole study.
Materials and methods
Reagents
Antibodies to CCL2, EGFR, PDGFRα, PDGFB and TNFRSF12a were purchased from Abcam China; CCR2, TGFβ and TNFSF12 were purchased from Affinity Biologicals, Inc; PDGFRβ and β-actin were purchased from Cell Signaling Technology; CD9 was purchased from Proteintech Group, Inc; whereas Trem2 was purchased from ABclonal Technology. Trizol reagents were purchased from Gibco BRL (NY, USA). Carbon tetrachloride, chloroform, methanol, and so forth were obtained from the Center of Equipment and Reagent (Chongqing Medical University). All the reagents were of analytical quality. FZHY was obtained from Shanghai HuangHai Pharmaceutical Co., Ltd (Shanghai, China).
Preparation of GXF
GXF formula was obtained from Chongqing Traditional Chinese Medicine Hospital. Specifically, GXF was composed of Huang Qin (Scutellariae Radix), Ku Shen (Sophorae Flavescentis Radix), Fang Feng (Saposhnikoviae Radix), Chen Pi (Citrus Reticulata) and Xi Yang Shen (Panacis Quinquefolii Radix), which were purchased from Beijing Kangrentang Pharmaceutical Company (China). According to the traditional water extraction method, the five herbs were soaked in water at a ratio of 6:4:2:1:1 for 60 min and then boiled. The extraction was filtered and concentrated for subsequent experiments.
UPLC-MS/MS analysis of GXF
GXF formula was analyzed using U3000 UPLC system equipped with Q Exactive Plus Orbitrap high-resolution mass spectrometer (Thermo Fisher Scientific). Chromatographic separation was performed on a ACQUITY UPLC HSS T3 column (2.1 × 100 mm, 1.8 μm, Waters). The column temperature was set at 35 ℃. The flow rate was 0.3 mL/min. The mobile phase consisted of deionized water with 0.1% formic acid (A) and acetonitrile with 0.1% formic acid (B). The mass spectrometer analysis was conducted in both positive and negative ion modes. The scan mass ratio was within the mass range of m/z 50–1500.
The analysis data were analyzed by Compound Discoverer software 3.2 (Thermo Fisher Scientific). The mass deviations of characteristic peak element matching, molecular formula prediction and isotope distribution matching were all set to within 5 ppm. Chemical identification was based on chromatographic elution behavior, mass spectrometry fragment patterns and mass spectrometry databases (MZ cloud, MZ vault).
Animals and experiment designs
Animals were purchased from the Laboratory Animal Center of Chongqing Medical University [SCXK(Yu)2018-0003]. The animal experiments were reviewed and approved by the Ethics Committee of Chongqing Traditional Chinese Medicine Hospital (2018-ky-48) and conformed to the Guidelines for the Care and Use of Laboratory Animals. 140 male Sprague–Dawley rats weighing 140 ~ 160 g were housed in a room at 22 ℃ with a 12-h light/dark cycle and free access to food and water. All rats acclimatized for 1 week before experiments. Liver fibrosis was induced by intraperitoneal (i.p.) injection of CCl4 dissolved in olive oil [1:1(v/v)] at a dose of 0.1 mL/100 g body weight (BW) twice-weekly for 8 weeks as shown in Fig. 2A. Rats were randomly divided into seven groups (n = 20 per group): (1–2) Control and GXF high-dose control groups were i.p. injected with olive oil and given 0.9% NaCl or 4.73 g/kg GXF daily respectively (p.o.); (3) CCl4 group; (4–6) CCl4 + GXF group (1.58, 3.15, 4.73 g/kg respectively, p.o.); (7) CCl4 + FZHY group (0.47 g/kg, p.o.). GXF and FZHY applied once daily during 5–8 weeks. Rats were anesthetized and sacrificed 48 h after the final CCl4 (or olive oil) injection. The livers were frozen and stored in liquid nitrogen or fixed in 4% buffered paraformaldehyde until further analysis. Serum was stored at − 80 ℃ for further analysis.
Detection of serum ALT, AST, Scr and CCL2 levels
Serum alanine aminotransferase (ALT), aspartate aminotransferase (AST) and Scr were measured by routine enzymatic assays kits (Maccura, Chengdu, China) using HITACHI Clinical Analyzer 7600 (Hitachi, Tokyo, Japan). The detection of these indicators was completed at the Clinical Laboratory Center of the Second Affiliated Hospital of Chongqing Medical University. The serum CCL2 level was measured by sandwich enzyme-linked immunosorbent assays (ELISA) using an ELISA kit according to the manufacturer’s instructions (Raybiotech, GA, USA).
Determination of hepatic hydroxyproline
Hydroxyproline (Hyp) contents of liver tissues were determined using an alkaline hydrolysis method described by the hydroxyproline test kit (Jiancheng Biotech, Jiangsu, China).
Fibrosis grade assessment
Grading of liver fibrosis was assessed by collagen stained with Sirius red, all rat livers treated for 8 weeks were evaluated. The fibrosis stage was evaluated according to the Ishak grading system on a scale from F0 to F6, higher score indicates greater fibrosis. F0: No fibrosis, F1: Fibrous expansion of some portal areas, F2: Fibrous expansion of most portal areas, F3: Fibrous expansion of most portal areas with portal to portal bridging, F4: Fibrous expansion of portal areas with marked bridging, F5: Marked bridging with occasional nodules (incomplete cirrhosis), F6: Cirrhosis, probable or definite [20]. Histological assessment was performed by two clinical pathologists without knowledge of the experimental design.
RNA-sequencing and bioinformatic analyses
RNA-Seq analysis was based on biological replicates with 8 rats in each experimental group. Library construction, RNA-seq and data analysis were performed by Novogene (Beijing, China). Total RNA was extracted from rat liver tissues with Trizol reagent (Invitrogen, USA). All samples passed RNA quality control examined by Agilent Bioanalyzer 2100, and mRNA libraries were sequenced using a 2 × 150 bp paired-end method by illumina novaseq 6000. The raw data were trimmed, filtered and qualified using FASTX (http://hannonlab.cshl.edu/fastx_toolkit/). Clean data was aligned to rat reference genome (rn5) using Hisat2 [21]. Feature Counts v1.5.0-p3 [22] was used to count the reads numbers mapped to each gene. And then FPKM [23] of each gene was calculated based on the length of the gene and reads count mapped to this gene. DESeq2 [24] was used to estimate the significance of differential genes (DEGs) between any two experimental groups according to the criteria of |log2(fold change) |> 1 and p-adjust < 0.05. KEGG is a database resource for understanding high-level functions and utilities of the biological system (http://www.genome.jp/kegg/). We used cluster Profiler R package to test the statistical enrichment of differential expression genes in KEGG pathways.
Immunohistochemistry
Immunohistochemistry for Trem2, CD9, PDGFRα and PDGFRβ were performed using commercial kit (ZSGB-BIO, China) according to the manufacturer’s instructions. In brief, deparaffinized and blocked 5 μM sections were incubated with anti-Trem2 antibody (1:100, A10482), anti-CD9 antibody (1:3000, 60232–1-Ig), anti-PDGFR alpha (PDGFRα) antibody (1:200, ab203491) and anti-PDGFR β antibody (1:100, 3169). The sections were then further incubated with the HRP-conjugated secondary antibody for 1 h at room temperature. After washing with PBS, 3,3'-adiaminobenzidine tetrahydrochloride (DAB) coloration signal was used. Images were acquired on a Slide scanning system Pannoramic MIDI II (3D HISTECH). Certified 3D HISTECH CaseViewer and imageJ software were applied to analyze the results.
RT-qPCR and western blot analysis
Total RNA was extracted from HSC-T6 cells or rat liver tissues with RNApure total RNA fast isolation kit (BioTeke, China) following the manufacturer’s instructions. The cDNAs were synthesized with the commercial kit (Takara, Japan). Gene expressions were measured by qPCR with CFX Connect™ Real-Time PCR System (Bio-Rad, USA). GAPDH was used as an internal control and the relative expression levels of mRNA were calculated using the 2−ΔΔCt method. The primer pairs used in the experiments were listed in Additional file 1: Table S1.
Total protein was extracted from pieces of liver tissues using whole cell lysis assay (KeyGen Biotech, Jiangsu, China). Protein concentration was determined by the BCA Protein Assay Kit (KeyGen Biotech, Jiangsu, China). Quantified proteins were separated on SDS-PAGE and transferred to PVDF membranes (Millipore Corporation, USA). After blocking, membranes were incubated with primary antibodies against CCL2 (1:1000, ab7202), anti-CCR2 (1:500, DF2711), TGFβ (1:1500, AF1027), PDGF-B (1:1500, ab178409), TNFSF12 (1:500, DF7444), EGFR (1:2000, ab52894), PDGFRα (1:1000, ab203491) and TNFRSF12a (1:1500, ab109365) at 4 °C overnight. The anti-β-actin (1:2000, 4970L) was probed as an internal control. Then, membranes were washed with TBST, and incubated with secondary antibodies for 2 h at room temperature. Protein bands were visualized by using ECL (Advansta, USA). Levels of target protein band densities were analyzed with ChemicDoc™ MP Imaging System (Bio-Rad, USA).
Cell culture, viability assay and treatment
HSC-T6 cells were cultured in DMEM (supplemented with 10% FBS and 1% P/S) at 37 ℃ in a 5% CO2 humidified atmosphere. Cell viability of HSC-T6 was assessed by cell counting kit-8 (CCK-8) assay (APExBIO Technology, USA). HSC-T6 cells (1 × 104 cells per well) were seeded in 96-well plates, then cells were treated with different doses of baicalin, matrine or hesperidin (6.25, 12.5, 25, 50, 100 and 200 µM) for 24 h. Then, a total 10 µL of CCK-8 reagent was added to each well. After incubation for 1 h, the absorbance at 450 nm of each well was measured by a microplate reader. For TGFβ1-induced HSCs activation, HSCs were stimulated with 10 ng/mL TGFβ1 (Peprotech, USA) for 24 h. For treatment groups, HSCs were treated with TGFβ1 (10 ng/mL) supplemented with 50 μM baicalin, matrine or hesperidin (MCE, USA) for 24 h.
Statistical analysis
GraphPad Prism 9.0 (GraphPad Software Inc., California, USA) and SPSS 25.0 (SPSS Inc., IL, USA) were used for statistical analysis. Continuous variables are presented as the mean ± SEM. The differences between two groups were assessed with Student's t-test. Statistical differences between groups were evaluated by one-way analysis of variance (ANOVA), and Ishak scores were performed with rank-sum test. A P-value of p < 0.05 was considered statistically significant.
Results
Identification of compounds in GXF
UPLC-MS/MS was employed to identify the main 10 compounds in GXF formula: baicalin, wogonoside, matrine, stachydrine, 5-O-Methylvisammioside, hesperidin, ursodeoxycholic acid, amygdalin, sophocarpine and paeoniflorin. The detailed results were presented in Table 1 (Additional file 2: Table S2).
GXF improved CCl4-induced liver function in rats
The process of animal experiment is shown in Fig. 2A. The body weight of rats was weighted once a week during the experiment process. The weight changes of the rats every 2 weeks were calculated, although the weight changes of the control group and the GXF control group were greater than the others, there was no statistical difference in the weight changes between the CCl4 group and the treatment groups (Fig. 2B). The body weight and liver weight of rats at the eighth week showed that GXF had no effect on the weight and liver weight/body weight ratio (Fig. 2C). Serum ALT, AST and Scr levels were used to assess liver and renal function. There was no significant difference in serum Scr level within each group, indicating that GXF had no effect on the renal function. Compared with control group, the serum ALT and AST levels of control drug group had no significant change, while CCl4 group had a significantly increase (p < 0.001). Further, the liver injury induced by CCl4 was reversed to varying degrees by GXF (1.58, 3.15, 4.73 g/kg/day) and FZHY (0.47 g/kg/day). Among the treatment groups, the serum ALT and AST levels of the GXF high-dose (4.73 g/kg/day) group had the most significant decrease compared with the CCl4 group (p < 0.001). In addition, the ALT and AST levels of the GXF high-dose group were lower than FZHY group (p < 0.05), indicating that its protective effect on liver is better than FZHY (Fig. 2D).
GXF alleviated liver histopathological changes
We observed that the liver surface of CCl4 and GXF low-dose groups were significantly rougher compared with the control group, whereas that of GXF high-dose and FZHY groups were significantly improved (Fig. 3A). HE staining showed that CCl4 caused collagen deposition, structural destruction and inflammatory cell infiltration in liver lobules and portal veins, while GXF high-dose and FZHY treatment alleviated these histopathological damages (Fig. 3B). Similarly, Sirius red staining showed that the collagen deposition in the CCl4 group was significantly increased compared with the control group, and collagen deposition could be reduced after treatment (Fig. 3C). The reduction of collagen in GXF high-dose and FZHY groups were the most significant, which could be more intuitively reflected by the panoramic scanning (Fig. 3D). According to the analysis of the ISHAK scoring standard, the incidence of liver cirrhosis in both GXF medium-dose and high-dose decreased from 76.5 to 55.6%. Although the F5/F6 scores of FZHY group accounted for the lowest proportion, the F6 scores of GXF high-dose group accounted for the lowest (Fig. 3E). The CPA and Hyp content of GXF high-dose group were the lowest and significantly lower than CCl4 group (p < 0.001 or p < 0.01) (Fig. 3F, G). Finally, the relative mRNA expression of collagen I and α-SMA in liver tissues of each treatment group were significantly lower than that of CCl4 group (p < 0.001 or p < 0.01) (Fig. 3H).
KEGG analysis revealed DEGs enriched in chemokine signaling pathway
To investigate the molecule mechanisms underlying the protective effect of GXF against CCl4-induced liver fibrosis in rats, hepatic gene expression profiles of control, CCl4 and GXF high-dose group were detected by RNA-seq (n = 8). 6636 DEGs were identified in CCl4 group compared with control group, 4394 genes of which were up-regulated and 2242 were down-regulated. Similarly, 2384 DEGs were identified in GXF high-dose group compared to CCl4 group, 829 genes of which were up-regulated and 2005 were down-regulated (|Foldchange|> 2, p-adjust < 0.05) (Fig. 4A). The KEGG pathway analyses of the difference genes between CCl4 and control group and that between GXF and CCl4 group were showed in Fig. 4B. 2392 genes were obtained by intersecting the DEGs between the CCl4 and control group and the DEGs between GXF and CCl4 group, and the KEGG analyses of these genes were showed in Fig. 4C (Additional file 3: Table S3). The results showed that these differential genes could be significantly enriched in the ECM-receptor interaction pathway, suggesting that GXF treatment could effectively inhibit collagen deposition, and the heat map results of this pathway revealing that GXF could significantly down-regulate the expression of Col1a1, Col1a2, Col4a1 and Col6a1 (Fig. 4D). In addition, the DEGs could also be significantly enriched in the chemokine signaling pathway. The heat map showed that GXF could significantly inhibit the expression of some chemokine-related genes such as CCL2, CXCL6 and CX3CL1 (Fig. 4E), suggesting that GXF may alleviate CCl4-induced liver fibrosis by inhibiting the recruitment of related immune cells.
GXF down-regulated CCL2 both in liver tissues and peripheral blood
According to the RNA-seq results, the expression levels of genes enriched in chemokine pathway were summarized, which showed in Fig. 5A. Between the CCl4 group and GXF group, there was no significant difference in expression levels of genes related to the recruitment of T cells and monocyte. However, compared with the CCl4 group, GXF treatment could significantly inhibit the expression of genes which related to macrophages such as CCL2, CX3CL1, CX3CR1, CXCL16 and CXCR6 (p < 0.05). The relative mRNA expression of CCL2, CCR2, CX3CL1, CX3CR1, CXCL16 and CXCR6 in liver tissues was detected by RT-qPCR. As shown in Fig. 5B, CCl4 intervention could significantly increase the expression of these genes, but it could be suppressed by GXF, of which CCL2 was the most significant (p < 0.05). Compared with GXF, FZHY had no obvious inhibitory effect on these genes. The western blot results showed that the protein expression of CCR2 did not differ between the groups, but the expression of CCL2 in the CCl4 group was significantly higher than that in the control group, which could be reduced after GXF treatment (p < 0.05) (Fig. 5C). In addition, the expression of CCL2 in serum was consistent with that in liver tissues (Fig. 5D). Therefore, GXF may delay the progression of liver fibrosis by inhibiting the recruitment of macrophages from peripheral blood to liver tissues.
GXF remodeled the fibrotic niche
Single-cell sequencing (scRNA-seq) studies have already provided a wealth of novel insights into cellular heterogeneity and the interactome present across the different cell types in liver fibrosis progression [18]. The interaction between macrophages, stellate cells and endothelial cells is closely related to liver fibrosis [19]. To interrogate how GXF attenuates fibrosis and to deepen our understanding of the mechanism, we integrated our total liver RNA-seq data with recent liver scRNA-seq studies. According to the RNA-seq results, the gene expression levels of genes related to scar-associated macrophages (SAMacs), scar-associated mesenchymal cells (SAMes) and scar-related endothelial (SAEndo) were summarized in Fig. 6A–C. Compared with the CCl4 group, RT-qPCR results showed that the unique markers of SAMacs such as Trem2 and CD9 were down-regulated by GXF (p < 0.05) (Fig. 6D). More intuitively, the number of Trem2+CD9+ macrophages were increased in the CCl4 and FZHY groups compared with the control group, while it was not observed in the GXF group (p < 0.05 or p < 0.01) (Fig. 6F). The distinct genes including PDGFRα and PDGFRβ characterized to SAMes were also obviously reduced by GXF (p < 0.01 or p < 0.001) (Fig. 6E), and this trend was highly consistent with the results of immunohistochemistry (p < 0.01) (Fig. 6G). Whereas, there was no obvious effect on the specific molecules of SAEndo subpopulations including Cd34, Ackr1, Csf1, Flt1 and Tek (Fig. 6C). Therefore, we speculated that GXF treatment could reduce the proliferation and activation of SAMacs and SAMes, but had little effect on SAEndo.
GXF interfered the interaction between macrophages and stellate cells
Ligands and receptors bridge the interaction between cells. TGFβ/EGFR, PDGFB/PDGFRα and TNFSF12/TNFRSF12a are important bonds between macrophages and stellate cells. We found that GXF could effectively inhibit the expression of PDGFB and TNFSF12 on SAMacs and its cognate receptor PDGFRα and TNFRSF12a expressed on SAMes (p < 0.05 or p < 0.01) (Fig. 7C–F). Although GXF had no effect on the expression of EGFR compared with the CCl4 group, it could significantly inhibit the expression of TGFβ in liver tissues (p < 0.05) (Fig. 7A, B). In addition, western blot results showed that FZHY also had a certain inhibitory effect on the expression of the above ligands and receptors, indicating that FZHY also inhibited the activation of stellate cells by interfering the interaction between macrophages and stellate cells. Overall, these data suggested that the anti-fibrosis effect of GXF is closely related to its inhibition expansion of scar-associated Trem2+CD9+ macrophage subpopulation and PDGFRα+ PDGFRβ+ scar-producing myofibroblasts in the damaged liver, and its remodeling of the fibrotic niche through regulation ligand-receptor interactions between SAMacs and SAMes including TGFβ/EGFR, PDGFB/PDGFRα and TNFSF12/TNFRSF12a signaling.
Verified potential pharmacodynamic substances of GXF
Based on the main compounds of GXF identified by UPLC-MS/MS, we assessed the potential pharmacodynamic substances on TGFβ1-activated HSC-T6 cells. Results of CCK-8 assay indicated that baicalin, matrine and hesperidin, with a concentration no higher than 50 μM, had no significant inhibitory effect on the viability of HSC-T6 cells (Fig. 8A–C). Next, activation of HSCs could be effectively reversed under the intervention of these three drugs at a concentration of 50 μM, as indicated by inhibited mRNA levels of TGFβ, α-SMA and collagen I (Fig. 8D–F). In vitro experiments demonstrated that baicalin, matrine and hesperidin inhibited the activation of HSCs, which deserves further attention.
Discussion
As there are no effective drugs so far, traditional Chinese medicine has advantages in the treatment of liver fibrosis [25, 26]. In this study, we used UPLC-MS/MS to identify the compound composition of GXF to deepen our understanding of the material basis of classic TCM formulas. Next, we have confirmed through animal experiments that the GXF high-dose has a good therapeutic effect on CCl4-induced liver fibrosis, which is no less effective than FZHY. Integration of our RNA-seq analysis and recent scRNA-seq data showed that GXF significantly downregulated CCL2 to inhibit the expansion of scar-associated Trem2+CD9+ macrophages subpopulation and PDGFRα+PDGFRβ+ scar-producing myofibroblasts in the damaged liver. Mechanisms of GXF inhibiting liver fibrosis are summarized in Fig. 9.
UPLC-MS/MS provides technical support for analyzing the compound composition of traditional Chinese medicine. In this study, we identified 10 main compounds of GXF, of which baicalin, matrine, hesperidin and paeoniflorin have been reported to have anti-inflammatory or antioxidant activity. Baicalin is a flavonoid compound with antibacterial and anti-inflammatory effects [27]. It has been reported that baicalin can effectively reduce the level of alanine aminotransferase during active hepatitis [28]. Matrine is also a compound with anti-inflammatory activity, and it has been reported to have therapeutic potential in chronic bronchitis and enteritis [29, 30]. In addition, studies have shown that hesperidin and paeoniflorin have certain antioxidant activities [31, 32]. In vitro experiments proved that baicalin, matrine and hesperidin inhibited the activation of HSCs in this study. In brief, the therapeutic effect of GXF on liver fibrosis is determined by multiple active compounds.
Liver pathology analysis is the gold standard for evaluating the efficacy of anti-fibrotic drugs. The qualitative analysis of histopathological results was performed by the ISHAK scoring standard. Although the proportion of F5/F6 in GXF high-dose (55.6%) was not lower than FZHY (47.4%), but the proportion of F6 in GXF high-dose (5.6%) was the lowest among the treatment groups. In addition, the F2/F3 scores of GXF high-dose accounted for the highest proportion (27.8%), revealing that GXF high-dose had the best reversal effect on CCl4-induced liver fibrosis. Hydroxyproline is also an indicator of collagen deposition [33]. The hydroxyproline content of GXF high-dose (0.39 μg/mg) was lower than the CCl4 group (0.53 μg/mg) and the FZHY group (0.43 μg/mg), which further proved the effectiveness of GXF high-dose from a quantitative perspective. We believe that GXF high-dose may have a better delaying effect on the progression of advanced cirrhosis.
As we know, hepatic inflammation is a major contributor to the pathogenesis of almost all liver diseases [34]. Low-molecular-weight proteins called chemokines such as CCL2, CX3CL1, CXCL16 and CCL5 are the main drivers of liver infiltration by immune cells such as macrophages, neutrophils and others during an inflammatory response [35, 36]. Chemokine CCL2/CCR2 signaling plays an important role in the process of fibrosis [37,38,39]. CCL2 is secreted by hepatic stellate cells, hepatocytes, biliary epithelial cells and Kupffer cells. CCR2, the only known receptor for CCL2, is expressed on monocytes and macrophages in the liver [37]. Following the liver injury induced by CCl4, there is recruitment of circulating LY6ChiCCR2+ monocytes that differentiate into liver monocyte-derived macrophages, resulting in a huge expansion of intrahepatic macrophages [40, 41]. CXCL16 is produced by macrophages, which is a survival and maturation factor for hepatic NKT cells, resulting in an accumulation of NKT cells at sites of injury [42]. CX3CL1 is widely expressed on immune and non-immune cells, and CX3CL1/CX3CR1 signaling has been shown to promote macrophage survival [43, 44]. In this study, compared with the CCl4 group and the FZHY group, GXF could inhibit the expression of CCL2 in peripheral blood and liver tissue, thereby inhibiting the homing of macrophages. Although FZHY had no effect on the expression of CCL2, it could inhibit the expression of CX3CL1 and CXCR6, indicating that FZHY acted on macrophages in a different mechanism.
Single-cell transcriptome technologies are transforming our understanding of cell diversity and disease pathogenesis. Liver fibrosis involves a complex interplay between multiple non-parenchymal cell (NPC) lineages including SAMacs, SAMes and SAEndo spatially located within areas of scarring, termed the fibrotic niche [18]. The monocyte-derived macrophages (MDMs), which recruited by CCL2/CCR2 signaling, have been showed to regulate a number of aspects of liver injury, including inflammation and fibrosis [45, 46]. A distinct subpopulation of Trem2+CD9+ scar-associated macrophages (SAMacs) were identified in these MDMs by scRNAseq. The SAMacs, which expand in fibrotic livers, are spatially localized to areas of scarring (termed the fibrotic niche) and promote the activation of quiescent HSCs into PDGFRα+ SAMes, which proliferate and produce fibrillar collagens within the fibrotic niche of diseased livers [19]. In addition, due to the complex disease process of liver fibrosis, the interaction between different cells may be an important driving force for disease progression. Specifically, SAMacs express ligands including TGFβ, TNFSF12 and PDGFB that could signal through their cognate receptors EGFR, TNFRSF12a and PDGFRα expressed on SAMes to potentially promote mesenchymal cell activation [47,48,49] and proliferation [50, 51]. In this study, both GXF and FZHY could inhibit the above-mentioned ligand-receptor interactions between SAMacs and SAMes. In general, the application of multiomic approaches helps to deepen our new understanding of the therapeutic targets of TCM formulas.
Conclusion
In this study, we identified the main compounds of GXF formula by UPLC-MS/MS. We verified that GXF had an anti-fibrosis effect by inhibiting the recruitment of macrophages and improving the microenvironment of hepatic stellate cells. In conclusion, this study reveals potential anti-fibrotic mechanism of GXF formula, and provides scientific basis for further studies and clinical applications.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Abbreviations
- GXF:
-
Ganxianfang
- TCM:
-
Traditional Chinese medicine
- CCl4 :
-
Carbon tetrachloride
- RNA-Seq:
-
RNA sequencing
- scRNA-seq:
-
Single-cell sequencing
- ALT:
-
Alanine aminotransferase
- AST:
-
Aspartate aminotransferase
- Scr:
-
Serum creatinine
- FZHY:
-
Fuzheng Huayu
- Hyp:
-
Hydroxyproline
- CPA:
-
Collagen proportion area
- α-SMA:
-
α-Smooth muscle actin
- KEGG:
-
Kyoto Encyclopedia of Genes and Genomes
- DEGs:
-
Differentially expressed genes
- ECM:
-
Extracellular matrix
- CCL2:
-
C–C motif chemokine ligand 2
- CCR2:
-
C–C motif chemokine receptor 2
- SAMacs:
-
Scar-associated macrophages
- SAMes:
-
Scar-associated mesenchymal cells
- SAEndo:
-
Scar-related endothelial
- HSCs:
-
Hepatic stellate cells
- Trem2:
-
Triggering receptor expressed on myeloid cells 2
- PDGFRα:
-
Platelet derived growth factor receptor alpha
- PDGFRβ:
-
Platelet derived growth factor receptor beta
- PDGFB:
-
Platelet derived growth factor subunit B
- TGFβ:
-
Transforming growth factor beta
- EGFR:
-
Epidermal growth factor receptor
- TNFSF12:
-
TNF superfamily member 12
- TNFRSF12a:
-
TNF receptor superfamily member 12a
References
Kisseleva T, Brenner D. Molecular and cellular mechanisms of liver fibrosis and its regression. Nat Rev Gastroenterol Hepatol. 2021;18:151–66.
Lee S-J. Mechanisms of fibrogenesis in liver cirrhosis: the molecular aspects of epithelial-mesenchymal transition. WJH. 2014;6:207.
Henderson NC, Rieder F, Wynn TA. Fibrosis: from mechanisms to medicines. Nature. 2020;587:555–66.
Schuppan D, Pinzani M. Anti-fibrotic therapy: lost in translation? J Hepatol. 2012;56:S66-74.
Dong S, Su S. Advances in mesenchymal stem cells combined with traditional Chinese medicine therapy for liver fibrosis. J Integr Med. 2014;12:147–55.
Zhang L, Schuppan D. Traditional Chinese Medicine (TCM) for fibrotic liver disease: hope and hype. J Hepatol. 2014;61:166–8.
Luk JM, Wang X, Liu P, Wong K-F, Chan K-L, Tong Y, et al. Traditional Chinese herbal medicines for treatment of liver fibrosis and cancer: from laboratory discovery to clinical evaluation. Liver Int. 2007;27:879–90.
Xu L, Liu P. Guidelines for diagnosis and treatment of hepatic fibrosis with integrated traditional Chinese and Western medicine (2019 edition). J Integr Med. 2020;18:203–13.
Liu H, Lv J, Zhao Z, Xiong A, Tan Y, Glenn JS, et al. Fuzhenghuayu decoction ameliorates hepatic fibrosis by attenuating experimental sinusoidal capillarization and liver angiogenesis. Sci Rep. 2019;9:18719.
Zhang M, Liu H, Huang K, Peng Y, Tao Y, Zhao C, et al. Fuzheng Huayu recipe prevented and treated CCl4-Induced mice liver fibrosis through regulating polarization and chemotaxis of intrahepatic macrophages via CCL2 and CX3CL1. Evid Based Complementary Altern Med. 2020;2020:8591892.
Fu Y, Xu B, Huang S, Luo X, Deng X, Luo S, et al. Baicalin prevents LPS-induced activation of TLR4/NF-κB p65 pathway and inflammation in mice via inhibiting the expression of CD14. Acta Pharmacol Sin. 2021;42:88–96.
Yao L, Wang S, Wei P, Bao K, Yuan W, Wang X, et al. Huangqi-Fangfeng protects against allergic airway remodeling through inhibiting epithelial–mesenchymal transition process in mice via regulating epithelial derived TGF-β1. Phytomedicine. 2019;64:153076.
Chen L, Lv D, Wang D, Chen X, Zhu Z, Cao Y, et al. A novel strategy of profiling the mechanism of herbal medicines by combining network pharmacology with plasma concentration determination and affinity constant measurement. Mol BioSyst. 2016;12:3347–56.
Kang J-W, Hong J-M, Lee S-M. Melatonin enhances mitophagy and mitochondrial biogenesis in rats with carbon tetrachloride-induced liver fibrosis. J Pineal Res. 2016;60:383–93.
Ogaly HA, Eltablawy NA, Abd-Elsalam RM. Antifibrogenic influence of Mentha piperita L. essential oil against CCl 4-induced liver fibrosis in rats. Oxid Med Cell Longev. 2018;2018:1–15.
Liu J, Kong D, Qiu J, Xie Y, Lu Z, Zhou C, et al. Praziquantel ameliorates CCl 4-induced liver fibrosis in mice by inhibiting TGF-β/Smad signalling via up-regulating Smad7 in hepatic stellate cells. Br J Pharmacol. 2019;176:4666–80.
Xiong X, Kuang H, Ansari S, Liu T, Gong J, Wang S, et al. Landscape of intercellular crosstalk in healthy and NASH liver revealed by single-cell secretome gene analysis. Mol Cell. 2019;75:644-660.e5.
Ramachandran P, Matchett KP, Dobie R, Wilson-Kanamori JR, Henderson NC. Single-cell technologies in hepatology: new insights into liver biology and disease pathogenesis. Nat Rev Gastroenterol Hepatol. 2020;17:457–72.
Ramachandran P, Dobie R, Wilson-Kanamori JR, Dora EF, Henderson BEP, Luu NT, et al. Resolving the fibrotic niche of human liver cirrhosis at single-cell level. Nature. 2019;575:512–8.
Ishak’ K, Baptista A, Grootes JD, Gudat F, Denk H, Desmet V, et al. Histological grading and staging of chronic hepatitis. J Hepatol. 1995;22:696–9.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5:621–8.
Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30:923–30.
Bray NL, Pimentel H, Melsted P, Pachter L. Near-optimal probabilistic RNA-seq quantification. Nat Biotechnol. 2016;34:525–7.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.
Li H. Advances in anti hepatic fibrotic therapy with Traditional Chinese Medicine herbal formula. J Ethnopharmacol. 2020;251:112442.
Zhao T, Mao L, Yu Z, Hui Y, Feng H, Wang X, et al. Therapeutic potential of bicyclol in liver diseases: lessons from a synthetic drug based on herbal derivative in traditional Chinese medicine. Int Immunopharmacol. 2021;91:107308.
Dong S, Zhong Y, Lu W, Li G, Jiang H, Mao B. Baicalin inhibits lipopolysaccharide-induced inflammation through signaling NF-κB pathway in HBE16 airway epithelial cells. Inflammation. 2015;38:1493–501.
Chen Y, Yang Y, Wang F, Yang X, Yao F, Ming K, et al. Antiviral effect of baicalin phospholipid complex against duck hepatitis A virus type 1. Poult Sci. 2018;97:2722–32.
Suo Z, Liu Y, Ferreri M, Zhang T, Liu Z, Mu X, et al. Impact of matrine on inflammation related factors in rat intestinal microvascular endothelial cells. J Ethnopharmacol. 2009;125:404–9.
Sun D, Wang J, Yang N, Ma H. Matrine suppresses airway inflammation by downregulating SOCS3 expression via inhibition of NF-κB signaling in airway epithelial cells and asthmatic mice. Biochem Biophys Res Commun. 2016;477:83–90.
Ali AM, Gabbar MA, Abdel-Twab SM, Fahmy EM, Ebaid H, Alhazza IM, et al. Antidiabetic potency, antioxidant effects, and mode of actions of Citrus reticulata fruit peel hydroethanolic extract, hesperidin, and quercetin in nicotinamide/streptozotocin-induced Wistar diabetic rats. Oxid Med Cell Longev. 2020;2020:1730492.
Mao Q-Q, Zhong X-M, Feng C-R, Pan A-J, Li Z-Y, Huang Z. Protective effects of paeoniflorin against glutamate-induced neurotoxicity in PC12 cells via antioxidant mechanisms and Ca2+ antagonism. Cell Mol Neurobiol. 2010;30:1059–66.
Zhou Y, Wu R, Cai F-F, Zhou W-J, Lu Y-Y, Zhang H, et al. Xiaoyaosan decoction alleviated rat liver fibrosis via the TGFβ/Smad and Akt/FoxO3 signaling pathways based on network pharmacology analysis. J Ethnopharmacol. 2021;264:113021.
Matyas C, Haskó G, Liaudet L, Trojnar E, Pacher P. Interplay of cardiovascular mediators, oxidative stress and inflammation in liver disease and its complications. Nat Rev Cardiol. 2021;18:117–35.
Marra F, Tacke F. Roles for chemokines in liver disease. Gastroenterology. 2014;147:577-594.e1.
Wen Y, Lambrecht J, Ju C, Tacke F. Hepatic macrophages in liver homeostasis and diseases-diversity, plasticity and therapeutic opportunities. Cell Mol Immunol. 2021;18:45–56.
Sahin H, Trautwein C, Wasmuth HE. Functional role of chemokines in liver disease models. Nat Rev Gastroenterol Hepatol. 2010;7:682–90.
Mandrekar P, Ambade A, Lim A, Szabo G, Catalano D. An essential role for monocyte chemoattractant protein-1 in alcoholic liver injury: regulation of proinflammatory cytokines and hepatic steatosis in mice. Hepatology. 2011;54:2185–97.
Baeck C, Wehr A, Karlmark KR, Heymann F, Vucur M, Gassler N, et al. Pharmacological inhibition of the chemokine CCL2 (MCP-1) diminishes liver macrophage infiltration and steatohepatitis in chronic hepatic injury. Gut. 2012;61:416–26.
Pellicoro A, Ramachandran P, Iredale JP, Fallowfield JA. Liver fibrosis and repair: immune regulation of wound healing in a solid organ. Nat Rev Immunol. 2014;14:181–94.
Krenkel O, Tacke F. Liver macrophages in tissue homeostasis and disease. Nat Rev Immunol. 2017;17:306–21.
Wehr A, Baeck C, Heymann F, Niemietz PM, Hammerich L, Martin C, et al. Chemokine receptor CXCR6-dependent hepatic NK T cell accumulation promotes inflammation and liver fibrosis. J Immunol. 2013;190:5226–36.
Castro-Sánchez S, García-Yagüe ÁJ, Kügler S, Lastres-Becker I. CX3CR1-deficient microglia shows impaired signalling of the transcription factor NRF2: implications in tauopathies. Redox Biol. 2019;22:101118.
Karlmark KR, Zimmermann HW, Roderburg C, Gassler N, Wasmuth HE, Luedde T, et al. The fractalkine receptor CX3CR1 protects against liver fibrosis by controlling differentiation and survival of infiltrating hepatic monocytes. Hepatology. 2010;52:1769–82.
Wu R, Nakatsu G, Zhang X, Yu J. Pathophysiological mechanisms and therapeutic potentials of macrophages in non-alcoholic steatohepatitis. Expert Opin Ther Targets. 2016;20:615–26.
Krenkel O, Puengel T, Govaere O, Abdallah AT, Mossanen JC, Kohlhepp M, et al. Therapeutic inhibition of inflammatory monocyte recruitment reduces steatohepatitis and liver fibrosis: metabolic liver disease. Hepatology. 2018;67:1270–83.
Minutti CM, Modak RV, Macdonald F, Li F, Smyth DJ, Dorward DA, et al. A macrophage-pericyte axis directs tissue restoration via amphiregulin-induced transforming growth factor beta activation. Immunity. 2019;50:645-654.e6.
McKee C, Sigala B, Soeda J, Mouralidarane A, Morgan M, Mazzoccoli G, et al. Amphiregulin activates human hepatic stellate cells and is upregulated in non alcoholic steatohepatitis. Sci Rep. 2015;5:8812.
Chen J, Chen J-K, Nagai K, Plieth D, Tan M, Lee T-C, et al. EGFR signaling promotes TGF β-dependent renal fibrosis. JASN. 2012;23:215–24.
Wilhelm A, Shepherd EL, Amatucci A, Munir M, Reynolds G, Humphreys E, et al. Interaction of TWEAK with Fn14 leads to the progression of fibrotic liver disease by directly modulating hepatic stellate cell proliferation: TWEAK as a regulator of liver fibrogenesis. J Pathol. 2016;239:109–21.
Makino K, Makino T, Stawski L, Mantero JC, Lafyatis R, Simms R, et al. Blockade of PDGF receptors by crenolanib has therapeutic effect in patient fibroblasts and in preclinical models of systemic sclerosis. J Investig Dermatol. 2017;137:1671–81.
Acknowledgements
Not applicable.
Funding
This work was supported by the National Science and Technology Major Project of China (Grant Numbers 2017ZX10202203-008, 2017ZX10202203-007, 2018ZX10302206-003), and a pilot project of clinical cooperation between traditional Chinese and western medicine for significant and complicated diseases of National Administration of Traditional Chinese Medicine: hepatic fibrosis.
Author information
Authors and Affiliations
Contributions
ZL and HX performed animal experiments and pathway analysis. DX and SX conducted the ELISA and qRT-PCR assay. JX and ZL performed western blot and immunohistochemical experiments. ZL analyzed the experimental data and wrote the manuscript. HR, PH, HL and MP applied for designed the experimental protocols and directed the manuscript writing. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
The animal study was reviewed and approved by the Ethics Committee of Chongqing Traditional Chinese Medicine Hospital (2018-ky-48).
Consent for publication
Not applicable.
Competing interests
The authors declare no conflict of interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Table S1.
Primer sequences of RT-PCR.
Additional file 2: Table S2.
Detailed information and detection map of the main compounds.
Additional file 3: Table S3.
KEGG pathways of DEGs with opposite trends in GXF and CCl4 group.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Liu, Z., Xiang, H., Xiang, D. et al. Revealing potential anti-fibrotic mechanism of Ganxianfang formula based on RNA sequence. Chin Med 17, 23 (2022). https://doi.org/10.1186/s13020-022-00579-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13020-022-00579-7