- Research article
- Open Access
- Open Peer Review
ZBP-89 and Sp1 contribute to Bak expression in hepatocellular carcinoma cells
BMC Cancer volume 18, Article number: 419 (2018)
Kruppel family member zinc binding protein 89 (ZBP-89), also known as ZNF148, regulates Bak expression via binding to GC-rich promoter domain. It is not clear if other GC-rich binding factors, such as Sp family members, can interact with ZBPp-89 on Bak expression. This study aims to elucidate the mechanism of Bak expression regulation by ZBP-89 and Sp proteins, based on in vitro experiment and The Cancer Genome Atlas (TCGA) hepatocellular carcinoma (HCC) data cohort.
We downloaded TCGA hepatocellular carcinoma (HCC) cohort data to analysis the association of Bak transcription level with ZBP-89 and Sp proteins transcription level. HCC cell lines and liver immortal non-tumour cell lines were used for mechanism study, including western blotting analysis, expression vector mediated gene expression and siRNA interference.
Results showed that cancer tissues have higher Bak transcription level compared with adjacent non-cancer tissues. Bak transcription level was correlated with Sp1 and Sp3 expression level, while no correlation was found in ZBP-89 and Bak, neither Sp2 nor Sp4. Mithramycin A (MMA) induced Bak expression in a dose-dependent manner. Western blotting results showed Sp1 overexpression increased Bak expression both in liver immortal non-tumour cells and HCC cells. Interference Sp1 expression could inhibit Bak expression alone. ZBP-89 siRNA suppressed Bak expression even in the presence of MMA treatment and S1 overexpression. Additionally, Bak and Sp1 level were associated with HCC patient survival.
Bak expression required ZBP-89 and Sp1 cooperative regulation simultaneously.
The Specificity protein (Sp) family members are GC boxes and CACCC boxes binding transcription factors, which are reported to participate in development, differentiation and tumorigenesis [1, 2]. Sp1 is also a C2H2-type zinc fingers Krüppel-like factor, and it is the first identified transcription factor . Numerous studies showed that Sp1 not only acted as a basal transcription factor, but also contributed to the regulation of many vital cellular genes [4, 5]. Overexpressed Sp1 differentially regulated a large number of genes which promoted cancer development, angiogenesis and metastasis . Sp3 is very similar to Sp1 both in structure and regulatory function. Sp3 and Sp1 showed equal affinity to GC box, indicating the potential competition in genes expression regulation . Reversely, Sp3 mediated transcriptional repression when it binded to GC-box, because Sp3 contained a suppressive domain . Additionally, Sp3 lacked the multi-merization domain, and thus Sp3 did not have the ability to super-activate promoters by bridging multiple Sp binding sites .
ZBP-89, also known as ZNF148, is a Krüppel-type zinc-finger transcription factor that binds to GC-rich sequences to activate or suppress gene transcription . ZBP-89 is also a member of the C2H2 zinc finger family subclass, as well as Sp family members. It is known that ZBP-89 participates in many genes transcription, such as Vimentin, Gastrin, p16, and Ornithine decarboxylase et al. These biological procedures are involved in cell growth, cell cycles, metabolism and T cell immunity. Study showed ZBP-89 competed with Sp1 for binding to induce Gastrin expression when ZBP-89 functions as a repressor . Interestingly, because Sp1 itself did not bind to the promoter element, ZBP-89 interacted with Sp1 to suppress Vimentin expression in vitro . In another case, ZBP-89 binded to the basal promoter of Pdcd4, also interacted with Sp family members to induce Pdcd4 protein expression .
Bak gene belongs to a large Bcl2 family, and acts as a pro-apoptosis protein to facilitate cellular apoptosis. It penetrates and makes pores in the mitochondrial membrane through forming heterodimer with Bax [14, 15]. The penetration releases cyto C into cytoplasm, and classic apoptosis signal pathway in the end [16, 17]. Currently, several cis-elements and associated factors are required to regulate expression of the bak gene. It was reported that Sp3 competed with Sp1 to transcriptional binding sites to increase Bak expression . Up-regulation of Bak by butyrate in the colon was associated with increased Sp3 binding which did not contain a multimerization domain . Instead, it contains an inhibitory domain located at N-terminal to its DNA-binding domain, which mediates transcriptional repression in some instances. Mithramycin A (MMA) is a GC-rich DNA binding chemical, which is also regarded as potent Sp family members inhibitor . In the Bak promoter analysis, there are many GC-rich sequence domain, indicating the possibility regulation by MMA,Sp family and ZBP-89 transcription factor. Previously, we reported that ZBP-89 bound to Bak gene promoter, and epigenetically regulated its transcription companied with histone deacetylase 3 and DNA methyltransferase 1 . Therefore, we speculate that ZBP-89 and Sp1/Sp3 are involved in hepatocellular carcinoma (HCC) and Bak induction. This study focuses on the elucidation of transcription factor Sp1/Sp3 and ZBP-89 involved Bak expression.
Patient cohort and mRNA quantification data
HCC patient cohort information and data of mRNA quantification data were downloaded from The Cancer Genome Atlas (TCGA) database (https://portal.gdc.cancer.gov/). Previously, we published the work how we got and analyzed the patient cohort and data . Totally 377 HCC patients with detailed mRNA expression level were collected. Additionally, 41 patients have non-tumor liver sample message which indicated 41 paired of tumor and non-tumor data were also included. Due to the missing follow-up survival data of some patients, we excluded those patients in the analysis of protein associated survival. Actually, 318 patients were collected for Kaplan–Meier survival analysis. Referring to published paper , cut-off point for gene grouping was set up and for survival analysis.
Chemicals and drugs
Antibodies against ZBP-89 (#sc-398148), Sp1 (#sc-420) and Sp3 (sc-28305) were obtained from Santa Cruz Biotechnology (Santa Cruz, CA). Antibodies against Bak (#12105) and Actin (#3700) were obtained from Cell Signaling Technology (Beverly, MA). Ad-ZBP-89 viral vector was a generous gift from Dr. JL Merchant (University of Michigan, MI) . Mithramycin A (MMA) was purchase from Cayman Inc. (#11434-1 mg). All other chemicals and reagents were purchased from Sigma-Aldrich (St. Louis, MO) unless otherwise specified.
HCC cell lines, HepG2 (ATCC® HB-8065™), Hep3B (ATCC® HB-8064™), PLC/PRF/5 (ATCC® CRL-8024™) and Huh-7 (ATCC® PTA-8561™) were purchased from ATCC, and kept in our department. LO2 (also called L-02, Accession number: CVCL_6926, # GNHu 6 of Cellbank of China) is a spontaneously immortalized cell line of hepatocyte, which was a gift from Dr. Zhou of Tongji Medical College, Huazhong University of Science and Technology. MIHA is also an immortalized cell line of hepatocyte and is a TGF-β-sensitive hepatocyte cell. Studies have revealed that the validation of MIHA as a valuable immortal hepatocyte in research by Prof. Ming-liang He (the Chinese University of Hong Kong) . MIHA cells were a generous gift from Prof. Ming-liang He. HKCI-3 and HKCI-4 were established by researchers of The Chinese University of Hong Kong affiliated Prince of Wales Hospital, and kept in the Dept. of Surgery lab . HepG2 and Huh-7 were cultured with MEM medium, supplemented with 10% FBS. Hep3B, PLC/PRF/5 and HKCI-4 were cultured with high glucose DMEM, supplemented with 10% FBS.
Western blotting was performed according to the previous publications with some modification [18, 19]. Briefly, cells were lysed by incubation on ice for 15 min in RIPA buffer to isolate total proteins. Proteins (30 μg)were separated on 10–15% SDS-PAGE gels and transferred to PVDF for immunoblotting. After membrane was blocked with 5% nonfat milk, the membrane was incubated with primary antibody for 2 h at room temperature, followed by incubation with a HRP-conjugated secondary antibody for 1 h at room temperature. The signal was visualized by ECL reagent kit (GE Healthcare) and analyzed by an Odyssey Scanner(Li-COR Biosciences Co.).
Small RNA interference
The siRNA interference experiments were performed according to our previous publications [20, 24]. ZBP-89 siRNA and Sp1/Sp3 siRNA were obtained from Santa Cruz Biotechnology. SiRNA oligonucleotides were transfected using Fugene lipid (Roche) according to the manufacturer’s protocol. Before transfection, old medium was replaced with fresh medium. The transfection volume was 1/10 of culture medium, such as 200 μL for 2 mL culture medium in six-well plate. We used 1 nM siRNA for one million cells, and make up the volume to 200 μL by adding OPTI-MEM.
All results are expressed as Mean ± SE of at least three independent experiments. Differences between groups were examined for statistical significance using Mann-Whitney (U-test) unless indicated specially. P<0.05 was regarded as statistically significant difference.
Bak transcription association with Sp1 and Sp3 expression
In the analysis of TCGA database, the comparison of bak transcription level between 40 cases of cancer tissue and the corresponding non-cancer tissues showed increased bak transcription level (Fig. 1a, p< 0.001). We analyzed the correlation of Bak with ZBP-89 and Sp1–4 using SPSS correlation method. Results showed only Sp1 and Sp3 transcription level were correlated with Bak transcription level (Fig. 1b and c, r=0.1889 and r = 0.1119, respectively). No correlation was found in ZBP-89 and Bak (Fig. 1d), neither Sp2 nor Sp4 (data not shown). The increased transcription level of Bak was contributed to the chemotherapy drug induction or anti-cancer related ROS apoptosis induction. The clinical treatment for patients was enclosed as Additional file 1: Table S1.
MMA increases Bak expression through blocking Sp-proteins binding activity
MMA is regarded as a potent Sp family inhibitor through displacing Sp1 binding to gene promoters. In the study, we used MMA as specific Sp proteins inhibitor. We found that increases Bak expression starting from the concentration of 5 nM in PLC/PRF/5 cells in 48 h incubation, and induces cleaved PARP level, which is a main mediator of apoptosis (Fig. 2a). Gray scale density quantification data indicated that Bak expression trend was accompanied with MMA increasing from 5 nM to 100 nM, with a peak value about 50% up-regulation compared with baseline (Fig. 2b, * p < 0.05 compared with control). The Bak induction rate was tested in liver immortal non-tumour cell lines and cancer cell lines. In the tested six HCC cell lines, MMA up-regulated Bak expression (C, column 3 to 8, * p < 0.05 compared with control), especially in PLC5/PRF/5 and Hep3B. In immortal cells, MIHA and LO2, which were less sensitive to MMA treatment, also showed increased Bak expression (C, column 1 and 2, * p < 0.05 compared with control). In addition, RT-PCR data showed Bak mRNA level increased significantly in HKCI-4, PLC/PRG/5 and Hep3B cells (Fig. 2d, * p < 0.05 compared with control). Therefore, we selected HKCI-4 and PLC/PRF/5 cancer cells as model in the subsequent study, as well as liver immortal non-tumour cells MIHA and LO2.
Sp1 increases Bak expression but not Sp3
We constructed Sp1 and Sp3 expression vectors by cloning the Sp1 and Sp3 coding sequences into pcDNA3.0 respectively. When MIHA and PLC/PRF/5 cells were in density of 80% confluent, recombinant gene plasmids was transfected into cell via Fugene lipid (Roche). Western blotting results showed Sp1 overexpression increased Bak expression both in MIHA and PLC/PRF/5 cells (Fig. 3a and b, lane 1, * p < 0.05 compared with control), as well as in the transfection of LO2 and HKCI-4 cells (Fig. 3c and d, lane 1, * p < 0.05 compared with control). We did not observe that Sp3 transfection induced Bak expression. In the co-transfection of Sp1 and Sp3 vectors, Sp1 and Sp3 do not compete with binding sites. Therefore, the outcome turns to up-regulate Bak level. The results were further confirmed by RT-PCR experiment, results showed that Sp1 overexpression induced Bak up-regulation but not Sp3 (Fig. 3e, * p < 0.05 compared with control).
siRNA interference Sp1 decreases Bak expression
We used Sp1 siRNA nucleotide to further explore its role in Bak expression. When Sp1 siRNA was transfected in MIHA, PLC/PRF/5 and HKCI-4 cells respectively, we found that Bak expression level reduction was accompanied with the siRNA interference of Sp1 (Fig. 4a and b, lane 1, * p < 0.05 compared with corresponding control). Results confirmed that Bak expression level was mainly associated with Sp1 protein level and transcriptive activity. RT-PCT experiment confirmed such findings (Fig. 4c, * p < 0.05 compared with corresponding control).
ZBP-89 cooperates with Sp1/Sp3 to facilitate Bak expression
Previously, our team elucidated the epigenetic mechanism of ZBP-89 mediated Bak increasing expression via HDACs and DNMT1 activity suppression [20, 24]. In the above experiments, we have already revealed that Bak expression could be up-regulated in the presence of MMA treatment, and Sp1 vector mediated overexpression. We suspected that Bak expression might be tightly related with ZBP-89 and Sp1/Sp3 protein level. Therefore, we designed the experiment to see if ZBP-89 could interact with Sp1/Sp3 to regulate Bak expression. Western blotting data showed that siRNA blocking ZBP-89 level decreased Bak expression even with treatment of MMA (Fig. 5a and b, lane 5). Additionally, ZBP-89 siRNA interference also turned down Bak expression even though Sp1 overexpressed (Fig. 5, lane 8). In the subsequent RT-PCR analysis, data also demonstrated that ZBB-89 siRNA suppressed Bak expression in the presence of MMA treatment or Sp1 overexpression (Fig. 5c, * p < 0.05 compared with control). Collectively, results revealed that Bak expression required ZBP-89 and Sp1 to work simultaneously.
Association of Bak, Sp1/Sp3 and ZBP-89 in patient survival
In order to explore if the expression of Bak, Sp1, Sp3 and ZBP-89 are associated with patients survival, we performed Kaplan–Meier survival analysis. Firstly, ROC curve was used to dichotomize these genes into high expression group and low expression group of all 318 HCC patients from TCGA database. Survival analysis was analyzed using Graphpad Prism 5.0 software, and p < 0.05 was set as statistic significant difference level. Results showed that only Bak and Sp1 expression were significantly correlated with patients survival (Fig. 6, p < 0.05, Log-rank test). Results indicated that high expression of Bak was associated with shorter living time. We speculated that high expression Bak might be induced by anti-cancer drug treatment, which damaged hepatocellular function and structure. In the end, patients died of liver dysfunction, multi-organs malfunction et al.
Bak is a major cell death initiator in the apoptotic signaling cascade . It plays a significant role in the therapy and survival of HCC. It has been reported that different agents induce apoptosis in HCC cells by inducing Bak expression, but the mechanism responsible for its regulation in HCC is largely unknown.
Sp1 is a ubiquitous transcription factor and its regulated genes include those involved in proliferation, apoptosis, senescence and angiogenesis. Sp3, which competes with Sp1 for GC-box binding sites, is normally a weak transcriptional activator. However, it could function as a transcriptional activator with similar potency to Sp1 in the absence of acetyltransferases . A similar switch of Sp1 for Sp3 has been observed at the Bak promoter following treatment with the HDAC inhibitors butyrate . However, the relationship between Sp1, Sp3 and Bak is unknown. After the analysis of TCGA HCC database, we observed that Bak transcription in cancer tissue was higher than that in non-cancer tissues, and its transcription was correlated with Sp1 and Sp3 transcription level. No correlation was found between Bak and ZBP-89, neither Sp2 nor Sp4 (data not shown). Further analysis found Bak and Sp1 transcription level were significantly correlated with patients survival. These results illustrated that Bak played an important role in the hepatocarcinogenesis or HCC therapy. In Western blotting analyses, we showed that all selected cell lines (HepG2, Hep3B, PLC/PRF/5,Huh-7, HKCI3, HKCI-4, LO2, and MIHA) were able to constitutively express endogenous Bak. The base level of Bak can be greatly enhanced by MMA (30 nM) treatment and the increase in Bak is functional in terms of apoptotic induction in PLC/PRF/5 cells.
MMA functions as GC-rich DNA sequence binding chemical, acting as a Sp family members’ inhibitor. It was shown to induce apoptosis through increase of Bak expression in colon cancer cells . We observed a regulation of Bak protein level after overexpression or knockdown of Sp1. Therefore we raised the question whether ZBP-89 could compete or cooperate with Sp1 to regulate the expression of the Bak gene. Indeed, Western blotting data showed MMA could not induce Bak expression if without ZBP-89 participation (Fig. 5). Additionally, when cells was treated with ZBP-89 siRNA, the consequence also turned down Bak expression even if Sp1 overexpression (Fig. 5, lane 1, 5th band). Results revealed that Bak expression required ZBP-89 and Sp1 to work simultaneously.
ZBP-89 and Sp family members are members of the zinc finger transcription factor and also subclasses of the C2H2 zinc finger family. ZBP-89 can efficiently inhibit cellular proliferation and induce apoptosis in human cancer cells . In our previous study, we have identified that ZBP-89 might target at − 457 to − 407 region of Bak promoter to induce its expression and subsequently apoptosis . Given the same bifunctional regulatory domains, ZBP-89 and Sp1 can interact to regulate a variety of genes. ZBP-89 appears to interact with Sp1, as binding to an overlapping site in DNA to regulate gastrin, two type I collagen genes , and ornithine decarboxylase . In another aspect, many findings also suggested that Sp1 and ZBP-89 regulated target genes in a cooperative manner, like vimentin , Pdcd4  and p21waf1/cip1 . In agreement with peers’ finding, we further showed the cooperative function of ZBP-89 and Sp1 on Bak transcription. Besides, the study also revealed their association in patient survival, as well as the Bak expression level.
Our findings suggested that ZBP-89 and Sp1 overexpression induced Bak expression in a genetic manner. Increased Bak level was associated with poor patient survival, whereas high level of Sp1 is a beneficial factor for patient survival.
Sp1 transcription factor
The Cancer Genome Atlas
Zinc binding protein 89
Shields JM, Yang VW. Identification of the DNA sequence that interacts with the gut-enriched Kruppel-like factor. Nucleic Acids Res. 1998;26(3):796–802.
Matsumoto N, Laub F, Aldabe R, Zhang W, Ramirez F, Yoshida T, Terada M. Cloning the cDNA for a new human zinc finger protein defines a group of closely related Kruppel-like transcription factors. J Biol Chem. 1998;273(43):28229–37.
Kaczynski J, Cook T, Urrutia R. Sp1- and Kruppel-like transcription factors. Genome Biol. 2003;4(2):206.
Oleaga C, Welten S, Belloc A, Sole A, Rodriguez L, Mencia N, Selga E, Tapias A, Noe V, Ciudad CJ. Identification of novel Sp1 targets involved in proliferation and cancer by functional genomics. Biochem Pharmacol. 2012;84(12):1581–91.
Gilmour J, Assi SA, Jaegle U, Kulu D, van de Werken H, Clarke D, Westhead DR, Philipsen S, Bonifer C. A crucial role for the ubiquitously expressed transcription factor Sp1 at early stages of hematopoietic specification. Development. 2014;141(12):2391–401.
Beishline K, Azizkhan-Clifford J. Sp1 and the 'hallmarks of cancer. FEBS J. 2015;282(2):224–58.
Wierstra I. Sp1: emerging roles--beyond constitutive activation of TATA-less housekeeping genes. Biochem Biophys Res Commun. 2008;372:1):1–13.
Dennig J, Beato M, Suske G. An inhibitor domain in Sp3 regulates its glutamine-rich activation domains. EMBO J. 1996;15(20):5659–67.
Mastrangelo IA, Courey AJ, Wall JS, Jackson SP, Hough PV. DNA looping and Sp1 multimer links: a mechanism for transcriptional synergism and enhancement. Proc Natl Acad Sci U S A. 1991;88(13):5670–4.
Merchant JL, Bai L, Okada M. ZBP-89 mediates butyrate regulation of gene expression. J Nutr. 2003;133(7 Suppl):2456S–60S.
Merchant JL, Iyer GR, Taylor BR, Kitchen JR, Mortensen ER, Wang Z, Flintoft RJ, Michel JB, Bassel-Duby R. ZBP-89, a Kruppel-like zinc finger protein, inhibits epidermal growth factor induction of the gastrin promoter. Mol Cell Biol. 1996;16(12):6644–53.
Wieczorek E, Lin Z, Perkins EB, Law DJ, Merchant JL, Zehner ZE. The zinc finger repressor, ZBP-89, binds to the silencer element of the human vimentin gene and complexes with the transcriptional activator, Sp1. J Biol Chem. 2000;275(17):12879–88.
Leupold JH, Asangani IA, Mudduluru G, Allgayer H. Promoter cloning and characterization of the human programmed cell death protein 4 (pdcd4) gene: evidence for ZBP-89 and Sp-binding motifs as essential Pdcd4 regulators. Biosci Rep. 2012;32(3):281–97.
Antignani A, Youle RJ. How do Bax and Bak lead to permeabilization of the outer mitochondrial membrane? Curr Opin Cell Biol. 2006;18(6):685–9.
Lomonosova E, Chinnadurai G. BH3-only proteins in apoptosis and beyond: an overview. Oncogene. 2008;27(Suppl 1):S2–19.
Degli Esposti M, Dive C. Mitochondrial membrane permeabilisation by Bax/Bak. Biochem Biophys Res Commun. 2003;304(3):455–61.
Westphal D, Dewson G, Czabotar PE, Kluck RM. Molecular biology of Bax and Bak activation and action. Biochim Biophys Acta. 2011;1813(4):521–31.
Chirakkal H, Leech SH, Brookes KE, Prais AL, Waby JS, Corfe BM. Upregulation of BAK by butyrate in the colon is associated with increased Sp3 binding. Oncogene. 2006;25(54):7192–200.
Sleiman SF, Berlin J, Basso M, Karuppagounder SS, Rohr J, Ratan RR. Histone deacetylase inhibitors and Mithramycin a impact a similar neuroprotective pathway at a crossroad between Cancer and neurodegeneration. Pharmaceuticals. 2011;4(8):1183–95.
Ye CG, Chen GG, Ho RL, Merchant JL, He ML, Lai PB. Epigenetic upregulation of Bak by ZBP-89 inhibits the growth of hepatocellular carcinoma. Biochim Biophys Acta. 2013;1833(12):2970–9.
Lu M, Kong X, Wang H, Huang G, Ye C, He Z. A novel microRNAs expression signature for hepatocellular carcinoma diagnosis and prognosis. Oncotarget. 2017;8(5):8775–84.
Zhang JF, He ML, Fu WM, Wang H, Chen LZ, Zhu X, Chen Y, Xie D, Lai P, Chen G, et al. Primate-specific microRNA-637 inhibits tumorigenesis in hepatocellular carcinoma by disrupting signal transducer and activator of transcription 3 signaling. Hepatology. 2011;54(6):2137–48.
Pang E, Wong N, Lai PB, To KF, Lau WY, Johnson PJ. Consistent chromosome 10 rearrangements in four newly established human hepatocellular carcinoma cell lines. Genes, chromosomes & cancer. 2002;33(2):150–9.
To AK, Chen GG, Chan UP, Ye C, Yun JP, Ho RL, Tessier A, Merchant JL, Lai PB. ZBP-89 enhances Bak expression and causes apoptosis in hepatocellular carcinoma cells. Biochim Biophys Acta. 2011;1813(1):222–30.
Cory S, Adams JM. The Bcl2 family: regulators of the cellular life-or-death switch. Nat Rev Cancer. 2002;2(9):647–56.
Steiner E, Holzmann K, Pirker C, Elbling L, Micksche M, Berger W. SP-transcription factors are involved in basal MVP promoter activity and its stimulation by HDAC inhibitors. Biochem Biophys Res Commun. 2004;317(1):235–43.
Gillissen B, Richter A, Richter A, Overkamp T, Essmann F, Hemmati PG, Preissner R, Belka C, Daniel PT. Targeted therapy of the XIAP/proteasome pathway overcomes TRAIL-resistance in carcinoma by switching apoptosis signaling to a Bax/Bak-independent 'type I' mode. Cell Death Dis. 2013;4:e643.
Law DJ, Labut EM, Merchant JL. Intestinal overexpression of ZNF148 suppresses ApcMin/+ neoplasia. Mammalian genome. 2006;17(10):999–1004.
Hasegawa T, Takeuchi A, Miyaishi O, Isobe K, de Crombrugghe B. Cloning and characterization of a transcription factor that binds to the proximal promoters of the two mouse type I collagen genes. J Biol Chem. 1997;272(8):4915–23.
Law GL, Itoh H, Law DJ, Mize GJ, Merchant JL, Morris DR. Transcription factor ZBP-89 regulates the activity of the ornithine decarboxylase promoter. J Biol Chem. 1998;273(32):19955–64.
Bai L, Merchant JL. Transcription factor ZBP-89 cooperates with histone acetyltransferase p300 during butyrate activation of p21waf1 transcription in human cells. J Biol Chem. 2000;275(39):30725–33.
This study, including design, collection data, analysis, experiments and interpretation of data, was supported by the National Natural Science Foundation of China (No.81572782 and No. 81502411). Writing of the manuscript was supported by Special Talent Supportive Grant (No. 4TZ160001G).
Availability of data and materials
Some of the datasets generated and analysis during the current study are available in the TCGA database https://portal.gdc.cancer.gov/, the rest can be obtained from corresponding author on reasonable request.
Ethics approval and consent to participate
Ethics approval and participant consent was not necessary as this study involved the use of a previously-published TCGA database, which is an open and free access repository according to the copyright and data disclaimers of National Cancer Institute of U.S.A. All usage of cell lines did not require Ethic approval.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Table summary of available patient and chemotherapy information. The supplemental table showed individual HCC patient and treatment one received, and also the therapy duration time since the initial treatment. (XLSX 12 kb)