key: cord-0050894-b5hda50d authors: Wu, Zhulin; He, Li; Wang, Lianan; Peng, Lisheng title: Systematically Exploring the Antitumor Mechanisms of Core Chinese Herbs on Hepatocellular Carcinoma: A Computational Study date: 2020-09-15 journal: Evid Based Complement Alternat Med DOI: 10.1155/2020/2396569 sha: a6e90a17e27806e91803cef8aec175eb02fc87ef doc_id: 50894 cord_uid: b5hda50d OBJECTIVE: Chinese herbs play a positive role in the management of hepatocellular carcinoma (HCC) in China. However, it is not clear which of Chinese herbs are critical for the treatment of HCC. Besides, mechanisms of CCHs in the treatment of HCC remain unclear. Hence, our goal is to identify the core Chinese herbs (CCHs) for treating HCC and explore their antitumor mechanism. METHODS: Firstly, clinical traditional Chinese medicine (TCM) prescriptions for HCC were collected from Chinese National Knowledge Infrastructure (CNKI) database, and then, data mining software was used to identify CCHs. After that, bioactive compounds and corresponding target genes of CCHs were obtained using three TCM databases, and target genes of HCC were acquired from MalaCards and OMIM. Subsequently, common target genes of CCHs and HCC were screened. Moreover, biological functions and pathways were analyzed, and Cytoscape plugin cytoHubba was used to identify hub genes. Finally, prognostic values of hub genes were verified by survival analysis, and the molecular docking approach was utilized to validate the interactions between targets and bioactive compounds of CCHs. RESULTS: Eight CCHs were determined from 630 prescriptions, and 100 bioactive compounds (e.g., quercetin and luteolin) and 126 common target genes were screened. Furthermore, common target genes of CCHs and HCC were mainly enriched in cancer-associated pathways, and six hub genes with statistical significance in survival analysis were selected as key target genes for molecular docking. Additionally, molecular docking showed that the bioactive compounds docked well with the protein receptors of key target genes. CONCLUSION: By combining data mining, network pharmacology, molecular docking, and survival analysis methods, we found that CCHs may play a therapeutic role in HCC through regulating the target genes and pathways related to cancer occurrence and development, angiogenesis, metastasis, and prognosis. As one of the most common cancers, hepatocellular carcinoma (HCC) is characterized by high morbidity and mortality rates. HCC is predicted to be the sixth most frequent cancer and the fourth most common cause of cancer death worldwide in 2018 [1] . In addition, advancedstage HCC has a 5-year survival rate of only 5-15% [2] . e chief etiologies for HCC include chronic infection with hepatitis B virus (HBV) or hepatitis C virus (HCV), alcoholic liver disease, and nonalcoholic fatty liver disease [3] . Currently, the therapeutic options for patients with HCC are poor and mainly include surgery, liver transplantation, chemoembolization, and molecularly targeted therapy [4, 5] . Moreover, most patients with late-stage HCC have lost the opportunity of surgical resection, and not all patients with advanced HCC are suitable for chemotherapy or targeted therapy. us, it is necessary to develop novel approaches for HCC control. More and more researchers are paying significant interest in traditional Chinese medicine (TCM) which has been used to treat cancer in China for a long time [6] . Furthermore, TCM is effective in improving symptoms, reducing side effects of chemotherapy, suppressing cancer cell growth, and regulating key intracellular signaling pathways [7] . As an important part of TCM, Chinese herbs play a positive role in the management of HCC in China [8] . Previous literature indicated that several TCM prescriptions consisting of Chinese herbs had anti-HCC effects in both basic and clinical research [9, 10] . Moreover, evidence demonstrated that TCM combined with chemotherapy showed significant efficacy and safety in improvement of life quality and reduction of chemotherapy side effects [11] . Nevertheless, TCM prescriptions for treating HCC are often based on the experience of TCM doctors, and it is unclear which core Chinese herbs (CCHs) are effective treatments for HCC. Additionally, the molecular mechanisms underlying the anti-HCC activity of CCHs are not completely clear. Currently, data mining analysis of TCM prescriptions has become a research focus in the TCM field [12] , which can be used to identify CCHs from a large number of clinical prescriptions. In the past, the complicated interactions between "multi-components" and "multi-targets" of Chinese herbs hindered the mechanism study of these herbs. At present, network pharmacology provides an effective solution to overcome these obstacles and can reveal the synergistic effects of complex Chinese herbs on human systems from a holistic view [13] . In the present study, data mining software was applied to determine CCHs for the treatment of HCC, and mechanisms of CCHs on HCC were analyzed by network pharmacology. Besides, survival analysis and molecular docking methods were used to validate the results of network pharmacology analysis. Clinical TCM prescriptions intended solely for HCC were collected from studies in the Chinese National Knowledge Infrastructure (CNKI) database. e search was conducted using the following search terms: "(Chinese medicine OR prescription OR decoction) AND (hepatocellular carcinoma OR liver cancer) AND clinical" (date: 1979 to 3 March 2020). e inclusion criteria for the studies included the following: (1) the first diagnosis of cancer patients being HCC; (2) clinical research on oral TCM prescriptions or oral TCM prescriptions combined with Western medicine in the treatment of HCC; (3) experience of TCM experts. Besides, exclusion criteria were as follows: (1) repeated literature and animal experiments; (2) TCM prescriptions in the literature being primarily for treating acute symptoms (e.g., cold and cough); (3) external prescriptions or Chinese patent medicine. Subsequently, eligible studies were screened, and TCM prescriptions were extracted from the included studies. Furthermore, Traditional Chinese Medicine Inheritance Support System (TCMISS, from the Institute of Chinese Materia Medica, China Academy of Chinese Medical Sciences, version 2.5) was employed to identify the CCHs from all prescriptions. TCMISS software has been widely used to analyze TCM prescriptions, and it has the function of text mining, association rules, and complex system entropy clustering methods [12] . In the present study, association rules in TCMISS were utilized to determine CCHs under the condition of support degree ≥126 (20%). Bioactive compounds and corresponding targets of CCHs were obtained using Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP, http://tcmspw.com/tcmsp.php) [14] , Integrative Pharmacology-Based Research Platform of TCM (TCMIP, http:// www.tcmip.cn/) [15] , and Bioinformatics Analysis Tool for Molecular Mechanism of TCM (BATMAN-TCM, http:// bionet.ncpsb.org/batman-tcm/) [16] . In this study, all TCM prescriptions composed of Chinese herbs were given orally, and the bioactive compounds were selected under the conditions of drug-likeness (DL) ≥0.18 (mean value for all molecules within the DrugBank database) and oral bioavailability (OB) ≥30% [17] . Additionally, all compounds are numbered using Mol ID from TCMSP, and all target names were converted to official gene names using the UniProt database (https://www.uniprot.org/). Moreover, potential therapeutic target genes associated with HCC were obtained from the MalaCards (https://www.malacards.org/) [18] and OMIM databases (https://omim.org/) [19] using "hepatocellular carcinoma" as the keyword, and known target genes for HCC were screened after removing duplicates. Venn diagrams were used to determine common target genes of CCHs and HCC by "VennDiagram" package in R 3.6.0 (https://www.r-project.org/). en, the interaction network of CCHs bioactive compounds and common target genes was built using Cytoscape 3.7.1 (https://cytoscape.org/), and the relationships between bioactive compounds of CCHs and common target genes were displayed in the interaction network. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomics (KEGG) pathway analyses for common target genes were done using the Database for Annotation, Visualization, and Integrated Discovery (DAVID, https:// david.ncifcrf.gov/) [20] , and false discovery rate (FDR) <0.05 was accepted as significant. In addition, GO analysis was performed according to three categories, namely, biological process (BP), cellular component (CC), and molecular function (MF), and the results were shown as a bar plot using the "ggplot2" package in R 3.6.0. Hub Genes. Protein-protein interaction (PPI) network analysis of the common target genes was completed using Search Tool for the Retrieval of Interacting Genes (STRING) database (https://string-db.org/) [21] with the highest confidence (score >0.9). en, hub genes were screened from the PPI network by degree, maximum neighborhood component (MNC), and maximal clique centrality (MCC) algorithms in the Cytoscape 3.7.1 plugin, cytoHubba [22, 23] . e overlap genes, predicted by all three algorithms, were selected as hub genes. e prognostic values of hub genes were assessed by survival analysis using the Kaplan-Meier plotter (http://kmplot.com/ analysis/) [24] . e Kaplan-Meier plotter includes data on HCC survival in the Cancer Genome Atlas (TCGA). In survival analysis, overall survival (OS) was analyzed by the Kaplan-Meier (KM) method (log-rank test), and a p value less than 0.05 was considered to be a significant difference. In this study, the hub genes that displayed statistically significant differences were considered to be the key target genes. Compounds. In the present study, molecular docking was used to verify the interactions between protein receptors of key target genes and bioactive compounds of CCHs. e corresponding protein receptors of key target genes were acquired using Protein Data Bank (PDB) database (https:// www.rcsb.org/) [25] , and we included protein receptors based on the following criteria: (1) e structure of protein receptors was determined by X-ray diffraction approach. (2) X-ray resolution <3Å was preferred. (3) Protein structure with initial ligand was also preferred. AutoDockTools (version 1.5.6, http://autodock.scripps.edu) was utilized to remove the original ligands (if any), excess protein chains, and water molecules of protein receptors [26] , and then hydrogens were added to the protein receptors. Furthermore, grid boxes in AutoDockTools were used to identify the docking coordinates. After that, bioactive compounds of CCHs corresponding to key targets were used as ligands, and the "mol2" files of compounds were obtained using TCMSP. Moreover, the files of protein receptors and ligands were converted into "PDBQT" format using AutoDockTools. Finally, AutoDock Vina program (http://vina.scripps.edu/) [27] was used to dock bioactive compounds into the corresponding protein receptors. Besides, docking results were analyzed and visualized using PyMOL (http://www.pymol. org/) and Discovery Studio 2016 (BIOVIA, San Diego, USA). e workflow of our study is displayed in Figure 1 . A total of 1472 studies were identified through CNKI database search. According to the criteria described above, 630 TCM prescriptions for the treatment of HCC were obtained from 527 studies. In addition, the result of TCMISS analysis showed that there were 180 different Chinese herbs used in all prescriptions, and a total of 19 combinations of Chinese herbs that were most commonly used in the TCM prescriptions were found (Table 1) . en, nine CCHs were identified by TCMISS, namely, Largehead Atractylodes Rhizome (Bai Zhu) , Poria (Fu Ling), Radix Bupleuri (Chai Hu), Radix Astragali (Huang Qi), Herba Hedyotis (Bai Hua She She Cao), Radix Codonopsis (Dang Shen), Radix Paeoniae Alba (Bai Shao), Radix Glycyrrhizae (Gan Cao), Herba Scutellariae Barbatae (Ban Zhi Lian), which were shown in a network (Figure 2(a) ). According to the theory of TCM, Radix Glycyrrhizae (Gan Cao) is often used as a harmonizing (Tiao He) drug, so it was not included in the following analysis. Finally, excluding Radix Glycyrrhizae (Gan Cao), eight CCHs were selected for further analysis. Bioactive compound counts of each core Chinese herb from the three databases are shown in Table 2 . After merging the data and removing duplicates, a total of 100 bioactive compounds were identified, and 840 target genes of eight CCHs were identified. Additionally, 191 and 745 HCC-related target genes were obtained from OMIM and MalaCards databases, respectively. Following the removal of duplicates, a total of 918 therapeutic target genes for HCC were found. A total of 126 common target genes of CCHs and HCC were identified using Venn diagram ( Figure 2 (b)), and the interaction network of CCHs bioactive compounds and common target genes was established, including 198 nodes (72 bioactive compounds and 126 genes) and 510 edges ( Figure 3 ). As shown in Figure 3 , eight CCHs were divided into 3 types: health-strengthening (Fu Zheng), heat-clearing and detoxicating (Qingre Jiedu), and relieving liver Qi stagnation (Shu Gan). Furthermore, the top 30 bioactive compounds are presented in Table 3 according to the gene count, and the results showed that quercetin, stigmasta-5,22dien-3-one, luteolin, wogonin, kaempferol, beta-sitosterol, baicalein, stigmasterol, pyrethrin II, etc. connected with most of the common target genes. e result of GO analysis revealed that common target genes were mainly enriched in negative regulation of apoptotic process, positive regulation of gene expression, response to drug, positive regulation of transcription from RNA polymerase II promoter, and cellular response to hypoxia (BP); enzyme binding, protein binding, transcription factor binding, identical protein binding, and protein heterodimerization activity (MF); cytosol, nucleus, cytoplasm, nucleoplasm, and phosphatidylinositol 3-kinase complex (CC). Top 5 GO terms of each category are shown in Figure 4 . Besides, KEGG pathway analysis demonstrated that most common target genes were related to pathways in cancer, proteoglycans in cancer, hepatitis B, PI3K-Akt, HIF-1 signaling pathway, etc., and the results of the top 20 signaling pathways are listed in Table 4 . e PPI network with 126 nodes and 698 edges was constructed by STRING and visualized by Cytoscape v. 3.7.1 ( Figure 5(a) ). en, hub genes were identified based on MNC, degree, and MCC methods ( Figure 5(b) ). e interaction network of eight hub genes is shown in Figure 5 Table 5 . In the present study, the binding between the receptor (protein receptor of key target gene) and the ligand (bioactive compound) was considered to be good if the affinity value <−5.0 kcal/mol, and the lower the affinity value was, the higher the affinity of the bioactive compound and the protein receptor was. As shown in Table 5 , stigmasta-5,22-dien-3-one, luteolin, quercetin, pyrethrin II, palbinone, and baicalein had strong binding affinities for the corresponding proteins, and molecular docking results are shown in Figures 7(a) -7(f ). According to the data mining results, eight Chinese herbs, including Largehead Atractylodes Rhizome, Poria, Radix Bupleuri, Radix Astragali, Herba Hedyotis, Radix Codonopsis, Radix Paeoniae Alba, and Herba Scutellariae Barbatae, were identified as the CCHs for the treatment of HCC. Based on the TCM theory, it is believed that the pathogenesis of HCC is due to "deficiency of healthy Qi," "liver Qi stagnation," "heat-toxicity," etc. In TCM, Largehead Atractylodes Rhizome, Poria, Radix Astragal, Radix Codonopsis, and Radix Paeoniae Alba can be used for "supporting the healthy energy"; Herba Hedyotis and Herba Scutellariae Barbatae have the effect of "clearing heat and removing toxicity"; Radix Bupleuri can be used for "relieving liver Qi stagnation." us, data mining results were highly consistent with the TCM theory. In addition, the bioactive compounds-common targets network showed that the critical bioactive compounds of CCHs could be quercetin, stigmasta-5,22-dien-3-one, luteolin, wogonin, kaempferol, betasitosterol, baicalein, Stigmasterol, Pyrethrin II, formononetin, (+)-catechin, and palbinone according to the gene counts (genes ≥10). Previous studies indicated that quercetin and luteolin could inhibit the growth of HCC cells and enhance chemotherapy efficacy [28, 29] . Also, wogonin and kaempferol could effectively suppress the proliferation and invasion of HCC cells by regulating EGFR signaling pathways and PI3K/AKT/mTOR pathway, respectively [30, 31] . Stigmasterol, beta-sitosterol, and baicalein have been found to induce apoptosis of HCC cells by upregulating proapoptotic gene (Bax) and downregulating antiapoptotic gene (Bcl-2) [32] [33] [34] . Experimental research showed that formononetin could impede the epithelial-mesenchymal transition (EMT) and malignant progression of HCC [35] . Besides, (+)-catechin could inhibit the proliferation of HCC cells via the caspase-dependent pathway [36] . Overall, previous studies strongly support our findings. e results of GO analysis showed that common target genes were significantly involved in cellular components and biological processes which are related to cell apoptosis, proliferation, differentiation, and various cellular functions. Moreover, the molecular functions of common target genes may be correlated with the physiological and metabolic processes of liver. KEGG analysis revealed that most of the common target genes were mainly enriched in cancerrelated signaling pathways. Previous literature reported that the PI3K/Akt pathway was activated in 30-50% of HCC and the upregulation of p-Akt was correlated with poor survival, metastasis, and vascular invasion in HCC patients [37, 38] . erefore, PI3K/Akt pathway could shed light on a novel strategy for drug development for HCC. Furthermore, a positive correlation between HBV infection and HCC was observed, and HBV infection can result in the activation of protooncogenes and inactivation of tumor suppressor genes [39, 40] . Proteoglycans are extracellular matrix components of liver microenvironment which play an important role in the progression of HCC and have the potential to be the HCC therapeutic target [41] . It is also reported that kindlin-2 is a member of the focal adhesion protein family which promotes HCC invasion, metastasis [42] . MicroRNA dysregulation has been found to be involved in all stages of HCC, and some microRNAs, such as miR-17-92, miR-21, and miR-221, are generally upregulated in HCC [43] . In addition, HIF-1 plays a critical role in immune escape and EMT of HCC [44] . Literature has shown that the interaction between thyroid hormone and its receptor plays an important role in the regulation of development and proliferation, and metastasis of HCC Evidence-Based Complementary and Alternative Medicine [45] . Research also showed that the occurrence of HBVrelated HCC activates the Ras/MAPK signaling pathway which is correlated with a poor prognosis [46, 47] . Besides, patients with cancer are more likely to suffer from various infections due to low immune function, which may be related to HTLV-I infection, Chagas disease, and Influenza A pathways. Taken together, our results are in line with previous studies, suggesting that CCHs may exert their antitumor effect in HCC by regulating the occurrence, progression, angiogenesis, and metastasis of HCC. Based on the results of the PPI network and survival curves, we identify six key target genes of CCHs in the treatment of HCC, namely, SRC, RHOA, PIK3R1, EGFR, VEGFA, and EGF. An experiment showed that HBV core protein could promote tumorigenesis of HCC cells by upregulating the expression of SRC protooncogene and then activating SRC/PI3K/Akt pathway [48] . A previous study demonstrated that RHOA (Ras homolog gene family member A) is commonly overexpressed in HCC, and its expression is associated with poor prognosis [49] . Additionally, evidence showed that knockdown of PIK3R1 promoted apoptosis of HCC cells and downregulated p-PI3K and p-AKT expressions in HCC cells [50] . As a growth factor, EGF plays a crucial role in cell proliferation and migration by binding to its receptor EGFR, and high expression of EGF could induce highly malignant HCC [51] . Previous studies suggested that EGFR is overexpressed or mutated in HCC and may be closely related to the formation, invasive growth, and clinical characteristics of HCC [52, 53] . Angiogenesis is closely related to tumor growth and invasion, and HCC is recognized as a typical angiogenic tumor [54] . VEGFA is an inducer of angiogenesis in HCC, and the expression of VEGFA in HCC was significantly higher than that in normal liver tissues [55] . Overall, previous studies made our results more reliable to some extent. e results of network pharmacology were also confirmed by molecular docking. Table 5 shows that the affinity values of all docking results were less than −5.0 kcal/mol, which indicated that the protein receptors of nine key target genes were docked well with the six different compounds of CCHs. As shown in the three-dimensional mode of Figure 7 , six active compounds were successfully docked to the active pocket of protein receptors (SRC, PIK3R1, EGFR, RHOA, VEGFA, EGF), and the two-dimensional diagram in Figure 7 also demonstrates the interactions between active compounds and protein receptors. According to our findings, van der Waals forces, hydrogen bonds, Alkyl, π-Alkyl, π-Cation, π-Sigma, etc. were shown to be involved in the interactions between receptors and compounds. For instance, Figure 7 (a) indicates that two hydrogen bonds between stigmasta-5,22-dien-3-one and amino acid residues (ARG-358 and SER-361) of PIK3R1 protein were generated. In addition, the van der Waals force, alkyl, and π-alkyl interactions between stigmasta-5,22-dien-3-one and other residues also play an important role. e interpretation of other molecular docking results could use this similar method. ese results have successfully validated the network pharmacology data from the perspective of molecular interactions. However, this paper lacks biological experimental confirmation. Both in vivo and in vitro experiments are required to validate our results, and we should take it into consideration in future work. Last but not least, there are few reports about some bioactive compounds (e.g., stigmasta-5,22-dien-3-one, pyrethrin II, palbinone), and further research of these compounds is still needed. To sum up, based on a combination of TCM prescription data mining, network pharmacology, KM survival analysis, and molecular docking, we identified eight CCHs for treating HCC and found that CCHs may play a therapeutic role in HCC by regulating the genes and pathways associated with cancer occurrence and development, angiogenesis, metastasis, and prognosis. Data Availability e data that support the findings of our study are openly available in http://tcmspw.com/tcmsp.php, http://www.tcmip. cn/, http://bionet.ncpsb.org/batman-tcm/, https://www. uniprot.org/, https://www.malacards.org/, https://omim.org/, https://david.ncifcrf.gov/, https://string-db.org/, http://kmplot. com/analysis/, and https://www.rcsb.org/. e data used and analyzed during the current study are available from the corresponding author on reasonable request. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries Feasibility of terahertz imaging for discrimination of human hepatocellular carcinoma Hepatocellular carcinoma Detection of doublecortin domain-containing 2 (DCDC2), a new candidate tumor suppressor gene of hepatocellular carcinoma, by triple combination array analysis Matrix metalloproteinase 2 knockdown suppresses the proliferation of HepG2 and Huh7 cells and enhances the cisplatin effect Comparison of efficacy and toxicity of traditional Chinese medicine (TCM) herbal mixture LQ and conventional chemotherapy on lung cancer metastasis and survival in mouse models Role of traditional Chinese medicine in the management of patients with hepatocellular carcinoma Preventive and therapeutic effects of Chinese herbal compounds against hepatocellular carcinoma Traditional herbal medicine: a review of potential of inhibitory hepatocellular carcinoma in basic research and clinical trial e Chinese medicine, jiedu recipe, inhibits the epithelial mesenchymal transition of hepatocellular carcinoma via the regulation of Smad2/3 dependent and independent pathways A systematic review and metaanalysis of traditional insect Chinese medicines combined chemotherapy for non-surgical hepatocellular carcinoma therapy Analysis on composition rules of Chinese patent drugs treating pain-related diseases based on data mining method Unraveling the action mechanism of polygonum cuspidatum by a network pharmacology approach TCMSP: a database of systems pharmacology for drug discovery from herbal medicines ETCM: an encyclopaedia of traditional Chinese medicine BATMAN-TCM: a Bioinformatics analysis Tool for molecular mechANism of traditional Chinese medicine Uncovering the mechanism of Maxing Ganshi Decoction on asthma from a systematic perspective: a network pharmacology study MalaCards: an amalgamated human disease compendium with diverse clinical and genetic annotation and structured search OMIM.org: leveraging knowledge across phenotype-gene relationships DAVID: database for annotation, visualization, and integrated discovery STRING v10: protein-protein interaction networks, integrated over the tree of life cytoHubba: identifying hub objects and sub-networks from complex interactome Study on the multitarget mechanism and key active ingredients of Herba siegesbeckiae and volatile oil against rheumatoid arthritis based on network pharmacology Online survival analysis software to assess the prognostic value of biomarkers using transcriptomic data in non-small-cell lung cancer e RCSB protein data bank: integrative view of protein, gene and 3D structural information Mechanism of Dayuanyin in the treatment of coronavirus disease 2019 based on network pharmacology and molecular docking AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading Quercetin induces apoptosis and enhances 5-FU therapeutic efficacy in hepatocellular carcinoma Luteolin synergizes the antitumor effects of 5-fluorouracil against human hepatocellular carcinoma cells through apoptosis induction and metabolism Wogonin inhibits the proliferation and invasion, and induces the apoptosis of HepG2 and Bel7402 HCC cells through NF-κB/Bcl-2, EGFR and EGFR downstream ERK/AKT signaling Kaempferol inhibits proliferation, migration, and invasion of liver cancer HepG2 cells by downregulation of microRNA-21 Stigmasterol isolated from marine microalgae Navicula incerta induces apoptosis in human hepatoma HepG2 cells Baicalein: a review of its anticancer effects and mechanisms in Hepatocellular Carcinoma Antiproliferative effects mechanism of beta-sitosterul in hepatoma HepG2 cells USP5 promotes epithelialmesenchymal transition by stabilizing SLUG in hepatocellular carcinoma A study on the role of (+)-catechin in suppression of HepG2 proliferation via caspase dependent pathway and enhancement of its in vitro and in vivo cytotoxic potential through liposomal formulation Pathways and targets in hepatocellular carcinoma Hepatocellular carcinoma: systemic therapies and future perspectives Hepatitis B e antigen and the risk of hepatocellular carcinoma Initial sites of hepadnavirus integration into host genome in human hepatocytes and in the woodchuck model of hepatitis B-associated hepatocellular carcinoma Proteoglycans are attractive biomarkers and therapeutic targets in hepatocellular carcinoma Kindlin-2 promotes hepatocellular carcinoma invasion and metastasis by increasing Wnt/ β-catenin signaling MicroRNAs as biomarkers for liver disease and hepatocellular carcinoma Development and validation of a prognostic classifier based on HIF-1 signaling for hepatocellular carcinoma yroid hormone in the regulation of hepatocellular carcinoma and its microenvironment Effects of hepatitis B virus X protein on the development of liver cancer e Ras/MAPK pathway and hepatocarcinoma: pathogenesis and therapeutic implications Hepatitis B virus core protein promotes hepatocarcinogenesis by enhancing Src expression and activating the Src/PI3K/Akt pathway e diagnostic and prognostic role of RhoA in hepatocellular carcinoma Overexpression of PIK3R1 promotes hepatocellular carcinoma progression EGF is highly expressed in hepatocellular carcinoma (HCC) and promotes motility of HCC cells via fibronectin Missense mutations in exons 18-24 of EGFR in hepatocellular carcinoma tissues EGFR signaling in liver diseases Long non-coding RNA UBE2CP3 enhances HCC cell secretion of VEGFA and promotes angiogenesis by activating ERK1/2/HIF-1α/VEGFA signalling in hepatocellular carcinoma MiR-497 suppresses angiogenesis and metastasis of hepatocellular carcinoma by inhibiting VEGFA and AEG-1