Abstract The Traditional Chinese Medicine compound preparation known as Diwu Yanggan capsule (DWYG) can effectively hinder the onset and progression of hepatocellular carcinoma (HCC), which is recognized worldwide as a significant contributor to fatalities associated with cancer. Nevertheless, the precise mechanisms implicated have remained ambiguous. In present study, the model of HCC was set up by the 2-acetylaminofluorene (2-AAF)/partial hepatectomy (PH) in rats. To confirm the differentially expressed genes (DEGs) identified in the microarray analysis, real-time quantitative reverse transcription PCR (qRT-PCR) was conducted. In the meantime, the liquid chromatography-quadrupole time of flight mass spectrometry (LC-QTOF-MS/MS) was employed to characterize the component profile of DWYG. Consequently, the DWYG treatment exhibited the ability to reverse 51 variation genes induced by 2-AAF/PH. Additionally, there was an overlap of 54 variation genes between the normal and model groups. Upon conducting RT-qPCR analysis, it was observed that the expression levels of all genes were increased by 2-AAF/PH and subsequently reversed after DWYG treatment. Notably, the fold change of expression levels for all genes was below 0.5, with 3 genes falling below 0.25. Moreover, an investigation was conducted to determine the signaling pathway that was activated/inhibited in the HCC group and subsequently reversed in the DWYG group. Moreover, the component profile of DWYG encompassed a comprehensive compilation of 206 compounds that were identified or characterized. The findings of this study elucidated the potential alleviative mechanisms of DWYG in the context of HCC, thereby holding significant implications for its future clinical utilization and widespread adoption. Keywords: Hepatocellular carcinoma, Diwu YangGan, Microarray analysis, LC-MS, Signaling pathways Graphical abstract [43]Image 1 [44]Open in a new tab 1. Introduction Globally, Hepatocellular carcinoma (HCC) is widely acknowledged as a prominent contributor to cancer cases in both males and females, with statistics showing that it is the fifth and ninth most common cancer cause in the former and latter respectively.[45]^1 The etiology and spread of HCC remain incompletely understood, impeding the development of early detection methods and effective treatment protocols.[46]^2^,[47]^3 Although surgical resection and liver transplantation show promise as primary treatment methods for HCC, their applications have not yet reached optimal levels and are associated with various side effects for patients.[48]^4^,[49]^5 Furthermore, the application of molecular techniques in HCC diagnosis has been reported to have limitations. Current therapeutic methods for liver cancer mainly include surgical resection, liver transplantation, local ablation, external radiation, trans-arterial therapies, chemotherapy, targeted therapy, and immunotherapy. Liver transplantation, curative surgical resection, percutaneous ethanol injection and radiofrequency ablation are considered curative treatments, while transarterial treatments, are currently recognized as treatments that may prolong survival time. Although the development of targeted therapy and immunotherapy for the treatment of HCC has brought new hope to patients with advanced HCC, the overall efficacy of these therapeutic methods remains dismal. Due to its growing significance, the mRNA-based transcriptomics technology is increasingly utilized for early disease diagnosis, disease progression assessment, and drug therapy.[50]6, [51]7, [52]8 In view of the afore mentioned, the search for precise prognosis of mRNA-based molecular signatures and treatment methods for controlling HCC have become crucial. The prevention and treatment of HCC by herbal medicine has been proven an effective therapeutic strategy.[53]^9 Traditional Chinese medicine is widely embraced worldwide as a form of complementary and alternative healthcare. Chinese herbal medicine has been widely used to treat cancer. In the realm of fundamental research, herbal medicine has exhibited efficacy in diminishing the proliferation, invasion, and metastasis of hepatocellular carcinoma (HCC). Moreover, empirical investigations have substantiated the advantageous impact of herbal medicine on survival rates and overall response rates among patients afflicted with HCC.[54]^10^,[55]^11 The investigation of bioactive compounds from traditional drug or food plants has always been a topic of interest for traditional and complementary Medicine. Over the past few years, researchers have shown interest in green chemistry to extract bioactive compound using environmentally benign agents such as plants, fruits, flowers, algae, yeasts, bacteria, fungi. Diwu Yanggan capsule (DWYG) is a patent herbal formulation (Patent No: 201210580999.2) permitted for use as a drug by the Hubei Food and Drug Administration after many years of intense research (Grant No. Z20113160).[56]^12 It contains five Chinese herbal medicines, including Artemisia scoparia Waldst. & Kit., Schisandra chinensis (Turcz.) Baill., Rehmannia glutinosa Libosch., Glycyrrhiza uralensis Fisch. and Curcuma longa L. Previous studies have demonstrated that the DWYG has exhibited hepatoprotective effects and has played significant roles in enhancing immune-modulatory responses, ameliorating liver lesions, promoting liver regeneration, exerting virus-killing activities, and inhibiting HCC growth.[57]12, [58]13, [59]14, [60]15 It has been shown that liver injury caused by an array of factors, can significantly hinder the regenerative effects of mature hepatocytes on liver parenchyma.[61]^16 2-acetylaminofluorene (2-AAF) exposure results in cell lysis and disrupts other liver functions. These have made the 2-AAF/partial hepatectomy (PH) animal model, the most used in assessing the mechanisms by which medicines exert therapeutic effects.[62]17, [63]18, [64]19 In our previous study, it has been demonstrated that DWYG in L-02 cell cultures can reduce alanine transaminase (ALT), aspartate aminotransferase (AST), and Bax levels and upregulate Bcl-2 expression. Additionally, DWYG may control HCC in 2-AAF/PH rats, possibly by modulating and restoring the liver regeneration microenvironment.[65]^12 However, the comprehensive mechanisms of 2-AAF/PH-induced rat pathogenesis and the pharmacological mechanisms of DWYG is still limited. Therefore, the present study aims to explore the mechanisms by which DWYG regulates hepatic transcriptome in 2-AAF/PH-induced hepatocellular carcinoma in view of enhancing its applications in liver treatment. In addition, a component profile of DWYG based on the high-resolution mass spectrometry was constructed and it might help to understand their mechanism on 2-AAF/PH-induced hepatocellular carcinoma. 2. Materials and methods 2.1. Materials and chemicals In accordance with previous descriptions,[66]^11 the Hubei Provincial Hospital of Traditional Chinese Medicine prepared the DWYG capsules. 2-AAF was obtained from Sigma-Aldrich (Saint Louis, MO, USA). Fisher Scientific (Fair Lawn, NJ, USA) supplied HPLC-grade methanol, while a Millipore water system (Bedford, MA, USA) was used to produce deionized water. An analytical grade formic acid was purchased from Sinopharm Chemical Reagent Co., Ltd (Shanghai, China). The remaining reagents used in this study were purchased from Sinopharm Chemical Reagent Co., Ltd (Shanghai, China). 2.2. Animals and experimental design 48 male Wistar rats (180 ± 10 g) that were Specific-pathogen-free were bought from Hubei Province Experimental Animal Research Center. The rodents were kept in regulated environmental settings, maintaining a temperature of 21 ± 2 °C and a light/dark cycle lasting for 12 h. The experimental protocol was approved by the Institutional Animal Care and Use Committee of the Hubei University of Chinese Medicine (No.202107002). The modeling of 2-AAF/PH was modified based on previous reports.[67]^20 After one week of acclimatization, the rats were divided into three groups: normal, model, and DWYG groups. Each group consisted of sixteen rats. For the model and DWYG groups, the rats underwent partial hepatectomy and then orally administered a daily dose of 20 mg kg^−1 2-AAF for one week. The DWYG capsules were suspended with distilled water in a solution with a concentration of 36 mg mL^−1. The rats of DWYG group were given 10 mL kg^−1·body weight^−1·day^−1 of DWYG solution orally. The rats in normal and model group were given distilled water orally at a dose of 10 mL kg^−1·body weight^−1·day^−1. Each group followed this regimen for one week. After receiving the last dose, the rats fasted for 12 h. On the 8th day, rats were anesthetized with a 2 % sodium pentobarbital intraperitoneal injection, and blood samples and liver tissues were gathered. The collected samples were then stored at −80 °C for microarray data analysis. 2.3. Assessment of liver function The serum levels of AST, ALT and gamma-glutamyl transpeptidase (γ-GT) were assessed as indicators of liver tumor progression[68]^16 and detected according to the kit instructions. The blood specimen underwent centrifugation at a rate of 1700×g utilizing a Beckman Coulter GS-6R centrifuge situated in Fullerton, USA, at a temperature of 4 °C for a period of 30 min. Afterwards, the levels of AST, ALT, and γ-GT were measured using a Roche Cobas 8000 modular analyzer, a Cobas-e601 automatic biochemical analyzer produced by Roche Diagnostics in Switzerland. The liver tissues collected from the animals were stored in a solution of formalin (10 % concentration) for further examination of histopathology. 2.4. Microarray hybridization and microarray data analysis According to kit instructions, total RNA was isolated from hepatic tissues using RNA isolation kit (Tiangen, China). RNA labeling and array hybridization were performed according to the One-Color Microarray-Based Gene Expression Analysis protocol (Agilent Technology) with minor modifications. An mRNA isolation kit (mRNA-ONLY Eukaryotic mRNA Isolation Kit, Epicentre) was used to remove rRNA from total RNA. After amplifying and trancribing each sample, fluorescent cRNA was generated by random priming methods (Arraystar Flash RNA Labeling Kit, Arraystar). The labeled cRNAs were purified by RNeasy Mini Kit (Qiagen). By using NanoDrop's ND-1000, the specific activity and concentration of the labeled cRNAs were determined. Each labeled cRNA of 1 μg was fragmented by adding 5 μL 10 × Blocking Agent and 1 μL of 25 × Fragmentation Buffer, then heated the mixture at 60 °C for 30 min. Finally, 25 μL 2 × GE Hybridization buffer was added to dilute the labeled cRNA. Assembled to the expression microarray slide, a hybridization solution of 50 μL were dispensed into the gasket slide. A Hybridization Oven (G2545a) from Agilent was used to incubate the slides at 65 °C for 17 h. DNA Microarray Scanner (G2505C) from Agilent was used to scan the hybridized arrays after washing, fixing, and scanning. Gene expression ratios were determined. Genes with p < 0.05 and |log2 (Fold Change) | >2 were assigned as differentially expressed genes (DEGs). KOBAS software was used to identify DEGs in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway by a two-tailed Fisher's exact test. The pathway with a corrected p value < 0.05 was considered significant. DEGs were visualized by a heat map using the heatmap2 function from the gplots R-package. 2.5. Real-time quantitative reverse transcription PCR To validate the DEGs from the microarray step, qRT-PCR analysis was performed. Several overlapped DEGs between the model and normal groups, DWYG and model groups comparisons were selected for analysis by qRT-PCR. Total RNA was extracted as described in section 2.3, and 500 ng RNA was reverse transcribed into first-strand cDNA using a Revert Aid First Strand cDNA Synthesis Kit (Thermo, USA). DEG amplifications were conducted by the FastStart Universal SYBR Green Master (Roche, USA) and specific primers (produced by Qing Ke, China) on a LightCycler480 II System (Roche, USA). Primers were designed using high-throughput primer design software according to the predicted sequence of the gene. The primer sequences were shown in [69]Table S1. Ct values were determined as previously described[70]^21 and the ΔΔCt method was used to determine the relative expression levels of target genes and the endogenous housekeeping gene GAPDH was used as an internal standard. 2.6. LC-QTOF-MS/MS analysis After adding 100 mg of DWYG powder to 5 mL of 75 % methanol in water, the mixture was vigorously vortexed for 5 min and then sonicated for 10 min. After being spun at a speed of 4000 revolutions per minute for a duration of 10 min, the solution was subjected to centrifugation and the resulting liquid, known as the extract, was meticulously extracted. With the ACQUITY UPLC M-Class system (Waters, Mass. USA), a Waters ACQUITY UPLC BEH C18 column (100 × 2.1 mm, 1.7 μm) from the USA was utilized with a flow rate of 0.3 mL min^−1. A total of 2.0 μL was injected. The water/formic acid (1000 1, V/V) was used as mobile phase A, while methanol served as mobile phase B. A binary gradient with linear interpolation was employed as follows: 0.01 min, 10 % B; 15min, 55 % B; 35min, 90 % B; 40min, 98 % B; 45min, 98 % B; 46min,10 % B; 50min, 10 % B. Mass spectrometric analyses were conducted using a Waters Xevo G2-XS QTof system (Waters, Mass., U.S.A.) equipped with an electrospray ionization source. Positive and negative ion electrospray data acquisition modes were utilized to conduct data analysis. It was determined that the following parameters were most optimal: the source temperature of 100 °C, the desolvation temperature of 500 °C, the cone gas flow of 50 L h^−1, the desolvation gas flow of 600 L h^−1, the cone voltage of 60 V, and the capillary voltage of 3 kV. For the full scan, the mass ranges were set at m/z 50–1500 Da with scan duration of 1 s. The LockSpray was used to ensure accuracy and reproducibility of all analyses Data were collected in the MS^E continuum mode and were analyzed by MassLynx v4.1 Software (Waters, Mass., U.S.A.). 3. Results 3.1. Liver function analysis As shown in [71]Table 1, significant (P﹤0.01) increases in AST, ALT and γ-GT were observed in serum from the 2-AAF/PH induced model group when compared with the control group. The rats treated with the DWYG resulted in 15.1 %, 22.5 % and 35.2 % reductions in the serum levels of AST, ALT and γ-GT activity, respectively, as compared with those from model group (all P﹤0.05). Table 1. Effect of DWYG on 2-AAF/PH induced hepatocellular carcinoma in activities of AST, ALT, and γ-GT in rats (U/L, mean ± S.D., n = 16). Group AST ALT γ-GT Control 121.80 ± 3.51 52.53 ± 0.81 7.05 ± 0.61 Model 175.31 ± 12.71** 89.10 ± 5.86** 29.50 ± 6.57** DWYG 148.84 ± 9.81^## 69.00 ± 6.76^## 19.10 ± 3.82^# [72]Open in a new tab **p﹤0.01, compared with the control group. ^#p﹤0.05. ^##p﹤0.01, compared with the model group. 3.2. Differential expression analysis To investigate the underlying mechanisms of DWYG on HCC in 2-AAF/PH rats, transcriptional changes in the liver were studied using microarray protocol. DEGs were characterized using two comparisons: normal and model groups, DWYG and model groups. As shown in [73]Fig. 1, there were 4935 DEGs between the normal and model groups, including 4166 up-regulation genes and 769 down-regulation genes. Moreover, a grand total of 827 DEGs were identified between the DWYG and model groups, with 255 genes exhibiting up-regulation and 472 genes showing down-regulation. Fig. 1. [74]Fig. 1 [75]Open in a new tab The number of up-regulated and down-regulated genes in the model/normal and the DWYG/model comparisons. 3.3. KEGG pathway enrichment analysis To analyze the DEGs from the biological pathway insight, DEGs were mapped into the KEGG database record. The KEGG pathway enrichment analysis and the top 20 up-regulated and down-regulated signaling pathways shared by DEGs were shown in [76]Fig. 2. Fig. 2. [77]Fig. 2 [78]Open in a new tab Bubble diagrams of KEGG pathway enrichment for DEGs in (A) model/normal (up-regulated), (B) model/normal (down-regulated), (C) DWYG/model (up-regulated) and (D) DWYG/model (down-regulated). In these pathways, chemokine signaling pathway, endocytosis, HTLV-1 infection, metabolic pathways, pathway in cancer, PI3K-Akt signaling pathway, Rap1 signaling pathway, and regulation of actin cytoskeleton were among the up-regulated pathways in the model/normal comparison and down-regulated in the DWYG/model comparison, suggesting that although these pathways were induced by 2-AAF/PH, they could be partially reversed by DWYG treatment. In addition, hepatitis B and pathways in cancer were among the down-regulated pathways in the model/normal comparison and up-regulated in the DWYG/model comparison. The pathways in cancer were enriched both in up-regulated and down-regulated of the model/normal comparison and DWYG/model comparison, indicating that there were a large number of targets were involved in this pathway by modeling and the treatment with DWYG. Similarly, the PI3K-Akt signaling pathway was considered as a critical pathway in HCC. 3.4. DEGs related to 2-AAF/PH and DWYG To identify DEGs whose expression levels were reversed or resistant to DWYG, the Venn diagrams between the model/normal and DWYG/model comparisons were both performed by selecting the genes of 4-fold change. After DWYG treatment, a total of 54 DEGs were overlapped between the normal/model group and the DWYG/model group ([79]Fig. 3A). As shown in [80]Figs. 3B and 48 of the 54 overlapped DEGs modified by DWYG were up-regulated in the model/normal comparison, whereas these genes were alleviated by the DWYG treatment. In addition, Ubd, Cyp26a1 and RT1-S2 genes were down-regulated in the model/normal comparison, while they were up-regulated in the DWYG/model comparison. Besides, 3 genes (Dock6, Per1 and Sertad4) showed no reversions in their expressions between the model/normal and DWYG/model comparisons. These findings implied that DWYG treatment could reverse 51 variation genes induced by 2-AAF/PH. Fig. 3. [81]Fig. 3 [82]Open in a new tab DWYG and 2-AAF/PH affect alterations in liver transcription. (A) The Venn diagram shows DEGs common to both the model/normal and DWYG/model genes. (B) Heat maps of the common DEGs between the model/normal and DWYG/model comparisons. In order to validate the accuracy of the microarray protocol, genes that were selected for further analysis by qRT-PCR included 2 HCC markers,[83]^22 Trem2 and Gpc3, and a calcium-binding protein and known to promote cell death via glycolysis,[84]^23 Efhd1, as well as a regulator of the BMP pathway,[85]^24 Grem2, and a therapeutic target of Notch signaling pathway,[86]^25 Itgb6, and a kinesin-encoding gene,[87]^26 Kif12, which were the overlapped DEGs. As shown in [88]Fig. 4, it was found that the fold change of expression levels among all genes were increased by 2-AAF/PH and reversed after DWYG treatment. In addition, all levels of the fold change of expression levels were below 0.5, among these were 3 genes which were below 0.25. Moreover, the results of microarray showed that these DEGs had similar expression trends when compared to the data of qRT-PCR, suggesting that the results of DEGs was reliable in this study. Fig. 4. [89]Fig. 4 [90]Open in a new tab qRT-PCR validation results of the DNA microarray. Blue line signifies DNA microarray data for the DEGs while the red bar stands for qRT-PCR results for the DEGs. The values were the means of three replicates of three different experiments. 3.5. Component profile of DWYG by LC-QTOF-MS/MS According to the LC retention time, accurate molecular weight and MS/MS behavior, a total of 210 compounds were identified or characterized. Chemical information of 210 compound characterized from DWYG was shown in [91]Table S2. The extracted ion chromatograms of each compound were depicted in [92]Fig. 5. The peaks in [93]Fig. 5 were extracted individually according to the different peak intensities to clearly display all the peaks. Their peak numbers of compounds are consistent with those in [94]Table S2. Fig. 5. [95]Fig. 5 [96]Open in a new tab Extracted ion chromatograms of the 206 compounds characterized from DWYG by LC-QTOF-MS/MS. The chemical constituents of turmeric mainly were curcuminoids. A total of 82 curcuminoids were identified from Turmeric in this study. The fragmentation patterns of curcumin were fully investigated in previous study.[97]^27^,[98]^28 One of the characteristic diagnostic ions that can be used to search for the same type of curcuminoids is the product ion at m/z 177.0552, which is combined with other product ions produced by various phenyl-substituted groups. Fifty-six peaks detected in DWYG were identified or characterized as compounds in Licorice. Most of them were flavonoids and six of them were triterpenoids. For instance, liquiritin (compound 137), showed a molecular ion [M+H]^+ (m/z 419.1310) at 11.60 min (t[R]) in positive ion mode. It yielded three prominent fragment ions at m/z 257.0797 ([M + H-Glc]^+), 229.0881 ([M + H-Glc-CO]^+), 137.0243 ([M + H-Glc-C[8]H[8]O]^+). Glycyrrhizin (compound 107) was also detected in positive ion mode at 22.65 min (t[R]). It showed a molecular ion [M+H]^+ at m/z 823.4088 (C[42]H[62]O[16]). The production of ions at m/z 647.3810 ([M + H-Glc]^+), 471.3462 ([M + H-2Glc]^+), 453.3364 ([M + H-2Glc-H[2]O]^+). These two compounds of retention time were confirmed by comparison with reference substances. A total of 15 compounds, including flavanones, coumarins and organic acid were detected in Artemisia capillaris. Arcapillin (compound 142), showed an accurate mass at m/z 361.0941 in positive ion mode, corresponding to the molecular formula C[18]H[16]O[8]. The ion at m/z 347.0714 was produced through losing the CH[3] group. Ion at m/z 333.0970 was formed by loss of CO from the parent ion. Ion at m/z 197.0476 was generated by elimination of C[9]H[8]O[3] from the parent ion. Forty-three compounds were detected in Schisandra chinensis. For example, Schisandrin (compound 154) was detected at 17.97 min (t[R]) in positive ion mode. It showed molecular ions [M+H]^+ at m/z 433.2193 (C[24]H[32]O[7]) and [M+Na]^+ at m/z 455.2076. The MS/MS fragmentation of m/z 433.2193 exhibited a series of major product ions at m/z 415.2108 [M + H–H[2]O]^+, 385.1989 [M + H–H[2]O–OMe]^+ and 361.1617 [M + H–H[2]O–C[4]H[6]]^+. In this study, 4 compounds were detected in Rehmanniae Radix Preparata as shown in [99]Table S2. All of them existed mainly in [M+Na]^+ form. For instance, Rehmannioside D (compound 209), showed a molecular ion [M+Na]^+ (m/z 709.2172) at 0.98 min (t[R]) in positive ion mode. The main fragment ions at m/z 669.2244, 527.1592, 365.1065, and 203.0505 of Rehmannioside D were observed. Furthermore, polysaccharides, the other major constituents of Rehmanniae Radix Preparata, were not detected in the existing detecting conditions because of their poor ionization efficiency in mass spectrometry. 4. Discussions The present study uses microarrays to investigates how DWYG treatment can significantly regulate HCC biomarkers in a 2-AAF/PH rat model. The five pathways and 54 overlapped DEGs were screened to be associated with the inhibitory mechanisms of DWYG on HCC. Meanwhile, the component profile of DWYG was characterized using liquid chromatography-quadrupole time of flight mass spectrometry (LC-QTOF-MS/MS). As a result, a total of 210 compounds were identified or characterized in DWYG, 10 of which were confirmed by corresponding standards, and the others were characterized by the comparing their MS/MS behavior with that of reported literatures. A further assignment was made to the individual herbs based on the compounds; 4 compounds were from Rehmanniae Radix Preparata, 15 were from Artemisia capillaris, 43 were from Schisandra chinensis, 82 were from Turmeric, and 56 were from Licorice. These findings could provide helpful chemical information for exploring DWYG efficacy and the mechanism of action. Curcumins and their derivatives have been used widely in preclinical in vitro and in vivo models of hepatocellular carcinoma.[100]^29 Besides the strong content in triterpenoids, such as glycyrrhizin and glycyrrhetinic acid, was demonstrated to suppress NF-κB activation in TNF-α-induced hepatocytes.[101]^30 Schisandrin A and Schisandrin B could induce Oatp1b1 expression and increase its transporter activity in the human hepatocellular liver carcinoma cell line.[102]^31 Already basic formulations likely contain ingredients with considerable biological activity. In addition, based on the component profile of DWYG capsule, it is more scientific and reasonable for the study of the spectrum–effect relationship between it and the mechanism of HCC. It was investigated that DWYG capsule has therapeutic effect on CCl4-induced liver injury in vivo and in L-02 cells by the regulation of Bax/Bcl-2 signaling pathway.[103]^32 In our previous study,[104]^15 the rat model was used to determine the effect of the DWYG capsule on the occurrence and development of liver cancer and to explore its mechanisms on the development and progression of liver cancer. The hepatoprotective effects of DWYG were investigated and compared with Sorafenib. The results showed that DWYG treatment could improve the liver regeneration microenvironment by regulating the RAF/MEK/ERK pathway, thereby inhibiting the occurrence and development of liver cancer. In another previous study,[105]^15 the clinical results suggest that DWYG could improve the liver histological response rate of patients with hepatitis B e antigen (HBeAg)-negative Chronic hepatitis B virus (CHB). Specifically, it has a stable long-term curative effect and can significantly reduce the incidence of liver cirrhosis. In addition to using antiviral drugs, treating HBeAg-negative CHB with DWYG is also an important way to improve the clinical effect by influencing host factors, which provides a reference for further research and clinical applications. Chemokine signaling pathway regulated HCC pro-inflammatory which could lead to hepatobiliary cancer, characterized by swollen hepatic tissues. Both past and recent research concur that the pathways of chemokines have a crucial impact on the creation of a tumor-promoting microenvironment and the advancement of hepatobiliary cancer.[106]^33^,[107]^34 Chemokines and their receptors play a crucial role in the migration of leukocytes and cancer cells, being linked to tumors at every stage of cancer growth and serving as important factors in tumor advancement. Chemokines play a crucial role in the immune response by interacting with different receptors to regulate the movement of immune cells.[108]^35 CX3CL1, a chemokine and adhesion molecule expressed by vascular endothelial cells.[109]^36 CX3CL1 could convert inflammatory Ly6C^high macrophages into reparative Ly6C^low macrophages. Restorative Ly6C^low not only display anti-inflammatory and tissue-protective characteristics but also facilitate liver tissue injury repair and ultimately impede the advancement of liver fibrosis.[110]^37 Licorice extract was discovered to inhibit the expression of chemokine mRNA induced by TNF-alpha, in a manner that depended on the dosage.[111]^38 TREM2 is a member of the receptor family with a single IgV domain, found on the surface of myeloid cells as transmembrane glycoproteins. Upon activation of TREM2, the cell surface exhibits expression of the chemokine receptor CCR7, thereby facilitating chemokine-directed migration towards chemokine receptor ligands. One further article has reported that the curcumin could regulate the TREM2, alleviating apoptosis.[112]^39 The secretory leukocyte protease inhibitor (SLPI) is a protein that belongs to the whey acidic protein four-disulfide core family. It is a non-glycosylated, single-chain protease inhibitor. SLPI's primary role is to hinder antimicrobial activity, regulate inflammation, suppress proteases, and regulate immune response.[113]^40 According to the reported article, SLPI hindered the activation of NF-κB, thereby suppressing the production of chemokines in response to microbial signals.[114]^41 During our investigation, we observed an up-regulation of the Cx3cl1, Trem2, and Slpi genes in the comparison between the model and normal conditions. This result indicated that DWYG could inhibit the progression of hepatobiliary cancer by Chemokine signaling pathway. The results of the current investigation validate that the PI3K/Akt pathway was up-regulated in the model group. There are several cellular processes that are influenced by the PI3K/Akt pathway, including growth and survival of tumors. The PI3K/Akt pathway has an impact on various cellular processes, such as the development and viability of tumors.[115]^42 Also, PI3K/Akt signaling contributes to tumorigenesis by working in conjunction with mTOR.[116]^43 Some literatures have reported that licorice, turmeric, Artemisia capillaris, Schisandra chinensis and Rehmanniae Radix Preparata could all inhibition of the PI3K-Akt signaling pathway.[117]44, [118]45, [119]46, [120]47, [121]48 The lignans of Schisandra chinensis has been reported to decrease the PI3K, Akt and Akt downstream signaling molecules Mcl-1 expression.[122]^49 Moreover, a recent study demonstrated that TREM2 functions as a cancer inhibitor in HCC by directing its attention towards the PI3K/Akt/β-catenin pathway.[123]^22 Through its ability to promote cell proliferation, angiogenesis, and inhibit apoptosis, CAPN6 promotes tumor development. Through the PI3K/AKT signaling pathway, CAPN6 is regulated in liver cancer, where it promotes proliferation and inhibits apoptosis.[124]^50 In order to effectively participate in cell migration and adhesion, SRPX2 is composed of cell secretory proteins and polysaccharides.[125]^43 In a prior investigation, it was demonstrated that suppression of SRPX2 effectively blocked the activation of the PI3K/Akt/mTOR pathway in prostate cancer cells.[126]^51 DPEP1, a metalloproteinase that relies on zinc, is commonly disrupted in numerous cancer types. By activating the PI3K/Akt/mTOR signaling pathways, it functioned as an oncogene in hepatoblastoma.[127]^52 However, Capn6, Srpx2 and Dpep1 were all down-regulated by DWYG in our study. The findings of the present study and indeed the theories put forward in available literature provide an update on the discovery and development of clinical substances that can disrupt HCC pro-inflammatory factors and pathways such as PI3K and Akt genes.[128]^53 Proteoglycans play a crucial role in the cellular and pericellular microenvironments, serving as important molecular agents in cancer biology.[129]^54 In our study, Gpc3, Cdcp1, Loxl1 and Ptprz1 genes were up-regulated in the model/normal comparison, while they were down-regulated in the DWYG/model comparison. Among of these genes, Gpc3 is a membrane heparan sulfate proteoglycan, is specifically expressed in HCC.[130]^55 The study conducted by Pan, Z. and colleagues has shown that increased expression of Gpc3, hinders the onset and progression of HCC.[131]^56 An expressed cell surface glycosylated transmembrane protein, Cdcp1 could be highly overexpressed in cancers such as HCC. Its overexpression was associated with more aggressive forms of these cancers.[132]^57 Inhibition of Loxl1 may delay liver fibrosis progression since Loxl1is involved in the activation of HSCs during fibrogenesis.[133]^58 Ptprz1 produces either chondroitin sulfate proteoglycan or a variant lacking proteoglycan in various types of cancer tissues. The various roles it plays in the development, advancement, and spread of cancer also have an impact on the treatment and outlook of the disease.[134]^59 Pathways in cancer is a series of signal pathways related to cancer, including Wnt and TGF-β etc. Dysregulation of the Wnt signaling pathway is linked to various types of cancer, such as breast cancer, prostate cancer, liver cancer, and more.[135]^60 PROM1, serving as a cancer stem cell indicator, is a subject of interest in Wnt signaling regenerative pathways in various types of cancer cells. The experiments demonstrated that GPC3 could promote the growth of HCC by stimulating canonical Wnt signaling.[136]^61 From the material basis of DWYG, we found that DWYG has great advantages in decrease the Wnt signaling pathway. For example, it has been reported that curcumin could inhibits HCC tumor growth by downregulated of Wnt/β-catenin signaling pathway.[137]^62^,[138]^63 Similarly, licorice flavonoids and catalpol could also decrease the Wnt signaling pathway. TGF-β is a crucial regulatory cytokine throughout the liver fibrosis progression.[139]^64 Throughout the progression of liver fibrosis, TGF-β plays a vital role as a regulatory cytokine. Currently, TGF-β is recognized as the most potent cytokine in stimulating hepatic fibrosis. It triggers HSC activation via the TGF-β-Smad signaling pathway, leading to enhanced synthesis of ECM in HSCs and modification of ECM metabolism, ultimately facilitating the progression of liver fibrosis. Research has indicated that blocking TGF-β signaling leads to numerous collaborative downstream consequences that are expected to enhance the clinical outcome in HCC.[140]^65 There are studies have demonstrated that the overexpression of Mfap4 could mitigate the TGF-β-induced enhanced expression of fibrosis-related proteins,[141]^66 and GPC3 could also affects HCC cell growth by TGF-β signaling pathway.[142]^56 Moreover, from the perspective of material fundamental, curcumin could significantly suppressed mRNA expression levels of TGF-β,[143]^67 and glycyrrhizin,[144]^68 catalpol[145]^69 and schisandrin B[146]^70 also have the same function. The movement of cells plays a crucial role in the invasion and spread of tumors, and controlling this mechanism will pave the way for effective cancer treatments. The actin cytoskeleton plays a crucial role in the movement of tumor cells and their ability to spread to other parts of the body.[147]^71 Recently, it has been reported that adenylate cyclase-associated protein 1, a protein responsible for regulating the actin cytoskeleton, is involved in cell motility and the development of pancreatic cancer.[148]^72 It has been reported that curcumin can inhibit cell migration and invasion, as well as downregulate the expression of the Cdc42 gene and Cdc42-related target genes. Additionally, it caused reorganizations of the actin cytoskeleton.[149]^73 Although the therapeutic effect of DWYG on HCC has been confirmed based on our previous and current research, its efficacy and safety of DWYG compared to other treatments for HCC remains to be clear. Additionally, the study did not investigate the potential side effects of DWYG, which is considered in our further research. In addition, based on the component profile of DWYG capsule, it is more scientific and reasonable for the study of the spectrum–effect relationship between it and the mechanism of HCC. The manuscript suggests that DWYG may be a promising area for further research and clinical development, and that it may have potential as a complementary or alternative treatment for HCC. 5. Conclusions The current investigation examined the transcriptomic consequences of DWYG treatment on the inhibition of HCC initiation and progression. Findings demonstrated that DWYG exhibited the ability to restore impaired organ and systemic functions through the modulation of various pathways, including chemokine signaling, endocytosis, HTLV-1 infection, metabolic pathways, pathway in cancer, PI3K-Akt signaling, Rap1 signaling, hepatitis B, and regulation of actin cytoskeleton. Furthermore, DWYG treatment effectively reversed the expression of 51 genes that were altered by 2-AAF/PH induction. Moreover, the component profile of DWYG provides evidence suggesting that these potential mechanisms may be associated with the active components present in the herbal formulation. Subsequent investigations will prioritize exploring the synergistic effects among the constituents in DWYG and their regulatory targets on HCC, potentially leading to advancements in the field of HCC therapy. Author contributions All authors contributed to the conception and design of the study. Writing – Original Draft Preparation and Writing – Review & Editing: Qingxin Shi and Jiangcheng He; Methodology and Software Guangya Chen and Jinlin Xu; Formal Analysis, Investigation, and Visualization: Zhaoxiang Zeng, Xueyan Zhao, Binbin Zhao, Xiang Gao and ZhihuaYe; Supervision, Project Administration, and Funding Acquisition: Mingzhong Xiao and Hanmin Li. All authors gave final approval and agreed to be accountable for the aspects of the work. Funding This research was funded by the National Natural Science Foundation of China with grant number [No. 81573815, No. 81973669, No. 81703912 and No. 81603484]. Data availability statement The original contributions presented in the study are included in the article/Supplementary Material; further inquiries could be directed to the corresponding authors. Declaration of generative AI in scientific writing The AI tools were not to the used to analyze and draw insights from data as part of the research process. Declaration of competing interest None. Footnotes Peer review under responsibility of The Center for Food and Biomolecules, National Taiwan University. ^Appendix A Supplementary data to this article can be found online at [150]https://doi.org/10.1016/j.jtcme.2023.12.002. Contributor Information Mingzhong Xiao, Email: xiaomingzhong@hbhtcm.com. Hanmin Li, Email: lihanmin69@126.com. Appendix A. Supplementary data The following are the Supplementary data to this article: Multimedia component 1 [151]mmc1.docx^ (30.2KB, docx) Multimedia component 2 [152]mmc2.docx^ (90.9KB, docx) References